{smcl}
{com}{sf}{ul off}{txt}{.-}
      name:  {res}<unnamed>
       {txt}log:  {res}/Users/stephenchang/Desktop/TICR Program/Spring 2017/EPI 265 Epidemiologic Methods III Research Methods in Chronic Disease Epidemiology/Lecture 8/HW8/assign8_stephenchang.smcl
  {txt}log type:  {res}smcl
 {txt}opened on:  {res}20 May 2017, 16:31:02
{txt}
{com}. 
. *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 {c -(}
. clear
. 
. set seed 98770
. set obs 500
{txt}{p}
number of observations (_N)  was 0,
now 500
{p_end}
{com}. 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 ///
> {c -(}
{txt}  2{com}.   foreach p2 in 1.0 0.1 ///
>   {c -(}
{txt}  3{com}.   
.     display "S1 param = `p2', S2 param = `p1'"
{txt}  4{com}. 
.     capture drop S1 S2 Xsel Ysel
{txt}  5{com}.   
.         gen S2 = runiform()>=`p1'*U
{txt}  6{com}.     gen S1 = runiform()<=`p2'*(X+Y)
{txt}  7{com}. 
.         gen Xsel = X if S1 == 1 & S2 == 1
{txt}  8{com}.         gen Ysel = Y if S1 == 1 & S2 == 1
{txt}  9{com}.         
.     * check & visualize data
.     summarize
{txt} 10{com}.     * scatter Xsel Ysel
.         
.         * simple linear regression of Y on X
.         regress Ysel Xsel
{txt} 11{com}.   {c )-}
{txt} 12{com}. {c )-}  
S1 param = 1.0, S2 param = 0.2
{txt}(189 missing values generated)
(189 missing values generated)

    Variable {c |}        Obs        Mean    Std. Dev.       Min        Max
{hline 13}{c +}{hline 57}
{space 3}errorTerm {c |}{res}        500   -.0069216    .5353013   -1.43531   1.540367
{txt}{space 11}U {c |}{res}        500    .4802203    .2819826   .0042005   .9938892
{txt}{space 11}X {c |}{res}        500    .6130487    .3009047   .0142714   1.172624
{txt}{space 11}Y {c |}{res}        500    .3053068    .6091797   -1.42284   2.335514
{txt}{space 10}S2 {c |}{res}        500        .894    .3081459          0          1
{txt}{hline 13}{c +}{hline 57}
{space 10}S1 {c |}{res}        500          .7    .4587165          0          1
{txt}{space 8}Xsel {c |}{res}        311    .6982614    .2802072   .0244896   1.172624
{txt}{space 8}Ysel {c |}{res}        311    .5316799    .5332607  -.7154276   2.335514

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

{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 1}        Ysel{col 14}{c |}      Coef.{col 26}   Std. Err.{col 38}      t{col 46}   P>|t|{col 54}     [95% Con{col 67}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 8}Xsel {c |}{col 14}{res}{space 2} .1693331{col 26}{space 2} .1078338{col 37}{space 1}    1.57{col 46}{space 3}0.117{col 54}{space 4}-.0428484{col 67}{space 3} .3815147
{txt}{space 7}_cons {c |}{col 14}{res}{space 2} .4134412{col 26}{space 2} .0811146{col 37}{space 1}    5.10{col 46}{space 3}0.000{col 54}{space 4} .2538343{col 67}{space 3}  .573048
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{res}S1 param = 0.1, S2 param = 0.2
{txt}(454 missing values generated)
(454 missing values generated)

    Variable {c |}        Obs        Mean    Std. Dev.       Min        Max
{hline 13}{c +}{hline 57}
{space 3}errorTerm {c |}{res}        500   -.0069216    .5353013   -1.43531   1.540367
{txt}{space 11}U {c |}{res}        500    .4802203    .2819826   .0042005   .9938892
{txt}{space 11}X {c |}{res}        500    .6130487    .3009047   .0142714   1.172624
{txt}{space 11}Y {c |}{res}        500    .3053068    .6091797   -1.42284   2.335514
{txt}{space 10}S2 {c |}{res}        500        .884    .3205456          0          1
{txt}{hline 13}{c +}{hline 57}
{space 10}S1 {c |}{res}        500        .104    .3055663          0          1
{txt}{space 8}Xsel {c |}{res}         46    .7491884    .2717626   .1120304   1.122617
{txt}{space 8}Ysel {c |}{res}         46    .7139529    .5350957  -.1785502   2.023463

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

{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 1}        Ysel{col 14}{c |}      Coef.{col 26}   Std. Err.{col 38}      t{col 46}   P>|t|{col 54}     [95% Con{col 67}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 8}Xsel {c |}{col 14}{res}{space 2}  .462005{col 26}{space 2} .2885482{col 37}{space 1}    1.60{col 46}{space 3}0.117{col 54}{space 4}-.1195256{col 67}{space 3} 1.043536
{txt}{space 7}_cons {c |}{col 14}{res}{space 2} .3678241{col 26}{space 2} .2296692{col 37}{space 1}    1.60{col 46}{space 3}0.116{col 54}{space 4}-.0950438{col 67}{space 3} .8306921
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{res}S1 param = 1.0, S2 param = 1.1
{txt}(344 missing values generated)
(344 missing values generated)

    Variable {c |}        Obs        Mean    Std. Dev.       Min        Max
{hline 13}{c +}{hline 57}
{space 3}errorTerm {c |}{res}        500   -.0069216    .5353013   -1.43531   1.540367
{txt}{space 11}U {c |}{res}        500    .4802203    .2819826   .0042005   .9938892
{txt}{space 11}X {c |}{res}        500    .6130487    .3009047   .0142714   1.172624
{txt}{space 11}Y {c |}{res}        500    .3053068    .6091797   -1.42284   2.335514
{txt}{space 10}S2 {c |}{res}        500        .506    .5004647          0          1
{txt}{hline 13}{c +}{hline 57}
{space 10}S1 {c |}{res}        500        .684    .4653785          0          1
{txt}{space 8}Xsel {c |}{res}        156    .6825249    .2661343   .0671131   1.098191
{txt}{space 8}Ysel {c |}{res}        156     .441816    .4958685  -.5124294   1.972102

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

{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 1}        Ysel{col 14}{c |}      Coef.{col 26}   Std. Err.{col 38}      t{col 46}   P>|t|{col 54}     [95% Con{col 67}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 8}Xsel {c |}{col 14}{res}{space 2}-.2610951{col 26}{space 2} .1486616{col 37}{space 1}   -1.76{col 46}{space 3}0.081{col 54}{space 4}-.5547743{col 67}{space 3} .0325841
{txt}{space 7}_cons {c |}{col 14}{res}{space 2} .6200199{col 26}{space 2} .1088599{col 37}{space 1}    5.70{col 46}{space 3}0.000{col 54}{space 4} .4049685{col 67}{space 3} .8350712
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{res}S1 param = 0.1, S2 param = 1.1
{txt}(488 missing values generated)
(488 missing values generated)

    Variable {c |}        Obs        Mean    Std. Dev.       Min        Max
{hline 13}{c +}{hline 57}
{space 3}errorTerm {c |}{res}        500   -.0069216    .5353013   -1.43531   1.540367
{txt}{space 11}U {c |}{res}        500    .4802203    .2819826   .0042005   .9938892
{txt}{space 11}X {c |}{res}        500    .6130487    .3009047   .0142714   1.172624
{txt}{space 11}Y {c |}{res}        500    .3053068    .6091797   -1.42284   2.335514
{txt}{space 10}S2 {c |}{res}        500        .462    .4990532          0          1
{txt}{hline 13}{c +}{hline 57}
{space 10}S1 {c |}{res}        500        .078    .2684402          0          1
{txt}{space 8}Xsel {c |}{res}         12    .7059129    .2996585   .1511894   1.054284
{txt}{space 8}Ysel {c |}{res}         12    .7253423    .6196588  -.2001904   1.972102

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

{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 1}        Ysel{col 14}{c |}      Coef.{col 26}   Std. Err.{col 38}      t{col 46}   P>|t|{col 54}     [95% Con{col 67}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 8}Xsel {c |}{col 14}{res}{space 2}-.2240837{col 26}{space 2} .6500714{col 37}{space 1}   -0.34{col 46}{space 3}0.737{col 54}{space 4}-1.672533{col 67}{space 3} 1.224366
{txt}{space 7}_cons {c |}{col 14}{res}{space 2} .8835259{col 26}{space 2} .4953464{col 37}{space 1}    1.78{col 46}{space 3}0.105{col 54}{space 4}-.2201747{col 67}{space 3} 1.987226
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{res}{com}. {c )-}
{txt}
{com}. 
. /* 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 {c -(}
. clear
. set obs 50000
{txt}{p}
number of observations (_N)  was 0,
now 50,000
{p_end}
{com}. 
. 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 ///
> {c -(}
{txt}  2{com}.   foreach p2 in 1.0 0.1 ///
>   {c -(}
{txt}  3{com}.   
.     display "S1 param = `p2', S2 param = `p1'"
{txt}  4{com}. 
.     capture drop S1 S2 Xsel Ysel
{txt}  5{com}.   
.         gen S2 = runiform()>=`p1'*U
{txt}  6{com}.     gen S1 = runiform()<=`p2'*(X+Y)
{txt}  7{com}. 
.         gen Xsel = X if S1 == 1 & S2 == 1
{txt}  8{com}.         gen Ysel = Y if S1 == 1 & S2 == 1
{txt}  9{com}.         
.         * simple linear regression of Y on X
.         bootstrap, reps(100) size(500): regress Ysel Xsel
{txt} 10{com}.   {c )-}
{txt} 11{com}. {c )-}
S1 param = 1.0, S2 param = 0.2
{txt}(18,997 missing values generated)
(18,997 missing values generated)
(running regress on estimation sample)

Bootstrap replications ({res}100{txt})
{hline 4}{c +}{hline 3} 1 {hline 3}{c +}{hline 3} 2 {hline 3}{c +}{hline 3} 3 {hline 3}{c +}{hline 3} 4 {hline 3}{c +}{hline 3} 5 
..................................................    50
..................................................   100
{res}
{txt}Linear regression{col 49}Number of obs{col 67}= {res}    31,003
{txt}{col 49}Replications{col 67}= {res}       100
{txt}{col 49}Wald chi2({res}1{txt}){col 67}= {res}      4.59
{txt}{col 49}Prob > chi2{col 67}= {res}    0.0322
{txt}{col 49}R-squared{col 67}= {res}    0.0099
{txt}{col 49}Adj R-squared{col 67}= {res}    0.0099
{txt}{col 49}Root MSE{col 67}= {res}    0.4577

{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 14}{c |}   Observed{col 26}   Bootstrap{col 54}         Norm{col 67}al-based
{col 1}        Ysel{col 14}{c |}      Coef.{col 26}   Std. Err.{col 38}      z{col 46}   P>|z|{col 54}     [95% Con{col 67}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 8}Xsel {c |}{col 14}{res}{space 2} .1680687{col 26}{space 2}  .078472{col 37}{space 1}    2.14{col 46}{space 3}0.032{col 54}{space 4} .0142664{col 67}{space 3}  .321871
{txt}{space 7}_cons {c |}{col 14}{res}{space 2} .4277872{col 26}{space 2} .0497253{col 37}{space 1}    8.60{col 46}{space 3}0.000{col 54}{space 4} .3303275{col 67}{space 3} .5252469
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
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 ({res}100{txt})
{hline 4}{c +}{hline 3} 1 {hline 3}{c +}{hline 3} 2 {hline 3}{c +}{hline 3} 3 {hline 3}{c +}{hline 3} 4 {hline 3}{c +}{hline 3} 5 
..................................................    50
..................................................   100
{res}
{txt}Linear regression{col 49}Number of obs{col 67}= {res}     4,290
{txt}{col 49}Replications{col 67}= {res}       100
{txt}{col 49}Wald chi2({res}1{txt}){col 67}= {res}     16.88
{txt}{col 49}Prob > chi2{col 67}= {res}    0.0000
{txt}{col 49}R-squared{col 67}= {res}    0.0353
{txt}{col 49}Adj R-squared{col 67}= {res}    0.0351
{txt}{col 49}Root MSE{col 67}= {res}    0.4942

{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 14}{c |}   Observed{col 26}   Bootstrap{col 54}         Norm{col 67}al-based
{col 1}        Ysel{col 14}{c |}      Coef.{col 26}   Std. Err.{col 38}      z{col 46}   P>|z|{col 54}     [95% Con{col 67}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 8}Xsel {c |}{col 14}{res}{space 2} .3485647{col 26}{space 2} .0848501{col 37}{space 1}    4.11{col 46}{space 3}0.000{col 54}{space 4} .1822615{col 67}{space 3} .5148679
{txt}{space 7}_cons {c |}{col 14}{res}{space 2} .4045199{col 26}{space 2} .0600283{col 37}{space 1}    6.74{col 46}{space 3}0.000{col 54}{space 4} .2868666{col 67}{space 3} .5221732
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
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 ({res}100{txt})
{hline 4}{c +}{hline 3} 1 {hline 3}{c +}{hline 3} 2 {hline 3}{c +}{hline 3} 3 {hline 3}{c +}{hline 3} 4 {hline 3}{c +}{hline 3} 5 
..................................................    50
..................................................   100
{res}
{txt}Linear regression{col 49}Number of obs{col 67}= {res}    14,087
{txt}{col 49}Replications{col 67}= {res}       100
{txt}{col 49}Wald chi2({res}1{txt}){col 67}= {res}      2.17
{txt}{col 49}Prob > chi2{col 67}= {res}    0.1410
{txt}{col 49}R-squared{col 67}= {res}    0.0035
{txt}{col 49}Adj R-squared{col 67}= {res}    0.0034
{txt}{col 49}Root MSE{col 67}= {res}    0.4254

{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 14}{c |}   Observed{col 26}   Bootstrap{col 54}         Norm{col 67}al-based
{col 1}        Ysel{col 14}{c |}      Coef.{col 26}   Std. Err.{col 38}      z{col 46}   P>|z|{col 54}     [95% Con{col 67}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 8}Xsel {c |}{col 14}{res}{space 2}-.0930817{col 26}{space 2} .0632368{col 37}{space 1}   -1.47{col 46}{space 3}0.141{col 54}{space 4}-.2170235{col 67}{space 3} .0308601
{txt}{space 7}_cons {c |}{col 14}{res}{space 2} .4974803{col 26}{space 2} .0399946{col 37}{space 1}   12.44{col 46}{space 3}0.000{col 54}{space 4} .4190923{col 67}{space 3} .5758683
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
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 ({res}100{txt})
{hline 4}{c +}{hline 3} 1 {hline 3}{c +}{hline 3} 2 {hline 3}{c +}{hline 3} 3 {hline 3}{c +}{hline 3} 4 {hline 3}{c +}{hline 3} 5 
..................................................    50
..................................................   100
{res}
{txt}Linear regression{col 49}Number of obs{col 67}= {res}     1,773
{txt}{col 49}Replications{col 67}= {res}       100
{txt}{col 49}Wald chi2({res}1{txt}){col 67}= {res}      0.09
{txt}{col 49}Prob > chi2{col 67}= {res}    0.7601
{txt}{col 49}R-squared{col 67}= {res}    0.0002
{txt}{col 49}Adj R-squared{col 67}= {res}   -0.0004
{txt}{col 49}Root MSE{col 67}= {res}    0.4484

{txt}{hline 13}{c TT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{col 14}{c |}   Observed{col 26}   Bootstrap{col 54}         Norm{col 67}al-based
{col 1}        Ysel{col 14}{c |}      Coef.{col 26}   Std. Err.{col 38}      z{col 46}   P>|z|{col 54}     [95% Con{col 67}f. Interval]
{hline 13}{c +}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{space 8}Xsel {c |}{col 14}{res}{space 2} .0221603{col 26}{space 2} .0725615{col 37}{space 1}    0.31{col 46}{space 3}0.760{col 54}{space 4}-.1200577{col 67}{space 3} .1643783
{txt}{space 7}_cons {c |}{col 14}{res}{space 2} .5035067{col 26}{space 2} .0501339{col 37}{space 1}   10.04{col 46}{space 3}0.000{col 54}{space 4} .4052461{col 67}{space 3} .6017674
{txt}{hline 13}{c BT}{hline 11}{hline 11}{hline 9}{hline 8}{hline 13}{hline 12}
{com}. {c )-}
{txt}
{com}.   
. /* 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 {c -(}   
. clear
. set seed 98770
. foreach p1 in 0.2 1.1 ///
> {c -(}
{txt}  2{com}.   foreach p2 in 1.0 0.1 ///
>   {c -(}
{txt}  3{com}. 
.     display "S1 param = `p2', S2 param = `p1'"
{txt}  4{com}.   
.         postfile results p using sim_`p1'_`p2', replace
{txt}  5{com}.         forvalues n=1/100 {c -(}
{txt}  6{com}.            set obs 500
{txt}  7{com}.            
.            capture drop errorTerm U X Y
{txt}  8{com}.        gen errorTerm = rnormal(0, 0.5)
{txt}  9{com}.        gen U = runiform(0, 1.)
{txt} 10{com}.        gen X = runiform(0, 1) + 0.2 * U
{txt} 11{com}.        gen Y = 0 + 0.5*X + 1.*X*(U-0.5) + errorTerm
{txt} 12{com}.   
.        capture drop S1 S2 Xsel Ysel
{txt} 13{com}.            gen S2 = runiform()>=`p1'*U
{txt} 14{com}.        gen S1 = runiform()<=`p2'*(X+Y)
{txt} 15{com}. 
.            gen Xsel = X if S1 == 1 & S2 == 1
{txt} 16{com}.            gen Ysel = Y if S1 == 1 & S2 == 1
{txt} 17{com}.         
.            quietly regress Ysel Xsel
{txt} 18{com}.            mat b = r(table)
{txt} 19{com}.            local p = b[4,1]
{txt} 20{com}.            post results (`p')
{txt} 21{com}.         {c )-}
{txt} 22{com}.         postclose results
{txt} 23{com}.   {c )-}
{txt} 24{com}. {c )-}
S1 param = 1.0, S2 param = 0.2
{txt}(note: file sim_0.2_1.0 not found)
{p}
number of observations (_N)  was 0,
now 500
{p_end}
(189 missing values generated)
(189 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(198 missing values generated)
(198 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(200 missing values generated)
(200 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(198 missing values generated)
(198 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(204 missing values generated)
(204 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(183 missing values generated)
(183 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(195 missing values generated)
(195 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(180 missing values generated)
(180 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(191 missing values generated)
(191 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(202 missing values generated)
(202 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(181 missing values generated)
(181 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(191 missing values generated)
(191 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(185 missing values generated)
(185 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(197 missing values generated)
(197 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(182 missing values generated)
(182 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(190 missing values generated)
(190 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(188 missing values generated)
(188 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(190 missing values generated)
(190 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(186 missing values generated)
(186 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(213 missing values generated)
(213 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(192 missing values generated)
(192 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(187 missing values generated)
(187 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(189 missing values generated)
(189 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(177 missing values generated)
(177 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(185 missing values generated)
(185 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(192 missing values generated)
(192 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(198 missing values generated)
(198 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(181 missing values generated)
(181 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(198 missing values generated)
(198 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(195 missing values generated)
(195 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(192 missing values generated)
(192 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(178 missing values generated)
(178 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(192 missing values generated)
(192 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(194 missing values generated)
(194 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(192 missing values generated)
(192 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(185 missing values generated)
(185 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(184 missing values generated)
(184 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(190 missing values generated)
(190 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(196 missing values generated)
(196 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(190 missing values generated)
(190 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(180 missing values generated)
(180 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(187 missing values generated)
(187 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(192 missing values generated)
(192 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(199 missing values generated)
(199 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(174 missing values generated)
(174 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(191 missing values generated)
(191 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(194 missing values generated)
(194 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(189 missing values generated)
(189 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(190 missing values generated)
(190 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(186 missing values generated)
(186 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(182 missing values generated)
(182 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(194 missing values generated)
(194 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(195 missing values generated)
(195 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(203 missing values generated)
(203 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(203 missing values generated)
(203 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(194 missing values generated)
(194 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(196 missing values generated)
(196 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(208 missing values generated)
(208 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(204 missing values generated)
(204 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(187 missing values generated)
(187 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(200 missing values generated)
(200 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(188 missing values generated)
(188 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(190 missing values generated)
(190 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(204 missing values generated)
(204 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(202 missing values generated)
(202 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(181 missing values generated)
(181 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(196 missing values generated)
(196 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(185 missing values generated)
(185 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(194 missing values generated)
(194 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(208 missing values generated)
(208 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(180 missing values generated)
(180 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(208 missing values generated)
(208 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(189 missing values generated)
(189 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(170 missing values generated)
(170 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(201 missing values generated)
(201 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(185 missing values generated)
(185 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(180 missing values generated)
(180 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(193 missing values generated)
(193 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(174 missing values generated)
(174 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(173 missing values generated)
(173 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(183 missing values generated)
(183 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(189 missing values generated)
(189 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(158 missing values generated)
(158 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(179 missing values generated)
(179 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(207 missing values generated)
(207 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(178 missing values generated)
(178 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(182 missing values generated)
(182 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(189 missing values generated)
(189 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(189 missing values generated)
(189 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(193 missing values generated)
(193 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(183 missing values generated)
(183 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(171 missing values generated)
(171 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(187 missing values generated)
(187 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(200 missing values generated)
(200 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(184 missing values generated)
(184 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(205 missing values generated)
(205 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(201 missing values generated)
(201 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(172 missing values generated)
(172 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(175 missing values generated)
(175 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(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)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(464 missing values generated)
(464 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(466 missing values generated)
(466 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(459 missing values generated)
(459 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(451 missing values generated)
(451 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(454 missing values generated)
(454 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(458 missing values generated)
(458 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(466 missing values generated)
(466 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(457 missing values generated)
(457 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(459 missing values generated)
(459 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(465 missing values generated)
(465 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(458 missing values generated)
(458 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(464 missing values generated)
(464 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(459 missing values generated)
(459 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(460 missing values generated)
(460 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(457 missing values generated)
(457 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(455 missing values generated)
(455 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(461 missing values generated)
(461 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(459 missing values generated)
(459 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(445 missing values generated)
(445 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(465 missing values generated)
(465 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(457 missing values generated)
(457 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(455 missing values generated)
(455 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(461 missing values generated)
(461 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(461 missing values generated)
(461 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(464 missing values generated)
(464 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(453 missing values generated)
(453 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(457 missing values generated)
(457 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(454 missing values generated)
(454 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(459 missing values generated)
(459 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(454 missing values generated)
(454 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(443 missing values generated)
(443 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(461 missing values generated)
(461 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(455 missing values generated)
(455 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(462 missing values generated)
(462 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(446 missing values generated)
(446 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(461 missing values generated)
(461 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(465 missing values generated)
(465 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(452 missing values generated)
(452 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(457 missing values generated)
(457 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(457 missing values generated)
(457 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(467 missing values generated)
(467 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(451 missing values generated)
(451 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(462 missing values generated)
(462 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(453 missing values generated)
(453 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(452 missing values generated)
(452 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(467 missing values generated)
(467 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(463 missing values generated)
(463 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(451 missing values generated)
(451 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(462 missing values generated)
(462 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(457 missing values generated)
(457 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(453 missing values generated)
(453 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(459 missing values generated)
(459 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(454 missing values generated)
(454 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(450 missing values generated)
(450 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(453 missing values generated)
(453 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(451 missing values generated)
(451 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(458 missing values generated)
(458 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(454 missing values generated)
(454 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(457 missing values generated)
(457 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(456 missing values generated)
(456 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(453 missing values generated)
(453 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(458 missing values generated)
(458 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(468 missing values generated)
(468 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(453 missing values generated)
(453 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(459 missing values generated)
(459 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(456 missing values generated)
(456 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(458 missing values generated)
(458 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(452 missing values generated)
(452 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(452 missing values generated)
(452 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(458 missing values generated)
(458 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(464 missing values generated)
(464 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(461 missing values generated)
(461 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(457 missing values generated)
(457 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(461 missing values generated)
(461 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(453 missing values generated)
(453 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(450 missing values generated)
(450 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(460 missing values generated)
(460 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(469 missing values generated)
(469 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(460 missing values generated)
(460 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(464 missing values generated)
(464 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(453 missing values generated)
(453 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(467 missing values generated)
(467 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(461 missing values generated)
(461 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(463 missing values generated)
(463 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(444 missing values generated)
(444 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(456 missing values generated)
(456 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(457 missing values generated)
(457 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(460 missing values generated)
(460 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(449 missing values generated)
(449 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(456 missing values generated)
(456 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(455 missing values generated)
(455 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(461 missing values generated)
(461 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(454 missing values generated)
(454 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(463 missing values generated)
(463 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(463 missing values generated)
(463 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(450 missing values generated)
(450 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(450 missing values generated)
(450 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(454 missing values generated)
(454 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(450 missing values generated)
(450 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(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)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(360 missing values generated)
(360 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(354 missing values generated)
(354 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(373 missing values generated)
(373 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(346 missing values generated)
(346 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(364 missing values generated)
(364 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(345 missing values generated)
(345 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(372 missing values generated)
(372 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(350 missing values generated)
(350 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(368 missing values generated)
(368 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(359 missing values generated)
(359 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(337 missing values generated)
(337 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(355 missing values generated)
(355 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(345 missing values generated)
(345 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(374 missing values generated)
(374 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(350 missing values generated)
(350 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(363 missing values generated)
(363 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(355 missing values generated)
(355 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(363 missing values generated)
(363 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(360 missing values generated)
(360 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(363 missing values generated)
(363 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(351 missing values generated)
(351 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(359 missing values generated)
(359 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(365 missing values generated)
(365 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(356 missing values generated)
(356 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(365 missing values generated)
(365 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(368 missing values generated)
(368 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(362 missing values generated)
(362 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(349 missing values generated)
(349 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(366 missing values generated)
(366 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(346 missing values generated)
(346 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(347 missing values generated)
(347 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(359 missing values generated)
(359 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(364 missing values generated)
(364 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(356 missing values generated)
(356 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(359 missing values generated)
(359 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(363 missing values generated)
(363 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(356 missing values generated)
(356 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(353 missing values generated)
(353 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(383 missing values generated)
(383 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(359 missing values generated)
(359 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(356 missing values generated)
(356 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(368 missing values generated)
(368 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(360 missing values generated)
(360 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(366 missing values generated)
(366 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(362 missing values generated)
(362 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(363 missing values generated)
(363 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(350 missing values generated)
(350 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(349 missing values generated)
(349 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(353 missing values generated)
(353 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(357 missing values generated)
(357 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(344 missing values generated)
(344 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(364 missing values generated)
(364 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(365 missing values generated)
(365 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(349 missing values generated)
(349 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(360 missing values generated)
(360 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(368 missing values generated)
(368 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(340 missing values generated)
(340 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(363 missing values generated)
(363 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(357 missing values generated)
(357 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(355 missing values generated)
(355 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(353 missing values generated)
(353 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(356 missing values generated)
(356 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(353 missing values generated)
(353 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(344 missing values generated)
(344 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(362 missing values generated)
(362 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(363 missing values generated)
(363 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(361 missing values generated)
(361 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(351 missing values generated)
(351 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(333 missing values generated)
(333 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(342 missing values generated)
(342 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(343 missing values generated)
(343 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(358 missing values generated)
(358 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(364 missing values generated)
(364 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(354 missing values generated)
(354 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(364 missing values generated)
(364 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(340 missing values generated)
(340 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(345 missing values generated)
(345 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(371 missing values generated)
(371 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(373 missing values generated)
(373 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(367 missing values generated)
(367 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(358 missing values generated)
(358 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(362 missing values generated)
(362 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(361 missing values generated)
(361 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(363 missing values generated)
(363 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(343 missing values generated)
(343 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(355 missing values generated)
(355 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(350 missing values generated)
(350 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(359 missing values generated)
(359 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(369 missing values generated)
(369 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(372 missing values generated)
(372 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(367 missing values generated)
(367 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(364 missing values generated)
(364 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(362 missing values generated)
(362 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(354 missing values generated)
(354 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(364 missing values generated)
(364 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(356 missing values generated)
(356 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(350 missing values generated)
(350 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(358 missing values generated)
(358 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(360 missing values generated)
(360 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(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)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(483 missing values generated)
(483 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(476 missing values generated)
(476 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(489 missing values generated)
(489 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(480 missing values generated)
(480 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(485 missing values generated)
(485 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(485 missing values generated)
(485 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(488 missing values generated)
(488 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(481 missing values generated)
(481 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(485 missing values generated)
(485 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(479 missing values generated)
(479 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(489 missing values generated)
(489 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(477 missing values generated)
(477 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(478 missing values generated)
(478 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(486 missing values generated)
(486 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(487 missing values generated)
(487 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(488 missing values generated)
(488 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(470 missing values generated)
(470 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(479 missing values generated)
(479 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(483 missing values generated)
(483 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(474 missing values generated)
(474 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(485 missing values generated)
(485 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(484 missing values generated)
(484 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(479 missing values generated)
(479 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(484 missing values generated)
(484 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(483 missing values generated)
(483 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(487 missing values generated)
(487 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(483 missing values generated)
(483 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(478 missing values generated)
(478 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(481 missing values generated)
(481 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(486 missing values generated)
(486 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(481 missing values generated)
(481 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(482 missing values generated)
(482 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(487 missing values generated)
(487 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(481 missing values generated)
(481 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(483 missing values generated)
(483 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(479 missing values generated)
(479 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(483 missing values generated)
(483 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(481 missing values generated)
(481 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(487 missing values generated)
(487 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(483 missing values generated)
(483 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(481 missing values generated)
(481 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(492 missing values generated)
(492 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(487 missing values generated)
(487 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(479 missing values generated)
(479 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(485 missing values generated)
(485 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(480 missing values generated)
(480 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(484 missing values generated)
(484 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(486 missing values generated)
(486 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(490 missing values generated)
(490 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(482 missing values generated)
(482 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(485 missing values generated)
(485 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(486 missing values generated)
(486 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(471 missing values generated)
(471 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(488 missing values generated)
(488 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(491 missing values generated)
(491 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(485 missing values generated)
(485 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(484 missing values generated)
(484 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(474 missing values generated)
(474 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(479 missing values generated)
(479 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(478 missing values generated)
(478 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(481 missing values generated)
(481 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(473 missing values generated)
(473 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(482 missing values generated)
(482 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(480 missing values generated)
(480 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(484 missing values generated)
(484 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(480 missing values generated)
(480 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(484 missing values generated)
(484 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(484 missing values generated)
(484 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(482 missing values generated)
(482 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(484 missing values generated)
(484 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(482 missing values generated)
(482 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(478 missing values generated)
(478 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(477 missing values generated)
(477 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(483 missing values generated)
(483 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(484 missing values generated)
(484 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(484 missing values generated)
(484 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(481 missing values generated)
(481 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(481 missing values generated)
(481 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(482 missing values generated)
(482 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(486 missing values generated)
(486 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(482 missing values generated)
(482 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(480 missing values generated)
(480 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(488 missing values generated)
(488 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(477 missing values generated)
(477 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(486 missing values generated)
(486 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(489 missing values generated)
(489 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(480 missing values generated)
(480 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(487 missing values generated)
(487 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(482 missing values generated)
(482 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(482 missing values generated)
(482 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(479 missing values generated)
(479 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(475 missing values generated)
(475 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(483 missing values generated)
(483 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(474 missing values generated)
(474 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(486 missing values generated)
(486 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(479 missing values generated)
(479 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(482 missing values generated)
(482 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(481 missing values generated)
(481 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(479 missing values generated)
(479 missing values generated)
{p}
number of observations (_N)  was 500,
now 500
{p_end}
(480 missing values generated)
(480 missing values generated)
{com}. 
. foreach p1 in 0.2 1.1 ///
> {c -(}
{txt}  2{com}.   foreach p2 in 1.0 0.1 ///
>   {c -(}
{txt}  3{com}.     clear
{txt}  4{com}.     use sim_`p1'_`p2'
{txt}  5{com}.         
.         display "S1 param = `p2', S2 param = `p1'"
{txt}  6{com}.         
.         * show total number of p values
.         summarize p
{txt}  7{com}.         
.         * show number of p values with p < 0.05
.         summarize p if p < 0.05
{txt}  8{com}.   {c )-}
{txt}  9{com}. {c )-}
S1 param = 1.0, S2 param = 0.2

{txt}    Variable {c |}        Obs        Mean    Std. Dev.       Min        Max
{hline 13}{c +}{hline 57}
{space 11}p {c |}{res}        100    .1966748     .261059   .0000864   .9832942

{txt}    Variable {c |}        Obs        Mean    Std. Dev.       Min        Max
{hline 13}{c +}{hline 57}
{space 11}p {c |}{res}         42    .0164047    .0158625   .0000864   .0474486
S1 param = 0.1, S2 param = 0.2

{txt}    Variable {c |}        Obs        Mean    Std. Dev.       Min        Max
{hline 13}{c +}{hline 57}
{space 11}p {c |}{res}        100     .328459     .291638   .0004549   .9922009

{txt}    Variable {c |}        Obs        Mean    Std. Dev.       Min        Max
{hline 13}{c +}{hline 57}
{space 11}p {c |}{res}         19    .0182594    .0130774   .0004549   .0458555
S1 param = 1.0, S2 param = 1.1

{txt}    Variable {c |}        Obs        Mean    Std. Dev.       Min        Max
{hline 13}{c +}{hline 57}
{space 11}p {c |}{res}        100    .4915957    .3285587   .0013235   .9996558

{txt}    Variable {c |}        Obs        Mean    Std. Dev.       Min        Max
{hline 13}{c +}{hline 57}
{space 11}p {c |}{res}         12    .0202482    .0146905   .0013235   .0381924
S1 param = 0.1, S2 param = 1.1

{txt}    Variable {c |}        Obs        Mean    Std. Dev.       Min        Max
{hline 13}{c +}{hline 57}
{space 11}p {c |}{res}        100    .5382526    .2797091   .0219545   .9840519

{txt}    Variable {c |}        Obs        Mean    Std. Dev.       Min        Max
{hline 13}{c +}{hline 57}
{space 11}p {c |}{res}          4    .0306168    .0081697   .0219545   .0408995
{com}. {c )-}
{txt}
{com}. 
. 
{txt}end of do-file

{com}. exit