R / agricolae / cv 찾기
#install.packages("agricolae") library(agricolae) tr=rep(c("A","S","U"), each=3) rp=rep(c(1:3),3) d=c(3.3,2.9,3.1, 3.1,2.7,3.3, 25.0,20.1,23.0) eff=data.frame(tr,rp,d) eff$tr=as.factor(eff$tr) eff$rp=as.factor(eff$rp) fix(eff) eff.aov=aov(d~tr+rp, data=eff) summary(eff.aov) dmrt=duncan.test(eff.aov, "tr") #print(dmrt) print(dmrt[["statistics"]]["CV"]) --------------------------------------------..
2020. 4. 16.