无法从R中的gamm模型提取logLik

Question

可能是一个愚蠢的问题，但是经过数小时的搜索，我仍然没有找到解决方案。 我已经运行了一个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时会返回相同的错误。

是因为我在伽马中有一个随机因素导致返回错误？

请帮忙！

干杯，泰森

Answer 1

根据文档 ， gamm调用的结果是一个列表，其中包含两个类型为gam和lme对象。

您可以获取lme模型的AIC和logLik，例如：

AIC(mod1$lme)

我对gamm模型不太了解，无法说出这是否等于您所追求的（您可能想在交叉验证中重新发布）