Dear Richard I think the results below are consistent
set.seed(1020) # created from EXP X1 <- rEXP(1000) Gexp <- gamlss(X1~1,family=EXP) #GAMLSS-RS iteration 1: Global Deviance = 1999.762 #GAMLSS-RS iteration 2: Global Deviance = 1999.762 Glno <- gamlss(X1~1,family=LOGNO) #GAMLSS-RS iteration 1: Global Deviance = 2213.252 #GAMLSS-RS iteration 2: Global Deviance = 2213.252 GAIC(Gexp, Glno) # df AIC #Gexp 1 2001.762 #Glno 2 2217.252 # EXP best # create from LOGNO X2 <- rLOGNO(1000) Aexp <- gamlss(X2~1,family=EXP) #GAMLSS-RS iteration 1: Global Deviance = 2954.628 #GAMLSS-RS iteration 2: Global Deviance = 2954.628 Alno <- gamlss(X2~1,family=LOGNO) #GAMLSS-RS iteration 1: Global Deviance = 2796.725 3GAMLSS-RS iteration 2: Global Deviance = 2796.725 GAIC(Aexp, Alno) # df AIC #Alno 2 2800.725 #Aexp 1 2956.628 # LOGNO best Mikis Prof Mikis Stasinopoulos [email protected] Companies Act 2006 : http://www.londonmet.ac.uk/companyinfo ______________________________________________ [email protected] mailing list https://stat.ethz.ch/mailman/listinfo/r-help PLEASE do read the posting guide http://www.R-project.org/posting-guide.html and provide commented, minimal, self-contained, reproducible code.
