作者arbiteron (..)
站內Statistics
標題Re: [程式] 用統計軟體R作ANOVA的分析
時間Sat Dec 12 17:31:44 2009
※ 引述《QQVila (雷曼×)》之銘言：
: 首先謝謝大家的建議:) 其實我也去google找了很久…
: 廢話不說請先看看哪裡有錯誤吧！ (詳情data請看上一篇！)
: > data <-read.table("r.txt",header=T)
: > data
: part r test tester
: 1 37 1 1 1
(略)
: 90 35 10 3 3
: > r.inter <- aov(part~test*tester,data=data)
: > summary(r.inter)
: Df Sum Sq Mean Sq F value Pr(>F)
: test 1 0.01667 0.01667 0.0004 0.9850
: tester 1 19.3 19.3 0.4109 0.5232
: test:tester 1 3.0 3.0 0.0645 0.8001
: Residuals 86 4032.1 46.9
: > anova(lm(part~test*tester,data=data))
: Analysis of Variance Table
: Response: part
: Df Sum Sq Mean Sq F value Pr(>F)
: test 1 0.01667 0.01667 0.0004 0.9850
: tester 1 19.3 19.3 0.4109 0.5232
: test:tester 1 3.0 3.0 0.0645 0.8001
: Residuals 86 4032.1 46.9
: 其實大家看就知道怪怪的…雖然兩張表寫出來的是一樣沒錯
: 但是Df的地方一看就知道是錯的，想當然爾後面的也是錯的！！
: 請問是哪裡出問題呢？
: 因為這題其實是課本上的，當然有一張他的ANOVA的表格。
: 感謝～
注意看aov裡頭的結構是類似回歸的語法
但是這個時候R將test及tester這兩個變數當成"數字"
而非"factor"
所以在執行aov之前需再加上
test <- factor(test)
tester <- factor(tester)
再執行即可
--
※ 發信站: 批踢踢實業坊(ptt.cc)
◆ From: 125.231.209.83
1F：推 QQVila:他說..找不到test這個物件..請問要怎麼解決= =" 12/12 17:42
2F：推 QQVila:我已經解決了@@!非常感謝喔~~有問題在請教各位~~^^ 12/12 17:52
3F：推 QQVila:抱歉..剛剛試過了還是不行說..囧? 12/12 18:09
4F：→ arbiteron:要先加上 attach(data) 12/12 19:01
5F：推 QQVila:抱歉剛剛離開了一下..還是不行說!?還是跟上面跑出來的一樣 12/12 19:18
6F：推 QQVila:成功了..謝謝!手邊其實沒有課本可以對答案,有誤在麻煩您^^" 12/12 19:55