## Not run:
# data(kidneySimTimeGroup)
# # running for samples in group 1
# G1 <- which(kidneySimTimeGroup$group=="G1")
# testLmmspline <- lmmSpline(data=kidneySimTimeGroup$data[G1,],
# time=kidneySimTimeGroup$time[G1],
# sampleID=kidneySimTimeGroup$sampleID[G1],keepModels=T)
#
#
# plot(testLmmspline, y=2)
# plot(testLmmspline, y=2, smooth=TRUE)
# # Don't keep the models to improve memory usage
# testLmmspline <- lmmSpline(data=kidneySimTimeGroup$data[G1,],
# time=kidneySimTimeGroup$time[G1],
# sampleID=kidneySimTimeGroup$sampleID[G1],keepModels=F)
#
# #plot only the fitted values
# plot(testLmmspline, y=2)
# #plot fitted values with original data
# plot(testLmmspline, y=2, data=kidneySimTimeGroup$data[G1,], time=kidneySimTimeGroup$time[G1])
#
# ## End(Not run)
Run the code above in your browser using DataLab