从 rptR 返回置信区间结果

Question

我正在使用 rptR 运行一个用于可重复性分析的 for 循环。我想将每次重复性分析运行的所有置信区间结果存储在数据框中。我在循环中使用以下代码返回 CI 值列表：

for(i in 1:P){
  m1<-rptR::rpt(Behaviour~Temperature+(1|ID), grname="ID", data=Data, datatype="Gaussian", nboot=10, 
  npermut=10)
  R_value[i] <- m1$R
  CI_value[i] <- m1$CI
}

但是，CI 值始终返回为 0.95 - 说明它是 95% 的置信区间，而不是为我提供 CI 值的实际范围。有谁知道如何返回实际的 CI 值？

这是仅运行一次 rptR 的模型输出：

Repeatability for ID
R  = 0.847
SE = 0.134
CI = [0.484, 0.982]
P  = 0.00152 [LRT]
     0.015 [Permutation]

Answer 1

我认为您希望提取 CI_emp 而不是 CI 并且因为将其存储在列表中会有多个值。试试看：

R_value <- numeric(P)
CI_value <- vector('list', P)

for(i in 1:P){
  m1<-rptR::rpt(Behaviour~Temperature+(1|ID), grname="ID", 
                data=Data, datatype="Gaussian", nboot=10, npermut=10)
  R_value[i] <- m1$R
  CI_value[[i]] <- m1$CI_emp
}