作者clever1510 (clever)
看板Statistics
標題[程式] R迴圈--簡單線性迴歸
時間Wed Mar 28 13:21:31 2012
[軟體程式類別]:
R
[程式問題]:
迴圈
[軟體熟悉度]:
請把以下不需要的部份刪除
新手(不到1個月)
[問題敘述]:
我想讓簡單線性迴歸重複執行很多次
將T-value 或 F-value 輸出成表格
例如我有一個10個column的表
第一個column 是Y 要"分別"對,其餘column 迴歸並檢測
我想的流程大概是這樣
C1 與 C2 做迴歸 --->T-value 1
C1 與 C3 做迴歸 --->T-value 2
.
.
.
C1 與 C9做迴歸 --->T-value 9
輸出 T-value1~9
可是不知道怎麼寫迴圈
[程式範例]:
這是目前我編輯好的部分
#open file
directory<-"C://";
filename<-"CP7_Ttest.txt";
#read file
a<-read.table(paste(directory, filename, sep=""), header=T, nrows=500, skip=0);
print(names(a));
Y<-a[,c(1)]; X<-a[,c(2:10)]; rm(a);
#?lm()
i=1
for(i in 1:9)
{
Xa = X[,i]
model = lm(formula = Y ~ Xa)
model
summary (model)
i=i+1
}
我參考R 大的意見修改了
可是對於結果輸出仍有疑問
不知道怎樣將不同次的結果輸出成一個文件
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 140.112.133.187
1F：→ rockken:用 for + if + paste 可以寫出來 03/28 15:09
※ 編輯: clever1510 來自: 140.112.133.187 (03/28 15:20)
※ 編輯: clever1510 來自: 140.112.133.187 (03/28 15:21)
※ 編輯: clever1510 來自: 140.112.133.187 (03/28 15:26)
※ 編輯: clever1510 來自: 140.112.133.187 (03/28 15:33)
※ 編輯: clever1510 來自: 140.112.133.187 (03/28 15:34)
2F：→ rockken:names(summary (model)) 看你要哪些 在宣告一矩陣 03/28 21:43
3F：→ rockken:每做完一次lm 存一次，最後在ptint 你存好矩陣即可 03/28 21:43