作者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/m.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