我有一些自相关数据,显示出正线性增长。当我使用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?
这个模型过度分散是不可接受的吗?我可以尝试其他适合更好的模型吗?