makeplot.m

cd 'C:\Users\jcohen2\Desktop\Bryan Graham Stata software\bandwidth_selection_600\dir_600_500reps_makegraph'

data = csvread('data_for_graph_600_1.csv',1,0);

number = data(:,1);
b_ate = data(:,2);
seate = data(:,3);
z = data(:,4);
id = data(:,5)
v1 = data(:,6);
yvlbe = data(:,7)
x = data(:,8);
y = data(:,9);
dtreat = data(:,10);
dtime = data(:,11);

scatter3(x,y,z);
xlabel('l1');
ylabel('l2');
zlabel('True ATEs for 600 obs');
title('Figure 3a: True ATEs for 600 Observations Simulations, Repetition #1');
saveas(gcf,'figure_trueate_600.pdf');

z = b_ate

scatter3(x,y,z);
xlabel('l1');
ylabel('l2');
zlabel('ATEs from GIPT for 600 obs');
title('Figure 3b: ATEs from GIPT for 600 Observations Simulation, Repetition #1');
saveas(gcf,'figure_GIPT_ate_600.pdf');

cd 'C:\Users\jcohen2\Desktop\Bryan Graham Stata software\bandwidth_selection_300\dir_300_500reps_makegraph'

data = csvread('data_for_graph_300_1.csv',1,0);

number = data(:,1);
b_ate = data(:,2);
seate = data(:,3);
z = data(:,4);
id = data(:,5)
v1 = data(:,6);
yvlbe = data(:,7)
x = data(:,8);
y = data(:,9);
dtreat = data(:,10);
dtime = data(:,11);

scatter3(x,y,z);
xlabel('l1');
ylabel('l2');
zlabel('True ATEs for 300 obs');
title('Figure 2a: True ATEs for 300 Observations Simulations, Repetition #1');
saveas(gcf,'figure_trueate_300.pdf');

z = b_ate

scatter3(x,y,z);
xlabel('l1');
ylabel('l2');
zlabel('ATEs from GIPT for 300 obs');
title('Figure 2b: ATEs from GIPT for 300 Observations Simulation, Repetition #1');
saveas(gcf,'figure_GIPT_ate_300.pdf');

cd 'C:\Users\jcohen2\Desktop\Bryan Graham Stata software'