Fatal编程技术网
  • r/
  • 混合模型残差的Durbin-Watson检验

混合模型残差的Durbin-Watson检验

r

混合模型残差的Durbin-Watson检验,r,mixed-models,temporal,autocorrelation,R,Mixed Models,Temporal,Autocorrelation,我试图计算一个有毒分布的混合模型的时间自相关，我想知道如何计算。我得到一个错误，上面写着“$operator not defined for this S4 class”，我可以在一个具有泊松分布的线性模型上成功地运行dwtest，但不是我真正想要的 成功的模型和代码： temp.nem.cuc.glm<-glm(AllDat$nem.cuc~ AllDat$year.collected, family=poisson(link="log")) summary(temp.nem.cuc.l

我试图计算一个有毒分布的混合模型的时间自相关，我想知道如何计算。我得到一个错误，上面写着“$operator not defined for this S4 class”，我可以在一个具有泊松分布的线性模型上成功地运行dwtest，但不是我真正想要的

成功的模型和代码：

temp.nem.cuc.glm<-glm(AllDat$nem.cuc~ AllDat$year.collected, family=poisson(link="log"))
summary(temp.nem.cuc.lm) 
time<-AllDat$year.collected
dwnem.cuc<-dwtest(temp.nem.cuc.lm, order.by = time, alternative = "two.sided", iterations = 50, exact = FALSE, tol = 1e-10)
dwnem.cuc

temp.nem.cuc.glmwhat package
dwtest
from？请给我们一个可复制的例子好吗？请不要从StackOverflow给我发邮件（除非特别要求）。相反，通过编辑您的问题发布您的可复制示例。dwtest是lmtest包的一部分。dwtest来自哪个包？请给我们一个可复制的例子好吗？请不要从StackOverflow给我发邮件（除非特别要求）。相反，通过编辑您的问题发布可复制的示例。dwtest是lmtest包的一部分。

#the model I am really interested in
nem.cuc.pois=glmer(nem.cuc~ I(year.collected-1930)+I(standard.length..mm./100) + (1|sites1), family = "poisson", data=AllDat)
time<-AllDat$year.collected
dwnemresid.cuc<-dwtest(nem.cuc.pois, order.by = time, alternative = "two.sided", iterations = 50, exact = FALSE, tol = 1e-10)
dwnem.cuc