为了包括多个因子,我们以长格式创建一个指标列,用于唯一标识项目所属的因子。...并且,因子间相关矩阵匹配来自多级的随机斜率相关。...----
例如,如果要对因子进行回归,可以使用相关矩阵的后验和solve()函数来得出回归中因子的系数。...在这里,我在因子2和3上回归因子1:
R <- extract(cfa.stan.fit, c("R[1, 2]", "R[1, 3]", "R[2, 3]"))
R <- cbind(R$`R[1,2...]`, R$`R[1,3]`, R$`R[2,3]`)
coefs <- matrix(NA, nrow(R), ncol(R) - 1)
for (i in 1:nrow(R)) {
m <- matrix