Content
-------------------------------------------------------------------------------
name:
log: /home/users/dlmiller.AD3/colinjhr/work2014/table2.log
log type: text
opened on: 26 Jun 2014, 16:55:37
. clear ;
. set seed 10101 ;
. set more off ;
. /* load up the programs. each do file has a program that will be used */
>
> do set_directory_macros ;
. #delimit ;
delimiter now ;
. local my_top_path = subinstr(c(pwd), "doug work\Colin\JHR cluster paper\work
> Summer 2014\empirical example code", "",.) ;
. global rawdat = "`my_top_path'data\CPS\July 2013 extract for cluster JHR" ;
. global madedat = "`my_top_path'madedata\Colin JHR paper" ;
. if "$S_OS" == "Unix" { ;
. global rawdat = "../data" ;
. global madedat = "../data" ;
. } ;
. di "$rawdat" ;
../data
. di "$madedat" ;
../data
. *global rawdat = "C:\data\CPS\July 2013 extract for cluster JHR" ;
. *global madedat = "C:\madedata\Colin JHR paper" ;
. *global rawdat = "C:\Users\dlmiller\data\CPS\July 2013 extract for cluster JH
> R" ;
. *global madedat = "C:\Users\dlmiller\madedata\Colin JHR paper" ;
.
end of do-file
. set matsize 10000 ;
. run t2_mcloop ;
. run CR23_IK_CSS ;
cr2cr3() not found
. run t2_display ;
. cap prog drop runme ;
. prog def runme ;
1. /* initialize the file names, lists, etc, that will be passed into the pgr
> ograms */
>
> local num_mc_reps = 100 ;
2. local mc_reps_smallG = 4000 ;
3. local mc_reps_bigG = 1000 ;
4. local num_bs_reps = 399 ;
5. *local num_states_list = "6 " ;
. local num_states_list = "6 10 20 30 50" ;
6. *local mysourcedata = "$madedat/CPS_2012_micro_small" ;
. local mysourcedata = "$madedat/CPS_2012_micro_medium" ;
7. /* ensure that the source data file is sorted on statefip */
> tempfile statelist ;
8. qui use "`mysourcedata'" ;
9. cap drop resid ;
10. sort statefip ;
11. qui save "`mysourcedata'" , replace ;
12. /* now, do the Monte Carlo exercise*/
> /* loop over number of states in the hum_states_list*/
>
> foreach ns in `num_states_list' { ;
13. /* do more MC reps for 6 and 10 states */
> local num_mc_reps = `mc_reps_bigG' ;
14. if `ns' < 11 { ;
15. local num_mc_reps = `mc_reps_smallG' ;
16. } ;
17. tempfile mcout_`ns' ;
18. di ;
19. di "$S_DATE $S_TIME" ;
20. do_the_monte_carlo , sourcedata("`mysourcedata'") savedata("`mcout
> _`ns''")
> mcreps(`num_mc_reps') bsreps(`num_bs_reps') numstates(`ns') ;
21. use `mcout_`ns'' , replace ;
22. qui gen datetime = "$S_DATE $S_TIME" ;
23. /* keeps track of when this batch of data was generated */
> qui save "`mcout_`ns''" , replace ;
24. qui save t2_mcout_G_`ns' , replace ;
25. } ;
26. /* this next loop will go through the "statelist" and swap-in the results
> into to_mcout_all
> this is to enable these results to be done in batches, in parallell.
> We just swap in those results in the list,
> leaving any results already in the file (and not in num_states_list)
> untouched. */
> tempfile toappend ;
27. foreach ns in `num_states_list' { ;
28. clear ;
29. cap erase `toappend' ;
30. cap use t2_mcout_all ;
31. cap drop if numstates == `ns' ;
32. cap save `toappend' , replace ;
33. use t2_mcout_G_`ns' ;
34. cap append using `toappend' ;
35. qui save t2_mcout_all , replace ;
36. } ;
37. display_results , results_data(t2_mcout_all)
> numstates_list("`num_states_list'") sourcedata("`mysourcedata'");
38. end ;
. runme ;
26 Jun 2014 16:55:37
MC Reps = 4000, BS Reps = 399, Number of states sampled = 6
Using data from ../data/CPS_2012_micro_medium
. 1 ...........................................................................
> ......................... 101 ...............................................
> ..................................................... 201 ...................
> .............................................................................
> .... 301 ....................................................................
> ................................ 401 ........................................
> ............................................................ 501 ............
> .............................................................................
> ........... 601 .............................................................
> ....................................... 701 .................................
> ................................................................... 801 .....
> .............................................................................
> .................. 901 ......................................................
> .............................................. 1001 .........................
> ........................................................................... 1
> 101 .........................................................................
> ........................... 1201 ............................................
> ........................................................ 1301 ...............
> .............................................................................
> ........ 1401 ...............................................................
> ..................................... 1501 ..................................
> .................................................................. 1601 .....
> .............................................................................
> .................. 1701 .....................................................
> ............................................... 1801 ........................
> ............................................................................
> 1901 ........................................................................
> ............................ 2001 ...........................................
> ......................................................... 2101 ..............
> .............................................................................
> ......... 2201 ..............................................................
> ...................................... 2301 .................................
> ................................................................... 2401 ....
> .............................................................................
> ................... 2501 ....................................................
> ................................................ 2601 .......................
> .............................................................................
> 2701 .......................................................................
> ............................. 2801 ..........................................
> .......................................................... 2901 .............
> .............................................................................
> .......... 3001 .............................................................
> ....................................... 3101 ................................
> .................................................................... 3201 ...
> .............................................................................
> .................... 3301 ...................................................
> ................................................. 3401 ......................
> .............................................................................
> . 3501 ......................................................................
> .............................. 3601 .........................................
> ........................................................... 3701 ............
> .............................................................................
> ........... 3801 ............................................................
> ........................................ 3901 ...............................
> ....................................................................
29 Jun 2014 21:26:57
MC Reps = 4000, BS Reps = 399, Number of states sampled = 10
Using data from ../data/CPS_2012_micro_medium
. 1 ...........................................................................
> ......................... 101 ...............................................
> ..................................................... 201 ...................
> .............................................................................
> .... 301 ....................................................................
> ................................ 401 ........................................
> ............................................................ 501 ............
> .............................................................................
> ........... 601 .............................................................
> ....................................... 701 .................................
> ................................................................... 801 .....
> .............................................................................
> .................. 901 ......................................................
> .............................................. 1001 .........................
> ........................................................................... 1
> 101 .........................................................................
> ........................... 1201 ............................................
> ........................................................ 1301 ...............
> .............................................................................
> ........ 1401 ...............................................................
> ..................................... 1501 ..................................
> .................................................................. 1601 .....
> .............................................................................
> .................. 1701 .....................................................
> ............................................... 1801 ........................
> ............................................................................
> 1901 ........................................................................
> ............................ 2001 ...........................................
> ......................................................... 2101 ..............
> .............................................................................
> ......... 2201 ..............................................................
> ...................................... 2301 .................................
> ................................................................... 2401 ....
> .............................................................................
> ................... 2501 ....................................................
> ................................................ 2601 .......................
> .............................................................................
> 2701 .......................................................................
> ............................. 2801 ..........................................
> .......................................................... 2901 .............
> .............................................................................
> .......... 3001 .............................................................
> ....................................... 3101 ................................
> .................................................................... 3201 ...
> .............................................................................
> .................... 3301 ...................................................
> ................................................. 3401 ......................
> .............................................................................
> . 3501 ......................................................................
> .............................. 3601 .........................................
> ........................................................... 3701 ............
> .............................................................................
> ........... 3801 ............................................................
> ........................................ 3901 ...............................
> ....................................................................
2 Jul 2014 12:02:35
MC Reps = 1000, BS Reps = 399, Number of states sampled = 20
Using data from ../data/CPS_2012_micro_medium
. 1 ...........................................................................
> ......................... 101 ...............................................
> ..................................................... 201 ...................
> .............................................................................
> .... 301 ....................................................................
> ................................ 401 ........................................
> ............................................................ 501 ............
> .............................................................................
> ........... 601 .............................................................
> ....................................... 701 .................................
> ................................................................... 801 .....
> .............................................................................
> .................. 901 ......................................................
> .............................................
3 Jul 2014 02:13:05
MC Reps = 1000, BS Reps = 399, Number of states sampled = 30
Using data from ../data/CPS_2012_micro_medium
. 1 ...........................................................................
> ......................... 101 ...............................................
> ..................................................... 201 ...................
> .............................................................................
> .... 301 ....................................................................
> ................................ 401 ........................................
> ............................................................ 501 ............
> .............................................................................
> ........... 601 .............................................................
> ....................................... 701 .................................
> ................................................................... 801 .....
> .............................................................................
> .................. 901 ......................................................
> .............................................
4 Jul 2014 01:45:35
MC Reps = 1000, BS Reps = 399, Number of states sampled = 50
Using data from ../data/CPS_2012_micro_medium
. 1 ...........................................................................
> ......................... 101 ...............................................
> ..................................................... 201 ...................
> .............................................................................
> .... 301 ....................................................................
> ................................ 401 ........................................
> ............................................................ 501 ............
> .............................................................................
> ........... 601 .............................................................
> ....................................... 701 .................................
> ................................................................... 801 .....
> .............................................................................
> .................. 901 ......................................................
> .............................................
Num Clusters, 6
Rej. rate rob_Tnk, 0.4390
Rej. rate clu_Tnk, 0.2147
Rej. rate clu_Tg1, 0.1253
Rej. rate clu_Tg2, 0.1052
Rej. rate CR2_Tg1, 0.0823
Rej. rate CR3_Tg1, 0.0483
Rej. rate CR2_Tik, 0.0515
Rej. rate clu_Tcss, 0.1143
Rej. rate bsse_Tg1, 0.0820
Rej. rate pairs_bs, 0.0088
Rej. rate rad_res_low, 0.0973
Rej. rate rad_res, 0.0683
Rej. rate rad_res_high, 0.0413
Rej. rate webb_res, 0.0788
Rej. rate rad_unres, 0.0855
I-K dof (mean), 3.3
I-K dof (p5), 2.7
I-K dof (p95), 3.8
CSS dof (mean), 4.7
Mean num. observations, 1554
Number MC replications, 4000
Number BS replications, 399
Avg miss prs BS reps, 12
Avg Rad Res p-avg, 0.4514
Avg Rad Res p-high, 0.4671
Avg Rad Res p-low, 0.4356
Created: 29 Jun 2014 21:26:57
Using: ../data/CPS_2012_micro_medium
-----------------------------------------
Num Clusters, 10
Rej. rate rob_Tnk, 0.4570
Rej. rate clu_Tnk, 0.1467
Rej. rate clu_Tg1, 0.1025
Rej. rate clu_Tg2, 0.0988
Rej. rate CR2_Tg1, 0.0698
Rej. rate CR3_Tg1, 0.0500
Rej. rate CR2_Tik, 0.0495
Rej. rate clu_Tcss, 0.0785
Rej. rate bsse_Tg1, 0.0722
Rej. rate pairs_bs, 0.0312
Rej. rate rad_res_low, 0.0653
Rej. rate rad_res, 0.0653
Rej. rate rad_res_high, 0.0643
Rej. rate webb_res, 0.0668
Rej. rate rad_unres, 0.0625
I-K dof (mean), 5.6
I-K dof (p5), 3.7
I-K dof (p95), 7.2
CSS dof (mean), 6.6
Mean num. observations, 2618
Number MC replications, 4000
Number BS replications, 399
Avg miss prs BS reps, 1
Avg Rad Res p-avg, 0.4833
Avg Rad Res p-high, 0.4844
Avg Rad Res p-low, 0.4823
Created: 2 Jul 2014 12:02:35
Using: ../data/CPS_2012_micro_medium
-----------------------------------------
Num Clusters, 20
Rej. rate rob_Tnk, 0.4710
Rej. rate clu_Tnk, 0.1040
Rej. rate clu_Tg1, 0.0820
Rej. rate clu_Tg2, 0.0760
Rej. rate CR2_Tg1, 0.0620
Rej. rate CR3_Tg1, 0.0500
Rej. rate CR2_Tik, 0.0470
Rej. rate clu_Tcss, 0.0570
Rej. rate bsse_Tg1, 0.0690
Rej. rate pairs_bs, 0.0460
Rej. rate rad_res_low, 0.0620
Rej. rate rad_res, 0.0620
Rej. rate rad_res_high, 0.0620
Rej. rate webb_res, 0.0610
Rej. rate rad_unres, 0.0500
I-K dof (mean), 9.4
I-K dof (p5), 4.9
I-K dof (p95), 14.5
CSS dof (mean), 9.9
Mean num. observations, 5210
Number MC replications, 1000
Number BS replications, 399
Avg miss prs BS reps, 0
Avg Rad Res p-avg, 0.4923
Avg Rad Res p-high, 0.4923
Avg Rad Res p-low, 0.4922
Created: 3 Jul 2014 02:13:05
Using: ../data/CPS_2012_micro_medium
-----------------------------------------
Num Clusters, 30
Rej. rate rob_Tnk, 0.4620
Rej. rate clu_Tnk, 0.0830
Rej. rate clu_Tg1, 0.0690
Rej. rate clu_Tg2, 0.0690
Rej. rate CR2_Tg1, 0.0600
Rej. rate CR3_Tg1, 0.0520
Rej. rate CR2_Tik, 0.0470
Rej. rate clu_Tcss, 0.0560
Rej. rate bsse_Tg1, 0.0670
Rej. rate pairs_bs, 0.0510
Rej. rate rad_res_low, 0.0510
Rej. rate rad_res, 0.0510
Rej. rate rad_res_high, 0.0510
Rej. rate webb_res, 0.0510
Rej. rate rad_unres, 0.0530
I-K dof (mean), 12.3
I-K dof (p5), 6.3
I-K dof (p95), 20.8
CSS dof (mean), 12.7
Mean num. observations, 7803
Number MC replications, 1000
Number BS replications, 399
Avg miss prs BS reps, 0
Avg Rad Res p-avg, 0.5034
Avg Rad Res p-high, 0.5035
Avg Rad Res p-low, 0.5034
Created: 4 Jul 2014 01:45:35
Using: ../data/CPS_2012_micro_medium
-----------------------------------------
Num Clusters, 50
Rej. rate rob_Tnk, 0.4980
Rej. rate clu_Tnk, 0.0780
Rej. rate clu_Tg1, 0.0750
Rej. rate clu_Tg2, 0.0750
Rej. rate CR2_Tg1, 0.0650
Rej. rate CR3_Tg1, 0.0610
Rej. rate CR2_Tik, 0.0540
Rej. rate clu_Tcss, 0.0610
Rej. rate bsse_Tg1, 0.0740
Rej. rate pairs_bs, 0.0610
Rej. rate rad_res_low, 0.0600
Rej. rate rad_res, 0.0600
Rej. rate rad_res_high, 0.0600
Rej. rate webb_res, 0.0610
Rej. rate rad_unres, 0.0560
I-K dof (mean), 16.9
I-K dof (p5), 9.6
I-K dof (p95), 29.5
CSS dof (mean), 17.0
Mean num. observations, 13055
Number MC replications, 1000
Number BS replications, 399
Avg miss prs BS reps, 0
Avg Rad Res p-avg, 0.4957
Avg Rad Res p-high, 0.4957
Avg Rad Res p-low, 0.4956
Created: 6 Jul 2014 01:53:55
Using: ../data/CPS_2012_micro_medium
-----------------------------------------
. log close ;
name:
log: /home/users/dlmiller.AD3/colinjhr/work2014/table2.log
log type: text
closed on: 6 Jul 2014, 01:53:56
-------------------------------------------------------------------------------