/* Table12_1997_2007_sim.sas

*/
libname public 'd:\a_data\proj\size_fun\exports\revision\wrap\public';run;run;


filename nu97 'd:\a_data\proj\size_fun\exports\revision\wrap\public\base97_nu_bsize.asc';


DATA nu_iter97; INFILE nu97;
input naicsindex naics iter_num boutiquesize beta1 beta2 beta3;
nuN=beta2;
nuT=beta3;
drop naics; *not character, so merge it in;

data nu_iter97; set nu_iter97;
if iter_num=10;
nuN97=nuN; nuT97=nuT;
keep naicsindex nuN97 nuT97;

*1997 naics_level;
*naics_level=naicsindex~naics~bni~sal_growth_chinaUS~sumc(est)~sumc(est_china)~alowerbar~(eta')~(gam');
filename n_lev97 'd:\a_data\proj\size_fun\exports\revision\wrap\public\base97_naics_level.asc'; 
DATA naics_level97; 
INFILE n_lev97;
input naicsindex naics iter_num  sal_growth_chinaUS est_modelUS est_chinaUS alowerbar eta1 eta2 gam1-gam177;
data naics_level97; 
set naics_level97;
if iter_num=10;
est_modelUS97=est_modelUS;
keep naicsindex est_modelUS97;

data nu_iter97;
merge nu_iter97 naics_level97; by naicsindex;
run;

*Now do 2007;
filename nu07 'd:\a_data\proj\size_fun\exports\revision\wrap\public\base07_nu_bsize.asc';

DATA nu_iter07; INFILE nu07;
input naicsindex naics iter_num boutiquesize nuN nuT beta1 beta2 beta3;

data nu_iter07; set nu_iter07;
if iter_num=10;
nuN07=nuN; nuT07=nuT;
keep naicsindex nuN07 nuT07;

*2007 naics_level
*naics_level=naicsindex~naics~bni~sumc(est)~alowerbar~(eta')~(gam');
filename n_lev07 'd:\a_data\proj\size_fun\exports\revision\wrap\public\base07_naics_level.asc';   
DATA naics_level07; 
INFILE n_lev07;
input naicsindex naics iter_num  est_modelUS alowerbar eta1 eta2 gam1-gam177;
data naics_level07; 
set naics_level07;
if iter_num=10;
est_modelUS07=est_modelUS;
keep naicsindex est_modelUS07;

data nu_iter07;
merge nu_iter07 naics_level07; by naicsindex;
run;




*grab totals for eachindustry form each file;
data est97;set public.ea_cm97;by naicsindex;if first.naicsindex;est97=estUS_LM; emp97=emphatUS_LM;keep naicsindex naics naicstext est97 emp97;
data est07;set public.ea_cbp07;by naicsindex;if first.naicsindex;est07=estUS_cbp; emp07=emphatUS_cbp;keep naicsindex est07 emp07;


data t;
merge nu_iter97(in=in1) nu_iter07 est97 est07 ; by naicsindex;
if in1;
share_N_97=100*(nuN97)/(nuN97 + nuT97*est_modelUS97 ); share_T_97 = 100-share_N_97;
share_N_07=100*(nuN07)/(nuN07 + nuT07*est_modelUS07 ); share_T_07 = 100-share_N_07;
estN97=.01*share_N_97*est97;
estT97=est97-estN97;
estN07=.01*share_N_07*est07;
estT07=est07-estN07;

run;

data tradedat;
set public.tradedat_forgauss;
length newChinacat $ 10;
     if newChina_shr07>=.50 then newChinacat='1 >=.5  ';
else if newChina_shr07>=.25 then newChinacat='2 .25-.5 ';
else if newChina_shr07>=.10 then newChinacat='3  .10-.25';
else if newChina_shr07>=.05 then newChinacat='4  .05-.10';
else if newChina_shr07>0    then newChinacat='5  .001-.05';
else if  newChina_shr07=0   then newChinacat='6  0';

     if newChina_shr07>=.25 then newChinacat2='1 >=.25  ';
else if newChina_shr07>=.05 then newChinacat2='2 .25-.05 ';
else if 0<=newChina_shr07<.05 then newChinacat2='3  <.05';
keep naics indcat newChina_shr07  newChinacat newChinacat2;

data t;
merge t(in=in1) tradedat;
by naics;
if in1;
if naics='339111' then delete;  *No trade data for 2010;

ods html file="d:\a_data\proj\size_fun\exports\revision\wrap\public\Table12_1997_2007_sim.html" style=minimal;

proc print data=t;
where naics='337122';
var naics naicstext est97 est07 estT97 estT07 estN97 estN07 share_N_97 share_N_07 ;
format share_N_97 share_N_07 estT97 estN97 estT07 estN07 7.1;
title 'where naics=337122';


proc sort data=t;
by newChinacat;

proc means data=t noprint;
by newChinacat;
var estT97 estN97 estT07 estN07 emp97 emp07 est97 est07 ;
output out=tsum sum=estT97 estN97 estT07 estN07 emp97 emp07 est97 est07 ;

data tsum; set tsum;
share_N_97=100*estN97/(estN97+estT97);
share_N_07=100*estN07/(estN07+estT07);
share_T_97=100-share_N_97;
share_T_07=100-share_N_07;
difN07=share_T_07-share_T_97;
est000_97=est97/1000;
est000_07=est07/1000;
format estT97 estN97 estT07 estN07 emp97 emp07 est97 est07 comma12.0 share_N_97 share_N_07 share_T_97 share_T_07  7.1;

proc print data=tsum;
var  newChinacat  est000_97 est000_07 share_T_97 share_T_07 difN07;
format est000_97 est000_07  7.0;
format share_T_97 share_T_07 difN07 7.0;
title 'by New China Share Category, 1997 to 2007';

run;

ods html close;
run;