> 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

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