115Thunderstorm or Raindrop Scatter Plots
call symputx(“Ref5”, refn5);
call symputx(“Ref6”, refn6);
call symputx(“Ref7”, refn7);
call symputx(“Ref8”, refn8);
call symputx(“Ref9”, refn9);
call symputx(“Ref10”, refn10);
run;
data mndiur_all;
merge iop_MnDiur subjid;
by siteid subjid;
run;
proc transpose data = MnDiur_all out = MnDiur_T prefix = vst_;
by siteid subjid subj_index;
var iop_MnDiur;
run;
** Mean Diurnal IOP at each visit;
data mndiur_vst;
set MnDiur_T;
drop _NAME_;
rename vst_1 = BSL vst_2 = WK4 vst_3 = WK8 vst_4 = WK12;
run;
proc format;
value indexid
1 = ‘ ’
OTHER = ‘ ’;
run;
** Annotate dataset to connect the 2 IOPs at baseline and;
** week12 of the same subject by a straight line;
** (forming a rain drop);
data anno_w12; ...