quietly {
  /*
     This program uses the "corr2data" command to create
       simulated data with specific parameters.  The
       paramters are passed from the parent program as
       arguments which are really local macro variables
       sequentially labelled `1', `2', etc.  The parameters
       are placed in matrices or vectors and then used by
       corr2data.
  */
  capture clear
  * First, input the means, sds, and correlations
  matrix input Mean = (`1' `2')
  matrix input SD   = (`3' `4')
  matrix input Corr = (1.00,`5' \ `5',1.00)

  * Now use corr2data to create a pseudo-simulation of the data
  corr2data age sbp, n(30) means(Mean) corr(Corr) sds(SD)

  * Confirm that all is well
  corr sbp age, means
  #delimit ;
  graph twoway (scatter sbp age),
    xtitle("Age")
    ytitle("SBP")
    xlabel(0(10)100, tposition(inside))
    ylabel(80(20)200, nogrid tposition(inside))
    text(199 1 "r= `5'" 
               "mean X= `1'; mean Y= `2'"
               "  sd X= `3';   sd Y= `4'", 
               placement(se) just(left) size(small))
    plotregion(lcolor(none))
    scheme(socy201)
    name("`6'", replace)
    nodraw
  ;
  #delimit cr
}