这似乎是一件简单的事情,但我无法让它发挥作用!
我有一个如下所示的循环:
PCPHENO<-merge(PCAEVEC,PHENOTYPE)
for (i in 1:5) {
summary(lm(Outcome ~ PCPHENO[,i+1], PCPHENO))
}
PCPHENO看起来像这样:
ID PC1 PC2 PC3 PC4 PC5 Outcome
Bob 0.1 0.4 -0.3 0.4 0.3 1
Carl -0.1 0.3 -0.2 0.2 0.2 1
Dan 0.6 0.2 0.3 -0.4 0.2 2
Ed 0.6 0.4 -0.3 0.4 0.3 2
Frank 0.2 0.4 -0.3 0.4 0.3 2
George 0.5 0.4 -0.3 0.4 0.3 1
当我单独运行命令 summary(...)时,它可以正常工作并输出数据。但是,当我运行这个循环时,看起来很好但是没有输出。
我确定我只是错过了一个命令,但我不知道它是什么!