50 Design & Analysis of Clinical Trials for Economic Evaluation
set lsmqaly;
run;
proc transpose data=lsmqaly2 out=lsmqaly3;
by bootno;
id trt;
var estimate;
run;
data bootdiffqaly;
set lsmqaly3;
incr_qaly_boot=a-b;**get mean incremental effects**;
run;
**merge bootstrap costs and qalys**;
data boot_all1;
merge bootdiffqaly bootdiffcosts;
by bootno;
icer=incr_costs_boot/incr_qaly_boot;**compute ICER**;
run;
**Sort ICER in order**;
proc sort data=boot_all1;
by icer;
run;
**statistics**;
proc univariate data=boot_all1;
var icer;
output out=final pctlpts=2.5, 97.5 pctlpre=ci;
run;
2A.2 STATA Command for Bootstrap
The STATA command for Bootstrap can also be accessed via http://www.
uphs.upenn.edu/dgimhsr/stat-cicer.htm.
First compute the NMB for each treatment. Assign a variable ...