r - 有没有办法解决 gls 模型中的过度分散问题?
问题描述
我有自相关数据显示正线性增加。当我使用 gls 对它们进行建模时,我认为摘要显示过度分散。
使用 GLMM 等时,我会更改错误结构,但我认为我无法更改 gls 中的家庭。有没有办法解决 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?
这个模型是否过度分散?有没有我可以尝试更合适的模型?
解决方案
推荐阅读
- python - 如何在 Pandas 中制作简单的频率表
- recursion - Rust 中的递归生成器
- r - 如何删除折线图(geom_line)R ggplot2上的一个点?
- java - 错误:找不到方法实现()
- python - 绘图线下方的所有网格线 - 带有双 x 轴 [Matplotlib]
- javascript - React js:图像上传在 CKEditor 中不起作用
- css - 动画反应模态
- arrays - 如何在MATLAB中为数组创建一个减法循环?
- sql - Where 子句 PL/SQL 上的标识符无效
- python - 为什么是“从 CAUSE 中引发异常”,而不是“从 CAUSE 中引发异常”?