Omega (hierarchichal)心理包置信区间

Question

我想知道是否有一种方法可以使用R中的心理包来计算omega (分层)的可靠性置信区间？

我知道在MBESS中使用ci.reliability可以做到这一点，但对于我正在测试的模型来说，心理更合适。谢谢。

Answer 1

令我惊讶的是，是的，您可以通过使用bootstrap来使用psych找到omega_h的置信区间。只需设置n.iter=some编号即可。如果关闭解决方案的打印，速度会更快。

test10 <- omega(HolzingerSwineford1939[,7:15],n.iter=10,plot=FALSE)

test10
 ...
Estimates and bootstrapped confidence intervals
          lower estimate upper
omega_h   0.086     0.45  0.64
alpha     0.567     0.76  0.80
omega_tot 0.789     0.85  0.87
G6        0.776     0.81  0.84
omega_lim 0.123     0.53  0.82`