Content
--------------------------------------------------------------------------------------------------------------------------------------------
name:
log: C:\doug work\Colin\JHR cluster paper\work Summer 2014\empirical example code\t2_pvalues_check.log
log type: text
opened on: 17 Jun 2014, 14:12:53
. clear ;
. set more off ;
. use t2_mcout_all ;
. foreach type in webb rad { ;
2. foreach gg in 6 10 20 50 { ;
3. quantile p_`type'_res if numstates == `gg' ,
> title(`type' G=`gg') name(g_`type'_`gg' , replace) ;
4. summ p_`type'_res if numstates == `gg' , detail ;
5. } ;
6. graph combine g_`type'_6 g_`type'_10 g_`type'_20 g_`type'_50 ;
7. graph export quantiles_`type'.png , replace ;
8. } ;
p_webb_res
-------------------------------------------------------------
Percentiles Smallest
1% .0075 0
5% .02125 .005
10% .035 .01 Obs 200
25% .1725 .01 Sum of Wgt. 200
50% .4675 Mean .460425
Largest Std. Dev. .310091
75% .73875 .985
90% .905 .99 Variance .0961564
95% .9475 .99 Skewness .0786398
99% .99 .9975 Kurtosis 1.704376
p_webb_res
-------------------------------------------------------------
Percentiles Smallest
1% .015 .01
5% .0425 .015
10% .0775 .015 Obs 200
25% .1925 .015 Sum of Wgt. 200
50% .50375 Mean .4918
Largest Std. Dev. .3101291
75% .76 .99
90% .925 .995 Variance .09618
95% .9675 1 Skewness .0391241
99% .9975 1 Kurtosis 1.639888
p_webb_res
-------------------------------------------------------------
Percentiles Smallest
1% .0075 0
5% .045 .005
10% .0825 .01 Obs 200
25% .21 .01 Sum of Wgt. 200
50% .4525 Mean .4797875
Largest Std. Dev. .3038135
75% .74 .98
90% .9175 .99 Variance .0923026
95% .9675 .99 Skewness .1561247
99% .99 .995 Kurtosis 1.732289
p_webb_res
-------------------------------------------------------------
Percentiles Smallest
1% .0075 0
5% .0475 .005
10% .0975 .01 Obs 200
25% .2625 .02 Sum of Wgt. 200
50% .55 Mean .5111
Largest Std. Dev. .2817223
75% .7325 .96
90% .885 .98 Variance .0793674
95% .935 .99 Skewness -.1614036
99% .985 1 Kurtosis 1.841761
(file quantiles_webb.png written in PNG format)
p_rad_res
-------------------------------------------------------------
Percentiles Smallest
1% .015 .01
5% .02375 .015
10% .04125 .015 Obs 200
25% .1825 .0175 Sum of Wgt. 200
50% .4675 Mean .459225
Largest Std. Dev. .3057698
75% .7225 .975
90% .88125 .98 Variance .0934952
95% .95 .9825 Skewness .0840705
99% .98125 .995 Kurtosis 1.727491
p_rad_res
-------------------------------------------------------------
Percentiles Smallest
1% .01375 .01
5% .04375 .01
10% .08 .0175 Obs 200
25% .205 .0225 Sum of Wgt. 200
50% .49125 Mean .488025
Largest Std. Dev. .3087339
75% .755 .9925
90% .92375 .995 Variance .0953166
95% .96375 1 Skewness .059974
99% .9975 1 Kurtosis 1.650329
p_rad_res
-------------------------------------------------------------
Percentiles Smallest
1% .005 0
5% .045 0
10% .0875 .01 Obs 200
25% .2025 .015 Sum of Wgt. 200
50% .45875 Mean .477475
Largest Std. Dev. .3018778
75% .73 .9875
90% .915 .99 Variance .0911302
95% .945 .995 Skewness .1547607
99% .9925 .995 Kurtosis 1.737922
p_rad_res
-------------------------------------------------------------
Percentiles Smallest
1% .005 0
5% .045 .005
10% .1025 .005 Obs 200
25% .2625 .025 Sum of Wgt. 200
50% .5425 Mean .50705
Largest Std. Dev. .2808794
75% .7325 .97
90% .865 .97 Variance .0788933
95% .935 .985 Skewness -.1481014
99% .9775 .99 Kurtosis 1.853053
(file quantiles_rad.png written in PNG format)
. foreach gg in 6 10 20 50 { ;
2. qqplot p_webb_res p_rad_res if numstates == `gg' ,
> title(G=`gg') name(g_qq_`gg' , replace) ;
3. graph twoway (scatter p_rad_res p_webb_res if numstates == `gg')
> (line p_webb_res p_webb_res if numstates == `gg') ,
> title(G=`gg') name(g_sc_`gg' , replace) ;
4. } ;
. graph combine g_qq_6 g_qq_10 g_qq_20 g_qq_50 ;
. graph export qq_plot.png , replace ;
(file qq_plot.png written in PNG format)
. graph combine g_sc_6 g_sc_10 g_sc_20 g_sc_50 ;
. graph export scatter_plot.png , replace ;
(file scatter_plot.png written in PNG format)
. log close ;
name:
log: C:\doug work\Colin\JHR cluster paper\work Summer 2014\empirical example code\t2_pvalues_check.log
log type: text
closed on: 17 Jun 2014, 14:13:02
--------------------------------------------------------------------------------------------------------------------------------------------