Displaying 1 result from an estimated 1 matches for "hhhdi1000".
Did you mean:
hhcdi1000
2010 Jun 02
1
Problems using gamlss to model zero-inflated and overdispersed count data: "the global deviance is increasing"
...s with a df=3. The best fitting model will then be choses on the basis of its AIC.
However, these models frequently fail to converge, and I am not sure why, and what to do about it. For example:
> model_Poisson <- 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= PO)
GAMLSS-RS iteration 1: Global Deviance = 1350.623
GAMLSS-RS iteration 2: Global Deviance = 1350.623
> model_ZIPoisson <- gamlss(duck ~ cs(HHCDI200,df=3) + cs(HHCDI1000,df=3) + cs(HHHDI200,df=3) + cs(HHHDI1000,df=3) + cs(LFAP200,df=3),data=data...