I think this should be a relatively general solution, doing some robustness checks
model_lin <- ' i =~ 1*t1 + 1*t2 + 1*t3 + 1*t4
s =~ 0*t1 + 1*t2 + 2*t3 + 3*t4 '
fitlin <- growth(model_lin, data=Demo.growth)
summary(fitlin)
mod_loadings_lin<-fitlin@Model@GLIST$lambda
factorscores_lin<-as.matrix(predict(fitlin))
new_growthdat_pred_linear<-t(mod_loadings%*%t(factorscores_lin))