-------------
      name:  <unnamed>
       log:  /Users/stephenchang/Desktop/TICR Program/Spring 2017/EPI 265 Epidemiologic Methods III Research Methods in Chronic Disease Epidemiology/Lecture 8/HW8/assign8_stephenchang.smcl
  log type:  smcl
 opened on:  20 May 2017, 16:31:02

.
. *Assignment8
.
. /*
> Modify your data generation code from last week to match the following causal
> DAG:
>
> Where S1 and S2 are indicators for whether an observation is selected into the
> analysis (ie, when S1=0 or S2=0, the observations are excluded).
>
> Write out data generating rules to match this causal DAG for at least 2 values
> S1 (e.g., S1 = 1 for 90% of observations vs S1=1 for 50% of observations) and 2
> values for S2.  Please be sure that your data generating model entails that U
> influences X, Y, and S2.
> For each, write a simulation model generating data with N=500 individuals
> following these data generating rules.  Report the structural causal model for
> each variable in the DAG, including the distributions of the error terms.  Use
> the statistical method of your choice to estimate the effect of X on Y in each
> scenario, assuming that U is unmeasured (ie, your estimate will likely be
> biased because you cannot account for U).
>
> Now assume the null, that X has no effect on Y.  Under each of your four
> scenarios above, what is:
> - The average estimate of the effect of X on Y (ie, if you repeat your
> chosen statistical analysis across repeated iterations, what is the average
> value of the parameter estimate for the association between X and Y)?
> - the type 1 error rate using an alpha threshold of 0.05 (ie, the percent of
> iterations under which you would find a statistically significant association
> between X and Y, using a p<.05 threshold)?*/
.
. if 1 {
. clear
.
. set seed 98770
. set obs 500


number of observations (_N)  was 0,
now 500


. gen errorTerm = rnormal(0, 0.5)
. gen U = runiform(0, 1.)
. gen X = runiform(0, 1) + 0.2 * U
. gen Y = 0 + 0.5*X + 1.*X*(U-0.5) + errorTerm
.
. * generate data
. foreach p1 in 0.2 1.1 ///
> {
  2.   foreach p2 in 1.0 0.1 ///
>   {
  3.
.     display "S1 param = `p2', S2 param = `p1'"
  4.
.     capture drop S1 S2 Xsel Ysel
  5.
.         gen S2 = runiform()>=`p1'*U
  6.     gen S1 = runiform()<=`p2'*(X+Y)
  7.
.         gen Xsel = X if S1 == 1 & S2 == 1
  8.         gen Ysel = Y if S1 == 1 & S2 == 1
  9.
.     * check & visualize data
.     summarize
 10.     * scatter Xsel Ysel
.
.         * simple linear regression of Y on X
.         regress Ysel Xsel
 11.   }
 12. }
S1 param = 1.0, S2 param = 0.2
(189 missing values generated)
(189 missing values generated)

    Variable |        Obs        Mean    Std. Dev.       Min        Max
-------------+---------------------------------------------------------
   errorTerm |        500   -.0069216    .5353013   -1.43531   1.540367
           U |        500    .4802203    .2819826   .0042005   .9938892
           X |        500    .6130487    .3009047   .0142714   1.172624
           Y |        500    .3053068    .6091797   -1.42284   2.335514
          S2 |        500        .894    .3081459          0          1
-------------+---------------------------------------------------------
          S1 |        500          .7    .4587165          0          1
        Xsel |        311    .6982614    .2802072   .0244896   1.172624
        Ysel |        311    .5316799    .5332607  -.7154276   2.335514

      Source |       SS           df       MS      Number of obs   =       311
-------------+----------------------------------   F(1, 309)       =      2.47
       Model |  .697917448         1  .697917448   Prob > F        =    0.1174
    Residual |  87.4558557       309  .283028659   R-squared       =    0.0079
-------------+----------------------------------   Adj R-squared   =    0.0047
       Total |  88.1537732       310   .28436701   Root MSE        =      .532

-------------+----------------------------------------------------------------
         Ysel |      Coef.    Std. Err.       t    P>|t|      [95% Con f. Interval]
-------------+----------------------------------------------------------------
        Xsel |    .1693331    .1078338      1.57    0.117     -.0428484     .3815147
       _cons |    .4134412    .0811146      5.10    0.000      .2538343      .573048
-------------+----------------------------------------------------------------
S1 param = 0.1, S2 param = 0.2
(454 missing values generated)
(454 missing values generated)

    Variable |        Obs        Mean    Std. Dev.       Min        Max
-------------+---------------------------------------------------------
   errorTerm |        500   -.0069216    .5353013   -1.43531   1.540367
           U |        500    .4802203    .2819826   .0042005   .9938892
           X |        500    .6130487    .3009047   .0142714   1.172624
           Y |        500    .3053068    .6091797   -1.42284   2.335514
          S2 |        500        .884    .3205456          0          1
-------------+---------------------------------------------------------
          S1 |        500        .104    .3055663          0          1
        Xsel |         46    .7491884    .2717626   .1120304   1.122617
        Ysel |         46    .7139529    .5350957  -.1785502   2.023463

      Source |       SS           df       MS      Number of obs   =        46
-------------+----------------------------------   F(1, 44)        =      2.56
       Model |  .709390187         1  .709390187   Prob > F        =    0.1165
    Residual |  12.1753422        44  .276712323   R-squared       =    0.0551
-------------+----------------------------------   Adj R-squared   =    0.0336
       Total |  12.8847324        45  .286327387   Root MSE        =    .52603

-------------+----------------------------------------------------------------
         Ysel |      Coef.    Std. Err.       t    P>|t|      [95% Con f. Interval]
-------------+----------------------------------------------------------------
        Xsel |     .462005    .2885482      1.60    0.117     -.1195256     1.043536
       _cons |    .3678241    .2296692      1.60    0.116     -.0950438     .8306921
-------------+----------------------------------------------------------------
S1 param = 1.0, S2 param = 1.1
(344 missing values generated)
(344 missing values generated)

    Variable |        Obs        Mean    Std. Dev.       Min        Max
-------------+---------------------------------------------------------
   errorTerm |        500   -.0069216    .5353013   -1.43531   1.540367
           U |        500    .4802203    .2819826   .0042005   .9938892
           X |        500    .6130487    .3009047   .0142714   1.172624
           Y |        500    .3053068    .6091797   -1.42284   2.335514
          S2 |        500        .506    .5004647          0          1
-------------+---------------------------------------------------------
          S1 |        500        .684    .4653785          0          1
        Xsel |        156    .6825249    .2661343   .0671131   1.098191
        Ysel |        156     .441816    .4958685  -.5124294   1.972102

      Source |       SS           df       MS      Number of obs   =       156
-------------+----------------------------------   F(1, 154)       =      3.08
       Model |  .748395057         1  .748395057   Prob > F        =    0.0810
    Residual |  37.3638668       154  .242622512   R-squared       =    0.0196
-------------+----------------------------------   Adj R-squared   =    0.0133
       Total |  38.1122619       155   .24588556   Root MSE        =    .49257

-------------+----------------------------------------------------------------
         Ysel |      Coef.    Std. Err.       t    P>|t|      [95% Con f. Interval]
-------------+----------------------------------------------------------------
        Xsel |   -.2610951    .1486616     -1.76    0.081     -.5547743     .0325841
       _cons |    .6200199    .1088599      5.70    0.000      .4049685     .8350712
-------------+----------------------------------------------------------------
S1 param = 0.1, S2 param = 1.1
(488 missing values generated)
(488 missing values generated)

    Variable |        Obs        Mean    Std. Dev.       Min        Max
-------------+---------------------------------------------------------
   errorTerm |        500   -.0069216    .5353013   -1.43531   1.540367
           U |        500    .4802203    .2819826   .0042005   .9938892
           X |        500    .6130487    .3009047   .0142714   1.172624
           Y |        500    .3053068    .6091797   -1.42284   2.335514
          S2 |        500        .462    .4990532          0          1
-------------+---------------------------------------------------------
          S1 |        500        .078    .2684402          0          1
        Xsel |         12    .7059129    .2996585   .1511894   1.054284
        Ysel |         12    .7253423    .6196588  -.2001904   1.972102

      Source |       SS           df       MS      Number of obs   =        12
-------------+----------------------------------   F(1, 10)        =      0.12
       Model |  .049598242         1  .049598242   Prob > F        =    0.7375
    Residual |  4.17414905        10  .417414905   R-squared       =    0.0117
-------------+----------------------------------   Adj R-squared   =   -0.0871
       Total |  4.22374729        11  .383977026   Root MSE        =    .64608

-------------+----------------------------------------------------------------
         Ysel |      Coef.    Std. Err.       t    P>|t|      [95% Con f. Interval]
-------------+----------------------------------------------------------------
        Xsel |   -.2240837    .6500714     -0.34    0.737     -1.672533     1.224366
       _cons |    .8835259    .4953464      1.78    0.105     -.2201747     1.987226
-------------+----------------------------------------------------------------
. }

.
. /* The average estimate of the effect of X on Y (ie, if you repeat your chosen
>    statistical analysis across repeated iterations, what is the average value
>    of the parameter estimate for the association between X and Y)? */
.
. if 1 {
. clear
. set obs 50000


number of observations (_N)  was 0,
now 50,000


.
. capture drop errorTerm U X Y
. gen errorTerm = rnormal(0, 0.5)
. gen U = runiform(0, 1.)
. gen X = runiform(0, 1) + 0.2 * U
. gen Y = 0 + 0.5*X + 1.*X*(U-0.5) + errorTerm
.
. * generate data
. foreach p1 in 0.2 1.1 ///
> {
  2.   foreach p2 in 1.0 0.1 ///
>   {
  3.
.     display "S1 param = `p2', S2 param = `p1'"
  4.
.     capture drop S1 S2 Xsel Ysel
  5.
.         gen S2 = runiform()>=`p1'*U
  6.     gen S1 = runiform()<=`p2'*(X+Y)
  7.
.         gen Xsel = X if S1 == 1 & S2 == 1
  8.         gen Ysel = Y if S1 == 1 & S2 == 1
  9.
.         * simple linear regression of Y on X
.         bootstrap, reps(100) size(500): regress Ysel Xsel
 10.   }
 11. }
S1 param = 1.0, S2 param = 0.2
(18,997 missing values generated)
(18,997 missing values generated)
(running regress on estimation sample)

Bootstrap replications (100)
----+--- 1 ---+--- 2 ---+--- 3 ---+--- 4 ---+--- 5
..................................................    50
..................................................   100

Linear regression                               Number of obs     =     31,003
                                                Replications      =        100
                                                Wald chi2(1)      =       4.59
                                                Prob > chi2       =     0.0322
                                                R-squared         =     0.0099
                                                Adj R-squared     =     0.0099
                                                Root MSE          =     0.4577

-------------+----------------------------------------------------------------
             |   Observed    Bootstrap                        Norm al-based
         Ysel |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
-------------+----------------------------------------------------------------
        Xsel |    .1680687     .078472      2.14    0.032      .0142664      .321871
       _cons |    .4277872    .0497253      8.60    0.000      .3303275     .5252469
-------------+----------------------------------------------------------------
S1 param = 0.1, S2 param = 0.2
(45,710 missing values generated)
(45,710 missing values generated)
(running regress on estimation sample)

Bootstrap replications (100)
----+--- 1 ---+--- 2 ---+--- 3 ---+--- 4 ---+--- 5
..................................................    50
..................................................   100

Linear regression                               Number of obs     =      4,290
                                                Replications      =        100
                                                Wald chi2(1)      =      16.88
                                                Prob > chi2       =     0.0000
                                                R-squared         =     0.0353
                                                Adj R-squared     =     0.0351
                                                Root MSE          =     0.4942

-------------+----------------------------------------------------------------
             |   Observed    Bootstrap                        Norm al-based
         Ysel |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
-------------+----------------------------------------------------------------
        Xsel |    .3485647    .0848501      4.11    0.000      .1822615     .5148679
       _cons |    .4045199    .0600283      6.74    0.000      .2868666     .5221732
-------------+----------------------------------------------------------------
S1 param = 1.0, S2 param = 1.1
(35,913 missing values generated)
(35,913 missing values generated)
(running regress on estimation sample)

Bootstrap replications (100)
----+--- 1 ---+--- 2 ---+--- 3 ---+--- 4 ---+--- 5
..................................................    50
..................................................   100

Linear regression                               Number of obs     =     14,087
                                                Replications      =        100
                                                Wald chi2(1)      =       2.17
                                                Prob > chi2       =     0.1410
                                                R-squared         =     0.0035
                                                Adj R-squared     =     0.0034
                                                Root MSE          =     0.4254

-------------+----------------------------------------------------------------
             |   Observed    Bootstrap                        Norm al-based
         Ysel |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
-------------+----------------------------------------------------------------
        Xsel |   -.0930817    .0632368     -1.47    0.141     -.2170235     .0308601
       _cons |    .4974803    .0399946     12.44    0.000      .4190923     .5758683
-------------+----------------------------------------------------------------
S1 param = 0.1, S2 param = 1.1
(48,227 missing values generated)
(48,227 missing values generated)
(running regress on estimation sample)

Bootstrap replications (100)
----+--- 1 ---+--- 2 ---+--- 3 ---+--- 4 ---+--- 5
..................................................    50
..................................................   100

Linear regression                               Number of obs     =      1,773
                                                Replications      =        100
                                                Wald chi2(1)      =       0.09
                                                Prob > chi2       =     0.7601
                                                R-squared         =     0.0002
                                                Adj R-squared     =    -0.0004
                                                Root MSE          =     0.4484

-------------+----------------------------------------------------------------
             |   Observed    Bootstrap                        Norm al-based
         Ysel |      Coef.    Std. Err.       z    P>|z|      [95% Con f. Interval]
-------------+----------------------------------------------------------------
        Xsel |    .0221603    .0725615      0.31    0.760     -.1200577     .1643783
       _cons |    .5035067    .0501339     10.04    0.000      .4052461     .6017674
-------------+----------------------------------------------------------------
. }

.
. /* the type 1 error rate using an alpha threshold of 0.05 (ie, the percent of
>    iterations under which you would find a statistically significant association
>    between X and Y, using a p<.05 threshold)? */
.
. if 1 {
. clear
. set seed 98770
. foreach p1 in 0.2 1.1 ///
> {
  2.   foreach p2 in 1.0 0.1 ///
>   {
  3.
.     display "S1 param = `p2', S2 param = `p1'"
  4.
.         postfile results p using sim_`p1'_`p2', replace
  5.         forvalues n=1/100 {
  6.            set obs 500
  7.
.            capture drop errorTerm U X Y
  8.        gen errorTerm = rnormal(0, 0.5)
  9.        gen U = runiform(0, 1.)
 10.        gen X = runiform(0, 1) + 0.2 * U
 11.        gen Y = 0 + 0.5*X + 1.*X*(U-0.5) + errorTerm
 12.
.        capture drop S1 S2 Xsel Ysel
 13.            gen S2 = runiform()>=`p1'*U
 14.        gen S1 = runiform()<=`p2'*(X+Y)
 15.
.            gen Xsel = X if S1 == 1 & S2 == 1
 16.            gen Ysel = Y if S1 == 1 & S2 == 1
 17.
.            quietly regress Ysel Xsel
 18.            mat b = r(table)
 19.            local p = b[4,1]
 20.            post results (`p')
 21.         }
 22.         postclose results
 23.   }
 24. }
S1 param = 1.0, S2 param = 0.2
(note: file sim_0.2_1.0 not found)


number of observations (_N)  was 0,
now 500


(189 missing values generated)
(189 missing values generated)


number of observations (_N)  was 500,
now 500


(198 missing values generated)
(198 missing values generated)


number of observations (_N)  was 500,
now 500


(200 missing values generated)
(200 missing values generated)


number of observations (_N)  was 500,
now 500


(198 missing values generated)
(198 missing values generated)


number of observations (_N)  was 500,
now 500


(204 missing values generated)
(204 missing values generated)


number of observations (_N)  was 500,
now 500


(183 missing values generated)
(183 missing values generated)


number of observations (_N)  was 500,
now 500


(195 missing values generated)
(195 missing values generated)


number of observations (_N)  was 500,
now 500


(180 missing values generated)
(180 missing values generated)


number of observations (_N)  was 500,
now 500


(191 missing values generated)
(191 missing values generated)


number of observations (_N)  was 500,
now 500


(202 missing values generated)
(202 missing values generated)


number of observations (_N)  was 500,
now 500


(181 missing values generated)
(181 missing values generated)


number of observations (_N)  was 500,
now 500


(191 missing values generated)
(191 missing values generated)


number of observations (_N)  was 500,
now 500


(185 missing values generated)
(185 missing values generated)


number of observations (_N)  was 500,
now 500


(197 missing values generated)
(197 missing values generated)


number of observations (_N)  was 500,
now 500


(182 missing values generated)
(182 missing values generated)


number of observations (_N)  was 500,
now 500


(190 missing values generated)
(190 missing values generated)


number of observations (_N)  was 500,
now 500


(188 missing values generated)
(188 missing values generated)


number of observations (_N)  was 500,
now 500


(190 missing values generated)
(190 missing values generated)


number of observations (_N)  was 500,
now 500


(186 missing values generated)
(186 missing values generated)


number of observations (_N)  was 500,
now 500


(213 missing values generated)
(213 missing values generated)


number of observations (_N)  was 500,
now 500


(192 missing values generated)
(192 missing values generated)


number of observations (_N)  was 500,
now 500


(187 missing values generated)
(187 missing values generated)


number of observations (_N)  was 500,
now 500


(189 missing values generated)
(189 missing values generated)


number of observations (_N)  was 500,
now 500


(177 missing values generated)
(177 missing values generated)


number of observations (_N)  was 500,
now 500


(185 missing values generated)
(185 missing values generated)


number of observations (_N)  was 500,
now 500


(192 missing values generated)
(192 missing values generated)


number of observations (_N)  was 500,
now 500


(198 missing values generated)
(198 missing values generated)


number of observations (_N)  was 500,
now 500


(181 missing values generated)
(181 missing values generated)


number of observations (_N)  was 500,
now 500


(198 missing values generated)
(198 missing values generated)


number of observations (_N)  was 500,
now 500


(195 missing values generated)
(195 missing values generated)


number of observations (_N)  was 500,
now 500


(192 missing values generated)
(192 missing values generated)


number of observations (_N)  was 500,
now 500


(178 missing values generated)
(178 missing values generated)


number of observations (_N)  was 500,
now 500


(192 missing values generated)
(192 missing values generated)


number of observations (_N)  was 500,
now 500


(194 missing values generated)
(194 missing values generated)


number of observations (_N)  was 500,
now 500


(192 missing values generated)
(192 missing values generated)


number of observations (_N)  was 500,
now 500


(185 missing values generated)
(185 missing values generated)


number of observations (_N)  was 500,
now 500


(184 missing values generated)
(184 missing values generated)


number of observations (_N)  was 500,
now 500


(190 missing values generated)
(190 missing values generated)


number of observations (_N)  was 500,
now 500


(196 missing values generated)
(196 missing values generated)


number of observations (_N)  was 500,
now 500


(190 missing values generated)
(190 missing values generated)


number of observations (_N)  was 500,
now 500


(180 missing values generated)
(180 missing values generated)


number of observations (_N)  was 500,
now 500


(187 missing values generated)
(187 missing values generated)


number of observations (_N)  was 500,
now 500


(192 missing values generated)
(192 missing values generated)


number of observations (_N)  was 500,
now 500


(199 missing values generated)
(199 missing values generated)


number of observations (_N)  was 500,
now 500


(174 missing values generated)
(174 missing values generated)


number of observations (_N)  was 500,
now 500


(191 missing values generated)
(191 missing values generated)


number of observations (_N)  was 500,
now 500


(194 missing values generated)
(194 missing values generated)


number of observations (_N)  was 500,
now 500


(189 missing values generated)
(189 missing values generated)


number of observations (_N)  was 500,
now 500


(190 missing values generated)
(190 missing values generated)


number of observations (_N)  was 500,
now 500


(186 missing values generated)
(186 missing values generated)


number of observations (_N)  was 500,
now 500


(182 missing values generated)
(182 missing values generated)


number of observations (_N)  was 500,
now 500


(194 missing values generated)
(194 missing values generated)


number of observations (_N)  was 500,
now 500


(195 missing values generated)
(195 missing values generated)


number of observations (_N)  was 500,
now 500


(203 missing values generated)
(203 missing values generated)


number of observations (_N)  was 500,
now 500


(203 missing values generated)
(203 missing values generated)


number of observations (_N)  was 500,
now 500


(194 missing values generated)
(194 missing values generated)


number of observations (_N)  was 500,
now 500


(196 missing values generated)
(196 missing values generated)


number of observations (_N)  was 500,
now 500


(208 missing values generated)
(208 missing values generated)


number of observations (_N)  was 500,
now 500


(204 missing values generated)
(204 missing values generated)


number of observations (_N)  was 500,
now 500


(187 missing values generated)
(187 missing values generated)


number of observations (_N)  was 500,
now 500


(200 missing values generated)
(200 missing values generated)


number of observations (_N)  was 500,
now 500


(188 missing values generated)
(188 missing values generated)


number of observations (_N)  was 500,
now 500


(190 missing values generated)
(190 missing values generated)


number of observations (_N)  was 500,
now 500


(204 missing values generated)
(204 missing values generated)


number of observations (_N)  was 500,
now 500


(202 missing values generated)
(202 missing values generated)


number of observations (_N)  was 500,
now 500


(181 missing values generated)
(181 missing values generated)


number of observations (_N)  was 500,
now 500


(196 missing values generated)
(196 missing values generated)


number of observations (_N)  was 500,
now 500


(185 missing values generated)
(185 missing values generated)


number of observations (_N)  was 500,
now 500


(194 missing values generated)
(194 missing values generated)


number of observations (_N)  was 500,
now 500


(208 missing values generated)
(208 missing values generated)


number of observations (_N)  was 500,
now 500


(180 missing values generated)
(180 missing values generated)


number of observations (_N)  was 500,
now 500


(208 missing values generated)
(208 missing values generated)


number of observations (_N)  was 500,
now 500


(189 missing values generated)
(189 missing values generated)


number of observations (_N)  was 500,
now 500


(170 missing values generated)
(170 missing values generated)


number of observations (_N)  was 500,
now 500


(201 missing values generated)
(201 missing values generated)


number of observations (_N)  was 500,
now 500


(185 missing values generated)
(185 missing values generated)


number of observations (_N)  was 500,
now 500


(180 missing values generated)
(180 missing values generated)


number of observations (_N)  was 500,
now 500


(193 missing values generated)
(193 missing values generated)


number of observations (_N)  was 500,
now 500


(174 missing values generated)
(174 missing values generated)


number of observations (_N)  was 500,
now 500


(173 missing values generated)
(173 missing values generated)


number of observations (_N)  was 500,
now 500


(183 missing values generated)
(183 missing values generated)


number of observations (_N)  was 500,
now 500


(189 missing values generated)
(189 missing values generated)


number of observations (_N)  was 500,
now 500


(158 missing values generated)
(158 missing values generated)


number of observations (_N)  was 500,
now 500


(179 missing values generated)
(179 missing values generated)


number of observations (_N)  was 500,
now 500


(207 missing values generated)
(207 missing values generated)


number of observations (_N)  was 500,
now 500


(178 missing values generated)
(178 missing values generated)


number of observations (_N)  was 500,
now 500


(182 missing values generated)
(182 missing values generated)


number of observations (_N)  was 500,
now 500


(189 missing values generated)
(189 missing values generated)


number of observations (_N)  was 500,
now 500


(189 missing values generated)
(189 missing values generated)


number of observations (_N)  was 500,
now 500


(193 missing values generated)
(193 missing values generated)


number of observations (_N)  was 500,
now 500


(183 missing values generated)
(183 missing values generated)


number of observations (_N)  was 500,
now 500


(171 missing values generated)
(171 missing values generated)


number of observations (_N)  was 500,
now 500


(187 missing values generated)
(187 missing values generated)


number of observations (_N)  was 500,
now 500


(200 missing values generated)
(200 missing values generated)


number of observations (_N)  was 500,
now 500


(184 missing values generated)
(184 missing values generated)


number of observations (_N)  was 500,
now 500


(205 missing values generated)
(205 missing values generated)


number of observations (_N)  was 500,
now 500


(201 missing values generated)
(201 missing values generated)


number of observations (_N)  was 500,
now 500


(172 missing values generated)
(172 missing values generated)


number of observations (_N)  was 500,
now 500


(175 missing values generated)
(175 missing values generated)


number of observations (_N)  was 500,
now 500


(193 missing values generated)
(193 missing values generated)
S1 param = 0.1, S2 param = 0.2
(note: file sim_0.2_0.1 not found)


number of observations (_N)  was 500,
now 500


(464 missing values generated)
(464 missing values generated)


number of observations (_N)  was 500,
now 500


(466 missing values generated)
(466 missing values generated)


number of observations (_N)  was 500,
now 500


(459 missing values generated)
(459 missing values generated)


number of observations (_N)  was 500,
now 500


(451 missing values generated)
(451 missing values generated)


number of observations (_N)  was 500,
now 500


(454 missing values generated)
(454 missing values generated)


number of observations (_N)  was 500,
now 500


(458 missing values generated)
(458 missing values generated)


number of observations (_N)  was 500,
now 500


(466 missing values generated)
(466 missing values generated)


number of observations (_N)  was 500,
now 500


(457 missing values generated)
(457 missing values generated)


number of observations (_N)  was 500,
now 500


(459 missing values generated)
(459 missing values generated)


number of observations (_N)  was 500,
now 500


(465 missing values generated)
(465 missing values generated)


number of observations (_N)  was 500,
now 500


(458 missing values generated)
(458 missing values generated)


number of observations (_N)  was 500,
now 500


(464 missing values generated)
(464 missing values generated)


number of observations (_N)  was 500,
now 500


(459 missing values generated)
(459 missing values generated)


number of observations (_N)  was 500,
now 500


(460 missing values generated)
(460 missing values generated)


number of observations (_N)  was 500,
now 500


(457 missing values generated)
(457 missing values generated)


number of observations (_N)  was 500,
now 500


(455 missing values generated)
(455 missing values generated)


number of observations (_N)  was 500,
now 500


(461 missing values generated)
(461 missing values generated)


number of observations (_N)  was 500,
now 500


(459 missing values generated)
(459 missing values generated)


number of observations (_N)  was 500,
now 500


(445 missing values generated)
(445 missing values generated)


number of observations (_N)  was 500,
now 500


(465 missing values generated)
(465 missing values generated)


number of observations (_N)  was 500,
now 500


(457 missing values generated)
(457 missing values generated)


number of observations (_N)  was 500,
now 500


(455 missing values generated)
(455 missing values generated)


number of observations (_N)  was 500,
now 500


(461 missing values generated)
(461 missing values generated)


number of observations (_N)  was 500,
now 500


(461 missing values generated)
(461 missing values generated)


number of observations (_N)  was 500,
now 500


(464 missing values generated)
(464 missing values generated)


number of observations (_N)  was 500,
now 500


(453 missing values generated)
(453 missing values generated)


number of observations (_N)  was 500,
now 500


(457 missing values generated)
(457 missing values generated)


number of observations (_N)  was 500,
now 500


(454 missing values generated)
(454 missing values generated)


number of observations (_N)  was 500,
now 500


(459 missing values generated)
(459 missing values generated)


number of observations (_N)  was 500,
now 500


(454 missing values generated)
(454 missing values generated)


number of observations (_N)  was 500,
now 500


(443 missing values generated)
(443 missing values generated)


number of observations (_N)  was 500,
now 500


(461 missing values generated)
(461 missing values generated)


number of observations (_N)  was 500,
now 500


(455 missing values generated)
(455 missing values generated)


number of observations (_N)  was 500,
now 500


(462 missing values generated)
(462 missing values generated)


number of observations (_N)  was 500,
now 500


(446 missing values generated)
(446 missing values generated)


number of observations (_N)  was 500,
now 500


(461 missing values generated)
(461 missing values generated)


number of observations (_N)  was 500,
now 500


(465 missing values generated)
(465 missing values generated)


number of observations (_N)  was 500,
now 500


(452 missing values generated)
(452 missing values generated)


number of observations (_N)  was 500,
now 500


(457 missing values generated)
(457 missing values generated)


number of observations (_N)  was 500,
now 500


(457 missing values generated)
(457 missing values generated)


number of observations (_N)  was 500,
now 500


(467 missing values generated)
(467 missing values generated)


number of observations (_N)  was 500,
now 500


(451 missing values generated)
(451 missing values generated)


number of observations (_N)  was 500,
now 500


(462 missing values generated)
(462 missing values generated)


number of observations (_N)  was 500,
now 500


(453 missing values generated)
(453 missing values generated)


number of observations (_N)  was 500,
now 500


(452 missing values generated)
(452 missing values generated)


number of observations (_N)  was 500,
now 500


(467 missing values generated)
(467 missing values generated)


number of observations (_N)  was 500,
now 500


(463 missing values generated)
(463 missing values generated)


number of observations (_N)  was 500,
now 500


(451 missing values generated)
(451 missing values generated)


number of observations (_N)  was 500,
now 500


(462 missing values generated)
(462 missing values generated)


number of observations (_N)  was 500,
now 500


(457 missing values generated)
(457 missing values generated)


number of observations (_N)  was 500,
now 500


(453 missing values generated)
(453 missing values generated)


number of observations (_N)  was 500,
now 500


(459 missing values generated)
(459 missing values generated)


number of observations (_N)  was 500,
now 500


(454 missing values generated)
(454 missing values generated)


number of observations (_N)  was 500,
now 500


(450 missing values generated)
(450 missing values generated)


number of observations (_N)  was 500,
now 500


(453 missing values generated)
(453 missing values generated)


number of observations (_N)  was 500,
now 500


(451 missing values generated)
(451 missing values generated)


number of observations (_N)  was 500,
now 500


(458 missing values generated)
(458 missing values generated)


number of observations (_N)  was 500,
now 500


(454 missing values generated)
(454 missing values generated)


number of observations (_N)  was 500,
now 500


(457 missing values generated)
(457 missing values generated)


number of observations (_N)  was 500,
now 500


(456 missing values generated)
(456 missing values generated)


number of observations (_N)  was 500,
now 500


(453 missing values generated)
(453 missing values generated)


number of observations (_N)  was 500,
now 500


(458 missing values generated)
(458 missing values generated)


number of observations (_N)  was 500,
now 500


(468 missing values generated)
(468 missing values generated)


number of observations (_N)  was 500,
now 500


(453 missing values generated)
(453 missing values generated)


number of observations (_N)  was 500,
now 500


(459 missing values generated)
(459 missing values generated)


number of observations (_N)  was 500,
now 500


(456 missing values generated)
(456 missing values generated)


number of observations (_N)  was 500,
now 500


(458 missing values generated)
(458 missing values generated)


number of observations (_N)  was 500,
now 500


(452 missing values generated)
(452 missing values generated)


number of observations (_N)  was 500,
now 500


(452 missing values generated)
(452 missing values generated)


number of observations (_N)  was 500,
now 500


(458 missing values generated)
(458 missing values generated)


number of observations (_N)  was 500,
now 500


(464 missing values generated)
(464 missing values generated)


number of observations (_N)  was 500,
now 500


(461 missing values generated)
(461 missing values generated)


number of observations (_N)  was 500,
now 500


(457 missing values generated)
(457 missing values generated)


number of observations (_N)  was 500,
now 500


(461 missing values generated)
(461 missing values generated)


number of observations (_N)  was 500,
now 500


(453 missing values generated)
(453 missing values generated)


number of observations (_N)  was 500,
now 500


(450 missing values generated)
(450 missing values generated)


number of observations (_N)  was 500,
now 500


(460 missing values generated)
(460 missing values generated)


number of observations (_N)  was 500,
now 500


(469 missing values generated)
(469 missing values generated)


number of observations (_N)  was 500,
now 500


(460 missing values generated)
(460 missing values generated)


number of observations (_N)  was 500,
now 500


(464 missing values generated)
(464 missing values generated)


number of observations (_N)  was 500,
now 500


(453 missing values generated)
(453 missing values generated)


number of observations (_N)  was 500,
now 500


(467 missing values generated)
(467 missing values generated)


number of observations (_N)  was 500,
now 500


(461 missing values generated)
(461 missing values generated)


number of observations (_N)  was 500,
now 500


(463 missing values generated)
(463 missing values generated)


number of observations (_N)  was 500,
now 500


(444 missing values generated)
(444 missing values generated)


number of observations (_N)  was 500,
now 500


(456 missing values generated)
(456 missing values generated)


number of observations (_N)  was 500,
now 500


(457 missing values generated)
(457 missing values generated)


number of observations (_N)  was 500,
now 500


(460 missing values generated)
(460 missing values generated)


number of observations (_N)  was 500,
now 500


(449 missing values generated)
(449 missing values generated)


number of observations (_N)  was 500,
now 500


(456 missing values generated)
(456 missing values generated)


number of observations (_N)  was 500,
now 500


(455 missing values generated)
(455 missing values generated)


number of observations (_N)  was 500,
now 500


(461 missing values generated)
(461 missing values generated)


number of observations (_N)  was 500,
now 500


(454 missing values generated)
(454 missing values generated)


number of observations (_N)  was 500,
now 500


(463 missing values generated)
(463 missing values generated)


number of observations (_N)  was 500,
now 500


(463 missing values generated)
(463 missing values generated)


number of observations (_N)  was 500,
now 500


(450 missing values generated)
(450 missing values generated)


number of observations (_N)  was 500,
now 500


(450 missing values generated)
(450 missing values generated)


number of observations (_N)  was 500,
now 500


(454 missing values generated)
(454 missing values generated)


number of observations (_N)  was 500,
now 500


(450 missing values generated)
(450 missing values generated)


number of observations (_N)  was 500,
now 500


(468 missing values generated)
(468 missing values generated)
S1 param = 1.0, S2 param = 1.1
(note: file sim_1.1_1.0 not found)


number of observations (_N)  was 500,
now 500


(360 missing values generated)
(360 missing values generated)


number of observations (_N)  was 500,
now 500


(354 missing values generated)
(354 missing values generated)


number of observations (_N)  was 500,
now 500


(373 missing values generated)
(373 missing values generated)


number of observations (_N)  was 500,
now 500


(346 missing values generated)
(346 missing values generated)


number of observations (_N)  was 500,
now 500


(364 missing values generated)
(364 missing values generated)


number of observations (_N)  was 500,
now 500


(345 missing values generated)
(345 missing values generated)


number of observations (_N)  was 500,
now 500


(372 missing values generated)
(372 missing values generated)


number of observations (_N)  was 500,
now 500


(350 missing values generated)
(350 missing values generated)


number of observations (_N)  was 500,
now 500


(368 missing values generated)
(368 missing values generated)


number of observations (_N)  was 500,
now 500


(359 missing values generated)
(359 missing values generated)


number of observations (_N)  was 500,
now 500


(337 missing values generated)
(337 missing values generated)


number of observations (_N)  was 500,
now 500


(355 missing values generated)
(355 missing values generated)


number of observations (_N)  was 500,
now 500


(345 missing values generated)
(345 missing values generated)


number of observations (_N)  was 500,
now 500


(374 missing values generated)
(374 missing values generated)


number of observations (_N)  was 500,
now 500


(350 missing values generated)
(350 missing values generated)


number of observations (_N)  was 500,
now 500


(363 missing values generated)
(363 missing values generated)


number of observations (_N)  was 500,
now 500


(355 missing values generated)
(355 missing values generated)


number of observations (_N)  was 500,
now 500


(363 missing values generated)
(363 missing values generated)


number of observations (_N)  was 500,
now 500


(360 missing values generated)
(360 missing values generated)


number of observations (_N)  was 500,
now 500


(363 missing values generated)
(363 missing values generated)


number of observations (_N)  was 500,
now 500


(351 missing values generated)
(351 missing values generated)


number of observations (_N)  was 500,
now 500


(359 missing values generated)
(359 missing values generated)


number of observations (_N)  was 500,
now 500


(365 missing values generated)
(365 missing values generated)


number of observations (_N)  was 500,
now 500


(356 missing values generated)
(356 missing values generated)


number of observations (_N)  was 500,
now 500


(365 missing values generated)
(365 missing values generated)


number of observations (_N)  was 500,
now 500


(368 missing values generated)
(368 missing values generated)


number of observations (_N)  was 500,
now 500


(362 missing values generated)
(362 missing values generated)


number of observations (_N)  was 500,
now 500


(349 missing values generated)
(349 missing values generated)


number of observations (_N)  was 500,
now 500


(366 missing values generated)
(366 missing values generated)


number of observations (_N)  was 500,
now 500


(346 missing values generated)
(346 missing values generated)


number of observations (_N)  was 500,
now 500


(347 missing values generated)
(347 missing values generated)


number of observations (_N)  was 500,
now 500


(359 missing values generated)
(359 missing values generated)


number of observations (_N)  was 500,
now 500


(364 missing values generated)
(364 missing values generated)


number of observations (_N)  was 500,
now 500


(356 missing values generated)
(356 missing values generated)


number of observations (_N)  was 500,
now 500


(359 missing values generated)
(359 missing values generated)


number of observations (_N)  was 500,
now 500


(363 missing values generated)
(363 missing values generated)


number of observations (_N)  was 500,
now 500


(356 missing values generated)
(356 missing values generated)


number of observations (_N)  was 500,
now 500


(353 missing values generated)
(353 missing values generated)


number of observations (_N)  was 500,
now 500


(383 missing values generated)
(383 missing values generated)


number of observations (_N)  was 500,
now 500


(359 missing values generated)
(359 missing values generated)


number of observations (_N)  was 500,
now 500


(356 missing values generated)
(356 missing values generated)


number of observations (_N)  was 500,
now 500


(368 missing values generated)
(368 missing values generated)


number of observations (_N)  was 500,
now 500


(360 missing values generated)
(360 missing values generated)


number of observations (_N)  was 500,
now 500


(366 missing values generated)
(366 missing values generated)


number of observations (_N)  was 500,
now 500


(362 missing values generated)
(362 missing values generated)


number of observations (_N)  was 500,
now 500


(363 missing values generated)
(363 missing values generated)


number of observations (_N)  was 500,
now 500


(350 missing values generated)
(350 missing values generated)


number of observations (_N)  was 500,
now 500


(349 missing values generated)
(349 missing values generated)


number of observations (_N)  was 500,
now 500


(353 missing values generated)
(353 missing values generated)


number of observations (_N)  was 500,
now 500


(357 missing values generated)
(357 missing values generated)


number of observations (_N)  was 500,
now 500


(344 missing values generated)
(344 missing values generated)


number of observations (_N)  was 500,
now 500


(364 missing values generated)
(364 missing values generated)


number of observations (_N)  was 500,
now 500


(365 missing values generated)
(365 missing values generated)


number of observations (_N)  was 500,
now 500


(349 missing values generated)
(349 missing values generated)


number of observations (_N)  was 500,
now 500


(360 missing values generated)
(360 missing values generated)


number of observations (_N)  was 500,
now 500


(368 missing values generated)
(368 missing values generated)


number of observations (_N)  was 500,
now 500


(340 missing values generated)
(340 missing values generated)


number of observations (_N)  was 500,
now 500


(363 missing values generated)
(363 missing values generated)


number of observations (_N)  was 500,
now 500


(357 missing values generated)
(357 missing values generated)


number of observations (_N)  was 500,
now 500


(355 missing values generated)
(355 missing values generated)


number of observations (_N)  was 500,
now 500


(353 missing values generated)
(353 missing values generated)


number of observations (_N)  was 500,
now 500


(356 missing values generated)
(356 missing values generated)


number of observations (_N)  was 500,
now 500


(353 missing values generated)
(353 missing values generated)


number of observations (_N)  was 500,
now 500


(344 missing values generated)
(344 missing values generated)


number of observations (_N)  was 500,
now 500


(362 missing values generated)
(362 missing values generated)


number of observations (_N)  was 500,
now 500


(363 missing values generated)
(363 missing values generated)


number of observations (_N)  was 500,
now 500


(361 missing values generated)
(361 missing values generated)


number of observations (_N)  was 500,
now 500


(351 missing values generated)
(351 missing values generated)


number of observations (_N)  was 500,
now 500


(333 missing values generated)
(333 missing values generated)


number of observations (_N)  was 500,
now 500


(342 missing values generated)
(342 missing values generated)


number of observations (_N)  was 500,
now 500


(343 missing values generated)
(343 missing values generated)


number of observations (_N)  was 500,
now 500


(358 missing values generated)
(358 missing values generated)


number of observations (_N)  was 500,
now 500


(364 missing values generated)
(364 missing values generated)


number of observations (_N)  was 500,
now 500


(354 missing values generated)
(354 missing values generated)


number of observations (_N)  was 500,
now 500


(364 missing values generated)
(364 missing values generated)


number of observations (_N)  was 500,
now 500


(340 missing values generated)
(340 missing values generated)


number of observations (_N)  was 500,
now 500


(345 missing values generated)
(345 missing values generated)


number of observations (_N)  was 500,
now 500


(371 missing values generated)
(371 missing values generated)


number of observations (_N)  was 500,
now 500


(373 missing values generated)
(373 missing values generated)


number of observations (_N)  was 500,
now 500


(367 missing values generated)
(367 missing values generated)


number of observations (_N)  was 500,
now 500


(358 missing values generated)
(358 missing values generated)


number of observations (_N)  was 500,
now 500


(362 missing values generated)
(362 missing values generated)


number of observations (_N)  was 500,
now 500


(361 missing values generated)
(361 missing values generated)


number of observations (_N)  was 500,
now 500


(363 missing values generated)
(363 missing values generated)


number of observations (_N)  was 500,
now 500


(343 missing values generated)
(343 missing values generated)


number of observations (_N)  was 500,
now 500


(355 missing values generated)
(355 missing values generated)


number of observations (_N)  was 500,
now 500


(350 missing values generated)
(350 missing values generated)


number of observations (_N)  was 500,
now 500


(359 missing values generated)
(359 missing values generated)


number of observations (_N)  was 500,
now 500


(369 missing values generated)
(369 missing values generated)


number of observations (_N)  was 500,
now 500


(372 missing values generated)
(372 missing values generated)


number of observations (_N)  was 500,
now 500


(367 missing values generated)
(367 missing values generated)


number of observations (_N)  was 500,
now 500


(364 missing values generated)
(364 missing values generated)


number of observations (_N)  was 500,
now 500


(362 missing values generated)
(362 missing values generated)


number of observations (_N)  was 500,
now 500


(354 missing values generated)
(354 missing values generated)


number of observations (_N)  was 500,
now 500


(364 missing values generated)
(364 missing values generated)


number of observations (_N)  was 500,
now 500


(356 missing values generated)
(356 missing values generated)


number of observations (_N)  was 500,
now 500


(350 missing values generated)
(350 missing values generated)


number of observations (_N)  was 500,
now 500


(358 missing values generated)
(358 missing values generated)


number of observations (_N)  was 500,
now 500


(360 missing values generated)
(360 missing values generated)


number of observations (_N)  was 500,
now 500


(356 missing values generated)
(356 missing values generated)
S1 param = 0.1, S2 param = 1.1
(note: file sim_1.1_0.1 not found)


number of observations (_N)  was 500,
now 500


(483 missing values generated)
(483 missing values generated)


number of observations (_N)  was 500,
now 500


(476 missing values generated)
(476 missing values generated)


number of observations (_N)  was 500,
now 500


(489 missing values generated)
(489 missing values generated)


number of observations (_N)  was 500,
now 500


(480 missing values generated)
(480 missing values generated)


number of observations (_N)  was 500,
now 500


(485 missing values generated)
(485 missing values generated)


number of observations (_N)  was 500,
now 500


(485 missing values generated)
(485 missing values generated)


number of observations (_N)  was 500,
now 500


(488 missing values generated)
(488 missing values generated)


number of observations (_N)  was 500,
now 500


(481 missing values generated)
(481 missing values generated)


number of observations (_N)  was 500,
now 500


(485 missing values generated)
(485 missing values generated)


number of observations (_N)  was 500,
now 500


(479 missing values generated)
(479 missing values generated)


number of observations (_N)  was 500,
now 500


(489 missing values generated)
(489 missing values generated)


number of observations (_N)  was 500,
now 500


(477 missing values generated)
(477 missing values generated)


number of observations (_N)  was 500,
now 500


(478 missing values generated)
(478 missing values generated)


number of observations (_N)  was 500,
now 500


(486 missing values generated)
(486 missing values generated)


number of observations (_N)  was 500,
now 500


(487 missing values generated)
(487 missing values generated)


number of observations (_N)  was 500,
now 500


(488 missing values generated)
(488 missing values generated)


number of observations (_N)  was 500,
now 500


(470 missing values generated)
(470 missing values generated)


number of observations (_N)  was 500,
now 500


(479 missing values generated)
(479 missing values generated)


number of observations (_N)  was 500,
now 500


(483 missing values generated)
(483 missing values generated)


number of observations (_N)  was 500,
now 500


(474 missing values generated)
(474 missing values generated)


number of observations (_N)  was 500,
now 500


(485 missing values generated)
(485 missing values generated)


number of observations (_N)  was 500,
now 500


(484 missing values generated)
(484 missing values generated)


number of observations (_N)  was 500,
now 500


(479 missing values generated)
(479 missing values generated)


number of observations (_N)  was 500,
now 500


(484 missing values generated)
(484 missing values generated)


number of observations (_N)  was 500,
now 500


(483 missing values generated)
(483 missing values generated)


number of observations (_N)  was 500,
now 500


(487 missing values generated)
(487 missing values generated)


number of observations (_N)  was 500,
now 500


(483 missing values generated)
(483 missing values generated)


number of observations (_N)  was 500,
now 500


(478 missing values generated)
(478 missing values generated)


number of observations (_N)  was 500,
now 500


(481 missing values generated)
(481 missing values generated)


number of observations (_N)  was 500,
now 500


(486 missing values generated)
(486 missing values generated)


number of observations (_N)  was 500,
now 500


(481 missing values generated)
(481 missing values generated)


number of observations (_N)  was 500,
now 500


(482 missing values generated)
(482 missing values generated)


number of observations (_N)  was 500,
now 500


(487 missing values generated)
(487 missing values generated)


number of observations (_N)  was 500,
now 500


(481 missing values generated)
(481 missing values generated)


number of observations (_N)  was 500,
now 500


(483 missing values generated)
(483 missing values generated)


number of observations (_N)  was 500,
now 500


(479 missing values generated)
(479 missing values generated)


number of observations (_N)  was 500,
now 500


(483 missing values generated)
(483 missing values generated)


number of observations (_N)  was 500,
now 500


(481 missing values generated)
(481 missing values generated)


number of observations (_N)  was 500,
now 500


(487 missing values generated)
(487 missing values generated)


number of observations (_N)  was 500,
now 500


(483 missing values generated)
(483 missing values generated)


number of observations (_N)  was 500,
now 500


(481 missing values generated)
(481 missing values generated)


number of observations (_N)  was 500,
now 500


(492 missing values generated)
(492 missing values generated)


number of observations (_N)  was 500,
now 500


(487 missing values generated)
(487 missing values generated)


number of observations (_N)  was 500,
now 500


(479 missing values generated)
(479 missing values generated)


number of observations (_N)  was 500,
now 500


(485 missing values generated)
(485 missing values generated)


number of observations (_N)  was 500,
now 500


(480 missing values generated)
(480 missing values generated)


number of observations (_N)  was 500,
now 500


(484 missing values generated)
(484 missing values generated)


number of observations (_N)  was 500,
now 500


(486 missing values generated)
(486 missing values generated)


number of observations (_N)  was 500,
now 500


(490 missing values generated)
(490 missing values generated)


number of observations (_N)  was 500,
now 500


(482 missing values generated)
(482 missing values generated)


number of observations (_N)  was 500,
now 500


(485 missing values generated)
(485 missing values generated)


number of observations (_N)  was 500,
now 500


(486 missing values generated)
(486 missing values generated)


number of observations (_N)  was 500,
now 500


(471 missing values generated)
(471 missing values generated)


number of observations (_N)  was 500,
now 500


(488 missing values generated)
(488 missing values generated)


number of observations (_N)  was 500,
now 500


(491 missing values generated)
(491 missing values generated)


number of observations (_N)  was 500,
now 500


(485 missing values generated)
(485 missing values generated)


number of observations (_N)  was 500,
now 500


(484 missing values generated)
(484 missing values generated)


number of observations (_N)  was 500,
now 500


(474 missing values generated)
(474 missing values generated)


number of observations (_N)  was 500,
now 500


(479 missing values generated)
(479 missing values generated)


number of observations (_N)  was 500,
now 500


(478 missing values generated)
(478 missing values generated)


number of observations (_N)  was 500,
now 500


(481 missing values generated)
(481 missing values generated)


number of observations (_N)  was 500,
now 500


(473 missing values generated)
(473 missing values generated)


number of observations (_N)  was 500,
now 500


(482 missing values generated)
(482 missing values generated)


number of observations (_N)  was 500,
now 500


(480 missing values generated)
(480 missing values generated)


number of observations (_N)  was 500,
now 500


(484 missing values generated)
(484 missing values generated)


number of observations (_N)  was 500,
now 500


(480 missing values generated)
(480 missing values generated)


number of observations (_N)  was 500,
now 500


(484 missing values generated)
(484 missing values generated)


number of observations (_N)  was 500,
now 500


(484 missing values generated)
(484 missing values generated)


number of observations (_N)  was 500,
now 500


(482 missing values generated)
(482 missing values generated)


number of observations (_N)  was 500,
now 500


(484 missing values generated)
(484 missing values generated)


number of observations (_N)  was 500,
now 500


(482 missing values generated)
(482 missing values generated)


number of observations (_N)  was 500,
now 500


(478 missing values generated)
(478 missing values generated)


number of observations (_N)  was 500,
now 500


(477 missing values generated)
(477 missing values generated)


number of observations (_N)  was 500,
now 500


(483 missing values generated)
(483 missing values generated)


number of observations (_N)  was 500,
now 500


(484 missing values generated)
(484 missing values generated)


number of observations (_N)  was 500,
now 500


(484 missing values generated)
(484 missing values generated)


number of observations (_N)  was 500,
now 500


(481 missing values generated)
(481 missing values generated)


number of observations (_N)  was 500,
now 500


(481 missing values generated)
(481 missing values generated)


number of observations (_N)  was 500,
now 500


(482 missing values generated)
(482 missing values generated)


number of observations (_N)  was 500,
now 500


(486 missing values generated)
(486 missing values generated)


number of observations (_N)  was 500,
now 500


(482 missing values generated)
(482 missing values generated)


number of observations (_N)  was 500,
now 500


(480 missing values generated)
(480 missing values generated)


number of observations (_N)  was 500,
now 500


(488 missing values generated)
(488 missing values generated)


number of observations (_N)  was 500,
now 500


(477 missing values generated)
(477 missing values generated)


number of observations (_N)  was 500,
now 500


(486 missing values generated)
(486 missing values generated)


number of observations (_N)  was 500,
now 500


(489 missing values generated)
(489 missing values generated)


number of observations (_N)  was 500,
now 500


(480 missing values generated)
(480 missing values generated)


number of observations (_N)  was 500,
now 500


(487 missing values generated)
(487 missing values generated)


number of observations (_N)  was 500,
now 500


(482 missing values generated)
(482 missing values generated)


number of observations (_N)  was 500,
now 500


(482 missing values generated)
(482 missing values generated)


number of observations (_N)  was 500,
now 500


(479 missing values generated)
(479 missing values generated)


number of observations (_N)  was 500,
now 500


(475 missing values generated)
(475 missing values generated)


number of observations (_N)  was 500,
now 500


(483 missing values generated)
(483 missing values generated)


number of observations (_N)  was 500,
now 500


(474 missing values generated)
(474 missing values generated)


number of observations (_N)  was 500,
now 500


(486 missing values generated)
(486 missing values generated)


number of observations (_N)  was 500,
now 500


(479 missing values generated)
(479 missing values generated)


number of observations (_N)  was 500,
now 500


(482 missing values generated)
(482 missing values generated)


number of observations (_N)  was 500,
now 500


(481 missing values generated)
(481 missing values generated)


number of observations (_N)  was 500,
now 500


(479 missing values generated)
(479 missing values generated)


number of observations (_N)  was 500,
now 500


(480 missing values generated)
(480 missing values generated)
.
. foreach p1 in 0.2 1.1 ///
> {
  2.   foreach p2 in 1.0 0.1 ///
>   {
  3.     clear
  4.     use sim_`p1'_`p2'
  5.
.         display "S1 param = `p2', S2 param = `p1'"
  6.
.         * show total number of p values
.         summarize p
  7.
.         * show number of p values with p < 0.05
.         summarize p if p < 0.05
  8.   }
  9. }
S1 param = 1.0, S2 param = 0.2

    Variable |        Obs        Mean    Std. Dev.       Min        Max
-------------+---------------------------------------------------------
           p |        100    .1966748     .261059   .0000864   .9832942

    Variable |        Obs        Mean    Std. Dev.       Min        Max
-------------+---------------------------------------------------------
           p |         42    .0164047    .0158625   .0000864   .0474486
S1 param = 0.1, S2 param = 0.2

    Variable |        Obs        Mean    Std. Dev.       Min        Max
-------------+---------------------------------------------------------
           p |        100     .328459     .291638   .0004549   .9922009

    Variable |        Obs        Mean    Std. Dev.       Min        Max
-------------+---------------------------------------------------------
           p |         19    .0182594    .0130774   .0004549   .0458555
S1 param = 1.0, S2 param = 1.1

    Variable |        Obs        Mean    Std. Dev.       Min        Max
-------------+---------------------------------------------------------
           p |        100    .4915957    .3285587   .0013235   .9996558

    Variable |        Obs        Mean    Std. Dev.       Min        Max
-------------+---------------------------------------------------------
           p |         12    .0202482    .0146905   .0013235   .0381924
S1 param = 0.1, S2 param = 1.1

    Variable |        Obs        Mean    Std. Dev.       Min        Max
-------------+---------------------------------------------------------
           p |        100    .5382526    .2797091   .0219545   .9840519

    Variable |        Obs        Mean    Std. Dev.       Min        Max
-------------+---------------------------------------------------------
           p |          4    .0306168    .0081697   .0219545   .0408995
. }

.
.
end of do-file

. exit