* proc genmod for Bernoulli (0/1) data with logit link; proc genmod data=heart descending; class id; model mi=oc age smoke / link=logit dist=bin; run; * column 1 subject id column 2 oral contraceptive indicator (0=no,1=yes) column 3 age (years) column 4 smoking indicator (0=no,1=yes) column 5 binary response -- whether MI has been suffered (0=no,1=yes); data heart; input id oc age smoke mi; datalines; 1 1 33 1 0 2 0 32 0 0 3 1 37 0 1 4 0 36 0 0 5 1 50 1 1 6 1 40 0 0 7 0 35 0 0 8 1 33 0 0 9 1 33 0 0 10 0 31 0 0 11 0 25 0 0 12 1 30 0 0 13 1 31 0 1 14 1 37 0 0 15 1 31 1 0 16 0 29 0 0 17 0 47 0 0 18 1 40 0 0 19 0 31 0 0 20 0 43 1 1 21 0 37 0 0 22 1 34 0 0 23 0 35 1 0 24 0 37 0 0 25 1 34 0 0 26 1 27 0 0 27 0 38 0 0 28 1 36 0 0 29 0 26 0 0 30 0 32 0 0 31 0 36 0 0 32 1 33 0 0 33 1 30 0 0 34 1 37 0 0 35 0 31 0 0 36 1 39 1 1 37 0 36 1 0 38 1 33 0 1 39 0 39 0 0 40 0 30 1 0 41 0 28 0 0 42 1 36 0 0 43 0 36 0 0 44 1 47 1 1 45 1 34 0 0 46 0 40 0 0 47 1 38 0 0 48 0 41 1 1 49 0 29 0 0 50 0 29 1 0 51 1 42 0 1 52 0 30 1 0 53 0 38 0 0 54 1 35 0 0 55 0 36 0 0 56 1 35 0 0 57 0 29 1 0 58 1 37 0 1 59 0 35 0 0 60 1 30 0 1 61 0 33 0 0 62 0 42 0 0 63 0 40 1 0 64 0 35 1 0 65 0 35 0 0 66 1 31 0 0 67 1 40 0 1 68 0 35 0 0 69 0 41 0 0 70 1 34 0 0 71 0 32 1 0 72 1 40 0 0 73 1 37 0 0 74 1 43 0 0 75 1 41 1 1 76 0 29 0 0 77 1 29 0 0 78 1 37 0 0 79 1 33 1 1 80 0 36 0 0 81 1 34 1 1 82 0 30 0 0 83 1 33 1 0 84 1 32 0 0 85 0 37 0 0 86 0 35 1 0 87 1 38 1 0 88 1 39 1 1 89 0 32 0 0 90 1 31 0 0 91 0 32 0 0 92 0 32 0 0 93 0 30 0 0 94 0 30 0 0 95 1 29 0 0 96 0 31 0 0 97 1 36 0 1 98 1 36 1 1 99 1 37 0 0 100 1 37 0 1 101 1 33 0 0 102 1 33 1 1 103 1 36 0 1 104 1 38 1 0 105 0 34 1 0 106 0 38 0 0 107 1 46 1 1 108 1 27 1 0 109 1 30 0 0 110 0 30 0 0 111 1 23 0 0 112 1 36 1 0 113 0 29 0 0 114 0 40 0 0 115 0 36 1 0 116 0 25 1 0 117 1 35 1 1 118 0 44 1 0 119 0 37 0 0 120 0 39 1 0 121 0 28 1 0 122 0 38 1 0 123 0 26 1 0 124 0 38 0 0 125 1 28 0 0 126 1 40 1 1 127 0 43 0 0 128 0 33 0 0 129 0 37 0 0 130 0 30 1 0 131 0 42 0 0 132 1 30 1 1 133 0 36 0 1 134 0 36 1 0 135 1 34 0 1 136 0 38 0 0 137 1 42 1 1 138 0 37 1 1 139 1 35 1 1 140 1 39 0 0 141 0 37 0 0 142 0 28 1 0 143 1 44 0 1 144 1 41 0 0 145 0 34 0 0 146 0 25 0 0 147 1 30 0 0 148 1 40 1 0 149 0 29 0 0 150 0 31 1 0 151 1 49 0 1 152 0 30 0 0 153 1 34 0 0 154 0 40 0 1 155 0 34 0 0 156 1 39 0 0 157 0 33 0 0 158 0 33 0 0 159 0 30 1 0 160 0 40 1 1 161 0 33 1 0 162 1 33 0 0 163 1 37 0 0 164 0 46 1 1 165 1 30 1 0 166 0 44 1 0 167 1 31 0 0 168 0 32 1 1 169 0 35 0 0 170 1 37 0 0 171 1 41 0 0 172 1 30 0 0 173 1 36 1 1 174 1 36 1 1 175 1 29 0 1 176 1 36 1 0 177 0 41 0 0 178 1 33 0 0 179 0 27 1 0 180 1 35 1 1 181 0 32 1 0 182 1 41 0 0 183 0 34 0 0 184 0 37 1 0 185 1 27 0 0 186 1 31 0 0 187 1 31 0 0 188 0 38 1 1 189 0 33 0 0 190 0 25 0 0 191 1 41 0 0 192 1 31 1 0 193 1 34 0 1 194 1 32 1 0 195 0 32 0 0 196 1 38 0 0 197 0 42 0 0 198 1 41 1 1 199 1 47 0 0 200 0 35 0 0 ;