Re: mvnrnd and positive semidefinite covariance matrices




"Frank Wu" <frank.z.wu@xxxxxxxxx> wrote in message
news:ef2c0c9.-1@xxxxxxxxxxxxxxxxxxx
Hi. Newbie question here.

Does anyone know why Matlab's mvnrnd function requires that the
covariance matrix be positive semidefinite? Is it possible to
generate random numbers from a multivariate normal distribution whose
covariance matrix is not positive semidefinite?

As context: I'm currently using the mvnrnd function to generate
random portfolio returns based off of random asset class
performances. Because the covariance matrix of the asset classes is
not positive semidefinite, I've had to use a conversion, which leads
to the portfolio not quite having the distribution that I expected.

Thanks!

A covariance matrix is always positive semidefinite (where this includes
positive definite as a subclass), by definition. If you think the
covariance matrix of your asset classes isn't, either you're thinking of
something other than covariance, or you're calculating it incorrectly. If
you're unsure of that, I suggest you consider the particular case of a
correlation matrix, and in particular examine the symmetry property of both
correlation and covariance matrices.


.



Relevant Pages

  • Re: mvnrnd and positive semidefinite covariance matrices
    ... covariance matrix is not positive semidefinite? ... to the portfolio not quite having the distribution that I expected. ... of real random variables must always be positive semidefinite! ...
    (comp.soft-sys.matlab)
  • mvnrnd and positive semidefinite covariance matrices
    ... Does anyone know why Matlab's mvnrnd function requires that the ... generate random numbers from a multivariate normal distribution whose ... covariance matrix is not positive semidefinite? ...
    (comp.soft-sys.matlab)
  • Re: mvnrnd and positive semidefinite covariance matrices
    ... Does anyone know why Matlab's mvnrnd function requires that the ... covariance matrix is not positive semidefinite? ... It's also sometimes possible have a true cov matrix that is nearly singular, but end up with an estimate that has eigenvalues that are slightly negative due to small numerical errors -- that's an easy fix. ...
    (comp.soft-sys.matlab)
  • Re: Q: properties of covariance matrices
    ... >Let x and y be n-tuplets of normally distributed random variables, ... > C = covbe the nxn covariance matrix of x and y. ... every positive semidefinite matrix is the covariance matrix of some ...
    (sci.math)
  • Re: Q: properties of covariance matrices
    ... >Let x and y be n-tuplets of normally distributed random variables, ... > C = covbe the nxn covariance matrix of x and y. ... every positive semidefinite matrix is the covariance matrix of some ...
    (sci.physics)