作者ndd2 (ndd2)
看板Statistics
标题[程式] Minitab、R 做2-way ANOVA分析结果？
时间Sun Sep 24 22:34:07 2017
[软体 ]
Minitab 及 R
[程式问题]:
用上述两个软体对自己设想的资料集做 ANOVA分析结果不一致，不知道问题在哪？
[软体熟悉度]:
新手
[问题叙述]:
资料集如下共12 observations,
其中money当response变数，team、sex当factors
ID money team sex
id1 123 A m
id2 432 B m
id3 122 C m
id4 134 D m
id5 45 A m
id6 34 B m
id7 626 C f
id8 266 D f
id9 664 A f
id10 345 B f
id11 452 C f
id12 453 D f
[程式范例]:
(1)用Minitab 18操作如下：
Stat>ANOVA>General Linear Model>Fit General Linear Model
结果得：
Analysis of Variance
Source DF Adj SS Adj MS F-Value P-Value
team 3 23887 7962 0.28 0.840
sex 1 296815 296815 10.34 0.032
team*sex 3 85640 28547 0.99 0.481
Error 4 114866 28717 Total 11 531472
(2)用R操作如下：
> nameslist<-read.table("name-list.csv",header=T,sep=",")
> lmod<-lm(nameslist$money~nameslist$team*nameslist$sex)
> anova(lmod)
结果得：
Analysis of Variance Table
Response: nameslist$money
Df Sum Sq Mean Sq F value Pr(>F)
nameslist$team 3 34150 11383 0.3964 0.76344
nameslist$sex 1 296815 296815 10.3360 0.03244 *
nameslist$team:nameslist$sex 3 85640 28547 0.9941 0.48099
Residuals 4 114867 28717
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
不知道为什麽结果不一致？
谢谢。
(另外我用书上的范例数子去跑程式，Minitab和R的结果竟然就一样!
，为何我自己的假想资料会出问题？
书上的范列数字如下：
MAJOR DEGREE OFFER
Acct Bach 26
Acct Bach 24.8
Acct Bach 29.7
Acct Bach 33.2
Acct Bach 31.8
Acct Bach 26.9
Acct Bach 30.9
Acct Bach 23.3
Acct Bach 29.1
Acct Bach 26.4
Acct Master 27.4
Acct Master 29.9
Acct Master 33.3
Acct Master 29.7
Acct Master 31.5
Acct Master 37.3
Acct Master 33.1
Acct Master 31.4
Acct Master 28.8
Acct Master 30.6
Mrkt Bach 25.5
Mrkt Bach 24.5
Mrkt Bach 24.2
Mrkt Bach 24.6
Mrkt Bach 21.8
Mrkt Bach 19.9
Mrkt Bach 23
Mrkt Bach 32
Mrkt Bach 23.8
Mrkt Bach 29.7
Mrkt Master 42.2
Mrkt Master 44.5
Mrkt Master 45
Mrkt Master 38.8
Mrkt Master 38.4
Mrkt Master 35.9
Mrkt Master 37.8
Mrkt Master 41.8
Mrkt Master 42.7
Mrkt Master 37.9
--
※ 发信站: 批踢踢实业坊(ptt.cc), 来自: 123.193.54.130
※ 文章网址: https://webptt.com/cn.aspx?n=bbs/Statistics/M.1506263650.A.19D.html
※ 编辑: ndd2 (123.193.54.130), 09/25/2017 00:24:04
1F：→ andrew43: 这组R code给你的是TYPE I SS，不是TYPE III 09/25 00:25
3F：→ andrew43: 重点在car::Anova(..., type = 3)及contrasts 09/25 00:28
4F：→ recorriendo: 每个cell有一样多观察个数十 结果才会完全相同 09/25 00:58