我正在尝试为广义添加剂混合模型进行模型选择(使用gamm4使用R中的MuMIn包进行模型选择 我基本上试图通过MuMIn和gamm4跟踪this文献中的模型选择。
我正在创建一个包含9个变量的模型,以及一个随机的单独效果 - 看起来像这样:
library(gamm4)
library(MuMIn)
SouthFull = gamm4(OtoWidth ~ s(Ages) +
LagFinFSLdat_AnnEsts +
FSL_months_Feb +
LagFSLdat_Annual +
LagFSLdat_Spring +
LagFSL_months_Oct+
FinFSLdat_SummerDat +
FinFSLdat_AutumnDat +
LagFSL_months_Nov , random = ~(1|FishName), data = South)
但是当我尝试使用dredge
函数时,它失败并显示错误消息:
(dd <- dredge(global.model=SouthFull))
Error in dredge(global.model = SouthFull) : result is empty
In addition: There were 50 or more warnings (use warnings() to see the first 50)
> warnings()
Warning messages:
1: In gamm4::gamm4(...) : family not recognized (model 1 skipped)
...
50: In gamm4::gamm4(...) : family not recognized (model 50 skipped)
这是由于MuMIn无法识别gamm4模型系列的问题吗?
完整代码&amp;数据可以在GitHub repository找到,具体问题在代码gamm_analysis.R中,第81行有特定警告。
提前谢谢
答案 0 :(得分:1)
我认为包装有问题。我有同样的错误。然后,我尝试运行在MumIn包的函数dredge的帮助下提供的示例,以查看我在做什么不同,即使使用以下示例,我也会遇到相同的错误:
fm1 <- lm(y ~ ., data = Cement)
dd <- dredge(fm1)
错误消息是:
Error in nobs.default(global.model) : no 'nobs' method is available