powered by
data(usair)
data(usair) str(usair) plot(usair) # a possible gamlss model # gamlss(library) #ap<-gamlss(y~cs(x1,2)+x2+x3+cs(x4,2)+x5+cs(x6,3)+x4:x5, # data=usair, family=GA(mu.link="inverse")) #
Run the code above in your browser using DataLab