195Sensitivity Analyses
var total_costs mean_eqos_pre mean_eqos_post os pfs pps;
run;
**USE IML to simulate multivariate data**;
proc iml ;
use pattern ;
read all var _num_ into F;
F=F`;
data=rannor(J(10000,5,0)); **10,000 simulations**;
data=data`;
Z=F*data;
Z=Z`;
x1=Z[,1]*9736+10290; ** linear transformation for a specified
mean and SD**;
x2=Z[,2]*SD1+mean1; **specify values for the mean and SD
here**;
x3=Z[,3]*SD2+mean2;
x4=Z[,4]*SD3+mean3;
x5=Z[,5]*SD4+mean4;
x6=Z[,6]*SD53+mean5;
Z=x1||x2||x3||x4||x5||x6;
create A from Z [COLNAME={x1 x2 x3 x4 x5 x6}];
append from Z;
quit;
proc means data=a n mean std skew kurt;
var x1 x2 x3 x4 x5 x6;
run;
proc means data=costs2 n mean std skew kurt;
class grash1_txt ;
var total_costs mean_eqos_pre mean_eqos_post os pfs pps;
run; ...