## 10000 variables, hence 10000 t-tests will be performed
set.seed(12345)
x1 <- rpois(10000, 5)
x2 <- rpois(10000, 5)
n1 <- rpois(10000, 40)
n2 <- rpois(10000, 40)
a <- proptests(x1, x2, n1, n2)
system.time( proptests(x1, x2, n1, n2) )
Run the code above in your browser using DataLab