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

也许是个愚蠢的问题，但经过几个小时的寻找，我仍然没有找到解决办法。 我运行了一个gamm模型，该模型似乎运行良好。现在我想使用在模型中计算的对数似然来计算模型的AIC（mod1）

模型运行如下所示

mod1<-gamm(Biomass..g. ~ Atoll + SizeCm + Depth + OuterOrLagoon,random=list(SiteID=~1),family=gaussian,data=RMI2014, method = "REML")

起初我认为这是错误的，因为我使用的是REML而不是最大似然，但在不使用REML时返回相同的错误

是因为gamma中有一个随机因子，所以返回了错误吗

请帮忙

干杯，泰森

根据，调用

gamm

的结果是一个列表，其中包含两个类型为

gam

和

lme

的对象

您可以获得lme车型的AIC和logLik，例如：

AIC(mod1$lme)

我对gamm模型了解不够，无法判断这是否等同于您所追求的（您可能希望重新发布）

AIC(mod1$lme)