→ 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