0

我正在尝试使用 dredge() 函数通过完成变量的每个组合（每个模型最多五个变量）并使用针对小样本量 (AICc) 校正的 AIC 比较模型来评估模型。

但是，我收到一个错误和两个警告消息，如下所示：

固定术语是“（拦截）”警告消息：1：在挖泥机中（MaxN.model，m.min = 2，m.max = 5）：比较由 REML 拟合的模型 2：在挖泥机中（MaxN.model，m.min = 2, m.max = 5) : 参数 'm.min' 和 'm.max' 已被弃用，使用 'm.lim' 代替

我尝试按照指定更改为“m.lim”，但出现错误：

疏浚错误（MaxN.model，m.lim = 5）：无效的“m.lim”值另外：警告消息：在疏浚中（MaxN.model，m.lim = 5）：比较由 REML 拟合的模型

我正在使用的代码是：

MaxN.model<-lme(T_MaxN~Seagrass.cover+composition.pca1+composition.pca2+Sg.Richness+traits.pca1+
              land.use.pc1+land.use.pc2+seascape.pc2+D.landing.site+T_Depth, 
                random=~1|site, data = sgdf, na.action = na.fail, method = "REML")
Dd_MaxN<-dredge(MaxN.model, m.min = 2 , m.max = 5)

我究竟做错了什么？

4

1 回答 1

0
  1. 你没有告诉我们你试图指定什么m.lim。?dredge说：

m.lim ...（可选）单个模型中项数的限制'c（lower，upper）'

所以你应该指定一个二元素数值（整数）向量。

  1. 您绝对应该使用method="ML"而不是method="REML". REML 的警告/错误非常严重；比较通过 REML 拟合的具有不同固定效应的模型将导致胡说八道。

所以你应该尝试：

MaxN.model <- lme(..., method = "ML")  ## where ... is the rest of your fit
Dd_MaxN <- dredge(MaxN.model, m.lim=c(2,5))
于 2019-09-30T21:34:23.567 回答