如何编写因子分析迭代程序

Question

目前，我正忙于一种新的因素分析。我创建了以下二进制运算符：

1. %M%从R数据中计算正定对称矩阵
2. %E%计算公共度估计，并将这些估计放在方阵M的主对角线上。
3. %F%一次从新的平方矩阵中提取一个因子。
4. %R%从%F%的因子负载中复制一个新的平方矩阵，并从原始矩阵中减去它。

这是一个周期。步骤2到4对残差重复，%F%提取第二个因子，%R%再现新的矩阵等。在这种方式下，完整的过程(步骤2-4，2-4)依次应用于越来越小的矩阵。

我的问题是:我如何在R中编程这些迭代？显然，在每个迭代中，步骤4是下一个迭代中从步骤2到步骤4的新周期的输入。

我用的一些R码：

迭代法

res=  a%M%lg3
est=  res%E%5
F1=   est%F%1

res=  F1%R%est
est=  res%E%5
F2=   est%F%1
res=   F2%R%est

est=  res%E%5
F3=   est%F%1
res=  F3%R%est

a=是一个R数据，lg3是新的二元log3共识系数。

我欢迎有此机会寻求帮助。

致以敬意，

Luc

Answer 1

您可以使用一个while循环，我认为：

while(..some stopping criteria..) {
  est <- res %E% 5
  res <- est %F% 1 %R% est
}

如果您只想运行四次，只需添加一个计数器即可。

for(i in 1:4) {
  est <- res %E% 5
  res <- est %F% 1 %R% est
}