R code :
library(MASS)
mu <- rep(0,4)
Sigma <- matrix(.7, nrow=4, ncol=4) + diag(4)*.3
rawvars <- mvrnorm(n=10000, mu=mu, Sigma=Sigma)
pvars <- pnorm(rawvars)
binomvars <- qpois(1-pvars, 3, .25)
I have not understood how does qpois(1-pvars, 3, .25) work ? What is the...