中的R

Question

運行鍼對同一主可變倍數，簡單線性迴歸我有一個data.frame（讓我們稱之爲DF）加載至R基本上如下所示：

  primary_variable  var1  var2  var3  var4...  var354 
sample1  5     1  4  3  2  1 
sample2  8     2  3  4  1  2 
sample3  7     3  2  1  4  3 
sample4  2     4  1  2  3  4 
. 
. 
. 
sample58  8     1  2  3  4 

基本上，我想運行簡單線性迴歸多次，比較所述主變量和所有其它的單獨如下：

reg <- lm(primary_variable~var1, data=DF) 
reg <- lm(primary_variable~var2, data=DF) 
reg <- lm(primary_variable~var3, data=DF) 
reg <- lm(primary_variable~var354, data=DF) 

並具備數據輸出到某種形式的表，這樣我可以爲他們所有讀出的P值。

這樣做最簡單的方法是什麼？

Answer 1

這裏是一個辦法做到這一點：

DF <- read.table(text="   primary_variable  var1  var2  var3  var4 
sample1  5     1  4  3  2 
sample2  8     2  3  4  1 
sample3  7     3  2  1  4 
sample4  2     4  1  2  3", header=TRUE) 

sapply(DF[,-1], function(x) summary(lm(DF[,1]~x))$coef[,"Pr(>|t|)"]) 

#     var1  var2  var3  var4 
# (Intercept) 0.1471971 0.477767 0.4023857 0.206388 
# x   0.5120500 0.512050 0.7072300 0.707230 

不過，我強烈建議你不要這樣做。請向統計員尋求幫助。