[程式] 在R裡面用多組資料作線性迴歸
------------------------------------------------------------------------
[軟體程式類別]: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)
討論串 (同標題文章)
以下文章回應了本文:
完整討論串 (本文為第 1 之 2 篇):