This program computes probability from a multivariate standard normal probability distribution given the z values and the correlations for up to 5 variables. The user can increase the number of variables by a slight modification of the program. A graphical representation of a bivariate standard normal distribution with no dependency between the two variables is shown in Figure 1.
In order for the program to work properly, the user need to place the input value in the same cell references as in Figure 2 or 3. Cell E14 of Figure 2 shows the probability from a bivariate standard normal distribution given z1 = 1 and z2 = 0 with correlation of 0.4. The probability in this case is 45.79% based on 100,000 iterations. Note that, unlike the probability derive from the Monte Carlo Integration program, which is computed based on the area under standard normal probability distribution from 0 to z. The probability computed in this program is based on the area under the probability distribution from negative infinity to z. Figure 3 shows the output when 5 variables are used.
The area of the cells, A16:F29, displays the standard deviation, the mean, and the correlation matrix for each variable based on the generated random numbers. The purpose of this section is to check how well the random numbers meet the specification. For example, the standard deviation and the mean from a standard normal probability distribution should be 1 and 0, respectively. The correlation matrix from the output should also match closely with the input.
Figure 1. Bivariate Standard Normal
Plot With Independency - Correlatoion = 0
Example of derivation of bivariate standard normal random numbers:
Figure 3. Multivariate Standard
Normal Computation With 5 Variables
Figure 2. Bivariate Standard Normal