cdf of multivariate normal random numbers

3 次查看（过去 30 天）
PChoppala
PChoppala 2013-8-8
Hi
I want to find the cdf of a set of 1000 bivariate normal random numbers that are generated by
X=mvnrnd([3,12],[1,.2;.2,1],1000);
To obtain the cdf for a 1d distribution, I can evaluate the cumsum of the probability of each sample after the samples have been sorted in ascending order. But I am confused on how to do this for a bivariate or (for greater dimension) case. I appreciate you help on this.

请先登录，再回答此问题。

回答（1 个）

Shashank Prasanna
  2 个评论
PChoppala
PChoppala 2013-8-8
To get the cdf I only have X. Is the below piece of code is a correct way of obtaining that (or what you are referring to)?
X=mvnrnd([3,12],[1,.3;.3,1],625);
mu=mean(X);
Sigma=corrcoef(X);
N=25;
[X1,X2] = meshgrid(linspace(1,5,N)', linspace(9,15,N)');
XX = [X1(:) X2(:)];
C = mvncdf(XX, mu, Sigma);
contour(X1,X2,reshape(C,N,N));
Shashank Prasanna
Shashank Prasanna 2013-8-8
It is correct for what it is doing. Also, you already have the means and covariance:
[3,12] and [1,0.3;0.3,1]
1) You don't don't have to re estimate them you can use them directly.
2) MVNCDF expects the Covariance and not correlation matrix. Use COV function instead.

请先登录，再进行评论。

类别

在 Help Center 和 File Exchange 中查找有关 Multivariate Normal Distribution 的更多信息

标签

产品

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by