proc import datafile = 'location of the excel spreadsheet'
	out = temp
	dbms = excelcs;
run;

%macro ttest(var);
proc ttest data=temp;
class tx_group;
var &var;
run;
%mend;

%MACRO BOOTSTRAP (DRP=,Y=,X=,INPUTDATA=,TITLE=);
TITLE1 'BOOTSTRAP CI';
TITLE2 "&TITLE";

PROC REG NOPRINT DATA=&INPUTDATA OUTEST=OUT;
	MODEL &Y = &X;

PROC PRINT DATA=OUT;

%INCLUDE 'location of jackboot sas program';
	%MACRO ANALYZE(DATA=,OUT=);
	PROC REG NOPRINT DATA= &DATA OUTEST=&OUT (DROP = &DRP _MODEL_ _TYPE_ _DEVAR_ _RMSE_);
		MODEL &Y = &X;
		%BYSTMT;
	%MEND ANALYZE;

TITLE1 'NORMAL CI WITH BIAS CORRECTION';
TITLE2 "&TITLE";
%BOOT (DATA= &INPUTDATA, SAMPLES=2000, RANDOM=12345, CHART=0);

TITLE1'BCA CI';
TITLE2"&TITLE";
%BOOTCI(BCA);

%MEND BOOTSTRAP;

title1'Day of surgery cost';
%ttest(cost_surg_norobot);
%bootstrap(drp=cost_surg_norobot, y=cost_surg_norobot, x=tx_group, inputdata=temp, title=day of surgery cost no robot);
%ttest(cost_surg_robot);
%bootstrap(drp=cost_surg_robot, y=cost_surg_robot, x=tx_group, inputdata=temp, title=day of surgery cost robot);

title1'6 week cost';
%ttest(cost_6wk_norobot);
%bootstrap(drp=cost_6wk_norobot, y=cost_6wk_norobot, x=tx_group, inputdata=temp, title=six wk total cost no robot);
%ttest(cost_6wk_robot);
%bootstrap(drp=cost_6wk_robot, y=cost_6wk_robot, x=tx_group, inputdata=temp, title=six wk total cost robot);

title1'Surgical time';
%ttest(surg_time);
%ttest(orsurg_time);