d <- read.table("rice112var.txt", header=T) 				#データの読み込み
Pd <- data.frame(Treat=factor(d$treat),DryWeight=d$data)		#フレームワーク作成


bymedian <- with(Pd, reorder(Treat, DryWeight, median))			#中央値の昇順に並べる場合
plot(DryWeight ~ bymedian,data=Pd, ylab="DryWeight(mg)", xlab="",btl="l",cex.axis=0.7,las=2) #箱ヒゲ図の描出


res1 <- aov(DryWeight~Treat,data=Pd)	#一元分散分析
summary(res1)				#分散分析表の出力


#有意差のアルファベット表示
library(agricolae)
res2 <- HSD.test(res1, "Treat", group=TRUE)
res2