1

当我将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并安装了相关软件包的最新版本。

4

1 回答 1

2

这是familyglmmTMB. 它现在在MuMInR-Forge 的 1.40.7 中修复:install.packages("MuMIn", repos="http://R-Forge.R-project.org")

于 2018-06-14T18:23:57.343 回答