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

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 和编码的新手，所以非常感谢任何帮助或建议

标签： rlme4mixed-modelsconfidence-intervalstatistics-bootstrap