(2) 帰無仮説からの２ｘ２表ランダムサンプリング

N1<-100
N2<-150
p1<-p2<-0.3
Niter<-10000
p.out<-rep(0,Niter)
v1<-c(rep(1,N1),rep(2,N2))
for(i in 1:Niter){
	v2<-sample(c(1,2),N1+N2,replace=TRUE,prob=c(p1,1-p1))
	#print(table(v1,v2))
	#p.out[i]<-chisq.test(c(rep(1,N1),rep(2,N2)),c(v1,v2))$p.value
	p.out[i]<-chisq.test(table(v1,v2),correct=FALSE)$p.value
}
plot(ppoints(Niter,a=0),sort(p.out))