s <- read.table("OXs.txt", header=T)			#ファイルの読み込み
Sd <- data.frame(LINES=factor(s$treat),INTENSITY=s$data) #フレームワーク作成
plot(INTENSITY ~LINES,data=Sd)　				#箱ヒゲ図
res1 <- aov(INTENSITY ~ LINES,data=Sd)			#一元分散分析
library(multcomp)					#パッケージ起動
res2 <- glht(res1, linfct = mcp(LINES = "Dunnett"))	#ダネット検定
summary(res2)						#結果を出力




#片側検定 コントロール（野生型）よりも形質値が高くなるはずの系統を検定する場合
res2 <- glht(res1, linfct = mcp(LINES = "Dunnett"), alternative = "greater")
confint(res2, level=0.95)
summary(res2)


#片側検定 コントロール（野生型）よりも形質値が低くなるはずの系統を検定する場合
res2 <- glht(res1, linfct = mcp(LINES = "Dunnett"), alternative = "less")
confint(res2, level=0.95)
summary(res2)