错误代码:Error in lm.wfit(X[onlydata, , drop = FALSE], y, w) :incompatible dimensions
原代码:library(gamlss)
library(readxl)
setwd("C:/Users/c000000wow/Desktop")
data=read_xlsx("early.xlsx")
newdata<-na.omit(data)
week<-data$week
long<-data$head hips long
gam <- gamlss(head hips long
cs(week,df=3),sigma.fo=cs(week),data=newdata, family=BCT)
summary(gam)
plot(gam)
centiles(gam,week,legend = FALSE,ylab = "Length",xlab = "Week",main = "",ylim=c(0,100))
result <- centiles.pred(gam,xname="week",xvalues=seq(8,15,1))