作者sseeaann (屌哥)
看板R_Language
標題[問題] 詢問有關迴歸分析 R的指令
時間Thu Aug 15 20:08:36 2013
我知道ANOVA這個指令可以叫出變異數分析表
不過這指令跑出來的是一個MODEL的各變數的變異數分析表
而小弟我要的則是一個MODEL的SST,SSR,SSE的變異數分析表來分析整個MODEL
知道指令的大大麻煩請解答一下吧
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 175.98.79.50
1F：→ celestialgod:summary(aov(lm(y~X))) 08/15 20:35
2F：→ diplazium:樓上的結果跟anova(lm(y~x)是一樣的，兩種都可以 08/15 21:59
3F：→ diplazium:不過我不懂原PO的提問,為何他用anova會得不到SSR,SSE呢? 08/15 22:00
4F：→ sseeaann:不知道ㄟ 就只有跑出各個變數的sum sqs、mean sq等等 08/15 23:38
5F：→ sseeaann:一樓的結果跟ANOVA一樣喔 08/15 23:40
6F：推 celestialgod:我看到的是不同的 明天把程式補上 08/16 00:54
7F：→ diplazium:好像都只能跑出各變數的SS,原PO不如自己做加總->SSR 08/16 12:35
8F：→ sseeaann:我想程式應該沒問題 只是普通的一次項回歸式 08/16 18:47
9F：→ sseeaann:7樓大可請教一下作加總得SSR? 直接加總還是有指令? 08/16 18:48
10F：→ celestialgod:我可能昨天眼殘，再跑一次真的一樣 08/16 19:24
11F：→ celestialgod:七樓大大說的應該是自己寫加總程式 08/16 19:25
12F：→ celestialgod:SS = anova(lm.out)$'Sum Sq'; SST = sum(SS) 08/16 19:27
13F：→ celestialgod:SSE = sum(SS[length(SS)]); SSR = SST - SSE 08/16 19:28
14F：→ sseeaann:感謝一、七樓大的幫忙^^ 08/16 22:29
15F：→ sseeaann:不過還是希望能知道如何跑出整體MODEL的ANOVA報表 08/16 22:30
16F：推 lin15:fit<-lm(y~x);fit1<-lm(y~1);anova(fit1,fit) 08/16 23:02