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