Displaying 1 result from an estimated 1 matches for "model_zinbi".
2010 Jun 02
1
Problems using gamlss to model zero-inflated and overdispersed count data: "the global deviance is increasing"
...df=3) + cs(HHHDI1000,df=3) + cs(LFAP200,df=3),data=data,family= NBII)
GAMLSS-RS iteration 1: Global Deviance = 284.5993
GAMLSS-RS iteration 2: Global Deviance = 281.9548
######......######
GAMLSS-RS iteration 5: Global Deviance = 280.7311
GAMLSS-RS iteration 15: Global Deviance = 280.6343
> model_ZINBI <- gamlss(duck ~ cs(HHCDI200,df=3) + cs(HHCDI1000,df=3) + cs(HHHDI200,df=3) + cs(HHHDI1000,df=3) + cs(LFAP200,df=3),data=data,family= ZINBI)
GAMLSS-RS iteration 1: Global Deviance = 1672.234
GAMLSS-RS iteration 2: Global Deviance = 544.742
GAMLSS-RS iteration 3: Global Deviance = 598.9939
Err...