作者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