本人有一系列的纵向数据,时间点(X)和每个人在对应时间点的测量值(Y)。文献方法描述:linear mixedmodels with a cubic spline mean structure and a random effects structure。本人采用SAS proc glimmix 和effect(时间点设置了3个knot points),因变量是测量值(Y),自变量是时间点(X),截距和斜率均random。
但是,问题:
1)本人没法进一步像文献一样做出10分位数,50分位数,90分位数的估计值,如表所示。
2)没法做出图。
图表均从文献截图
表注释:表中时间点列对应图横轴,10th,50th, 90th分别代表分位数,对应表中三条红线。(图和表数据不完全重合)
图注释:横轴代表时间点,纵轴代表测量值,三条红线有低到高分别代表10分位/50分位/90分位估计值,忽略黑色虚线。
请求大神指导,十分感谢。