[程式] 在R裡面用多組資料作線性迴歸

看板Statistics作者 (十字路口)時間13年前 (2013/01/09 13:37), 編輯推噓0(000)
留言0則, 0人參與, 最新討論串1/2 (看更多)
------------------------------------------------------------------------ [軟體程式類別]: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)
文章代碼(AID): #1GxG8IVE (Statistics)
文章代碼(AID): #1GxG8IVE (Statistics)