Univariate and Bivariate Normal Distributions
The following was implemented in Maple by Marcus Davidsson (2008) davidsson_marcus@hotmail.com
1) A Univariate Normal Distribution
A univariate normal distribution has a probability density function equal to
![p(x) = `+`(`/`(`*`(`/`(1, 2), `*`(`^`(2, `/`(1, 2)), `*`(exp(`+`(`-`(`/`(`*`(`/`(1, 2), `*`(`^`(`+`(x, `-`(mu[x])), 2))), `*`(`^`(sigma[x], 2))))))))), `*`(`^`(Pi, `/`(1, 2)), `*`(sigma[x]))))](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_3.gif) |
(1) |
with mean
![mu[x]](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_5.gif) |
(2) |
and standard deviation
![sigma[x]](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_7.gif) |
(3) |
We can plot a univariate normal distribution as follows
![restart; -1; with(stats); -1; plot(statevalf[pdf, normald[0, 2]], -8 .. 8, color = black, thickness = 3, labels = [x, p(x)])](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_8.gif)
![restart; -1; with(stats); -1; plot(statevalf[pdf, normald[0, 2]], -8 .. 8, color = black, thickness = 3, labels = [x, p(x)])](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_9.gif)
![restart; -1; with(stats); -1; plot(statevalf[pdf, normald[0, 2]], -8 .. 8, color = black, thickness = 3, labels = [x, p(x)])](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_10.gif)
We should also note that changing the values of the mean and standard deviation results in different shapes and of
the Probability Density Function (PDF) as seen in the figures below where normal [mean, standard deviation].
Different values for the mean give us:
![restart; -1; with(stats); -1; with(plots); -1; `:=`(x1, plot(statevalf[pdf, normald[0, 1]], -4 .. 4, color = black, thickness = 3)); -1; `:=`(x2, plot(statevalf[pdf, normald[2, 1]], -6 .. 6, color = r...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_13.gif)
![restart; -1; with(stats); -1; with(plots); -1; `:=`(x1, plot(statevalf[pdf, normald[0, 1]], -4 .. 4, color = black, thickness = 3)); -1; `:=`(x2, plot(statevalf[pdf, normald[2, 1]], -6 .. 6, color = r...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_14.gif)
![restart; -1; with(stats); -1; with(plots); -1; `:=`(x1, plot(statevalf[pdf, normald[0, 1]], -4 .. 4, color = black, thickness = 3)); -1; `:=`(x2, plot(statevalf[pdf, normald[2, 1]], -6 .. 6, color = r...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_15.gif)
![restart; -1; with(stats); -1; with(plots); -1; `:=`(x1, plot(statevalf[pdf, normald[0, 1]], -4 .. 4, color = black, thickness = 3)); -1; `:=`(x2, plot(statevalf[pdf, normald[2, 1]], -6 .. 6, color = r...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_16.gif)
![restart; -1; with(stats); -1; with(plots); -1; `:=`(x1, plot(statevalf[pdf, normald[0, 1]], -4 .. 4, color = black, thickness = 3)); -1; `:=`(x2, plot(statevalf[pdf, normald[2, 1]], -6 .. 6, color = r...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_17.gif)
![restart; -1; with(stats); -1; with(plots); -1; `:=`(x1, plot(statevalf[pdf, normald[0, 1]], -4 .. 4, color = black, thickness = 3)); -1; `:=`(x2, plot(statevalf[pdf, normald[2, 1]], -6 .. 6, color = r...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_18.gif)
Different values for the standard deviation
![restart; -1; with(stats); -1; with(plots); -1; `:=`(x1, plot(statevalf[pdf, normald[0, 1]], -4 .. 4, color = black, thickness = 3)); -1; `:=`(x2, plot(statevalf[pdf, normald[0, 2]], -6 .. 6, color = r...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_21.gif)
![restart; -1; with(stats); -1; with(plots); -1; `:=`(x1, plot(statevalf[pdf, normald[0, 1]], -4 .. 4, color = black, thickness = 3)); -1; `:=`(x2, plot(statevalf[pdf, normald[0, 2]], -6 .. 6, color = r...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_22.gif)
![restart; -1; with(stats); -1; with(plots); -1; `:=`(x1, plot(statevalf[pdf, normald[0, 1]], -4 .. 4, color = black, thickness = 3)); -1; `:=`(x2, plot(statevalf[pdf, normald[0, 2]], -6 .. 6, color = r...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_23.gif)
![restart; -1; with(stats); -1; with(plots); -1; `:=`(x1, plot(statevalf[pdf, normald[0, 1]], -4 .. 4, color = black, thickness = 3)); -1; `:=`(x2, plot(statevalf[pdf, normald[0, 2]], -6 .. 6, color = r...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_24.gif)
![restart; -1; with(stats); -1; with(plots); -1; `:=`(x1, plot(statevalf[pdf, normald[0, 1]], -4 .. 4, color = black, thickness = 3)); -1; `:=`(x2, plot(statevalf[pdf, normald[0, 2]], -6 .. 6, color = r...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_25.gif)
![restart; -1; with(stats); -1; with(plots); -1; `:=`(x1, plot(statevalf[pdf, normald[0, 1]], -4 .. 4, color = black, thickness = 3)); -1; `:=`(x2, plot(statevalf[pdf, normald[0, 2]], -6 .. 6, color = r...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_26.gif)
A standard normal distribution has a probability density function equal to
 |
(4) |
with mean
 |
(5) |
and standard deviation
 |
(6) |
We can plot a standard normal distribution as follows
![restart; -1; with(stats); -1; plot(statevalf[pdf, normald[0, 1]], -4 .. 4, color = black, thickness = 3, labels = [x, p(x)])](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_35.gif)
![restart; -1; with(stats); -1; plot(statevalf[pdf, normald[0, 1]], -4 .. 4, color = black, thickness = 3, labels = [x, p(x)])](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_36.gif)
![restart; -1; with(stats); -1; plot(statevalf[pdf, normald[0, 1]], -4 .. 4, color = black, thickness = 3, labels = [x, p(x)])](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_37.gif)
We should also note that a standard normal distribution has an kurtosis equal to 3.
 |
(7) |
2) A Bivariate Normal Distribution
A bivariate normal distribution has a probability density function equal to
![P(x, y) = `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(sqrt(`+`(1, `-`(`*`(`^`(rho, 2)))))))))))), `*`(exp(`+`(`-`(`/`(`*`(`+`(`/`(`*`(`^`(`+`(x, `-`(mu[x])), 2)), `*`(`^`(sigma[x], 2)))...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_43.gif)
![P(x, y) = `+`(`/`(`*`(`/`(1, 2), `*`(exp(`+`(`-`(`/`(`*`(`+`(`/`(`*`(`^`(`+`(x, `-`(mu[x])), 2)), `*`(`^`(sigma[x], 2))), `-`(`/`(`*`(2, `*`(rho, `*`(`+`(x, `-`(mu[x])), `*`(`+`(y, `-`(mu[y])))))), `*...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_45.gif) |
(8) |
We can plot the bivariate normal distribution if we assume different values of
![restart; 1; with(plots); -1; `:=`(rho, 0); -1; `:=`(mu[x], 0); -1; `:=`(mu[y], 0); -1; `:=`(sigma[x], 1); -1; `:=`(sigma[y], 1); -1; `:=`(f, `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_47.gif)
![restart; 1; with(plots); -1; `:=`(rho, 0); -1; `:=`(mu[x], 0); -1; `:=`(mu[y], 0); -1; `:=`(sigma[x], 1); -1; `:=`(sigma[y], 1); -1; `:=`(f, `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_48.gif)
![restart; 1; with(plots); -1; `:=`(rho, 0); -1; `:=`(mu[x], 0); -1; `:=`(mu[y], 0); -1; `:=`(sigma[x], 1); -1; `:=`(sigma[y], 1); -1; `:=`(f, `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_49.gif)
![restart; 1; with(plots); -1; `:=`(rho, 0); -1; `:=`(mu[x], 0); -1; `:=`(mu[y], 0); -1; `:=`(sigma[x], 1); -1; `:=`(sigma[y], 1); -1; `:=`(f, `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_50.gif)
![restart; 1; with(plots); -1; `:=`(rho, 0); -1; `:=`(mu[x], 0); -1; `:=`(mu[y], 0); -1; `:=`(sigma[x], 1); -1; `:=`(sigma[y], 1); -1; `:=`(f, `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_51.gif)
![restart; 1; with(plots); -1; `:=`(rho, 0); -1; `:=`(mu[x], 0); -1; `:=`(mu[y], 0); -1; `:=`(sigma[x], 1); -1; `:=`(sigma[y], 1); -1; `:=`(f, `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_52.gif)
![restart; 1; with(plots); -1; `:=`(rho, 0); -1; `:=`(mu[x], 0); -1; `:=`(mu[y], 0); -1; `:=`(sigma[x], 1); -1; `:=`(sigma[y], 1); -1; `:=`(f, `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_53.gif)
![restart; 1; with(plots); -1; `:=`(rho, 0); -1; `:=`(mu[x], 0); -1; `:=`(mu[y], 0); -1; `:=`(sigma[x], 1); -1; `:=`(sigma[y], 1); -1; `:=`(f, `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_54.gif)
![restart; 1; with(plots); -1; `:=`(rho, 0); -1; `:=`(mu[x], 0); -1; `:=`(mu[y], 0); -1; `:=`(sigma[x], 1); -1; `:=`(sigma[y], 1); -1; `:=`(f, `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_55.gif)
![restart; 1; with(plots); -1; `:=`(rho, 0); -1; `:=`(mu[x], 0); -1; `:=`(mu[y], 0); -1; `:=`(sigma[x], 1); -1; `:=`(sigma[y], 1); -1; `:=`(f, `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_56.gif)
![restart; 1; with(plots); -1; `:=`(rho, 0); -1; `:=`(mu[x], 0); -1; `:=`(mu[y], 0); -1; `:=`(sigma[x], 1); -1; `:=`(sigma[y], 1); -1; `:=`(f, `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_57.gif)
![restart; 1; with(plots); -1; `:=`(rho, 0); -1; `:=`(mu[x], 0); -1; `:=`(mu[y], 0); -1; `:=`(sigma[x], 1); -1; `:=`(sigma[y], 1); -1; `:=`(f, `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_58.gif)
![restart; 1; with(plots); -1; `:=`(rho, 0); -1; `:=`(mu[x], 0); -1; `:=`(mu[y], 0); -1; `:=`(sigma[x], 1); -1; `:=`(sigma[y], 1); -1; `:=`(f, `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_61.gif)
![restart; 1; with(plots); -1; `:=`(rho, 0); -1; `:=`(mu[x], 0); -1; `:=`(mu[y], 0); -1; `:=`(sigma[x], 1); -1; `:=`(sigma[y], 1); -1; `:=`(f, `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_62.gif)
![restart; 1; with(plots); -1; `:=`(rho, 0); -1; `:=`(mu[x], 0); -1; `:=`(mu[y], 0); -1; `:=`(sigma[x], 1); -1; `:=`(sigma[y], 1); -1; `:=`(f, `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_63.gif)
![restart; 1; with(plots); -1; `:=`(rho, 0); -1; `:=`(mu[x], 0); -1; `:=`(mu[y], 0); -1; `:=`(sigma[x], 1); -1; `:=`(sigma[y], 1); -1; `:=`(f, `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_64.gif)
![restart; 1; with(plots); -1; `:=`(rho, 0); -1; `:=`(mu[x], 0); -1; `:=`(mu[y], 0); -1; `:=`(sigma[x], 1); -1; `:=`(sigma[y], 1); -1; `:=`(f, `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_65.gif)
![restart; 1; with(plots); -1; `:=`(rho, 0); -1; `:=`(mu[x], 0); -1; `:=`(mu[y], 0); -1; `:=`(sigma[x], 1); -1; `:=`(sigma[y], 1); -1; `:=`(f, `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_66.gif)
![restart; 1; with(plots); -1; `:=`(rho, 0); -1; `:=`(mu[x], 0); -1; `:=`(mu[y], 0); -1; `:=`(sigma[x], 1); -1; `:=`(sigma[y], 1); -1; `:=`(f, `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_67.gif)
![restart; 1; with(plots); -1; `:=`(rho, 0); -1; `:=`(mu[x], 0); -1; `:=`(mu[y], 0); -1; `:=`(sigma[x], 1); -1; `:=`(sigma[y], 1); -1; `:=`(f, `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_68.gif)
![restart; 1; with(plots); -1; `:=`(rho, 0); -1; `:=`(mu[x], 0); -1; `:=`(mu[y], 0); -1; `:=`(sigma[x], 1); -1; `:=`(sigma[y], 1); -1; `:=`(f, `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_69.gif)
![restart; 1; with(plots); -1; `:=`(rho, 0); -1; `:=`(mu[x], 0); -1; `:=`(mu[y], 0); -1; `:=`(sigma[x], 1); -1; `:=`(sigma[y], 1); -1; `:=`(f, `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_70.gif)
![restart; 1; with(plots); -1; `:=`(rho, 0); -1; `:=`(mu[x], 0); -1; `:=`(mu[y], 0); -1; `:=`(sigma[x], 1); -1; `:=`(sigma[y], 1); -1; `:=`(f, `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_71.gif)
![restart; 1; with(plots); -1; `:=`(rho, 0); -1; `:=`(mu[x], 0); -1; `:=`(mu[y], 0); -1; `:=`(sigma[x], 1); -1; `:=`(sigma[y], 1); -1; `:=`(f, `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_72.gif)
We can now simulate how the bivariate normal distribution would look with different mean values of x
(just click on the chart and click play)
![restart; 1; with(plots); -1; `:=`(rho, 0); -1; `:=`(mu[y], 0); -1; `:=`(sigma[x], 1); -1; `:=`(sigma[y], 1); -1; `:=`(f, `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(sqrt(`+`(1, `-`(`*`(...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_76.gif)
![restart; 1; with(plots); -1; `:=`(rho, 0); -1; `:=`(mu[y], 0); -1; `:=`(sigma[x], 1); -1; `:=`(sigma[y], 1); -1; `:=`(f, `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(sqrt(`+`(1, `-`(`*`(...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_77.gif)
![restart; 1; with(plots); -1; `:=`(rho, 0); -1; `:=`(mu[y], 0); -1; `:=`(sigma[x], 1); -1; `:=`(sigma[y], 1); -1; `:=`(f, `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(sqrt(`+`(1, `-`(`*`(...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_78.gif)
![restart; 1; with(plots); -1; `:=`(rho, 0); -1; `:=`(mu[y], 0); -1; `:=`(sigma[x], 1); -1; `:=`(sigma[y], 1); -1; `:=`(f, `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(sqrt(`+`(1, `-`(`*`(...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_79.gif)
![restart; 1; with(plots); -1; `:=`(rho, 0); -1; `:=`(mu[y], 0); -1; `:=`(sigma[x], 1); -1; `:=`(sigma[y], 1); -1; `:=`(f, `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(sqrt(`+`(1, `-`(`*`(...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_80.gif)
![restart; 1; with(plots); -1; `:=`(rho, 0); -1; `:=`(mu[y], 0); -1; `:=`(sigma[x], 1); -1; `:=`(sigma[y], 1); -1; `:=`(f, `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(sqrt(`+`(1, `-`(`*`(...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_81.gif)
![restart; 1; with(plots); -1; `:=`(rho, 0); -1; `:=`(mu[y], 0); -1; `:=`(sigma[x], 1); -1; `:=`(sigma[y], 1); -1; `:=`(f, `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(sqrt(`+`(1, `-`(`*`(...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_82.gif)
![restart; 1; with(plots); -1; `:=`(rho, 0); -1; `:=`(mu[y], 0); -1; `:=`(sigma[x], 1); -1; `:=`(sigma[y], 1); -1; `:=`(f, `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(sqrt(`+`(1, `-`(`*`(...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_83.gif)
![restart; 1; with(plots); -1; `:=`(rho, 0); -1; `:=`(mu[y], 0); -1; `:=`(sigma[x], 1); -1; `:=`(sigma[y], 1); -1; `:=`(f, `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(sqrt(`+`(1, `-`(`*`(...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_84.gif)
![restart; 1; with(plots); -1; `:=`(rho, 0); -1; `:=`(mu[y], 0); -1; `:=`(sigma[x], 1); -1; `:=`(sigma[y], 1); -1; `:=`(f, `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(sqrt(`+`(1, `-`(`*`(...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_85.gif)
![restart; 1; with(plots); -1; `:=`(rho, 0); -1; `:=`(mu[y], 0); -1; `:=`(sigma[x], 1); -1; `:=`(sigma[y], 1); -1; `:=`(f, `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(sqrt(`+`(1, `-`(`*`(...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_86.gif)
![restart; 1; with(plots); -1; `:=`(rho, 0); -1; `:=`(mu[y], 0); -1; `:=`(sigma[x], 1); -1; `:=`(sigma[y], 1); -1; `:=`(f, `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(sqrt(`+`(1, `-`(`*`(...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_87.gif)
We can now simulate how the distribution would look with different mean values of y
(just click on the chart and click play)
![restart; 1; with(plots); -1; `:=`(rho, 0); -1; `:=`(mu[x], 0); -1; `:=`(sigma[x], 1); -1; `:=`(sigma[y], 1); -1; `:=`(f, `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(sqrt(`+`(1, `-`(`*`(...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_91.gif)
![restart; 1; with(plots); -1; `:=`(rho, 0); -1; `:=`(mu[x], 0); -1; `:=`(sigma[x], 1); -1; `:=`(sigma[y], 1); -1; `:=`(f, `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(sqrt(`+`(1, `-`(`*`(...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_92.gif)
![restart; 1; with(plots); -1; `:=`(rho, 0); -1; `:=`(mu[x], 0); -1; `:=`(sigma[x], 1); -1; `:=`(sigma[y], 1); -1; `:=`(f, `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(sqrt(`+`(1, `-`(`*`(...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_93.gif)
![restart; 1; with(plots); -1; `:=`(rho, 0); -1; `:=`(mu[x], 0); -1; `:=`(sigma[x], 1); -1; `:=`(sigma[y], 1); -1; `:=`(f, `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(sqrt(`+`(1, `-`(`*`(...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_94.gif)
![restart; 1; with(plots); -1; `:=`(rho, 0); -1; `:=`(mu[x], 0); -1; `:=`(sigma[x], 1); -1; `:=`(sigma[y], 1); -1; `:=`(f, `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(sqrt(`+`(1, `-`(`*`(...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_95.gif)
![restart; 1; with(plots); -1; `:=`(rho, 0); -1; `:=`(mu[x], 0); -1; `:=`(sigma[x], 1); -1; `:=`(sigma[y], 1); -1; `:=`(f, `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(sqrt(`+`(1, `-`(`*`(...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_96.gif)
![restart; 1; with(plots); -1; `:=`(rho, 0); -1; `:=`(mu[x], 0); -1; `:=`(sigma[x], 1); -1; `:=`(sigma[y], 1); -1; `:=`(f, `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(sqrt(`+`(1, `-`(`*`(...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_97.gif)
![restart; 1; with(plots); -1; `:=`(rho, 0); -1; `:=`(mu[x], 0); -1; `:=`(sigma[x], 1); -1; `:=`(sigma[y], 1); -1; `:=`(f, `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(sqrt(`+`(1, `-`(`*`(...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_98.gif)
![restart; 1; with(plots); -1; `:=`(rho, 0); -1; `:=`(mu[x], 0); -1; `:=`(sigma[x], 1); -1; `:=`(sigma[y], 1); -1; `:=`(f, `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(sqrt(`+`(1, `-`(`*`(...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_99.gif)
![restart; 1; with(plots); -1; `:=`(rho, 0); -1; `:=`(mu[x], 0); -1; `:=`(sigma[x], 1); -1; `:=`(sigma[y], 1); -1; `:=`(f, `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(sqrt(`+`(1, `-`(`*`(...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_100.gif)
![restart; 1; with(plots); -1; `:=`(rho, 0); -1; `:=`(mu[x], 0); -1; `:=`(sigma[x], 1); -1; `:=`(sigma[y], 1); -1; `:=`(f, `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(sqrt(`+`(1, `-`(`*`(...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_101.gif)
![restart; 1; with(plots); -1; `:=`(rho, 0); -1; `:=`(mu[x], 0); -1; `:=`(sigma[x], 1); -1; `:=`(sigma[y], 1); -1; `:=`(f, `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(sqrt(`+`(1, `-`(`*`(...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_102.gif)
We can now simulate how the bivariate normal distribution would look with different standard deviation values of x
(just click on the chart and click play)
![restart; 1; with(plots); -1; `:=`(rho, 0); -1; `:=`(mu[x], 0); -1; `:=`(mu[y], 0); -1; `:=`(sigma[y], 1); -1; `:=`(f, `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(sqrt(`+`(1, `-`(`*`(`^`...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_106.gif)
![restart; 1; with(plots); -1; `:=`(rho, 0); -1; `:=`(mu[x], 0); -1; `:=`(mu[y], 0); -1; `:=`(sigma[y], 1); -1; `:=`(f, `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(sqrt(`+`(1, `-`(`*`(`^`...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_107.gif)
![restart; 1; with(plots); -1; `:=`(rho, 0); -1; `:=`(mu[x], 0); -1; `:=`(mu[y], 0); -1; `:=`(sigma[y], 1); -1; `:=`(f, `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(sqrt(`+`(1, `-`(`*`(`^`...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_108.gif)
![restart; 1; with(plots); -1; `:=`(rho, 0); -1; `:=`(mu[x], 0); -1; `:=`(mu[y], 0); -1; `:=`(sigma[y], 1); -1; `:=`(f, `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(sqrt(`+`(1, `-`(`*`(`^`...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_109.gif)
![restart; 1; with(plots); -1; `:=`(rho, 0); -1; `:=`(mu[x], 0); -1; `:=`(mu[y], 0); -1; `:=`(sigma[y], 1); -1; `:=`(f, `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(sqrt(`+`(1, `-`(`*`(`^`...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_110.gif)
![restart; 1; with(plots); -1; `:=`(rho, 0); -1; `:=`(mu[x], 0); -1; `:=`(mu[y], 0); -1; `:=`(sigma[y], 1); -1; `:=`(f, `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(sqrt(`+`(1, `-`(`*`(`^`...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_111.gif)
![restart; 1; with(plots); -1; `:=`(rho, 0); -1; `:=`(mu[x], 0); -1; `:=`(mu[y], 0); -1; `:=`(sigma[y], 1); -1; `:=`(f, `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(sqrt(`+`(1, `-`(`*`(`^`...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_112.gif)
![restart; 1; with(plots); -1; `:=`(rho, 0); -1; `:=`(mu[x], 0); -1; `:=`(mu[y], 0); -1; `:=`(sigma[y], 1); -1; `:=`(f, `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(sqrt(`+`(1, `-`(`*`(`^`...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_113.gif)
![restart; 1; with(plots); -1; `:=`(rho, 0); -1; `:=`(mu[x], 0); -1; `:=`(mu[y], 0); -1; `:=`(sigma[y], 1); -1; `:=`(f, `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(sqrt(`+`(1, `-`(`*`(`^`...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_114.gif)
![restart; 1; with(plots); -1; `:=`(rho, 0); -1; `:=`(mu[x], 0); -1; `:=`(mu[y], 0); -1; `:=`(sigma[y], 1); -1; `:=`(f, `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(sqrt(`+`(1, `-`(`*`(`^`...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_115.gif)
![restart; 1; with(plots); -1; `:=`(rho, 0); -1; `:=`(mu[x], 0); -1; `:=`(mu[y], 0); -1; `:=`(sigma[y], 1); -1; `:=`(f, `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(sqrt(`+`(1, `-`(`*`(`^`...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_116.gif)
![restart; 1; with(plots); -1; `:=`(rho, 0); -1; `:=`(mu[x], 0); -1; `:=`(mu[y], 0); -1; `:=`(sigma[y], 1); -1; `:=`(f, `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(sqrt(`+`(1, `-`(`*`(`^`...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_117.gif)
We can now simulate how the bivariate normal distribution would look with different standard deviation values of y
(just click on the chart and click play)
![restart; 1; with(plots); -1; `:=`(rho, 0); -1; `:=`(mu[x], 0); -1; `:=`(mu[y], 0); -1; `:=`(sigma[x], 1); -1; `:=`(f, `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(sqrt(`+`(1, `-`(`*`(`^`...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_121.gif)
![restart; 1; with(plots); -1; `:=`(rho, 0); -1; `:=`(mu[x], 0); -1; `:=`(mu[y], 0); -1; `:=`(sigma[x], 1); -1; `:=`(f, `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(sqrt(`+`(1, `-`(`*`(`^`...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_122.gif)
![restart; 1; with(plots); -1; `:=`(rho, 0); -1; `:=`(mu[x], 0); -1; `:=`(mu[y], 0); -1; `:=`(sigma[x], 1); -1; `:=`(f, `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(sqrt(`+`(1, `-`(`*`(`^`...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_123.gif)
![restart; 1; with(plots); -1; `:=`(rho, 0); -1; `:=`(mu[x], 0); -1; `:=`(mu[y], 0); -1; `:=`(sigma[x], 1); -1; `:=`(f, `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(sqrt(`+`(1, `-`(`*`(`^`...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_124.gif)
![restart; 1; with(plots); -1; `:=`(rho, 0); -1; `:=`(mu[x], 0); -1; `:=`(mu[y], 0); -1; `:=`(sigma[x], 1); -1; `:=`(f, `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(sqrt(`+`(1, `-`(`*`(`^`...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_125.gif)
![restart; 1; with(plots); -1; `:=`(rho, 0); -1; `:=`(mu[x], 0); -1; `:=`(mu[y], 0); -1; `:=`(sigma[x], 1); -1; `:=`(f, `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(sqrt(`+`(1, `-`(`*`(`^`...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_126.gif)
![restart; 1; with(plots); -1; `:=`(rho, 0); -1; `:=`(mu[x], 0); -1; `:=`(mu[y], 0); -1; `:=`(sigma[x], 1); -1; `:=`(f, `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(sqrt(`+`(1, `-`(`*`(`^`...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_127.gif)
![restart; 1; with(plots); -1; `:=`(rho, 0); -1; `:=`(mu[x], 0); -1; `:=`(mu[y], 0); -1; `:=`(sigma[x], 1); -1; `:=`(f, `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(sqrt(`+`(1, `-`(`*`(`^`...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_128.gif)
![restart; 1; with(plots); -1; `:=`(rho, 0); -1; `:=`(mu[x], 0); -1; `:=`(mu[y], 0); -1; `:=`(sigma[x], 1); -1; `:=`(f, `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(sqrt(`+`(1, `-`(`*`(`^`...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_129.gif)
![restart; 1; with(plots); -1; `:=`(rho, 0); -1; `:=`(mu[x], 0); -1; `:=`(mu[y], 0); -1; `:=`(sigma[x], 1); -1; `:=`(f, `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(sqrt(`+`(1, `-`(`*`(`^`...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_130.gif)
![restart; 1; with(plots); -1; `:=`(rho, 0); -1; `:=`(mu[x], 0); -1; `:=`(mu[y], 0); -1; `:=`(sigma[x], 1); -1; `:=`(f, `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(sqrt(`+`(1, `-`(`*`(`^`...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_131.gif)
![restart; 1; with(plots); -1; `:=`(rho, 0); -1; `:=`(mu[x], 0); -1; `:=`(mu[y], 0); -1; `:=`(sigma[x], 1); -1; `:=`(f, `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(sqrt(`+`(1, `-`(`*`(`^`...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_132.gif)
For a standard normal distribution (mean=0 and standrad deviation=1) so we get:
![restart; 1; `:=`(mu[x], 0); -1; `:=`(mu[y], 0); -1; `:=`(sigma[x], 1); -1; `:=`(sigma[y], 1); -1; P(x, y) = `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(sqrt(`+`(1, `-`(`*`(`^`(rho, 2)))...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_135.gif)
![restart; 1; `:=`(mu[x], 0); -1; `:=`(mu[y], 0); -1; `:=`(sigma[x], 1); -1; `:=`(sigma[y], 1); -1; P(x, y) = `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(sqrt(`+`(1, `-`(`*`(`^`(rho, 2)))...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_136.gif)
![restart; 1; `:=`(mu[x], 0); -1; `:=`(mu[y], 0); -1; `:=`(sigma[x], 1); -1; `:=`(sigma[y], 1); -1; P(x, y) = `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(sqrt(`+`(1, `-`(`*`(`^`(rho, 2)))...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_137.gif)
![restart; 1; `:=`(mu[x], 0); -1; `:=`(mu[y], 0); -1; `:=`(sigma[x], 1); -1; `:=`(sigma[y], 1); -1; P(x, y) = `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(sqrt(`+`(1, `-`(`*`(`^`(rho, 2)))...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_138.gif)
![restart; 1; `:=`(mu[x], 0); -1; `:=`(mu[y], 0); -1; `:=`(sigma[x], 1); -1; `:=`(sigma[y], 1); -1; P(x, y) = `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(sqrt(`+`(1, `-`(`*`(`^`(rho, 2)))...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_139.gif)
![restart; 1; `:=`(mu[x], 0); -1; `:=`(mu[y], 0); -1; `:=`(sigma[x], 1); -1; `:=`(sigma[y], 1); -1; P(x, y) = `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(sqrt(`+`(1, `-`(`*`(`^`(rho, 2)))...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_140.gif)
![restart; 1; `:=`(mu[x], 0); -1; `:=`(mu[y], 0); -1; `:=`(sigma[x], 1); -1; `:=`(sigma[y], 1); -1; P(x, y) = `*`(`/`(`+`(`*`(2, `*`(Pi, `*`(sigma[x], `*`(sigma[y], `*`(sqrt(`+`(1, `-`(`*`(`^`(rho, 2)))...](/view.aspx?SI=7169/Univariate_Bivariate_Normal_Distributions_141.gif)
 |
(9) |
We can plot the standard normal distribution if we for example assume that the correlation coefficient ρ is 0






We can now simulate how the distribution would look with different values of
(just click on the chart and click play)







We can now plot the all distributions for for different values of







We can also plot the contour plots








The End !
Legal Notice: � Maplesoft, a division of Waterloo Maple Inc. 2008. Maplesoft and Maple are trademarks of Waterloo Maple Inc. Neither Maplesoft nor the authors are responsible for any errors contained within and are not liable for any damages resulting from the use of this material. This application is intended for non-commercial, non-profit use only. Contact the authors for permission if you wish to use this application in for-profit activities.