如何获得lavaan cfa中每个项目的误差方差？

Question

我想从 lavaan CFA 中提取每个项目（即指标）的误差方差。

library(lavaan)

HS.model <- ' visual  =~ x1 + x2 + x3
              textual =~ x4 + x5 + x6
              speed   =~ x7 + x8 + x9 '

fit <- cfa(HS.model, data = HolzingerSwineford1939)
summary(fit, fit.measures = TRUE)

我知道我可以用 inspect(fit, what = "est")$lambda 提取 lambda。如何获得每个项目的误差方差？

Answer 1

lavInspect(fit, what = "est")$theta

这是残差的整个协方差矩阵。对角线有残差。