我正在使用gamlss包中的示例代码绘制百分位曲线:
library(gamlss)
data(abdom)
lms(y,x , data=abdom, n.cyc=30)
它绘制了自己的百分位曲线集。如何选择仅绘制第10,第50和第90百分位曲线?此外,我想避免绘制点,以便只绘制曲线。 谢谢你的帮助。
答案 0 :(得分:4)
阅读帮助页面总是一个好主意:
> centiles(h,xvar=abdom$x, cent=c(10,50,90), points=FALSE)
% of cases below 10 centile is 8.688525
% of cases below 50 centile is 50.16393
% of cases below 90 centile is 90