Question

我有一个使用glmmTMB构建的glmm模型：

M2 <- glmmTMB(TotalSG_18_beta ~ Sim_mean_CI + (1|Site) + (1|Transect), data = CI_simulations_18_beta_transformed, family=list(family="beta",link="logit"))

我已经转换了数据，并且正在使用Beta分布式模型，因为我的数据不正常。我想使用lmer软件包中的bootMer来使用引导程序验证该模型的适当性。我想将数据集拆分为80:20，因此请使用80％的值来预测20％，然后与模型输出进行比较。但是我似乎无法从自举模型中提取置信区间。

我尝试了以下代码：

b_test_1 <- bootMer(M2,
                      FUN = function(x) confint(x),
                      nsim=10, 
                      use.u = TRUE, 
                      type = "parametric", 
                      .progress="txt")

我看到以下输出：

Bootstrap Statistics :
             original       bias    std. error
t1*     -0.8895147064 -0.493229558   0.7950155
t2*     -0.8897969448 -0.493181702   0.7950244
t3*     -0.8953021883 -0.492248240   0.7951980
t4*     -0.8919079383 -0.492823765   0.7950909
t5*     -0.8954317661 -0.492226269   0.7952021

以下是摘要的输出：

summary(b_test_1)
          Length Class      Mode     
t0         15    -none-     numeric  
t         120    -none-     numeric  
R           1    -none-     numeric  
data        4    data.frame list     
seed      626    -none-     numeric  
statistic   1    -none-     function 
sim         1    -none-     character
call        7    -none-     call     
ran.gen     1    -none-     character
mle         2    -none-     list

我真的不明白输出。我期望看到我的上下CI。我想我在FUN = function（x）的行上出错了，但是我似乎无法找到我应该放在这里的东西。当我尝试以下功能时：

boot.ci(b_test_1, level = 0.95)

我收到以下错误消息：

Error in boot.ci(b_test_1, length = 0.95) : 't' must of length 10

有人可以帮助我理解此输出的含义以及为什么我收到此错误吗？我是R和编码方面的新手，因此非常感谢任何帮助或建议

我可以在lmer中使用bootMer获得置信区间吗？

0 个答案: