data(heights)
## fit model
m1 <- sitar(x=age, y=height, id=id, data=heights, df=5)
## predictions at level 0
predict(m1, newdata=data.frame(x=9:16), level=0)
## predictions at level 1 for subject 5
predict(m1, newdata=data.frame(x=9:16, id=5), level=1)Run the code above in your browser using DataLab