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这个问题解决了吗？我也遇到了同样的问题，求解答