作者swwf (十字路口)
看板Statistics
标题[程式] 在R里面用多组资料作线性回归
时间Wed Jan 9 13:37:19 2013
------------------------------------------------------------------------
[软体程式类别]:R
[程式问题]: 资料处理、回归
[软体熟悉度]:
低(1~3个月)
[问题叙述]:
是这样,我有多组数据需要跑线性回归,y~ ax+b
如果只是单纯跑一组的话我还会弄,但如果是很多分组分层的话就不知道怎麽写回圈
原始数据是X(时间),和Y(存活比率)作简单直线回归,
现在如果我数据有多项不同的处理,整理结果如下样式
处理A 处理B 处理C 重复 时间 存活比率
7个层级 2个层级 6个 4个 X1~Xn Y1~ Yn
要注意的是,每个组合的时间n是不一样的观察数量,有的组合观察数少(早早全死光),
有的高(等很久才死完)
现在想要跑这7*2*6*4个直线回归,把作回归得来的值给输出,
并把这些直线回归作图,请问我的回圈程式要怎麽写?
搞不清楚要用apply还是tapply...
请详尽叙述遭遇到的问题,可能的话,分点叙述你要处理的流程
[程式范例]:
虽然张贴程式很可怕,但基本上有些程式还是要张贴才能解决
for (i in treat1)
{for (j in treat2)
{for (k in treat3)
{for (l in rep)
{lm.reg=lm(x~y)
}
}
}
}
summary[i,j,k,l](lm.reg)
plot[i,j,k,l](x,y)
-----------------------------------------------------------------------------
不好意思,我是R和写程式的新手,真的不太知道要怎麽样跑这种回圈,另外想要一页
多图把结果的图都画一下,想请教各位大大指导我该如何写这回圈,谢谢。
打扰各位不好意思 >_<
--
※ 发信站: 批踢踢实业坊(ptt.cc)
◆ From: 86.166.252.36
※ 编辑: swwf 来自: 86.166.252.36 (01/09 13:38)