Как извлечь факторные нагрузки из лавы?

Question

Как я могу получить таблицу со всеми скрытыми факторами и загрузкой каждого элемента измерения по всем факторам? Я не могу найти способ вытащить это из подходящей модели лавы. Вот общий код, который я использую для подбора модели.

library(lavaan)
fit <- sem(mySemModel, data=df, std.ov=TRUE, std.lv=TRUE)
summary(fit, fit.measures=TRUE, rsq=TRUE, standardized=TRUE)

Я ищу тот же вывод, который вы получите от EFA. Например, если я запустил код:

library(psych)
myFA <- fa(tpblatentData, 2)
print(myFA)

Я бы получил что-то вроде этого:

               PA1   PA2
Qitem1              0.74
Qitem2              0.82
Qitem3              0.87
Qitem4        0.98      
Qitem5        0.94      
Qitem6        0.89      

Answer 1

Вы можете получить стандартизированные загрузки модели в матричной форме, используя inspect функция от lavaan пакет. Следующий код вернет матрицы лямбда (коэффициент загрузки), тета (ковариационная матрица наблюдаемой ошибки), psi (матрица скрытой ковариации) и бета (скрытые пути).

inspect(fit,what="std")

Из вашего примера видно, что вы ищете коэффициенты загрузки, которые находятся в лямбда-матрице:

inspect(fit,what="std")$lambda

Аналогичным образом вы можете извлечь нестандартные параметры, указав "est" вместо "std".

Answer 2

Я нашел ответ Джоэла полезным. Еще одна вещь, которая может помочь, - это приписывание результата нагрузок переменной.

Поскольку функция inspect() возвращает список, для меня это было полезно:

model_loadings <- inspect(model_fit, what = "std")[["lambda"]]

Теперь я могу использовать эти значения для вычисления других интересных вещей.