Dear R experts, I am tring to run GAMLSS to model the association between the two variables (x,y) and I am getting the following error message: *** Initial fit*** GAMLSS-RS iteration 1: Global Deviance = -277.9826 GAMLSS-RS iteration 2: Global Deviance = -277.9826 *** Fitting BCCGo *** BCCGo failed *** Fitting BCPEo *** BCPEo failed *** Fitting BCTo *** BCTo failed *** Refitting NO *** GAMLSS-RS iteration 1: Global Deviance = -289.9546 GAMLSS-RS iteration 2: Global Deviance = -297.5721 GAMLSS-RS iteration 3: Global Deviance = -300.6583 GAMLSS-RS iteration 4: Global Deviance = -302.7118 GAMLSS-RS iteration 5: Global Deviance = -304.5811 GAMLSS-RS iteration 6: Global Deviance = -306.7035 GAMLSS-RS iteration 7: Global Deviance = -309.4601 GAMLSS-RS iteration 8: Global Deviance = -313.0182 GAMLSS-RS iteration 9: Global Deviance = -317.0507 GAMLSS-RS iteration 10: Global Deviance = -321.1582 GAMLSS-RS iteration 11: Global Deviance = -325.0426 GAMLSS-RS iteration 12: Global Deviance = -328.6318 GAMLSS-RS iteration 13: Global Deviance = -334.9106 GAMLSS-RS iteration 14: Global Deviance = -352.5725 Error in glim.fit(f = sigma.object, X = sigma.X, y = y, w = w, fv = sigma, : object 'fit' not found In addition: Warning messages: 1: In is.na(data) : is.na() applied to non-(list or vector) of type 'NULL' 2: In is.na(data) : is.na() applied to non-(list or vector) of type 'NULL' This is my code which generated the error message. library(gamlss) x<-c(13.25,8.91,11.32,8.04,9.23,11.69,8.19,13.28,10.84,13.63,14.65,15.71,15.28,10.862,6.533,9.231,7.109,16.31,11.24,17.7,14.2,10.19,12.25,9.27,16.93,8.28,8.45,16.08,9.92,11.91,14.79,18.77,19.78,10.75,7.68,11.92,15.27,24.82,11.58,15.21,10.74,8.191,8.91,7.17,12.59,11.56,11.31,9.35,10.9,8.3,10.82,10.21,14.18,16.28,8.56,8.591,11.4,22.76,10.52,25.34,10.541,14.43,19.23,14.36,11.34,8.164,7.958,8.15,8.887,10.68,11.32,9.529,7.604,21.15,23.22,23.35,20.04,8.5,8.52,10.73,13.18,7.26,12.1,12.23,11.26,8.91,9.001,9,12.9,9.54,11.2,9.48,11.09) y<-c(0.27052,0.34393,0.40514,0.34321,0.30607,0.28391,0.37087,0.24221,0.24616,0.34823,0.39945,0.37309,0.30546,0.39562,0.27128,0.29521,0.26446,0.32227,0.3859,0.28562,0.29736,0.28727,0.34267,0.23611,0.4089,0.32197,0.3153,0.2651,0.31764,0.5713,0.37181,0.4212,0.29231,0.34072,0.28202,0.26942,0.46582,0.34381,0.30618,0.29451,0.32669,0.3671,0.38266,0.27047,0.345,0.38236,0.39403,0.38572,0.34314,0.3133,0.31352,0.32085,0.38288,0.35723,0.32862,0.28227,0.3791,0.35457,0.27427,0.33545,0.29433,0.27307,0.26075,0.33414,0.25645,0.32292,0.23796,0.37676,0.27803,0.4421,0.32455,0.31943,0.27655,0.32069,0.35513,0.36519,0.4359,0.30588,0.31528,0.33549,0.37415,0.27627,0.27067,0.26746,0.29562,0.3008,0.32832,0.2972,0.32896,0.36565,0.28812,0.34541,0.24058) out <- lms(y, x, data=NULL,family = "BCCG") Can you please help me in resolving the issue. Thank you very much for the help in advance. -- Thanks Warm Regards, Tanmay [[alternative HTML version deleted]]