qqplot p_webb_res p_rad_res if numstates == `gg' , 
		title(G=`gg') name(g_qq_`gg' , replace)

		di _column(50) "G(`clusvar')" _column(68) "=    " _result(18)

graph export G6_Wild_CDFs.png , replace

local G_ar1 = r(r)

local G_ols = e(N_clust) 					

graph combine g_qq_6 g_qq_10 g_qq_20 g_qq_50

graph combine g_qq_6 g_qq_10 g_qq_20 g_qq_50

graph combine g_sc_6 g_sc_10 g_sc_20 g_sc_50

graph combine g_sc_6 g_sc_10 g_sc_20 g_sc_50

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)

graph combine g_`type'_6 g_`type'_10 g_`type'_20 g_`type'_50

graph combine g_`type'_6 g_`type'_10 g_`type'_20 g_`type'_50

		local Gclusvar`c' = _result(18)

		    di "          See Cameron, Gelbach & Miller, "

local options "BY(string) GENerate(string) Detail"

foreach gg in 6 10 20 50 {

di in gr "Number of records is  "in ye "`N'"

/* partial out differences in demographics across state years.  So, get regression on */
egen styr = group(statefip year)

egen statefip2 = group(statefip id)

			local grouplist

                sub_robust `if' `in' `weight', groupvar(`groupvar') xxinv(`xxinv') plusminus(`plusminus') resid(`resid') running_sum(`running_sum') touse(`touse')

	syntax [if] [in] [aweight fweight iweight pweight /] , groupvar(string) xxinv(string) plusminus(string) resid(string) running_sum(string) touse(string)