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;

if chd=1 then chdny=.;
if chd=0 then chdny=0;
if chd>1 then chdny=1;

if age<=48 then ageny=0;
if 48<age<=52 then ageny=1;
if 52<age<=56 then ageny=2;
if 56<age then ageny=3;

run;

/* Resultater, s.2-6 */

proc genmod descending;
class ageny;
model chdny=ageny/dist=bin type3;
estimate 'alder 1 vs. 0' ageny -1 1 0 0/exp;
estimate 'alder 2 vs. 0' ageny -1 0 1 0/exp;
estimate 'alder 3 vs. 0' ageny -1 0 0 1/exp;
run;    

proc genmod descending;
model chdny=ageny/dist=bin type3;
estimate 'age odds ratio' ageny 1/exp;
run;    

/* Testet for linearitet (departures from trend) fås ved at trække
   deviance for de to seneste modeller fra hinanden */  

/* Resultater, s.8-9 */

proc genmod descending;
model chdny=age/dist=bin type3;
estimate 'age odds ratio pr år.' age 1/exp;
run;    

data framing; set framing;

agesq=age*age;

run;

proc genmod descending;
model chdny=age agesq/dist=bin type3;
run;    

/* Resultater, s.12-13 */

/* Først ageny som kategorisk (class) variabel */    

/* For at TESTE for ingen interaktion, tages både sex ageny og sex*ageny med */
    
proc genmod descending;
class ageny sex;
model chdny=ageny sex sex*ageny/dist=bin type3;
run;    

/* Dernæst ageny som kvantitativ variabel */    

/* For at TESTE for ingen interaktion, tages både sex ageny og sex*ageny med */
    
proc genmod descending;
class sex;
model chdny=ageny sex sex*ageny/dist=bin type3;
run;    

/* Dernæst age som kvantitativ variabel */    

/* For at TESTE for ingen interaktion, tages både sex age og sex*age med */
    
proc genmod descending;
class sex;
model chdny=age sex sex*age/dist=bin type3;
run;    

/* Resultater, s.14-17 */

/* For at ESTIMERE -OR'r separat for M og K udelades leddet ageny */

proc genmod descending;
class sex ageny;
model chdny=sex ageny*sex/dist=bin type3;
estimate 'alder 1 vs. 0: M' ageny*sex -1 1 0 0 0 0 0 0/exp;
estimate 'alder 2 vs. 0: M' ageny*sex -1 0 1 0 0 0 0 0/exp;
estimate 'alder 3 vs. 0: M' ageny*sex -1 0 0 1 0 0 0 0/exp;
estimate 'alder 1 vs. 0: K' ageny*sex 0 0 0 0 -1 1 0 0/exp;
estimate 'alder 2 vs. 0: K' ageny*sex 0 0 0 0 -1 0 1 0/exp;
estimate 'alder 3 vs. 0: K' ageny*sex 0 0 0 0 -1 0 0 1/exp;
run;    



proc genmod descending;
class sex;
model chdny=sex sex*ageny/dist=bin type3;
estimate 'alder, OR pr. gruppe: M' sex*ageny 1 0/exp;
estimate 'alder, OR pr. gruppe: K' sex*ageny 0 1/exp;
run;