alpha1 <- .95 alpha2 <- .99 var1 <- -qnorm(1-alpha1, 20, 25) var2 <- -qnorm(1-alpha2, 20, 25) var12 <- -25*qnorm(1-alpha1)-20 var22 <- -25*qnorm(1-alpha2)-20