data zadatak7;
input spol zaposlen dijete dob BMI BMIkat;
cards;
2 0 0 23 23 1
2 0 0 21 26 2
2 0 0 33 27 2
2 0 0 28 18 1
2 0 0 33 24 1
2 0 0 21 35 2
2 0 0 34 21 1
2 0 0 26 20 1
2 0 0 20 34 2
2 0 0 33 37 3
2 0 0 35 27 2
2 0 0 33 29 2
2 0 0 27 35 2
2 0 0 26 34 2
2 0 1 20 28 2
2 0 1 18 36 3
2 0 1 23 32 2
2 1 0 30 37 3
2 1 0 30 18 1
2 1 0 22 22 1
2 1 0 34 33 2
2 1 0 31 32 2
2 1 0 35 18 1
2 1 0 29 32 2
2 1 0 35 31 2
2 1 0 30 28 2
2 1 1 19 21 1
2 1 1 33 24 1
2 1 1 34 35 2
2 1 1 35 31 2
2 1 1 32 38 3
2 1 1 34 36 3
2 1 1 34 22 1
2 1 1 24 30 2
2 1 1 32 25 2
2 1 1 31 34 2
2 1 1 30 19 1
2 1 1 32 20 1
2 1 1 22 27 2
2 1 1 20 20 1
2 1 1 21 35 2
2 1 1 24 19 1
2 1 1 20 37 3
2 1 1 24 36 3
1 1 0 21 33 2
1 1 0 19 18 1
1 1 0 19 22 1
1 1 0 27 30 2
1 1 0 32 25 2
1 1 0 21 31 2
1 1 0 23 26 2
1 1 0 27 26 2
1 1 0 25 28 2
1 1 0 32 38 3
1 1 1 20 25 2
1 1 1 32 32 2
1 1 1 27 25 2
1 1 1 29 24 1
1 1 1 22 27 2
1 1 1 24 20 1
1 1 1 20 26 2
1 1 1 28 30 2
1 1 1 31 38 3
1 1 1 23 28 2
1 1 1 23 22 1
1 1 1 32 28 2
1 1 1 25 28 2
1 1 1 20 36 3
1 1 1 29 28 2
1 1 1 30 20 1
1 1 1 29 30 2
1 0 1 23 38 3
1 0 1 31 23 1
1 0 1 31 34 2
1 0 1 19 23 1
1 0 1 25 38 3
1 0 1 21 20 1
1 0 1 34 35 2
1 0 1 30 32 2
1 0 1 28 19 1
1 0 1 23 35 2
1 0 1 32 38 3
1 0 1 34 37 3
1 0 1 29 26 2
1 0 1 22 35 2
1 0 1 21 26 2
1 0 0 22 32 2
1 0 0 35 35 2
1 0 0 31 22 1
1 0 0 22 18 1
1 0 0 21 21 1
1 0 0 19 33 2
1 0 0 18 36 3
1 0 0 24 25 2
1 0 0 25 33 2
1 0 0 34 27 2
1 0 0 26 31 2
1 0 0 22 29 2
1 0 0 26 31 2
1 0 0 20 28 2
;

title"Modeliranje BMI pomocu dummy varijabli baseline 25-30";
data dummy; set zadatak7;
if BMI < 25 then do; BMI0=1;BMI2=0;end;
if BMI <= 30 && BMI >= 25 then do; BMI0=0;BMI2=0;end;
if BMI > 30 then do; BMI0=0;BMI2=1;end;


title" Univarijatna COX regresija";
proc phreg data=dummy ;
model dob*dijete(0)=BMI0 BMI2 / rl;
run;

proc phreg data=zadatak7 ;
model dob*dijete(0)= spol/rl ;
run;

proc phreg data=zadatak7 ;
model dob*dijete(0)= zaposlen /rl;
run;

title"Multipla COX regresija";
proc phreg data=dummy ;
model dob*dijete(0)=spol zaposlen BMI0 BMI2 / rl;
run;

title"Multipla cox regresija-stepwise";
proc phreg data=dummy ;
model dob*dijete(0)=spol zaposlen BMI0 BMI2/selection=stepwise rl details  ;
run;

title "interakcija";
data zadatak7; set zadatak7;
zs = zaposlen*spol;
proc phreg data=zadatak7;
model dob*dijete(0)= zaposlen spol zs / rl ;
run;