dat <- read.delim("http://phi.med.gunma-u.ac.jp/msb/data/p06.txt") attach(dat) layout(matrix(c(1,3,2,4),2,2)) matplot(rbind(rep(1,10),rep(2,10)),rbind(NDS,NDE),type="l", ylim=c(20,35),ylab="‘Ìd(g)",main="•’ʉaŒQ",xaxt="n",xlab="") axis(1,1:2,c("ŠJŽnŽž","I—¹Žž")) matplot(rbind(rep(1,10),rep(2,10)),rbind(HFDS,HFDE),type="l", ylim=c(20,35),ylab="‘Ìd(g)",main="‚Ž‰–b‰aŒQ",xaxt="n",xlab="") axis(1,1:2,c("ŠJŽnŽž","I—¹Žž")) NDD <- NDE-NDS HFDD <- HFDE-HFDS WeightGain <- c(NDD,HFDD) Diet <- as.factor(c(rep("•’ʉa",10),rep("‚Ž‰–b‰a",10))) IX <- c(1.1,2.1) MWG <- tapply(WeightGain,Diet,mean) SDWG <- tapply(WeightGain,Diet,sd) stripchart(WeightGain~Diet,method="jitter",vert=T,ylim=c(0,3), ylab="‚PTŠÔ‚Ì‘Ìd‘‰Á(g)") points(IX,MWG,pch=18) arrows(IX,MWG-SDWG,IX,MWG+SDWG,angle=90,code=3) boxplot(WeightGain~Diet,ylim=c(0,3),ylab="‚PTŠÔ‚Ì‘Ìd‘‰Á(g)") print(res<-var.test(NDD,HFDD)) VAREQ <- ifelse(res$p.value < 0.05, FALSE, TRUE) t.test(NDD,HFDD,var.equal=VAREQ) detach(dat)