data zadatak5;
input dob tip trajanje CRP FIB HCY ADN bolest BMI BMIkat spol RRS RRD HbA1c GUKb LEUK Ac_UR AER_0 BMI3$;
cards;
69 2 17 12 2.8 8 10.1 0 29 2 0 116 99 5.8 12.5 3.5 239 27.8 25-30
30 2 13 13 3 14 15.4 0 29 2 0 111 77 5.5 17.7 6.9 478 96.3 25-30
78 2 8 0 5.8 12 23.5 1 19 1 0 128 97 9.5 2.9 6.9 330 28.8 <25
76 1 20 12 3.2 13 18.4 1 40 3 1 115 108 8.4 15.5 4.7 234 77.1 >30
70 2 17 3 6.5 15 5.6 0 40 3 0 147 110 7.6 5.5 12 213 16.8 >30
52 0 7 7 2.8 17 3.8 0 38 3 0 126 78 5.4 7 8.5 368 112 >30
37 0 14 4 2.8 7 16.5 1 39 3 0 136 108 7.1 5.3 10.5 505 32.4 >30
66 1 10 6 4.6 21 19.2 1 25 2 1 152 94 5.5 15.8 10.1 469 82.5 25-30
68 1 16 6 5.6 17 19.2 1 27 2 0 135 75 7.5 13.2 9.1 476 77.2 25-30
29 2 17 2 2.3 11 10.4 0 35 3 0 139 118 4.8 7.4 5.5 498 118.6 >30
69 2 10 7 3.5 6 12.7 0 32 3 1 147 120 10.6 7.2 8.5 433 51.6 >30
22 2 1 4 6.1 19 7 0 40 3 0 135 115 8.1 8.5 12.7 318 80.3 >30
65 2 3 14 2.7 11 15.5 0 19 1 1 136 75 9.1 3.9 10.8 309 50.4 <25
62 2 14 2 3.3 9 2.5 0 39 3 1 132 106 6.3 17.2 9.2 183 21.8 >30
56 2 11 12 3.5 17 17.5 1 21 1 1 120 95 5.7 15.1 7.7 447 97.6 <25
35 1 10 7 6 16 24.3 1 20 1 1 138 76 6.9 15.6 12.5 136 81 <25
27 0 1 4 3.7 9 17.4 1 36 3 0 130 91 6.7 7 12.1 275 40.7 >30
60 0 14 11 6.2 17 8.5 0 20 1 0 107 80 10.7 7 4.6 445 10.3 <25
41 2 2 12 1.8 9 2.6 0 33 3 0 150 93 5.3 17.7 11.9 346 109.3 >30
57 0 10 2 5.9 10 11.9 0 27 2 1 131 85 8.3 7.5 10.2 538 36.2 25-30
68 1 6 10 5.3 14 9.6 0 23 1 0 156 111 10.8 6 4.6 418 5.4 <25
40 2 3 11 4.1 17 8.4 0 40 3 0 105 90 4.8 14 8.9 263 60.6 >30
37 1 6 5 3.5 15 15.2 0 31 3 0 142 94 10 11 8.1 149 87 >30
46 2 16 10 2.1 12 22.2 1 29 2 0 142 98 5.2 5.9 10.3 128 60 25-30
51 0 5 9 2.9 11 9.6 0 26 2 1 105 80 8.1 9 11.8 132 115.1 25-30
67 2 17 7 5.1 15 31 1 18 1 0 160 115 5.4 17.5 10.9 424 87.2 <25
28 0 18 4 4.8 14 26.7 1 29 2 1 125 98 10.2 9 9.3 310 16.3 25-30
52 2 18 7 4.2 17 26.6 1 26 2 0 151 80 10.6 2.5 9.8 373 45.8 25-30
23 2 15 14 2.6 6 27 1 39 3 0 130 109 5.7 7.7 12.8 239 33 >30
37 2 4 8 2.8 9 25.7 1 25 2 1 153 82 9.9 7.8 3.5 541 33.8 25-30
51 2 16 14 2.1 23 23.3 1 24 1 1 129 78 6.2 7.8 12.1 230 47.4 <25
72 2 15 8 2.5 22 30.4 1 40 3 1 138 94 9.4 6.5 7.3 320 115.7 >30
43 2 1 4 4.2 22 26.8 1 33 3 0 123 77 5.6 17.9 12.5 149 18.6 >30
38 2 4 0 3.4 21 4.1 0 27 2 0 119 109 5.8 17.2 9.9 392 91.9 25-30
18 2 14 0 4.9 15 15 0 36 3 0 119 80 4.9 7.3 7.2 127 21.8 >30
60 2 4 10 4.9 9 15.4 0 21 1 1 121 114 7.8 8.3 10.8 521 104.5 <25
43 2 2 3 2.5 7 13.5 0 30 2 0 130 80 5.5 12.8 11.6 276 6.4 25-30
39 0 13 6 4.5 23 26 1 18 1 1 128 88 7.9 5 3.8 392 30.8 <25
20 2 12 4 1.9 11 13.5 0 21 1 0 135 90 6.6 13.9 3.5 314 90.2 <25
55 1 19 2 5.2 7 24.8 1 39 3 0 130 90 6.7 6.9 11.1 422 49 >30
42 0 17 12 2 9 2.9 0 19 1 1 117 98 7.7 2.2 11.5 161 42.6 <25
44 0 8 4 1.9 6 6.6 0 39 3 0 147 82 5.7 7 11.8 425 54.9 >30
18 1 6 10 4.2 6 5.5 0 20 1 1 111 107 6.6 16.9 7.6 495 9.9 <25
64 2 15 6 5.5 18 30 1 35 3 0 145 110 5.5 9.4 11.3 362 70.1 >30
77 2 4 8 5.1 17 15.8 0 27 2 1 100 75 9.9 9 7.3 261 67 25-30
25 1 2 12 1.9 14 13.1 0 26 2 1 115 107 7.4 17.8 3.7 425 113.1 25-30
79 2 2 0 5.6 8 12.5 0 32 3 0 122 87 7.7 17.6 5.9 545 30.9 >30
27 0 12 2 4.6 19 1.7 0 18 1 0 117 77 10.2 2.1 10.9 501 66.2 <25
64 1 16 10 5.3 21 21.5 1 32 3 0 138 83 5.6 17.7 6.3 453 18.7 >30
41 0 1 11 4.5 9 6.1 0 20 1 1 126 101 6.2 6.1 10.1 414 29.2 <25
58 2 8 6 5.6 19 4.1 0 30 2 1 143 94 4.8 12.8 11.1 429 75.7 25-30
77 2 15 10 2.3 21 22.9 1 30 2 0 132 86 5.2 14.3 10.9 400 22.5 25-30
19 2 11 5 2.2 19 3 0 40 3 0 150 119 5.3 4.4 8.9 358 22.4 >30
49 0 13 9 5 13 12.4 0 25 2 0 153 78 5.3 2.8 7.7 367 35.2 25-30
26 0 5 4 5.1 22 7.2 0 26 2 1 103 101 7.6 5.6 9.4 155 39.1 25-30
69 0 8 6 3 18 29.8 1 33 3 0 154 81 9.3 6 5.6 164 91.3 >30
65 2 6 8 4.3 12 28.7 1 21 1 1 133 106 9.2 17.4 11.4 245 101.7 <25
47 2 2 9 1.9 14 11.4 0 38 3 0 110 108 8 11.3 5.6 477 60.7 >30
52 0 10 8 4 23 24.3 1 28 2 1 129 98 9.9 7 6 198 76.3 25-30
23 0 19 10 5.5 18 25.6 1 35 3 0 156 117 9.5 4 3.6 291 36.9 >30
37 2 13 11 4.7 7 7.6 0 19 1 1 129 96 6.7 15.1 12.6 474 52.6 <25
45 2 14 14 6.1 14 25.3 1 39 3 0 118 100 5.5 13.9 12 411 103.5 >30
49 0 2 6 4.6 16 5.9 0 26 2 0 139 80 10.7 3 10.1 279 78.1 25-30
26 2 7 5 3.7 17 30.9 1 22 1 0 139 116 7.2 11.8 12 204 55.5 <25
69 2 11 7 4.3 15 9.3 0 27 2 1 135 110 5.8 11.8 6.9 312 33.4 25-30
41 2 1 12 5.5 20 21.3 1 35 3 0 106 83 9.7 6.4 10.6 331 95.2 >30
44 0 20 14 4.1 11 14.1 0 38 3 0 148 105 7.9 4 11.9 164 74.5 >30
72 2 4 11 4.7 7 15.3 0 32 3 0 119 105 10.3 6 4.4 476 12.1 >30
31 1 19 11 5.2 6 30 1 33 3 1 146 82 10.7 5.4 4.6 359 71.8 >30
35 1 6 6 4 16 9.8 0 20 1 0 125 77 5.4 11.9 6.1 174 95.8 <25
49 2 9 7 6.1 6 1.8 0 21 1 1 111 99 6.3 8.7 7 336 97.5 <25
49 2 6 6 4.3 20 28 1 29 2 0 107 92 5.3 12.5 12.9 547 15.6 25-30
64 1 9 2 5.5 19 20.6 1 38 3 1 118 90 5 17.8 11.9 121 110.8 >30
61 2 16 5 6.3 12 8 0 35 3 0 138 102 10.5 16.2 11.8 214 40.7 >30
68 2 15 13 3.9 6 25.4 1 42 3 0 102 80 5.8 5 6.9 502 4.3 >30
25 2 16 2 5.3 19 8.8 0 35 3 1 123 118 8 13.6 4.8 257 22.4 >30
66 2 15 0 4.3 12 28.5 1 30 2 0 104 76 5 6.4 5.9 175 92 25-30
67 2 2 6 6.3 10 20.4 1 21 1 0 103 83 6.4 10.9 4.8 498 92.6 <25
43 1 5 3 3.9 7 18.8 1 20 1 0 145 114 9.4 5 6.1 473 52.4 <25
55 2 17 6 3.1 9 11.4 0 29 2 1 157 116 6.5 17.5 4.1 277 103 25-30
;

data dummy; set zadatak5;
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;
run;


proc logistic data=dummy descending;
model bolest=BMI0 BMI2/lackfit rsq outroc=rocgraf;
run;

title" Univarijatna logisticka regresija CRP";
proc logistic data = zadatak5 descending;
model bolest = CRP / lackfit rsq outroc=rocgraf;
run;

title" Univarijatna logisticka regresija HCY";
proc logistic data = zadatak5 descending;
model bolest = HCY / lackfit rsq outroc=rocgraf;
run;


title"Multipla logisticka regresija";
proc logistic data = dummy descending;
model bolest = CRP HCY BMI0 BMI2;
run;

title"Multipla logisticka regresija-stepwise";
proc logistic data = dummy descending;
model bolest = CRP HCY BMI0 BMI2 / selection=stepwise;
run;