大家好 我想用R语言gamlss做个参考值范围 我师兄给我分享了一段程序 有些不懂 麻烦解

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

gamlss(grip~cs(log(age),df=0)+height,sigma.fo=~1,data=mp,family=NO)这是原程序代码
grip代表握力的意思 age代表年龄 height代表身高
麻烦解释下cs（log（age），df=0）什么意思 sigma.fo=~1什么意思，family=no什么意思

扫码加我 拉你入群

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

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

分享0 收藏2 回帖

关键词：参考值范围 参考值 R语言 height family

https://cran.r-project.org/web/packages/gamlss/gamlss.pdf

page 39、40

xwan288 发表于 2018-5-19 16:26
https://cran.r-project.org/web/packages/gamlss/gamlss.pdf

page 39、40

> hannanBUA<-gamlss(BUAnormal~cs(log(age1),df=2),sigma.fo=~1,data=na.omit(hannan),family=BCCG)
GAMLSS-RS iteration 1: Global Deviance = 11445.14
GAMLSS-RS iteration 2: Global Deviance = 11445.12
GAMLSS-RS iteration 3: Global Deviance = 11445.12
> 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))
Error in if (any(mu < 0)) stop(paste("mu must be positive", "\n", "")) :
  missing value where TRUE/FALSE needed

我想做个参考值范围 程序提示错误 这是什么原因呢

can't comment without data

楼主你好，我跟你一样想做一个医学参考值。不知道你的问题解决了吗？能否请教你的分布模型是怎么确定的？各参数的自由度应该怎么得出？麻烦了，非常感谢！

楼主您好！Error in if (any(mu < 0)) stop(paste("mu must be positive", "\n", "")) :
  missing value where TRUE/FALSE needed这个问题解决了吗？我也遇到了同样的问题，求解答