作者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