(3) 対立仮説からの２ｘ２表ランダムサンプリング

N1<-100
N2<-150
p1<-0.3
p2<-0.4
Niter<-10000
p.out<-rep(0,Niter)
v1<-c(rep(1,N1),rep(2,N2))
for(i in 1:Niter){
	v2<-c(sample(c(1,2),N1,replace=TRUE,prob=c(p1,1-p1)),sample(c(1,2),N2,replace=TRUE,prob=c(p2,1-p2)))
	#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))