y <- NULL
for(i in 1:1000) {
x <- matrix(round(runif(4)*1000),2,2)
a <- x[1,1]; b <- x[1,2]
c <- x[2,1]; d <- x[2,2]
COR <- (a*d - b*c) / sqrt(
sum(x[1,]) * sum(x[2,]) *
sum(x[,1])* sum(x[,2]) )
or <- c(oddsratio(x, log=F))
YULE <- (or-1)/(1+or)
# P <- round(100*x[1,1]/sum(x[1,]),1)
y <- rbind(y, c(COR,YULE,min(x)) )
}
plot(y[,1],y[,2],type="n",ylim=c(0,1),xlim=c(0,1))
text(y[,1],y[,2],y[,3],cex=.6)
