#การเปรียบเทียบค่าเฉลี่ย
เริ่มต้นด้วยการเรียก library(agricolae) แล้วสั่งด้วยคำสั่งข้างล่างนี้ได้เลยครับlibrary(agricolae)
#LSD
LSD.test(yield,trt,21,2468,group=T,console=T)LSD.test(md,'trt',alpha=0.05,console=T)
tapply(yield,trt,sd)
tapply(yield,trt,mean)
tapply(yield,trt,var)
tapply(yield,trt,sd)/sqrt(4)
#duncan
duncan.test(yield,trt,21,2468,group=T,console=T)duncan.test(md,'trt',alpha=0.05,console=T)
#qtukey
qtukey(0.05,7,21,lower.tail=F)qtukey(0.95,7,21)
#scheffe
scheffe.test(yield,trt,21,2468.4,group=T,console=T)
scheffe.test(md,'trt',alpha=0.05,console=T)
#การเปรียบเทีบบแบบออธาโกนอล
d<-read.table('clipboard',header=TRUE)d
attach(d)
names (d)
blk<-factor(blk)
trt<-factor(trt)
coef<-matrix(c(-3,1,1,1,0,1,1,-2,0,-1,1,0),nrow=4,ncol=3,byrow=F)
coef
contrasts(trt)<-coef
contrasts(trt)
md<-aov(yield~blk+trt)
summary(md)
summary(md,split=list(trt=list(1,2,3)))
summary(md,split=list(trt=list('T1vs(T2,T3,T4)'=1,'(T2,T3)vsT4'=2,'T2vsT3'=3)))
#การเปรียบเทียบแนวโน้ม
d<-read.table('clipboard',header=TRUE)
d
attach(d)
names (d)
blk<-factor(blk)
trt<-factor(trt)
coef<-matrix(c(-3,-1,1,3,1,-1,-1,1,-1,3,-3,1),nrow=4,ncol=3,byrow=F)
coef
contrasts(trt)<-coef
contrasts(trt)
md<-aov(yield~blk+trt)
summary(md)
summary(md,split=list(trt=list('Li'=1,'Qu'=2,'Cu'=3)))
#หรือ
summary(md,split=list(trt=list('Linear'=1,'Quadratic'=2,'Cubic'=3)))
y<-tapply(yield,trt,sum)
y
x<-c(5,10,15,20)
x
md<-lm(y~x+I(x^2))
summary(md)
detach(d)
#การเปรียบเทียบแนวโน้ม แบบลัดไม่ต้องเปิดตาราง
d<-read.table('clipboard',header=TRUE)d
attach(d)
names (d)
blk<-factor(blk)
trt<-factor(trt)
coef<-contr.poly(4)
coef
contrasts(trt)<-coef
trt
md<-aov(yield~blk+trt)
summary(md)
summary(md,split=list(trt=list('Li'=1,'Qu'=2,'Cu'=3)))
#การเปรียบเทียบระหว่างกลุ่มร่วมกับการเปรียบเทียบแนวโน้ม
d<-read.table('clipboard',header=TRUE)d
attach(d)
names (d)
blk<-factor(blk)
trt<-factor(trt)
d1<-read.table('clipboard',header=TRUE)
d1
coef<-as.matrix(d1)
coef
is.matrix(coef)
md<-aov(yield~blk+trt)
summary(md)
detach(d)
rm(list=ls(all=TRUE))
