как программировать итерации факториала

Question

В настоящее время я очень занят новым аналитическим фактором. Я создал следующие бинарные операторы:

1. % M% вычисляет квадрат положительно определенную симметричную матрицу из R dataframe

2. % E% вычисляет оценку соборности и размещает их на главной диагонали квадратной матрицы М

3. % F% извлекает один фактор за один раз из новой квадратной матрицы

4. % R% воспроизводят новую квадратную матрицу из факторных нагрузок% F% и вычитает I t от исходной матрицы

Это один цикл. Шаги 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 

а = является R dataframe, 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 
}