本人有一系列的纵向数据，时间点（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分位估计值，忽略黑色虚线。

请求大神指导，十分感谢。