data framing;

infile 't:\framing.txt' firstobs=2;
input id sex age frw sbp sbp10 dbp chol cig chd yrs_chd death yrs_dth cause;

logsbp=log(sbp);
agegrp=(age>=49)+(age>=53)+(age>=57);

diffsbp=sbp10-sbp;
difflog=log(sbp10)-logsbp;

run;

proc sort; by sex; run;

/* spm. 1 */
    
proc univariate;
    var sbp logsbp;
    by sex;
run;    

/* spm. 2 */

proc ttest;
class sex;
var sbp;
run;

proc ttest;
class sex;
var logsbp;
run;

proc genmod;
class sex;
model sbp=sex/type3;
run;

proc genmod;
class sex;
model logsbp=sex/type3;
run;

/* spm. 3 */

proc univariate;
    var diffsbp difflog;
run;    

/* spm. 4 */

proc genmod;
class agegrp;
model logsbp=agegrp/type3;
run;

proc npar1way;
    var logsbp;
    class agegrp;
run;

/* spm. 5 */

proc genmod;
model logsbp=frw/type3;
run;

proc genmod;
class sex agegrp;
model logsbp=frw sex agegrp/type3;
run;

/* spm. 6 */

proc lifetest notable plots=(s);
    time yrs_dth*cause(0);
    strata sex;
run;

/* spm. 7 */

proc lifetest notable plots=(s);
    time yrs_dth*cause(0);
    strata agegrp;
run;