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