libname heartlib "/home/u63742108/heartlib";
run;

proc means data=heartlib.myheart nway;
class weight_statusnew;
var chd;
output out=Out1 Mean=CHD_Mean;
run;

data out2;
set out1;
Logit = log( CHD_Mean / (1 - CHD_Mean) );
run;

proc print data=out2;
run;

proc sgplot data=out2;
series y=logit x=weight_statusnew / markers;
label logit="Logit of Developing CHD" weight_statusnew="Weight_Status";
yaxis values=(-2 to 0.25 by .25);
run;

proc rank data=heartlib.myheart out=out3 groups=20;
var mrw height weight;
ranks mrw_Group height_Group weight_Group;
run;

proc means data=out3;
class mrw_group;
var mrw;
run;

proc means data=out3;
class height_group;
var height;
run;

proc means data=out3;
class weight_group;
var weight;
run;

proc means data=out3 nway;
class mrw_group;
var mrw;
output out=Out4a Mean=CHD_Mean mean(mrw)=MRW_Mean;
run;

proc print data=out4a noobs;
run;

proc means data=out3 nway;
class height_group;
var height;
output out=Out4b Mean=CHD_Mean mean(height)=height_Mean;
run;

proc print data=out4b noobs;
run;

proc means data=out3 nway;
class weight_group;
var weight;
output out=Out4c Mean=CHD_Mean mean(weight)=weight_Mean;
run;

proc print data=out4c noobs;
run;

data out5a;
set out4a;
Logit_Dia = log( CHD_Mean / (1 - CHD_Mean) );
run;

proc print data=out5a noobs;
run;

data out5b;
set out4b;
Logit_Dia = log( CHD_Mean / (1 - CHD_Mean) );
run;

proc print data=out5b noobs;
run;

data out5c;
set out4c;
Logit_Dia = log( CHD_Mean / (1 - CHD_Mean) );
run;

proc print data=out5c noobs;
run;