BookmarkSubscribeRSS Feed
JennyQLi
Calcite | Level 5


I'm trying to simulate data from a multivariate normal distribution with mean=MU and variance=MVNCOV:

For example, if we want to run 1000 simulations from a multivariate normal with MU={6.3 1.8 3.7} and

MVNCOV={2, 3, 1,

                 3, 8, 5,

                 1, 5, 4}

Are the following codes correct? Thanks!

proc iml;

  nobs=1000;

  MU={6.3 1.8 3.7);

  MVNCOV={2 3 1,

                   3 8 5,

                   1 5 4};

  MVS=ROOT(MVNCOV);

  RandSeed=floor(RANUNI(round(time(),16.0))*10000000); /* simulate random seed from system time */

  x=J(nobs,1,1)*MU + NORMAL(J(nobs,3,RandSeed))*MVS;

  create _simparm0_ from x(| COLNAME={P1 P2 P3} |); append from x;

quit;

run;

2 REPLIES 2
PGStats
Opal | Level 21

I can't help you with proc IML. But you might want to look at proc SIMNORMAL. Its sole purpose is the generation of correlated normal random numbers.

PG

PG

sas-innovate-white.png

Our biggest data and AI event of the year.

Don’t miss the livestream kicking off May 7. It’s free. It’s easy. And it’s the best seat in the house.

Join us virtually with our complimentary SAS Innovate Digital Pass. Watch live or on-demand in multiple languages, with translations available to help you get the most out of every session.

 

Register now!

From The DO Loop
Want more? Visit our blog for more articles like these.
Discussion stats
  • 2 replies
  • 2125 views
  • 6 likes
  • 3 in conversation