gls模型中是否有解决过度分散的方法?

时间:2019-04-23 15:04:35

标签: r nlme

我有一些自相关数据,显示出正线性增长。当我使用gls对其建模时,我认为摘要显示出过度分散。

当使用GLMM等时,我会更改错误结构,但我认为我不能一眼就能改变家庭。 gls模型中是否有解决过度分散的方法?

library(nlme)    

eggs<-c(4087,2761,3807,4158,2046,4757,2984,3316,3143,
3042,4429,3335,5124,2464,3713,3028,5739,4671,3799,6167,2937,5031)

year<-seq(1997,2018,1)

m1<-glm(nest~y)

summary(m1)

acf(resid(m1),type="p")

残差在一年后显示自相关,所以我需要使用gls

m1.gls<-gls(nest~y,correlation=corARMA(p=1), method="ML")

summary(m1.gls)

acf(resid(m1.gls),type="p")

m0.gls<-gls(nest~y,correlation=NULL,method="ML")

AIC(m0.gls,m1.gls)

anova(m1.gls,m0.gls)

#L.Ratio = chisqu = 7.382355; p = 0.0066
#autocorrelation definitely significant

null.gls<-gls(nest~1,correlation=corARMA(p=1), method="ML")

anova(m1.gls,null.gls)

#L.Ratio = chisq = 7.956113; p = 0.0048
#linear trend is significant

summary(m1.gls)

#residual se = 970 over 22 df
#does this show overdispersion?

这个模型过度分散是不可接受的吗?我可以尝试其他适合更好的模型吗?

0 个答案:

没有答案