Drawing Values From The Distribution
A widely used method for drawing a random vector x from the N-dimensional multivariate normal distribution with mean vector μ and covariance matrix Σ works as follows:
- Find any real matrix A such that A AT = Σ. When Σ is positive-definite, the Cholesky decomposition is typically used, and the extended form of this decomposition can be always be used (as the covariance matrix may be only positive semi-definite) in both cases a suitable matrix A is obtained. An alternative is to use the matrix A = UΛ½ obtained from a spectral decomposition Σ = UΛUT of Σ. The former approach is more computationally straightforward but the matrices A change for different orderings of the elements of the random vector, while the latter approach gives matrices that are related by simple re-orderings. In theory both approaches give equally good ways of determining a suitable matrix A, but there are differences in compuation time.
- Let z = (z1, …, zN)T be a vector whose components are N independent standard normal variates (which can be generated, for example, by using the Box–Muller transform).
- Let x be μ + Az. This has the desired distribution due to the affine transformation property.
Read more about this topic: Multivariate Normal Distribution
Famous quotes containing the words drawing, values and/or distribution:
“Wild as it was, it was hard for me to get rid of the associations of the settlements. Any steady and monotonous sound, to which I did not distinctly attend, passed for a sound of human industry.... Our minds anywhere, when left to themselves, are always thus busily drawing conclusions from false premises.”
—Henry David Thoreau (1817–1862)
“... the loss of belief in future states is politically, though certainly not spiritually, the most significant distinction between our present period and the centuries before. And this loss is definite. For no matter how religious our world may turn again, or how much authentic faith still exists in it, or how deeply our moral values may be rooted in our religious systems, the fear of hell is no longer among the motives which would prevent or stimulate the actions of a majority.”
—Hannah Arendt (1906–1975)
“The man who pretends that the distribution of income in this country reflects the distribution of ability or character is an ignoramus. The man who says that it could by any possible political device be made to do so is an unpractical visionary. But the man who says that it ought to do so is something worse than an ignoramous and more disastrous than a visionary: he is, in the profoundest Scriptural sense of the word, a fool.”
—George Bernard Shaw (1856–1950)