The SIMNORMAL Procedure

Getting Started: SIMNORMAL Procedure

The following example illustrates the use of PROC SIMNORMAL to generate two normal random variates that have specified means and covariance.

In this example, the means and covariances are given; these might have come from previous experiments, observational studies, or other considerations.

First you create a _TYPE_=COV data set as the input data set, and then you run PROC SIMNORM with NUMREAL=5000, creating a sample that contains 5,000 observations. The simple statistics of this sample are checked using PROC CORR. The results are shown in Figure 90.1.

data scov(type=COV) ;
   input _TYPE_ $ 1-4 _NAME_ $ 9-10 S1 S2 ;
   datalines ;
COV     S1      1.915  0.3873
COV     S2      0.3873 4.321
MEAN            1.305  2.003
run;
proc simnorm data=scov outsim=ssim
               numreal = 5000
               seed = 54321 ;
   var s1 s2 ;
run;
proc corr data=ssim cov ;
   var s1 s2 ;
   title "Statistics for PROC SIMNORM Sample Using NUMREAL=5000" ;
run;

Figure 90.1: Statistics for PROC SIMNORM Sample Using NUMREAL=5000

Statistics for PROC SIMNORM Sample Using NUMREAL=5000

The CORR Procedure

2 Variables: S1 S2

Covariance Matrix, DF = 4999
  S1 S2
S1 1.895805499 0.424837163
S2 0.424837163 4.132974275

Simple Statistics
Variable N Mean Std Dev Sum Minimum Maximum
S1 5000 1.30254 1.37688 6513 -3.90682 6.49864
S2 5000 1.98790 2.03297 9940 -5.69812 9.42833