R语言GAMLSS问题

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

本人想做个医学参考值范围应用了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)，谢谢

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏1 回帖

关键词：GAM R语言 Hannan iGroup normal