本人想做个医学参考值范围应用了GAMLSS模型，现在已经选好分布和模型，但是想知道p5 p25 p50 p75 p95百分位数，
我只会写画图的程序centiles(hannanBUA,xvar=hannan$age1,cent=c(2.5,5,25,50,75,95,97.5),col.cent=c(2,3,4,5,4,3,2))，请问算百分位数的程序是什么？
这是我选择的模型，agegroup和bmigroup都是自变量，BUAnormal是因变量，hannanBUA<-gamlss(BUAnormal~cs(log(agegroup),df=0)+cs(log(bmigroup),df=0),sigma.fo=~cs(agegroup,df=0)+cs(log(bmigroup),df=0),nu.formula=~cs(agegroup,df=0),tau.formula=~1,data=hannan,family=BCT)，谢谢