当我将model.sel
函数应用于g0
下面的模型时,它运行良好。但是,当我使用我的数据创建类似的模型时,我收到如下所示的错误。为什么我可能会收到此错误的任何线索?
data(sleepstudy,package="lme4")
g0 <- glmmTMB(Reaction~Days+(Days|Subject),sleepstudy)
model.sel(g0)
上面运行良好,但是,这不是:
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
并安装了相关软件包的最新版本。