#test.csv読み込み
test=read.csv("test.csv")
test
summary(test)

#様々なグラフ作成
boxplot(test)
boxplot(test,col="white")
install.packages("beeswarm")
library(beeswarm)
beeswarm(test)
beeswarm(test,method="hex",pch=20)
boxplot(test)
beeswarm(test,method="hex",pch=20,add=T)

#二元配置分散分析
twoway=read.csv("twoway.csv")
twoway
str(twoway)
twoway$dose=factor(twoway$dose,levels=c ("Control","1","10","100"))
twoway$add =factor(twoway$add,levels=c("none","A","B"))
str(twoway)
res=aov(live~dose+add+dose*add, data=twoway)
summary(res)