可能是一个愚蠢的问题,但经过数小时寻找答案后,我仍然没有找到解决方案。我已经运行了一个 gamm 模型,该模型似乎工作正常。我现在想使用模型(mod1)中计算的对数似然来计算模型的 AIC。
模型运行如下
mod1<-gamm(Biomass..g. ~ Atoll + SizeCm + Depth + OuterOrLagoon,random=list(SiteID=~1),family=gaussian,data=RMI2014, method = "REML")
然后在模型运行后,我尝试使用以下代码获取 AIC,该代码返回后续错误
> AIC(mod1, k=2)
Error in UseMethod("logLik") :
no applicable method for 'logLik' applied to an object of class "c('gamm', 'list')"
最初我认为这是错误的,因为我使用的是 REML 而不是 Max Likelihood,但不使用 REML 时会返回相同的错误。
是因为我在返回错误的伽玛中有一个随机因素吗?
请帮忙!
干杯,泰森