Question

当我将model.sel函数应用于下面的g0模型时，它运行正常。但是，当我使用我的数据创建类似的模型时，我得到如下所示的错误。为什么我会收到这个错误的任何线索？

data(sleepstudy,package="lme4")
g0 <- glmmTMB(Reaction~Days+(Days|Subject),sleepstudy)
model.sel(g0)

上面运行正常，但是，这不是：

数据位置：https://mqoutlook-my.sharepoint.com/:f:/g/personal/joseph_mbui_mq_edu_au/EjFIajIJxxxBsIBcMi3GyIMBc6FyN4a5y-39cDev2Aoyng?e=tUTj2W

datT<-read.csv('datT.csv')
myModel<-glmmTMB(y~x1+x2+ (1|randomEffect),list(family="beta",link="logit"), data=datT)
model.sel(myModel)
Error in vapply(unique(f), function(x) if (is.na(x)) NA_character_ else formals(get(x))$link,  : 
values must be length 1,
but FUN(X[[1]]) result is length 0

两种配置之间的唯一区别是我适合beta家庭

我尝试切换na.options并安装了相关软件包的最新版本。

Answer 1

这是由family的异常glmmTMB表述引起的。现在已在R-Forge的MuMIn 1.40.7中修复：install.packages("MuMIn", repos="http://R-Forge.R-project.org")

使用glmmTMB模型的MuMIn model.sel错误

1 个答案: