***************************************************************
***							    			
*** DO-file for:					    			
***							    			
***	Gretchen Casper and Claudiu Tufis		    		
***							    		
***	Correlation versus Interchangeability: The Limited  	
***	Robustness of Empirical Findings on Democracy	    	
*** 	using Highly Correlated Datasets		    	
***							    	
***	Political Analysis, Vol.11 No.2			    	
***************************************************************

set more off

#delimit;

log using "C:\temp\PA_LOG_file.log", replace;

set mem 80m;

set matsize 200;

use "C:\temp\PA_DTA_file.dta", clear;

summarize _all;

tis year;

iis id;

sort id year;

tsset id year, yearly;

***************************
***					
*** Table 1 - Column 1  
***				
***************************;

xtpcse polityiv L1.pcaplog L1.rgdppcgr L1.open L1.cpi L1.prime 
	L1.second L1.presiden L1.parliamn L1.bksfrac if dv1==1 & 
	dv2==1, pairwise c(a);

***************************
***				
*** Table 1 - Column 2  
***				
***************************;

xtpcse poly12 L1.pcaplog L1.rgdppcgr L1.open L1.cpi L1.prime 
	L1.second L1.presiden L1.parliamn L1.bksfrac if dv1==1 & 
	dv2==1, pairwise c(a);

***************************
***				
*** Table 1 - Column 3  
***				
***************************;

drop if year == 1950;
drop if year == 1974;
drop if year == 1975;
drop if year == 1976;
drop if year == 1977;
drop if year == 1978;
drop if year == 1979;
drop if year == 1980;
drop if year == 1981;
drop if year == 1982;
drop if year == 1983;
drop if year == 1984;
drop if year == 1985;
drop if year == 1986;
drop if year == 1987;
drop if year == 1988;
drop if year == 1989;
drop if year == 1990;
drop if year == 1991;
drop if year == 1992;

tsset id year, yearly;

xtpcse polityiv L1.pcaplog L1.rgdppcgr L1.open L1.cpi L1.prime 
	L1.second L1.presiden L1.parliamn L1.bksfrac if dv1==1 & 
	dv2==1, pairwise c(a);

***************************
***				
*** Table 1 - Column 4  
***				
***************************;

xtpcse poly12 L1.pcaplog L1.rgdppcgr L1.open L1.cpi L1.prime 
	L1.second L1.presiden L1.parliamn L1.bksfrac if dv1==1 & 
	dv2==1, pairwise c(a);

clear;

use "C:\temp\PA_DTA_file.dta", clear;

***************************
***				
*** Table 1 - Column 5  
***				
***************************;

drop if year == 1950;
drop if year == 1951;
drop if year == 1952;
drop if year == 1953;
drop if year == 1954;
drop if year == 1955;
drop if year == 1956;
drop if year == 1957;
drop if year == 1958;
drop if year == 1959;
drop if year == 1960;
drop if year == 1961;
drop if year == 1962;
drop if year == 1963;
drop if year == 1964;
drop if year == 1965;
drop if year == 1966;
drop if year == 1967;
drop if year == 1968;
drop if year == 1969;
drop if year == 1970;
drop if year == 1971;
drop if year == 1972;
drop if year == 1973;
drop if year == 1974;

tsset id year, yearly;

xtpcse polityiv L1.pcaplog L1.rgdppcgr L1.open L1.cpi L1.prime 
	L1.second L1.presiden L1.parliamn L1.bksfrac if dv1==1 & 
	dv2==1 & dv3==1, pairwise c(a);

***************************
***				
*** Table 1 - Column 6  
***				
***************************;

xtpcse poly12 L1.pcaplog L1.rgdppcgr L1.open L1.cpi L1.prime 
	L1.second L1.presiden L1.parliamn L1.bksfrac if dv1==1 & 
	dv2==1 & dv3==1, pairwise c(a);

***************************
***				
*** Table 1 - Column 7  
***				
***************************;

xtpcse fhscore L1.pcaplog L1.rgdppcgr L1.open L1.cpi L1.prime 
	L1.second L1.presiden L1.parliamn L1.bksfrac if dv1==1 & 
	dv2==1 & dv3==1, pairwise c(a);

clear;

use "C:\temp\PA_DTA_file.dta", clear;

***************************
***				
*** Table 2 - Column 1  
***				
***************************;

xtpcse polityiv L1.pcaplog L1.rgdppcgr L1.open L1.cpi L1.prime 
	L1.second L1.presiden L1.parliamn L1.bksfrac if dv1==1 & 
	dv2==1 & dv3==1 & oecd==1, pairwise c(a);

***************************
***				
*** Table 2 - Column 2  
***				
***************************;

xtpcse poly12 L1.pcaplog L1.rgdppcgr L1.open L1.cpi L1.prime 
	L1.second L1.presiden L1.parliamn L1.bksfrac if dv1==1 & 
	dv2==1 & dv3==1 & oecd==1, pairwise c(a);

***************************
***				
*** Table 2 - Column 3  
***				
***************************;

xtpcse fhscore L1.pcaplog L1.rgdppcgr L1.open L1.cpi L1.prime 
	L1.second L1.presiden L1.parliamn L1.bksfrac if dv1==1 & 
	dv2==1 & dv3==1 & oecd==1, pairwise c(a);

***************************
***				
*** Table 2 - Column 4  
***				
***************************;

xtpcse polityiv L1.pcaplog L1.rgdppcgr L1.open L1.cpi L1.prime 
	L1.second L1.presiden L1.parliamn L1.bksfrac if dv1==1 & 
	dv2==1 & dv3==1 & oecd==0, pairwise c(a);

***************************
***				
*** Table 2 - Column 5  
***				
***************************;

xtpcse poly12 L1.pcaplog L1.rgdppcgr L1.open L1.cpi L1.prime 
	L1.second L1.presiden L1.parliamn L1.bksfrac if dv1==1 & 
	dv2==1 & dv3==1 & oecd==0, pairwise c(a);

***************************
***				
*** Table 2 - Column 6  
***				
***************************;

xtpcse fhscore L1.pcaplog L1.rgdppcgr L1.open L1.cpi L1.prime 
	L1.second L1.presiden L1.parliamn L1.bksfrac if dv1==1 & 
	dv2==1 & dv3==1 & oecd==0, pairwise c(a);

log close;

clear;

exit;