Re: [程式] 用統計軟體R作ANOVA的分析
※ 引述《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
推
12/12 17:42, , 1F
12/12 17:42, 1F
推
12/12 17:52, , 2F
12/12 17:52, 2F
推
12/12 18:09, , 3F
12/12 18:09, 3F
→
12/12 19:01, , 4F
12/12 19:01, 4F
推
12/12 19:18, , 5F
12/12 19:18, 5F
推
12/12 19:55, , 6F
12/12 19:55, 6F
討論串 (同標題文章)