メモ

• var.equal=TRUEの影響を見る

ma <- 120
mb <- 120
sa <- 10
sb <- 10 + (0:(-9))
na <- 15
nb <- 10:15
sb.nb <- expand.grid(sb,nb)
n.iter <- 10000
ps <- matrix(0,length(sb.nb[,1]),n.iter)
for(i in 1:length(sb.nb[,1])){
	for(j in 1:n.iter){
		x <- rnorm(na,ma,sa)
		y <- rnorm(sb.nb[i,2],mb,sb.nb[i,1])
		ps[i,j] <- t.test(x,y,var.equal=TRUE)$p.value
	}
}
matplot((apply(ps,1,sort)),type="l")
image(matrix(apply(ps,1,quantile,0.5),length(sb),length(nb)))
#plot(apply(ps,1,quantile,0.5))