2023 R배우기(반복없는 이원 분산분석1)_1031

[storyteller]

2023 R배우기(반복없는 이원 분산분석1)_1031

특성치문제(난괴법)이원배치_1101.xlsx

9.57KB

---

2023 R배우기(반복없는 이원 분산분석1)_1031.pdf

571.10KB

---

< 반복없는 이원분산 1101>

trt<-c(1,1,1,1, 2,2,2,2, 3,3,3,3 )
block<-c(1,2,3,4, 1,2,3,4, 1,2,3,4)
y<-c(20, 22, 18, 25,
16, 18, 17, 19,
30, 34, 29, 27)
trt<-as.factor(trt)
block<-as.factor(block)

mean(y)
tapply(y,trt,mean)
tapply(y,trt,var)
tapply(y,block,mean)
tapply(y,block,var)
tapply(y,trt,mean)-mean(y)
tapply(y,block,mean)-mean(y)

aov.out<-aov(y~trt+ block)
summary(aov.out)
aov.out$fitted
aov.out$residuals
shapiro.test(aov.out$residuals)

par(mfrow=c(1,2))
plot(y~trt)
plot(y~block)


pairwise.t.test(y,trt,p.adjust="none",pool.sd=TRUE)
pairwise.t.test(y,trt,p.adjust="none")
pairwise.t.test(y,trt,p.adjust="bonferroni",pool.sd=TRUE)
a.tukey<-TukeyHSD(aov.out, ordered=TRUE)
a.tukey


par(mfrow=c(2,2))
plot(aov.out)
 plot(a.tukey)