ID <- NULL
for(i in 1:50){ID <- c(ID,rep(i,10))}

time <- rep(seq(1,10),50)
time2 <- time^2

Y <- rnorm(500)


A <- NULL
A <- as.data.frame(cbind(ID,time1,time2,Y))

gmm <- hlme(Y ~ time + time2, ng = 3, mixture = Y ~ time + time2, subject = 'ID',data  = A)

summary(gmm)