Пример
OUTCOME=0:9
X1=1:10
X2=5:14
Z1=9:18
Z2=21:30
Z3=41:50
Z4=52:61
Z5=20:29
DATA=data.frame(OUTCOME,X1,X2,Z1,Z2,Z3,Z4,Z5)
ADJUSTMENTS=c("X1","X2")
Z=c(Z1,Z2,Z3,Z4,Z5)
for (i in Z){
coef(lm(OUTCOME~ADJUSTMENTS + Z));
lm(OUTCOME~ADJUSTMENTS + Z)$r.squared;
}
У меня есть данные, и я надеюсь сделать это: DV не изменяется, и все COVS включают X1, НО я хочу зациклить для VARS Z1, Z2, Z3, Z4, Z5, Z6, Z7, Z8 и так далее. Затем я собираюсь сохранить этот коэффициент только для Z var, чтобы у меня был новый фрейм данных:
Z COEF RSQUA
Z1 VALUE VALUE
Z2 VALUE VALUE
Z3 VALUE VALUE
Z4 VALUE VALUE
Z5 VALUE VALUE
Z6 VALUE VALUE
Z7 VALUE VALUE
Z8 VALUE VALUE