search for: mvrs

...hermore, the overall probability to get either outcome of the binary variable should be 0.5. Below you can see the R code (I use for simplicity a diagonal matrix in rmvnorm even if it produces no correlated sample): "sampleCop" <- function(n = 1000, rho = 0.2) { require(splus2R) mvrs <- rmvnorm(n + 1, mean = rep(0, 3), cov = diag(3)) pmvrs <- pnorm(mvrs, 0, 1) var1 <- matrix(0, nrow = n + 1, ncol = 1) var1[1] <- qbinom(pmvrs[1, 1], 1, 0.5) if(var1[1] == 0) var1[nrow(mvrs)] <- -1 for(i in 1:(nrow(pmvrs) - 1)) { if(pmvrs[i + 1, 1] <= rho) var1[i + 1] &lt...

...nclude -DNDEBUG -I/usr/local/include -fpic -g -O2 -c me_ols.c -o me_ols.o gcc -std=gnu99 -I/ebio/abt6/jmao/R/R-2.15.0/include -DNDEBUG -I/usr/local/include -fpic -g -O2 -c mvr.c -o mvr.o gcc -std=gnu99 -I/ebio/abt6/jmao/R/R-2.15.0/include -DNDEBUG -I/usr/local/include -fpic -g -O2 -c mvrs.c -o mvrs.o gcc -std=gnu99 -I/ebio/abt6/jmao/R/R-2.15.0/include -DNDEBUG -I/usr/local/include -fpic -g -O2 -c newick.c -o newick.o gcc -std=gnu99 -I/ebio/abt6/jmao/R/R-2.15.0/include -DNDEBUG -I/usr/local/include -fpic -g -O2 -c nj.c -o nj.o gcc -std=gnu99 -I/ebio/abt6/jmao/R/R-2.15.0/inc...