This page is maintained by Alan Agresti. Last Updated: April 8, 2010.
(Our Department Home Page)
This site contains data sets that are not shown completely in the text. All but the first and third one are in the form of SAS programs for performing analyses such as shown in the text. The third data set can also be copied in a SAS data file.
A.1 Mental Impairment by Life Events and SES, in Table 3.3
mental ses life
1 1 1 1 1 9 1 1 4 1 1 3 1 0 2 1 1 0 1 0 1 1 1 3 1 1 3 1 1 7 1 0 1 1 0 2 2 1 5 2 0 6 2 1 3 2 0 1 2 1 8 2 1 2 2 0 5 2 1 5 2 1 9 2 0 3 2 1 3 2 1 1 3 0 0 3 1 4 3 0 3 3 0 9 3 1 6 3 0 4 3 0 3 4 1 8 4 1 2 4 1 7 4 0 5 4 0 4 4 0 4 4 1 8 4 0 8 4 0 9
A.2 SAS file for Example in Sections 9.2.3 and 10.2.2 on Arthritis Clinical Trial (data from Stu Lipsitz)
options ps=70 ls=70;
DATA arthritis;
input id sex ara age trt baseline extra y1 y2 y3 ;
if age=999 then delete;
if y1 = 1 or y1 = 2 then y1 = 1;
if y1 = 3 then y1 = 2;
if y1 = 4 or y1 = 5 then y1 = 3;
if y2 = 1 or y2 = 2 then y2 = 1;
if y2 = 3 then y2 = 2;
if y2 = 4 or y2 = 5 then y2 = 3;
if y3 = 1 or y3 = 2 then y3 = 1;
if y3 = 3 then y3 = 2;
if y3 = 4 or y3 = 5 then y3 = 3;
datalines;
1 2 3 54 2 2 3 4 5 5 2 1 3 41 1 3 4 4 4 4 3 2 2 48 2 3 3 3 4 4 4 2 2 40 1 3 4 4 3 4 5 2 1 29 2 3 4 4 3 4 6 2 3 43 2 2 3 3 2 3 7 2 2 55 2 4 4 3 2 3 8 1 2 49 1 3 3 3 2 2 9 2 3 64 1 2 2 2 2 1 10 2 2 51 2 3 4 3 3 4 11 2 3 46 2 4 4 4 4 4 12 2 1 59 1 4 4 4 5 4 13 1 3 58 2 3 2 3 2 3 14 2 3 49 1 4 4 5 5 5 15 2 1 48 2 3 3 3 3 3 16 2 3 60 1 4 3 4 3 3 17 2 3 62 1 3 3 4 3 4 18 1 3 65 2 2 2 2 2 2 19 1 2 56 1 4 3 4 3 3 20 2 2 64 1 3 3 4 3 3 22 2 3 42 1 3 3 3 3 4 23 2 3 59 2 4 3 3 1 2 24 2 3 61 2 4 4 4 4 4 25 2 3 59 1 4 4 4 5 5 27 2 2 29 1 3 4 3 3 4 28 1 2 41 2 3 3 4 4 4 29 2 3 65 1 3 3 4 3 4 30 1 2 22 1 3 5 3 3 4 31 2 3 55 2 3 4 4 3 4 32 2 3 62 2 3 3 3 4 4 33 2 2 31 2 3 2 2 3 5 34 2 1 43 1 2 3 4 4 3 35 2 2 59 2 2 3 4 4 4 36 2 2 57 1 3 2 3 3 3 37 1 3 48 2 3 5 4 5 4 58 2 1 61 1 3 2 1 3 4 59 2 2 35 2 3 3 3 4 4 60 2 2 61 1 2 3 4 4 3 61 2 1 56 2 4 3 3 3 3 62 2 3 56 2 3 3 3 4 2 63 2 2 60 1 4 3 3 3 4 64 2 1 52 2 3 3 4 3 4 65 2 1 55 1 3 3 4 3 4 66 2 2 64 2 3 2 3 3 4 68 1 1 50 1 2 2 3 2 5 77 1 3 57 1 4 3 4 4 4 78 2 3 50 2 3 3 3 4 4 79 2 1 56 1 3 2 3 2 3 81 2 3 56 2 2 3 2 4 3 82 2 3 64 2 4 3 4 3 4 83 2 3 64 1 3 2 3 3 2 85 2 3 33 1 3 3 4 4 4 87 1 3 48 2 3 3 2 4 3 88 1 3 35 2 5 5 5 5 5 89 2 3 55 1 2 3 2 1 2 90 2 1 50 2 4 4 3 3 4 92 2 1 45 2 3 2 3 3 2 93 1 1 46 2 4 4 3 5 5 94 2 3 55 1 4 4 4 2 4 96 2 3 51 2 4 2 2 2 2 99 1 1 55 1 5 5 5 5 4 100 2 3 34 2 2 2 4 5 5 101 2 1 53 2 3 2 2 2 3 103 2 3 61 1 3 4 2 3 3 105 2 2 51 1 4 4 2 5 4 107 2 3 57 1 4 1 3 3 4 108 2 3 42 1 5 5 5 5 5 110 1 2 52 2 4 3 3 4 5 112 1 3 44 2 3 3 4 4 3 113 1 3 31 1 4 4 4 4 4 116 2 2 64 1 3 3 4 4 4 133 2 2 58 1 1 1 2 1 2 134 2 3 50 1 2 2 2 3 2 136 2 2 47 1 4 3 4 2 3 137 1 3 57 2 3 3 3 3 3 138 2 3 46 1 3 1 3 2 2 139 2 2 39 1 4 3 1 3 1 140 2 2 64 2 3 3 4 3 3 142 2 2 57 2 1 1 3 3 5 144 2 2 31 1 2 2 2 3 3 145 2 2 58 1 2 3 4 4 4 146 2 2 60 2 4 4 3 4 5 147 2 2 36 2 4 4 4 4 5 148 1 3 60 2 2 2 2 4 3 149 2 3 63 1 3 3 4 4 4 150 2 2 58 1 1 3 1 3 3 151 2 2 31 2 3 2 4 2 3 153 2 2 61 1 2 3 3 4 3 155 2 2 49 2 2 2 3 3 4 158 1 2 44 1 3 4 4 3 4 159 1 2 64 2 3 3 3 3 3 160 2 3 41 1 3 3 4 4 4 161 1 2 31 2 3 3 3 3 3 163 2 2 63 1 3 3 4 3 4 167 2 2 41 2 2 2 2 2 2 169 2 2 62 1 3 3 3 1 3 170 2 2 62 2 3 3 3 3 4 171 2 3 54 2 2 3 4 3 3 172 2 2 42 1 3 3 3 3 4 175 2 2 33 2 1 5 5 5 5 176 2 2 40 1 2 3 4 3 4 177 2 2 50 2 3 3 3 4 4 179 2 2 25 2 3 2 2 3 3 180 2 2 46 2 3 3 4 4 4 181 2 2 66 1 1 3 3 3 3 184 2 2 32 2 2 3 4 4 5 186 2 2 54 1 4 4 4 4 4 187 2 2 58 1 3 3 3 4 3 188 2 2 33 2 3 3 4 4 4 189 1 2 62 1 2 2 3 2 2 190 2 3 62 2 3 3 3 2 3 192 2 2 26 2 3 3 3 4 3 193 2 2 61 1 3 4 3 4 4 208 2 2 46 1 3 2 3 3 3 209 2 1 64 2 3 3 3 3 3 211 2 2 21 1 4 3 4 5 4 213 1 1 50 2 3 3 3 4 5 214 2 2 45 1 2 2 2 3 3 215 2 1 50 2 2 3 5 5 4 218 2 3 60 1 1 2 5 3 4 219 1 1 52 1 4 5 4 4 3 220 2 3 42 2 3 3 3 3 2 221 2 2 65 2 2 3 3 2 3 223 2 3 48 2 5 4 4 5 5 229 1 1 60 2 3 3 2 3 2 230 2 2 62 1 3 4 3 4 4 231 2 2 63 1 3 3 3 3 3 233 2 1 62 2 3 3 4 3 3 234 2 3 54 2 2 2 3 3 2 236 2 1 37 1 2 4 2 3 3 237 2 2 53 2 2 3 4 5 5 238 2 2 63 1 3 3 3 1 2 239 2 3 46 2 4 4 4 4 5 240 2 2 55 2 3 4 5 5 5 241 2 1 36 1 3 3 2 3 3 243 2 2 55 1 3 3 3 3 3 244 2 1 33 1 4 4 5 4 5 245 2 2 51 1 4 3 3 2 3 246 1 3 60 2 2 3 3 3 3 247 2 3 61 2 2 3 4 3 5 248 1 3 41 1 4 3 5 5 5 249 2 1 63 2 3 3 3 2 3 250 2 2 55 2 2 2 4 4 4 265 2 2 34 1 3 3 3 2 2 266 2 2 62 1 3 2 3 2 2 267 2 3 53 2 2 3 3 3 3 268 2 3 24 1 3 3 3 3 3 269 2 3 48 1 2 1 2 2 3 271 1 3 52 1 2 3 3 3 3 272 2 2 55 2 1 2 4 3 4 273 2 1 60 2 2 3 3 3 2 275 2 2 36 2 1 1 3 2 3 276 2 1 28 2 3 3 2 4 2 277 2 2 28 1 3 3 3 2 2 278 2 3 34 2 2 2 3 3 3 279 1 2 26 1 3 2 2 3 3 281 2 2 56 2 4 5 4 4 5 282 2 2 57 1 3 3 3 3 3 287 2 3 55 1 1 1 2 2 1 288 2 2 55 2 1 1 3 3 3 289 1 2 46 1 4 3 4 4 4 290 2 3 59 1 3 3 3 3 3 291 2 3 51 1 3 4 4 4 4 292 2 3 54 1 2 2 4 3 3 294 2 3 58 1 3 3 1 4 3 298 1 3 47 1 1 4 4 2 3 299 2 2 54 2 3 3 3 3 3 300 2 3 60 2 2 2 4 4 5 301 2 2 45 2 3 2 3 2 3 302 2 2 33 1 4 3 5 2 5 303 2 2 57 1 2 1 1 2 1 306 2 3 38 1 2 3 4 4 4 308 2 3 51 2 4 3 4 4 3 309 2 3 55 2 2 2 3 3 3 310 1 3 33 1 1 1 2 1 1 311 1 3 57 1 2 2 3 3 3 313 2 3 37 1 2 2 2 4 2 315 1 2 55 2 2 2 3 3 4 317 2 3 63 2 2 1 3 2 3 318 2 3 45 2 5 4 4 4 5 319 2 3 54 1 3 3 3 3 4 320 1 3 63 1 2 1 2 2 1 321 1 3 53 2 4 4 4 4 4 322 2 2 49 2 3 4 4 1 3 325 2 3 64 2 3 3 4 4 3 326 2 3 45 1 3 2 4 3 4 328 2 3 55 1 3 3 3 4 4 331 2 9 65 2 3 4 3 5 1 333 2 3 58 2 4 4 4 3 3 335 2 3 49 2 4 4 4 4 3 381 2 3 38 2 3 2 2 1 4 383 2 2 57 2 2 3 3 4 4 384 2 3 57 1 2 3 3 2 2 385 1 2 51 2 2 2 3 3 2 386 2 3 58 2 2 2 3 4 3 387 1 2 65 1 3 3 2 3 2 388 1 2 50 2 3 3 3 2 4 389 2 2 52 1 5 4 4 4 4 390 2 2 65 2 5 4 4 4 5 391 2 3 60 1 5 3 5 4 3 392 1 3 58 1 1 2 2 1 2 393 1 1 58 1 3 3 2 2 2 394 2 1 55 2 4 4 4 3 2 395 1 1 59 1 3 2 2 4 5 396 2 9 54 2 4 3 4 4 4 397 1 2 44 2 4 4 4 4 4 398 1 1 44 1 2 3 4 4 3 399 1 2 47 2 3 4 3 3 4 400 2 2 55 1 4 4 4 3 3 401 2 1 53 1 4 4 4 4 4 402 2 3 999 2 1 1 2 1 2 403 2 2 63 2 3 2 4 4 4 404 2 2 57 1 3 3 4 4 4 405 2 3 55 1 4 1 4 4 5 406 2 1 58 2 1 1 1 2 4 407 2 3 62 1 1 3 2 1 1 408 2 3 53 2 2 1 2 2 3 409 2 2 25 1 3 2 2 3 2 410 2 2 58 2 2 2 3 2 2 421 1 3 62 1 1 2 3 2 2 422 1 3 32 2 3 2 3 3 3 423 1 2 63 2 5 4 4 4 4 424 1 3 39 2 4 5 4 5 5 426 1 1 53 1 3 2 2 3 2 427 1 2 62 2 2 3 3 4 4 428 1 1 61 1 3 3 4 3 4 429 1 3 64 2 5 3 4 4 5 430 1 2 57 2 3 2 3 3 4 431 1 2 35 2 3 4 4 5 4 432 1 1 54 1 3 4 3 3 3 433 1 1 63 2 3 3 3 4 3 434 1 2 50 1 2 2 3 2 2 435 1 2 35 1 4 4 2 4 3 436 1 2 35 2 3 3 4 4 4 437 1 2 63 1 3 3 2 3 3 438 1 1 55 1 3 2 3 2 3 457 2 1 38 1 3 4 4 4 4 458 2 1 22 2 1 2 3 3 5 459 2 1 28 2 4 4 5 4 5 460 2 2 51 1 3 3 3 3 4 463 1 2 62 1 2 1 1 2 2 465 2 1 33 1 4 4 4 2 3 466 2 2 36 2 4 4 4 5 5 467 1 1 31 1 2 2 2 3 3 468 2 1 39 2 3 3 3 2 2 469 2 1 33 1 2 3 1 2 2 470 2 1 52 2 3 2 3 4 3 471 1 2 27 2 3 3 3 4 4 472 2 2 41 2 2 2 2 2 2 473 2 1 31 1 3 2 2 2 2 474 2 2 48 2 2 2 3 2 4 487 2 1 48 2 3 3 1 4 3 488 2 2 42 2 3 3 3 3 4 490 2 2 50 1 2 2 3 3 2 491 2 2 53 2 4 3 3 5 4 497 1 2 63 2 3 3 3 3 4 498 1 2 55 1 5 5 5 3 4 499 2 2 37 1 4 4 4 3 5 500 2 1 39 2 3 3 3 3 3 501 1 2 35 2 3 3 3 4 3 502 1 3 35 2 3 3 3 3 3 503 1 2 65 1 3 2 2 3 3 504 2 1 58 2 4 4 4 3 5 505 2 3 40 1 3 2 4 1 3 506 2 3 33 1 2 1 2 2 1 508 2 3 33 2 4 3 4 5 5 509 2 1 40 2 3 2 3 2 2 510 1 2 60 1 2 2 2 2 2 511 2 3 59 1 3 2 3 4 4 512 2 3 62 2 3 1 3 3 3 516 2 1 54 2 3 2 3 2 4 517 2 2 55 2 1 1 4 4 3 518 2 2 45 1 3 4 2 4 4 519 1 3 65 2 2 3 3 3 3 520 2 2 51 1 3 2 3 4 3 522 1 2 56 1 3 2 3 3 3 523 1 2 64 2 2 2 3 3 2 524 1 2 36 2 3 1 3 3 3 525 2 3 61 1 2 3 2 3 3 526 1 2 46 1 2 2 3 3 2 527 1 2 60 1 3 3 4 4 4 528 2 2 56 1 2 2 3 3 3 529 1 2 63 1 3 4 3 3 4 530 2 3 56 1 3 2 3 3 3 532 2 3 60 2 1 2 2 2 3 533 2 3 58 2 1 2 3 2 1 534 2 2 61 1 1 1 2 1 2 535 2 3 64 1 2 1 1 2 2 536 2 2 55 1 2 2 2 3 3 109 2 3 57 2 4 4 4 4 . 135 2 2 63 2 2 2 3 . . 143 1 2 42 2 1 1 3 . 2 157 2 2 37 1 3 2 3 . 3 165 2 2 52 2 2 3 4 3 . 270 2 1 52 1 3 3 . 3 3 280 1 2 47 2 3 3 . . . 312 2 3 63 2 1 2 3 2 . 330 2 3 60 1 3 3 2 2 2 334 2 2 63 2 2 3 . 2 2 382 1 3 33 1 3 2 1 3 . 425 1 2 57 1 4 3 3 2 . 489 2 1 41 2 3 3 3 . . 521 2 3 29 2 3 3 2 . .;
A.3 Data file for Example in Section 9.1 on Eye Disease Risk Factors (from John Williamson). The variables include DIAB for "Duration of diabetes" (in years), GH for "glycosylated haemoglobin level", DBP for "Diastolic blood pressure", SEX for "Gender" (male = 0, female = 1), PROT for "Proteinuria" (absent = 0, present = 1), and RME (right eye) and LME (left eye) for "Macular oedema" (absent = 0, present = 1).
Obs rrl lrl rl rme lme rre lre riop liop age diab gh 1 4.0 2 28 0 0 -0.250 -1.125 14 14 29.9 10.3 13.7 2 3.0 4 27 0 0 -2.500 -3.125 16 19 18.6 9.9 13.5 3 4.0 4 26 0 0 -0.750 -1.500 10 11 28.7 15.6 13.8 4 1.0 1 40 0 0 1.000 1.125 12 10 20.1 34.5 8.4 5 4.0 3 27 0 0 -1.000 -1.250 17 17 29.3 22.9 12.8 6 4.0 4 26 0 1 0.750 0.750 8 8 21.8 26.0 13.0 7 2.0 3 32 0 0 0.000 0.000 20 15 21.3 13.8 11.1 8 4.0 3 27 1 0 -0.750 -1.375 16 19 27.5 31.1 11.3 9 1.0 1 40 0 0 0.125 0.000 15 12 28.7 2.6 15.1 10 4.0 3 27 0 0 1.500 1.750 14 16 27.1 25.8 17.2 11 4.0 3 27 1 0 -1.000 -0.750 19 18 18.0 39.6 13.4 12 6.0 6 13 0 0 -4.125 -4.500 17 16 5.8 31.4 9.2 13 1.0 1 40 0 0 0.250 0.625 20 17 11.2 7.3 14.6 14 1.0 1 40 0 0 0.000 0.000 15 12 22.9 8.3 15.4 15 6.0 7 6 1 1 0.500 10.000 18 18 29.6 27.6 12.8 16 4.0 3 27 1 0 -1.250 -1.000 16 16 23.7 18.4 14.0 17 1.0 1 40 0 0 -4.125 -3.875 16 18 17.1 4.4 11.7 18 1.0 1 40 0 0 -0.500 0.000 20 20 24.9 15.8 11.4 19 4.0 4 26 0 0 -0.250 -0.250 15 15 3.7 25.5 12.3 20 1.0 1 40 0 0 -0.625 -0.125 23 19 19.0 4.5 12.6 21 1.0 1 40 0 0 -1.375 -1.250 15 15 25.3 5.3 9.9 22 4.0 4 26 0 0 -7.875 -7.125 15 15 10.8 13.5 10.7 23 2.0 2 35 0 0 -1.000 -0.875 13 14 3.3 17.4 12.8 24 1.0 1 40 0 0 -0.250 -3.250 20 20 26.8 3.3 10.7 25 3.0 3 31 0 0 0.250 -0.500 16 15 21.4 12.2 15.2 26 6.0 6 13 1 0 -0.375 -0.125 17 16 4.7 19.5 9.8 27 1.0 1 40 0 0 0.000 0.000 16 18 14.3 4.0 8.5 28 2.0 1 37 0 0 0.000 -0.125 15 14 14.9 16.4 10.2 29 6.0 6 13 0 0 0.000 -0.250 19 19 12.2 22.4 12.3 30 6.0 6 13 1 0 -2.250 -1.750 16 13 10.7 30.5 14.3 31 3.0 2 32 0 0 -5.375 -6.375 18 15 7.8 15.5 11.8 Obs sbp dbp bmi pr sex prot dose resid rerl lerl rl2 id 1 123 89 3.391 43 1 0 2 1 3 2 5 1 2 122 90 3.336 43 1 0 2 2 2 3 5 2 3 118 85 3.529 51 1 1 1 1 3 3 6 3 4 186 99 4.455 42 1 0 2 1 1 1 1 4 5 135 84 4.486 45 1 0 1 2 3 2 5 5 6 128 70 3.067 47 1 0 1 2 3 3 6 6 7 106 71 3.506 43 1 0 1 2 2 2 3 7 8 123 74 3.502 37 1 0 1 1 3 2 5 8 9 108 73 2.739 32 1 0 1 2 1 1 1 9 10 122 74 4.885 45 1 0 1 2 3 2 5 10 11 142 80 4.325 45 2 0 1 1 3 2 5 11 12 162 66 3.156 44 2 0 2 1 4 4 10 12 13 121 77 3.531 49 2 0 1 1 1 1 1 13 14 127 80 2.944 40 2 0 1 1 1 1 1 14 15 205 83 3.516 55 2 1 2 2 4 4 10 15 16 127 85 3.003 55 1 0 1 1 3 2 5 16 17 121 73 3.240 30 1 0 1 1 1 1 1 17 18 120 87 4.380 49 1 0 1 1 1 1 1 18 19 148 103 3.477 42 1 0 1 1 3 3 6 19 20 113 73 3.477 36 1 0 1 1 1 1 1 20 21 114 78 3.052 44 1 0 1 1 1 1 1 21 22 110 71 3.218 35 2 0 2 1 3 3 6 22 23 127 66 3.443 40 1 0 2 1 2 2 3 23 24 148 108 6.397 36 2 0 2 2 1 1 1 24 25 125 82 3.806 39 1 0 1 1 2 2 3 25 26 125 84 3.491 41 2 0 2 1 4 4 10 26 27 113 75 3.180 36 1 0 2 1 1 1 1 27 28 130 78 2.952 39 2 0 2 2 2 1 2 28 29 124 79 2.380 53 2 0 1 1 4 4 10 29 30 164 97 3.487 43 2 1 1 1 4 4 10 30 31 102 65 2.866 37 2 0 2 1 2 2 3 31 Obs rrl lrl rl rme lme rre lre riop liop age diab gh 32 4.0 4 26 0 0 -0.125 -1.125 22 21 28.1 29.5 11.7 33 4.0 3 27 0 0 -1.750 -1.250 15 19 18.3 15.6 11.9 34 3.0 3 31 0 0 0.250 0.750 22 21 22.4 5.9 13.4 35 4.0 4 26 0 0 -1.250 -0.750 15 15 14.1 16.6 14.7 36 2.0 2 35 0 0 -4.625 -4.375 14 17 9.8 11.0 14.6 37 3.0 3 31 0 0 -3.625 -3.750 13 13 17.8 27.4 13.7 38 1.0 1 40 0 0 -1.000 -1.250 20 20 13.9 16.5 9.6 39 4.0 4 26 0 0 -0.875 -0.625 18 18 16.1 29.5 9.4 40 3.0 4 27 0 0 -1.750 -2.375 24 18 20.2 15.7 16.7 41 1.0 1 40 0 0 -3.750 -3.625 16 16 28.8 1.9 7.3 42 4.0 2 28 0 0 12.625 13.625 14 15 26.8 29.6 14.3 43 3.0 4 27 0 0 -3.250 -3.500 19 19 16.1 7.4 12.0 44 3.0 2 32 0 0 -2.750 -4.125 10 12 29.1 19.8 10.6 45 2.0 2 35 0 0 -0.250 -0.250 13 15 17.9 7.2 14.8 46 2.0 2 35 0 0 -2.000 -1.875 9 10 7.2 23.0 13.7 47 6.0 7 6 0 0 -2.000 -2.250 20 20 12.1 14.4 23.3 48 1.0 2 37 0 0 -0.250 -0.500 15 13 28.3 5.9 14.4 49 3.0 4 27 0 0 -4.750 -4.250 17 16 18.5 17.6 11.9 50 6.0 6 13 0 0 -0.750 0.125 15 15 27.2 6.3 11.8 51 4.0 3 27 0 0 0.375 0.500 20 20 17.1 26.4 7.7 52 1.0 2 37 0 0 1.000 0.000 14 14 10.7 2.5 14.3 53 3.0 4 27 0 0 0.000 0.000 15 15 5.7 14.6 10.1 54 1.5 1 39 0 0 0.000 -0.250 17 13 14.1 3.3 8.6 55 1.0 1 40 0 0 0.000 0.000 18 19 12.7 6.1 12.6 56 1.0 1 40 0 0 0.125 0.125 20 20 8.3 2.0 14.2 57 1.0 1 40 0 0 0.125 -0.250 17 16 12.4 3.1 11.5 58 1.0 1 40 0 0 0.000 0.250 21 21 10.8 6.1 15.0 59 1.0 1 40 0 0 -1.875 -1.375 13 11 7.6 5.1 13.3 60 1.5 2 36 0 0 -0.375 0.375 10 10 12.0 3.3 21.5 61 2.0 2 35 0 0 4.375 2.750 12 13 11.0 7.6 18.7 62 1.5 1 39 0 0 0.000 0.000 13 13 6.9 8.0 17.3 Obs sbp dbp bmi pr sex prot dose resid rerl lerl rl2 id 32 153 89 4.376 50 1 0 2 1 3 3 6 32 33 126 80 3.561 47 1 0 2 1 3 2 5 33 34 115 91 3.640 50 2 0 2 2 2 2 3 34 35 114 75 3.016 40 1 0 1 1 3 3 6 35 36 113 82 3.490 51 2 0 2 2 2 2 3 36 37 165 85 3.358 49 2 0 1 2 2 2 3 37 38 129 80 3.398 31 1 0 1 1 1 1 1 38 39 152 73 3.241 42 2 0 1 1 3 3 6 39 40 141 88 3.349 48 2 1 1 1 2 3 5 40 41 112 75 3.663 41 1 0 1 1 1 1 1 41 42 158 92 3.285 52 2 0 1 1 3 2 5 42 43 106 87 4.185 51 2 0 1 1 2 3 5 43 44 122 60 3.842 37 2 0 1 2 2 2 3 44 45 120 78 3.616 45 1 0 2 1 2 2 3 45 46 93 78 3.267 39 2 0 1 1 2 2 3 46 47 118 89 3.032 43 2 1 1 1 4 4 10 47 48 103 73 2.801 48 2 0 1 1 1 2 2 48 49 132 87 3.036 32 2 0 1 1 2 3 5 49 50 125 84 3.987 39 1 0 1 1 4 4 10 50 51 139 93 3.459 36 1 1 2 1 3 2 5 51 52 108 78 2.542 32 2 0 2 1 1 2 2 52 53 116 81 3.013 32 1 1 2 1 2 3 5 53 54 96 82 2.699 40 2 0 1 1 2 1 2 54 55 106 53 3.606 38 2 0 2 2 1 1 1 55 56 95 63 2.446 37 1 0 2 1 1 1 1 56 57 124 78 2.605 42 2 0 2 1 1 1 1 57 58 115 81 2.738 39 1 0 2 2 1 1 1 58 59 100 75 3.210 52 2 0 2 2 1 1 1 59 60 121 87 2.671 46 2 0 1 2 2 2 3 60 61 109 83 3.298 39 2 0 2 2 2 2 3 61 62 104 74 2.388 38 1 0 2 1 2 1 2 62 Obs rrl lrl rl rme lme rre lre riop liop age diab gh 63 1.0 1 40 0 0 0.000 -0.375 17 16 2.9 12.4 11.7 64 1.0 1 40 0 0 0.000 0.000 16 13 8.3 2.2 17.2 65 3.0 4 27 0 0 -1.250 -1.250 17 17 6.0 14.4 9.1 66 1.0 1 40 0 0 0.250 0.250 18 18 6.6 6.8 11.3 67 1.0 1 40 0 0 0.000 0.000 19 19 13.6 1.5 12.1 68 2.0 1 37 0 0 0.125 0.125 21 20 9.2 2.7 16.0 69 1.0 1 40 0 0 -2.750 -2.875 20 20 9.2 4.1 12.3 70 1.0 2 37 0 0 0.000 0.000 19 18 1.7 8.7 12.2 71 4.0 3 27 0 0 -1.250 -1.375 15 14 27.1 13.5 14.9 72 2.0 2 35 0 0 0.000 0.000 20 22 16.8 4.0 13.5 73 4.0 5 21 0 1 -1.375 -1.125 24 26 12.0 13.6 12.4 74 1.0 1 40 0 0 -6.000 -4.750 14 14 25.0 2.6 18.2 75 4.0 4 26 0 0 0.000 -0.250 13 14 17.5 11.4 13.2 76 1.0 1 40 0 0 0.000 0.000 17 16 21.5 5.9 13.4 77 1.0 1 40 0 0 -7.750 -7.625 21 21 17.3 6.9 9.6 78 1.0 2 37 0 0 -2.500 -3.000 18 21 28.0 9.5 12.2 79 3.0 3 31 0 0 3.750 3.250 16 16 28.7 33.5 10.3 80 2.0 2 35 0 0 -0.625 -0.750 18 20 13.0 13.5 12.8 81 4.0 4 26 0 0 -1.875 -0.875 15 12 4.8 20.5 13.7 82 4.0 3 27 0 0 -2.750 -4.750 14 15 19.2 19.4 11.2 83 6.0 4 15 1 0 -4.500 -5.125 18 20 21.6 26.6 13.5 84 3.0 2 32 0 0 0.375 0.375 24 26 17.1 29.6 12.6 85 8.0 6 3 1 1 0.750 -0.250 5 16 18.8 12.1 9.6 86 4.0 4 26 0 0 -1.500 -0.500 22 22 21.3 16.4 9.7 87 2.0 4 28 0 0 0.500 0.625 12 12 28.5 21.5 11.1 88 4.0 4 26 0 0 -0.250 0.250 27 30 12.2 11.5 18.5 89 4.0 4 26 0 0 -0.375 -0.250 18 18 21.2 11.5 15.8 90 6.0 6 13 0 0 0.000 0.250 15 18 4.6 28.7 11.0 91 3.0 6 16 0 0 -3.750 -3.625 14 15 12.8 33.9 10.1 92 2.0 2 35 0 0 0.250 0.125 15 14 28.0 5.6 13.7 93 4.0 4 26 0 0 -2.125 -2.375 14 12 3.7 19.6 16.8 Obs sbp dbp bmi pr sex prot dose resid rerl lerl rl2 id 63 127 79 3.510 38 2 0 2 1 1 1 1 63 64 95 75 2.319 39 2 0 2 2 1 1 1 64 65 131 83 3.291 39 1 0 2 1 2 3 5 65 66 101 84 2.539 39 2 0 1 1 1 1 1 66 67 114 68 3.477 33 1 0 2 1 1 1 1 67 68 125 85 3.121 64 2 0 2 1 2 1 2 68 69 100 82 3.008 44 2 0 2 1 1 1 1 69 70 108 81 2.438 39 1 0 1 1 1 2 2 70 71 124 74 3.158 32 2 0 1 1 3 2 5 71 72 115 80 3.496 45 1 0 1 1 2 2 3 72 73 149 104 3.483 46 1 1 1 1 3 3 6 73 74 128 79 3.646 45 2 1 1 1 1 1 1 74 75 111 94 2.873 49 1 1 2 1 3 3 6 75 76 113 78 3.585 39 1 0 1 1 1 1 1 76 77 134 70 3.035 36 2 0 1 1 1 1 1 77 78 110 64 3.715 38 1 0 1 1 1 2 2 78 79 140 68 3.371 42 2 0 1 1 2 2 3 79 80 115 81 3.211 40 2 0 2 1 2 2 3 80 81 153 101 3.526 62 2 0 2 1 3 3 6 81 82 121 66 3.237 32 1 0 1 1 3 2 5 82 83 161 94 3.702 41 2 1 1 1 4 3 9 83 84 134 74 2.863 33 2 0 1 1 2 2 3 84 85 204 100 4.868 40 2 1 1 2 4 4 10 85 86 133 92 3.764 41 1 0 2 1 3 3 6 86 87 130 77 4.111 51 2 1 2 1 2 3 5 87 88 100 80 2.042 51 2 1 2 1 3 3 6 88 89 113 89 3.142 45 2 0 1 1 3 3 6 89 90 115 70 3.126 43 1 1 1 1 4 4 10 90 91 90 69 3.250 41 2 0 2 1 2 4 8 91 92 105 77 3.599 44 2 0 1 1 2 2 3 92 93 143 104 4.123 38 2 1 1 1 3 3 6 93 Obs rrl lrl rl rme lme rre lre riop liop age diab gh 94 1.0 2 37 0 0 0.000 0.000 20 20 24.5 7.5 15.0 95 1.0 1 40 0 0 0.000 0.125 17 17 8.5 3.2 11.3 96 4.0 4 26 0 0 -2.625 -1.500 22 22 13.0 16.5 9.3 97 4.0 4 26 0 0 -1.000 -1.000 14 13 8.1 16.8 13.2 98 1.0 1 40 0 0 -1.750 -2.000 14 14 25.4 5.9 14.0 99 1.5 1 39 0 0 -0.125 0.000 18 18 26.0 3.4 8.6 100 2.0 2 35 0 0 -0.125 0.000 17 18 16.1 23.5 8.9 101 4.0 3 27 0 0 0.125 0.625 12 12 5.0 26.4 9.1 102 1.0 1 40 0 0 -0.875 -1.000 21 17 12.3 12.5 13.9 103 2.0 1 37 0 0 -0.750 -0.625 15 13 11.0 7.5 7.8 104 3.0 3 31 0 0 -2.500 -2.500 16 15 15.0 8.6 13.0 105 1.0 1 40 0 0 -0.250 -0.125 20 20 18.8 7.1 9.5 106 4.0 4 26 0 0 -2.375 -2.625 16 16 1.7 21.3 9.1 107 4.0 2 28 0 0 -1.500 -1.500 18 18 11.8 14.4 15.0 108 2.0 2 35 0 0 -3.625 -2.750 21 22 12.0 8.2 15.6 109 2.0 2 35 0 0 -0.125 0.000 20 19 9.9 7.4 12.5 110 1.0 1 40 0 0 0.125 -0.125 12 12 8.4 5.9 9.8 111 1.0 1 40 0 0 0.000 0.000 16 14 14.0 3.2 12.0 112 1.0 3 34 0 0 -0.500 -0.500 13 13 3.7 10.5 12.3 113 3.0 3 31 0 0 -0.250 0.000 21 20 11.5 10.4 10.5 114 4.0 3 27 0 0 -1.625 -1.750 19 21 6.1 10.7 14.2 115 3.0 3 31 0 0 1.250 1.375 18 17 6.3 12.5 10.3 116 1.0 1 40 0 0 -2.625 -2.375 19 19 14.0 5.8 10.5 117 1.0 1 40 0 0 0.000 0.000 14 12 14.8 2.6 12.9 118 1.0 1 40 0 0 -3.625 0.625 16 15 9.4 4.4 12.6 119 1.0 1 40 0 0 0.000 -0.500 13 13 5.4 7.8 16.5 120 1.0 1 40 0 0 0.000 0.000 18 18 6.5 4.5 13.2 121 2.0 3 32 0 0 0.000 0.000 15 15 4.4 12.4 12.5 122 3.0 2 32 0 0 -3.250 -2.875 22 21 9.6 10.6 15.7 123 1.0 1 40 0 0 0.000 0.000 15 14 15.8 3.8 11.6 124 4.0 6 15 0 0 12.000 0.750 12 12 19.9 30.8 12.7 Obs sbp dbp bmi pr sex prot dose resid rerl lerl rl2 id 94 113 72 3.050 43 2 0 1 1 1 2 2 94 95 100 60 2.814 37 2 0 2 1 1 1 1 95 96 128 84 3.412 41 1 1 2 1 3 3 6 96 97 118 70 3.115 41 2 0 2 2 3 3 6 97 98 119 83 3.346 38 2 0 1 1 1 1 1 98 99 103 69 3.263 41 2 0 2 1 2 1 2 99 100 100 72 2.785 48 2 0 2 1 2 2 3 100 101 109 65 3.370 45 2 0 2 2 3 2 5 101 102 103 81 3.253 47 2 0 1 1 1 1 1 102 103 114 70 3.469 47 2 0 2 1 2 1 2 103 104 122 78 3.502 42 2 1 2 1 2 2 3 104 105 109 69 3.195 50 1 0 2 1 1 1 1 105 106 115 75 3.356 36 2 0 2 1 3 3 6 106 107 139 90 3.287 33 2 1 1 1 3 2 5 107 108 119 87 3.418 41 2 0 2 1 2 2 3 108 109 122 73 3.366 39 1 0 1 1 2 2 3 109 110 106 81 3.123 44 1 0 1 2 1 1 1 110 111 117 79 3.138 40 1 0 1 1 1 1 1 111 112 109 60 3.119 38 2 0 1 1 1 2 2 112 113 121 76 3.252 60 2 0 1 1 2 2 3 113 114 166 108 3.601 44 2 1 2 2 3 2 5 114 115 124 92 3.385 48 1 1 2 1 2 2 3 115 116 131 88 3.529 43 1 1 1 1 1 1 1 116 117 124 78 3.448 33 1 0 2 1 1 1 1 117 118 107 82 2.832 50 2 0 2 1 1 1 1 118 119 109 61 2.405 44 1 0 2 2 1 1 1 119 120 109 89 2.381 61 2 0 2 1 1 1 1 120 121 109 74 2.599 37 1 0 2 1 2 2 3 121 122 130 93 3.766 47 2 0 1 1 2 2 3 122 123 121 71 3.312 44 1 0 2 1 1 1 1 123 124 102 66 2.683 40 2 0 1 1 3 4 9 124 Obs rrl lrl rl rme lme rre lre riop liop age diab gh 125 3 4 27 0 0 -4.250 -5.250 15 15 11.0 16.1 9.4 126 3 3 31 0 0 0.000 0.125 16 17 15.6 9.5 13.3 127 3 3 31 0 0 0.000 0.000 15 16 20.4 9.9 13.4 128 4 3 27 0 0 0.125 -0.125 16 13 25.1 20.4 13.9 129 3 4 27 0 0 -7.500 -8.125 12 12 2.5 24.4 13.9 130 4 4 26 0 0 0.000 0.000 16 15 12.6 11.6 17.9 131 3 3 31 0 0 -3.750 -3.375 16 16 14.8 16.4 12.5 132 1 1 40 0 0 0.500 0.500 18 19 26.6 4.5 15.7 133 2 4 28 0 0 -3.250 -3.500 15 15 5.5 21.6 19.7 134 3 3 31 0 0 -0.875 -3.125 14 16 17.7 10.6 10.4 135 3 4 27 0 0 -5.750 -5.250 18 16 12.9 10.4 9.8 136 3 4 27 0 0 -4.500 -4.750 18 20 27.4 31.6 12.1 137 1 1 40 0 0 -3.125 -3.750 14 14 9.9 18.1 12.5 138 2 2 35 0 0 -1.125 -1.750 17 17 12.7 11.6 13.2 139 2 2 35 0 0 0.000 0.000 16 14 17.6 14.7 12.5 140 6 6 13 1 1 -8.500 -2.875 13 15 11.5 24.4 10.7 141 2 2 35 0 0 -1.750 -2.000 18 18 15.0 20.6 13.2 142 3 6 16 0 0 -1.750 -1.250 28 31 11.8 31.6 11.8 143 2 2 35 0 0 -0.500 0.250 10 10 22.3 9.0 17.7 144 2 2 35 0 0 -2.000 -1.000 19 17 16.9 9.5 9.7 145 1 1 40 0 0 -1.875 0.125 13 14 25.2 10.6 12.5 146 1 1 40 0 0 0.250 0.125 18 17 29.7 3.2 19.9 147 3 2 32 0 0 -0.875 -1.125 12 12 16.1 10.5 13.6 148 2 2 35 0 0 -1.500 -1.125 15 12 25.7 7.2 12.7 149 2 2 35 0 0 -2.500 -2.875 15 15 7.5 26.2 11.6 150 2 2 35 0 0 -2.500 -1.625 14 14 12.9 19.3 10.6 151 3 3 31 0 1 2.500 3.125 15 15 15.2 31.5 11.3 152 1 1 40 0 0 0.000 -0.500 18 18 23.7 8.4 10.7 153 3 3 31 0 0 0.125 0.000 20 20 4.9 13.7 12.2 154 1 1 40 0 0 9.125 9.750 11 13 15.0 14.4 9.6 155 4 3 27 0 0 -3.000 -3.375 18 16 17.1 12.5 14.1 Obs sbp dbp bmi pr sex prot dose resid rerl lerl rl2 id 125 121 75 4.363 36 2 0 1 1 2 3 5 125 126 124 96 4.037 41 2 0 1 1 2 2 3 126 127 125 76 3.419 39 1 0 1 2 2 2 3 127 128 117 74 3.169 46 2 0 1 1 3 2 5 128 129 116 72 2.632 40 2 0 1 1 2 3 5 129 130 103 75 3.539 44 2 0 1 1 3 3 6 130 131 115 89 3.118 47 1 0 1 1 2 2 3 131 132 131 95 2.877 44 1 0 1 1 1 1 1 132 133 108 82 3.039 41 2 0 1 1 2 3 5 133 134 104 72 2.802 30 2 0 2 1 2 2 3 134 135 125 91 4.025 42 1 0 2 1 2 3 5 135 136 131 77 3.148 42 2 0 2 1 2 3 5 136 137 103 63 3.020 30 2 0 1 1 1 1 1 137 138 130 95 3.113 33 2 0 1 1 2 2 3 138 139 125 89 5.382 59 2 0 1 1 2 2 3 139 140 135 91 2.942 49 2 1 2 1 4 4 10 140 141 137 89 3.011 43 2 0 1 1 2 2 3 141 142 110 77 3.312 46 2 0 2 1 2 4 8 142 143 102 67 3.215 48 1 0 1 2 2 2 3 143 144 112 90 3.419 36 2 0 1 1 2 2 3 144 145 109 71 3.212 43 2 0 1 1 1 1 1 145 146 131 74 3.346 38 1 0 0 1 1 1 1 146 147 90 56 3.210 46 2 0 1 1 2 2 3 147 148 137 107 3.773 60 2 0 1 1 2 2 3 148 149 91 71 3.174 38 2 1 1 1 2 2 3 149 150 111 77 3.756 34 2 0 2 1 2 2 3 150 151 116 77 3.363 38 2 0 1 1 2 2 3 151 152 105 74 3.829 35 1 0 1 1 1 1 1 152 153 102 67 3.209 35 1 0 1 1 2 2 3 153 154 149 90 3.438 48 2 0 1 1 1 1 1 154 155 117 89 4.113 45 2 0 1 1 3 2 5 155 Obs rrl lrl rl rme lme rre lre riop liop age diab gh 156 3 4 27 0 0 -0.375 -0.250 14 12 10.5 10.6 11.0 157 1 1 40 0 0 0.500 -0.375 19 18 9.0 4.2 18.8 158 1 1 40 0 0 0.000 0.000 20 16 12.7 9.5 9.8 159 3 3 31 0 0 -1.000 -1.625 11 15 12.5 12.7 9.4 160 6 3 16 0 0 0.625 1.125 15 13 22.6 35.4 11.3 161 3 2 32 0 0 0.000 0.000 16 16 12.1 17.1 13.0 162 1 1 40 0 0 -0.750 -0.875 20 20 22.5 9.5 9.1 163 2 2 35 0 0 -2.500 -1.500 18 17 14.7 26.3 11.9 164 1 1 40 0 0 -1.625 -2.125 20 20 20.7 8.4 10.8 165 4 6 15 1 1 0.000 -0.625 16 17 25.8 17.8 13.7 166 4 6 15 0 0 -7.125 -6.875 14 14 19.9 8.3 13.7 167 4 4 26 0 0 1.375 14.500 14 17 28.2 24.0 11.6 168 4 4 26 0 0 3.000 3.750 12 14 17.1 21.3 10.0 169 2 2 35 0 0 0.000 0.000 13 11 15.4 6.6 11.3 170 4 4 26 0 0 -3.875 -4.125 11 13 17.4 7.6 17.2 171 3 3 31 0 0 -0.250 -0.375 20 19 13.1 9.6 11.4 172 1 1 40 0 0 0.125 -0.500 18 17 21.5 2.9 10.6 173 1 1 40 0 0 0.000 0.000 14 14 5.9 4.2 13.0 174 2 2 35 0 0 -5.875 -4.625 15 15 13.8 11.6 12.5 175 4 4 26 0 1 -5.000 -5.750 24 22 14.8 25.3 10.0 176 2 2 35 0 0 -0.625 0.000 16 16 20.6 6.4 11.3 177 1 1 40 0 0 -3.375 -2.750 16 16 20.1 2.0 9.6 178 4 4 26 0 0 -4.000 -4.375 14 15 14.9 16.6 16.2 179 7 7 5 0 1 -3.750 -4.000 16 15 11.3 15.5 20.8 180 2 2 35 0 0 -1.375 -0.750 22 22 18.5 3.3 9.7 181 4 4 26 0 0 -2.875 -4.250 17 15 9.6 20.5 15.4 182 3 3 31 0 0 0.000 0.000 16 16 9.8 15.6 18.6 183 4 4 26 0 0 -0.500 0.375 18 20 20.4 18.6 12.1 184 2 1 37 0 0 0.000 0.000 17 16 24.6 6.9 12.1 185 4 4 26 0 0 -4.500 -3.500 16 15 14.9 14.4 9.8 186 4 4 26 1 0 -1.875 -1.500 12 13 14.5 10.4 14.0 Obs sbp dbp bmi pr sex prot dose resid rerl lerl rl2 id 156 128 93 3.687 49 2 1 2 1 2 3 5 156 157 107 71 3.335 40 1 0 2 1 1 1 1 157 158 113 69 4.058 47 2 1 1 1 1 1 1 158 159 130 82 3.621 32 1 1 1 1 2 2 3 159 160 141 77 3.049 39 2 0 1 1 4 2 8 160 161 121 70 3.205 37 2 0 1 1 2 2 3 161 162 124 92 3.300 40 1 0 1 1 1 1 1 162 163 101 67 3.439 44 2 0 1 1 2 2 3 163 164 128 81 3.380 50 2 0 2 1 1 1 1 164 165 201 99 4.837 38 2 1 2 2 3 4 9 165 166 141 88 3.119 50 1 0 2 1 3 4 9 166 167 134 82 3.428 39 1 0 2 2 3 3 6 167 168 124 75 3.089 50 2 0 2 1 3 3 6 168 169 110 84 3.327 33 2 0 1 1 2 2 3 169 170 120 85 3.127 35 1 0 1 2 3 3 6 170 171 154 86 3.074 58 1 1 1 1 2 2 3 171 172 124 83 3.013 38 1 1 2 1 1 1 1 172 173 100 73 2.502 44 2 0 2 1 1 1 1 173 174 115 86 3.118 50 2 0 1 1 2 2 3 174 175 125 83 3.217 48 1 1 1 1 3 3 6 175 176 147 84 3.109 44 1 0 2 1 2 2 3 176 177 124 92 3.088 38 2 0 1 1 1 1 1 177 178 122 87 3.546 42 1 1 1 1 3 3 6 178 179 159 108 2.846 46 1 1 1 1 4 4 10 179 180 120 82 3.794 48 2 0 1 2 2 2 3 180 181 123 76 3.092 56 1 1 1 1 3 3 6 181 182 134 95 2.733 49 2 0 2 1 2 2 3 182 183 119 91 2.931 42 1 0 1 1 3 3 6 183 184 115 71 3.060 34 1 0 1 1 2 1 2 184 185 131 63 3.651 34 1 0 2 1 3 3 6 185 186 114 65 2.826 44 1 0 1 1 3 3 6 186 Obs rrl lrl rl rme lme rre lre riop liop age diab gh 187 1 1.0 40 0 0 0.000 -0.750 11 13 17.9 3.0 11.6 188 1 1.0 40 0 0 0.000 0.000 14 12 23.1 10.5 14.3 189 1 2.0 37 0 0 0.500 0.375 15 15 29.1 8.0 13.5 190 2 2.0 35 0 0 -4.125 -4.750 18 18 20.0 7.0 11.4 191 1 1.0 40 0 0 0.125 0.000 15 15 29.8 2.2 10.0 192 1 1.0 40 0 0 -0.625 0.000 17 17 17.7 4.2 12.9 193 2 3.0 32 0 0 0.000 -0.125 13 15 1.7 22.5 11.2 194 4 4.0 26 0 0 -2.375 -2.500 19 20 19.7 24.4 10.5 195 4 3.0 27 0 0 -1.875 -1.750 16 16 3.7 15.5 10.4 196 1 1.0 40 0 0 -2.125 -2.375 11 13 8.7 4.5 13.8 197 1 1.0 40 0 0 -0.750 -0.500 21 21 9.9 9.0 12.2 198 2 2.0 35 0 0 -1.375 -1.125 19 18 10.5 6.0 18.0 199 1 1.0 40 0 0 -0.125 -0.125 16 16 5.1 2.8 12.5 200 1 1.0 40 0 0 -1.750 -1.500 16 16 10.0 1.5 15.3 201 1 1.0 40 0 0 0.125 0.000 15 15 12.9 4.0 18.5 202 1 1.0 40 0 0 -1.000 -1.000 16 16 11.9 5.2 14.7 203 2 1.0 37 0 0 0.250 -0.750 16 13 3.7 5.7 14.7 204 3 3.0 31 0 0 -0.125 -0.375 22 22 5.7 10.3 20.0 205 4 4.0 26 0 0 -1.000 -1.000 6 8 16.0 18.2 12.8 206 1 1.0 40 0 0 -4.125 -4.125 17 16 27.5 9.6 14.9 207 1 1.0 40 0 0 -1.500 -1.625 13 14 22.5 6.3 11.7 208 6 6.0 13 1 1 -0.625 0.000 15 13 8.8 18.7 8.1 209 3 4.0 27 0 0 -0.125 0.125 13 12 8.8 9.2 15.0 210 3 2.0 32 0 0 0.125 0.500 20 19 13.3 15.6 11.4 211 1 1.0 40 0 0 -15.125 -15.000 21 22 25.0 12.2 14.8 212 3 4.0 27 0 0 -2.125 -2.000 20 20 7.9 18.6 13.8 213 4 2.0 28 0 0 -2.625 -1.250 15 14 15.8 17.3 10.3 214 4 4.0 26 0 0 -0.250 0.375 18 18 12.8 11.3 9.6 215 2 2.0 35 0 0 -1.000 -0.625 26 24 11.8 11.3 9.3 216 2 3.0 32 0 0 -0.250 -0.375 16 16 11.8 12.3 6.0 217 1 1.0 40 0 0 -1.000 -0.750 16 15 23.7 9.1 9.2 Obs sbp dbp bmi pr sex prot dose resid rerl lerl rl2 id 187 118 73 3.578 37 1 0 1 1 1 1 1 187 188 108 70 2.910 49 2 0 2 1 1 1 1 188 189 144 84 3.270 39 2 0 1 1 1 2 2 189 190 128 94 3.987 48 1 0 2 1 2 2 3 190 191 134 75 3.348 51 2 0 1 1 1 1 1 191 192 130 50 2.924 44 1 0 1 1 1 1 1 192 193 107 72 3.280 40 2 0 1 1 2 2 3 193 194 147 88 3.952 44 2 0 1 1 3 3 6 194 195 132 77 3.255 37 1 1 2 1 3 2 5 195 196 122 80 2.636 36 1 0 1 1 1 1 1 196 197 134 75 3.585 33 1 0 1 1 1 1 1 197 198 120 87 2.939 41 1 0 1 1 2 2 3 198 199 107 70 2.297 43 1 0 1 1 1 1 1 199 200 94 70 2.192 47 1 0 1 1 1 1 1 200 201 119 77 3.229 39 2 0 1 2 1 1 1 201 202 98 60 2.827 43 1 0 1 1 1 1 1 202 203 113 82 2.746 41 2 0 1 2 2 1 2 203 204 116 87 2.649 38 1 0 1 1 2 2 3 204 205 144 89 3.269 52 2 1 2 1 3 3 6 205 206 104 69 3.384 32 1 0 1 1 1 1 1 206 207 88 61 2.863 25 2 0 1 1 1 1 1 207 208 149 96 3.281 42 1 1 2 1 4 4 10 208 209 129 82 3.768 43 2 0 2 1 2 3 5 209 210 132 84 3.784 32 1 0 2 1 2 2 3 210 211 122 79 4.308 46 1 0 2 2 1 1 1 211 212 98 73 3.252 45 2 0 2 1 2 3 5 212 213 117 84 3.541 33 1 0 2 1 3 2 5 213 214 101 65 2.654 30 2 0 2 1 3 3 6 214 215 119 89 3.279 42 2 0 2 1 2 2 3 215 216 112 70 3.377 39 2 0 1 1 2 2 3 216 217 108 76 2.834 41 2 0 2 1 1 1 1 217 Obs rrl lrl rl rme lme rre lre riop liop age diab gh 218 1 1.0 40 0 0 -8.750 -3.875 12 13 22.8 11.3 9.3 219 4 4.0 26 0 0 0.000 0.125 13 14 9.9 13.6 11.9 220 2 1.0 37 0 0 -4.250 -4.500 18 15 23.2 4.8 14.4 221 2 2.0 35 0 0 -1.250 -1.625 16 19 14.7 6.7 14.2 222 2 1.5 36 0 0 0.000 0.125 18 18 6.9 11.4 14.2 223 1 1.5 39 0 0 0.000 0.000 16 13 3.0 15.8 12.2 224 1 1.0 40 0 0 -2.125 -1.875 15 15 9.9 3.9 17.4 225 1 1.0 40 0 0 0.375 -0.125 15 15 8.8 4.8 9.6 226 1 1.0 40 0 0 -0.250 -0.125 12 14 6.9 7.8 17.3 227 1 1.0 40 0 0 0.125 0.000 9 10 2.9 9.5 11.0 228 2 2.0 35 0 0 -4.375 -4.750 21 21 14.3 6.5 13.3 229 2 1.5 36 0 0 0.125 0.125 18 18 5.1 8.0 14.7 230 2 2.0 35 0 0 0.250 -0.375 14 16 10.5 7.3 10.9 231 2 2.0 35 0 0 -0.125 -0.250 13 15 10.6 9.2 13.6 232 3 3.0 31 0 0 -2.375 -2.375 17 15 11.6 11.0 12.1 233 1 1.0 40 0 0 0.250 -0.500 20 19 13.6 5.0 12.9 234 4 4.0 26 0 0 0.000 0.000 17 17 9.6 11.9 14.3 235 3 2.0 32 0 0 -0.125 0.000 17 18 2.8 14.8 15.2 236 1 1.0 40 0 0 0.500 0.000 18 18 3.5 17.3 10.0 237 1 1.0 40 0 0 0.000 -0.250 14 14 1.2 10.1 9.4 238 1 1.0 40 0 0 0.000 0.000 15 15 17.3 4.0 22.6 239 1 1.0 40 0 0 0.000 0.000 12 12 5.6 3.0 12.6 240 1 1.0 40 0 0 -2.125 -2.250 19 20 14.1 5.0 17.3 241 1 1.0 40 0 0 0.000 0.000 11 9 4.7 4.2 9.1 242 1 1.0 40 0 0 -2.875 -3.375 14 14 14.8 2.2 10.8 243 2 2.0 35 0 0 -2.250 -3.875 18 18 11.2 6.8 10.3 244 2 2.0 35 0 0 -2.500 -2.625 18 15 11.5 10.2 9.2 245 3 2.0 32 0 0 0.000 0.125 14 16 11.8 4.1 14.0 246 1 1.0 40 0 0 0.000 0.000 20 19 10.1 6.1 10.0 247 1 1.0 40 0 0 0.000 0.000 13 12 5.3 6.8 11.0 248 2 1.0 37 0 0 0.000 0.250 17 17 4.1 8.9 13.0 Obs sbp dbp bmi pr sex prot dose resid rerl lerl rl2 id 218 160 96 3.826 47 2 0 1 1 1 1 1 218 219 119 78 3.192 45 2 0 2 1 3 3 6 219 220 113 76 3.796 37 1 0 2 2 2 1 2 220 221 117 60 2.792 47 1 0 1 1 2 2 3 221 222 127 86 3.270 38 2 1 1 1 2 2 3 222 223 114 73 3.363 34 1 0 2 2 1 2 2 223 224 119 76 2.185 51 2 0 2 2 1 1 1 224 225 115 69 3.142 34 1 0 2 2 1 1 1 225 226 111 69 2.578 41 2 0 2 2 1 1 1 226 227 112 71 2.605 39 1 0 2 1 1 1 1 227 228 101 72 3.039 34 2 0 2 2 2 2 3 228 229 104 82 2.153 58 2 0 2 1 2 2 3 229 230 121 50 3.318 39 1 0 2 1 2 2 3 230 231 99 56 3.797 34 1 0 1 1 2 2 3 231 232 115 75 3.393 41 2 0 2 2 2 2 3 232 233 111 73 3.855 39 2 0 2 1 1 1 1 233 234 143 86 4.230 33 1 1 2 2 3 3 6 234 235 123 89 3.263 37 1 0 2 2 2 2 3 235 236 104 76 3.256 47 2 1 1 1 1 1 1 236 237 118 89 2.646 42 2 0 2 1 1 1 1 237 238 143 90 2.863 50 1 0 1 1 1 1 1 238 239 112 85 2.348 52 1 0 2 1 1 1 1 239 240 104 78 2.591 36 2 0 2 1 1 1 1 240 241 96 70 2.874 50 1 0 1 1 1 1 1 241 242 118 78 3.068 48 1 0 2 2 1 1 1 242 243 117 84 3.308 39 1 0 2 2 2 2 3 243 244 123 64 3.504 41 1 0 1 1 2 2 3 244 245 111 74 3.121 34 2 0 2 1 2 2 3 245 246 132 62 3.618 32 1 0 2 1 1 1 1 246 247 87 59 2.555 32 1 0 2 1 1 1 1 247 248 100 81 3.413 45 2 0 2 2 2 1 2 248 Obs rrl lrl rl rme lme rre lre riop liop age diab gh 249 2.0 2 35 0 0 0.000 0.000 12 13 10.1 9.1 18.6 250 1.0 1 40 0 0 0.000 0.000 19 17 6.6 2.9 10.7 251 1.0 1 40 0 0 0.000 0.125 15 17 2.8 7.8 11.8 252 4.0 4 26 0 0 0.000 -0.625 10 9 8.1 16.9 12.9 253 2.0 3 32 0 0 0.000 0.000 15 15 4.3 17.9 14.0 254 1.0 1 40 0 0 0.000 0.000 12 13 10.0 2.2 11.2 255 1.0 1 40 0 0 -0.625 -0.375 19 17 6.3 8.3 11.3 256 1.0 1 40 0 0 -1.125 -1.375 18 20 15.5 3.1 13.0 257 1.0 1 40 0 0 -3.125 -2.750 13 15 5.8 7.8 13.8 258 6.0 7 6 0 0 -5.375 -4.625 21 21 17.1 16.6 12.4 259 1.0 2 37 0 0 -3.000 -4.500 18 16 27.6 3.2 9.5 260 1.0 2 37 0 0 0.000 0.250 18 18 25.6 5.6 10.5 261 4.0 3 27 0 0 0.750 -0.125 20 20 23.2 12.5 12.2 262 2.0 2 35 0 0 -0.250 -0.875 19 18 12.2 10.4 14.1 263 1.0 1 40 0 0 0.000 0.000 20 17 25.3 9.5 13.5 264 3.0 3 31 0 0 -4.250 -4.125 16 16 16.7 12.3 9.6 265 6.0 4 15 1 0 -0.375 -0.375 14 13 14.5 19.4 11.6 266 1.0 1 40 0 0 0.125 -0.375 16 17 13.8 2.1 11.5 267 1.0 1 40 0 0 -0.250 0.250 20 21 12.4 1.8 13.7 268 1.0 1 40 0 0 -6.375 -6.125 24 25 22.8 9.4 11.0 269 1.0 1 40 0 0 -5.250 -5.125 18 19 14.4 5.9 11.7 270 1.0 1 40 0 0 -2.500 -3.500 12 15 26.3 6.3 12.4 271 1.0 2 37 0 0 0.250 0.000 17 16 27.2 3.3 15.8 272 4.0 3 27 0 0 0.000 0.000 19 19 17.9 21.6 10.5 273 1.0 2 37 0 0 0.000 0.000 20 18 12.0 9.2 14.7 274 2.0 2 35 0 0 -3.000 -3.000 21 20 8.7 6.6 13.4 275 4.0 4 26 1 0 0.500 1.250 23 22 22.4 33.7 16.7 276 1.0 2 37 0 0 -0.250 0.000 14 14 10.2 7.6 15.2 277 2.0 2 35 0 0 -3.250 -4.500 10 10 16.2 7.8 14.3 278 2.0 2 35 0 0 -0.625 -1.375 17 15 14.3 12.6 9.4 279 6.0 6 13 0 0 -0.875 -0.500 14 15 14.7 34.6 11.7 Obs sbp dbp bmi pr sex prot dose resid rerl lerl rl2 id 249 106 73 3.195 41 2 0 2 2 2 2 3 249 250 94 66 2.402 34 2 0 2 1 1 1 1 250 251 103 80 2.405 36 2 1 2 1 1 1 1 251 252 110 71 3.126 33 1 0 1 2 3 3 6 252 253 103 75 3.106 40 1 0 1 1 2 2 3 253 254 117 66 2.648 34 1 0 2 2 1 1 1 254 255 120 76 3.396 36 1 0 2 1 1 1 1 255 256 120 93 4.122 54 2 0 2 1 1 1 1 256 257 105 61 2.610 43 1 0 1 1 1 1 1 257 258 156 110 4.772 50 1 1 2 1 4 4 10 258 259 97 62 3.528 36 2 0 2 1 1 2 2 259 260 132 74 4.526 43 2 0 1 1 1 2 2 260 261 119 84 3.538 39 1 0 2 1 3 2 5 261 262 121 83 3.461 41 2 0 2 1 2 2 3 262 263 141 85 4.362 36 2 0 1 2 1 1 1 263 264 117 81 3.089 40 1 1 1 1 2 2 3 264 265 117 86 4.170 42 2 0 1 1 4 3 9 265 266 112 92 2.467 44 1 1 1 1 1 1 1 266 267 125 75 3.288 38 1 0 2 1 1 1 1 267 268 152 92 3.729 59 1 0 1 1 1 1 1 268 269 126 69 3.675 32 1 0 1 1 1 1 1 269 270 123 85 2.972 49 1 0 1 1 1 1 1 270 271 109 84 3.159 46 1 0 1 1 1 2 2 271 272 114 70 2.653 44 2 0 1 1 3 2 5 272 273 136 48 3.076 41 1 0 1 1 1 2 2 273 274 117 62 2.401 40 1 0 2 1 2 2 3 274 275 146 75 2.845 47 2 0 1 1 3 3 6 275 276 124 78 3.041 34 1 0 1 1 1 2 2 276 277 110 64 2.866 49 2 1 2 1 2 2 3 277 278 127 95 3.370 43 2 0 1 1 2 2 3 278 279 175 86 3.335 40 1 1 2 1 4 4 10 279 Obs rrl lrl rl rme lme rre lre riop liop age diab gh 280 1.0 1 40 0 0 -1.500 -1.750 18 16 9.4 5.6 13.8 281 1.0 1 40 0 0 -0.125 0.125 12 13 3.3 7.3 14.7 282 2.0 1 37 0 0 0.000 0.000 14 13 15.7 3.7 8.6 283 3.0 3 31 0 0 -3.500 -4.250 17 20 11.9 21.7 14.1 284 2.0 1 37 0 0 0.000 0.000 17 19 7.4 8.5 8.0 285 2.0 2 35 0 0 0.000 0.000 19 21 12.8 6.9 18.4 286 4.0 2 28 0 0 12.875 0.250 22 22 28.5 30.7 14.6 287 3.0 2 32 0 0 0.000 0.000 12 15 24.7 9.7 14.1 288 3.0 3 31 1 1 -1.000 -1.125 18 18 26.0 24.7 13.8 289 2.0 2 35 0 0 1.125 0.500 14 15 19.6 29.7 10.1 290 2.0 4 28 0 0 -3.750 -2.500 14 14 13.4 36.1 7.9 291 1.0 1 40 0 0 -4.250 -3.375 18 16 27.4 5.0 11.5 292 1.0 1 40 0 0 0.375 -1.500 17 15 29.8 4.6 13.4 293 3.0 2 32 0 0 -0.125 -0.375 11 9 9.5 8.2 17.7 294 1.0 1 40 0 0 0.000 0.000 18 20 20.0 1.2 10.4 295 5.0 4 21 0 0 -3.250 -3.625 17 17 18.2 23.7 16.5 296 4.0 4 26 0 0 -3.875 -5.375 14 14 19.2 21.8 16.0 297 1.0 1 40 0 0 0.000 -0.250 15 16 3.7 5.1 9.7 298 3.0 3 31 0 0 0.125 -0.375 25 23 22.7 8.6 11.7 299 1.0 1 40 0 0 0.000 0.000 10 10 4.8 1.7 11.7 300 2.0 2 35 0 0 0.000 0.000 17 16 26.8 10.3 12.0 301 2.0 2 35 0 0 -1.375 -1.500 21 25 14.5 1.6 15.0 302 4.0 4 26 1 0 -2.000 -1.875 15 16 15.9 10.9 15.9 303 1.0 1 40 0 0 -0.250 -0.750 12 13 14.9 11.3 13.7 304 1.0 1 40 0 0 -1.000 -1.375 17 15 29.4 7.6 11.6 305 1.0 1 40 0 0 -1.250 -1.250 15 15 8.4 7.4 7.8 306 1.0 1 40 0 0 0.125 0.125 19 20 9.7 3.2 14.5 307 1.0 1 40 0 0 -0.625 0.375 15 16 11.5 1.3 12.6 308 1.5 3 33 0 0 -2.750 -3.500 15 15 7.0 9.8 16.7 309 1.0 1 40 0 0 -0.250 -0.250 13 13 14.4 1.7 8.3 310 2.0 1 37 0 0 -0.125 -0.125 22 22 5.6 10.3 13.4 Obs sbp dbp bmi pr sex prot dose resid rerl lerl rl2 id 280 110 80 3.070 59 2 0 2 1 1 1 1 280 281 107 79 2.276 39 2 0 1 1 1 1 1 281 282 129 74 3.629 38 1 0 2 1 2 1 2 282 283 99 77 2.657 51 2 1 2 1 2 2 3 283 284 104 66 2.646 44 1 0 2 1 2 1 2 284 285 106 73 3.418 32 1 0 1 1 2 2 3 285 286 214 70 4.206 54 2 0 2 1 3 2 5 286 287 115 90 3.832 45 1 0 1 1 2 2 3 287 288 136 89 3.382 38 1 0 2 1 2 2 3 288 289 111 75 2.502 40 2 0 1 1 2 2 3 289 290 105 69 3.601 30 1 0 2 1 2 3 5 290 291 116 83 3.179 38 2 0 2 1 1 1 1 291 292 126 60 6.341 38 1 0 1 1 1 1 1 292 293 126 87 3.809 40 2 0 2 1 2 2 3 293 294 110 78 3.067 43 2 0 1 1 1 1 1 294 295 132 82 5.553 41 2 0 1 1 3 3 6 295 296 131 70 2.928 50 1 0 1 1 3 3 6 296 297 94 70 2.563 40 2 0 2 1 1 1 1 297 298 125 84 3.871 40 1 0 2 1 2 2 3 298 299 96 70 2.362 47 1 0 2 2 1 1 1 299 300 119 83 3.502 37 1 0 1 1 2 2 3 300 301 109 89 3.072 47 2 0 2 1 2 2 3 301 302 115 59 2.976 45 2 0 1 2 3 3 6 302 303 111 64 2.991 40 2 0 1 1 1 1 1 303 304 125 84 4.511 43 2 1 2 1 1 1 1 304 305 114 66 3.037 30 1 0 2 1 1 1 1 305 306 101 74 2.711 39 2 0 2 1 1 1 1 306 307 109 79 2.223 43 2 0 2 1 1 1 1 307 308 121 84 2.890 52 2 0 2 1 2 2 3 308 309 115 70 3.018 41 2 0 2 1 1 1 1 309 310 114 81 2.917 38 1 1 2 2 2 1 2 310 Obs rrl lrl rl rme lme rre lre riop liop age diab gh 311 2 2 35 0 0 -2.875 -2.750 14 13 16.2 2.4 13.7 312 1 1 40 0 0 0.000 0.000 16 17 10.2 5.9 9.9 313 2 3 32 0 0 -1.250 -1.375 30 23 9.9 13.8 13.7 314 1 1 40 0 0 -1.500 3.000 19 18 26.9 3.0 12.7 315 3 3 31 0 0 -1.250 -1.500 19 17 23.1 9.4 18.1 316 1 2 37 0 0 0.375 0.625 16 15 19.3 23.7 12.2 317 6 6 13 0 0 0.625 0.000 10 14 16.9 23.7 12.7 318 2 2 35 0 0 -3.500 -3.000 16 18 17.3 36.8 11.2 319 2 2 35 0 0 1.250 1.125 13 10 12.7 25.8 13.6 320 3 4 27 0 0 -0.125 0.000 17 18 12.1 16.7 13.9 321 4 4 26 0 0 0.000 0.000 28 24 13.3 18.2 10.7 322 1 1 40 0 0 0.000 0.000 12 12 20.8 1.2 8.0 323 1 1 40 0 0 0.000 0.125 16 13 12.8 3.0 13.3 324 1 1 40 0 0 0.875 0.125 20 23 6.3 2.5 8.9 325 4 4 26 0 0 -0.875 -0.875 16 16 22.7 14.7 11.7 326 2 3 32 0 0 -0.375 -0.250 20 18 24.1 16.7 12.3 327 2 1 37 0 0 0.375 0.500 17 20 5.4 23.1 11.0 328 1 1 40 0 0 -1.000 -0.875 17 17 13.2 3.0 16.0 329 2 2 35 0 0 0.000 -0.250 23 18 9.2 10.2 14.9 330 2 2 35 0 0 -0.250 -0.625 20 20 13.2 11.0 14.4 331 2 2 35 0 0 -1.000 -1.000 22 22 11.9 17.6 10.7 332 1 2 37 0 0 -0.375 -0.875 17 17 12.8 3.9 12.7 333 1 1 40 0 0 0.000 0.000 17 17 10.6 1.7 15.5 334 5 5 20 0 0 -1.000 -1.250 12 11 18.8 20.7 16.5 335 1 1 40 0 0 0.000 0.000 16 13 8.9 1.5 11.8 336 2 1 37 0 0 2.000 1.750 20 19 12.5 47.8 13.8 337 1 1 40 0 0 -0.375 -0.750 18 19 6.1 3.0 10.5 338 1 1 40 0 0 0.000 -1.000 20 22 5.1 5.0 10.9 339 3 4 27 0 0 -2.250 -2.750 16 17 26.2 6.9 11.3 340 4 4 26 0 0 1.250 2.000 20 17 29.4 23.8 11.1 341 2 3 32 0 0 -2.250 -3.875 14 14 25.6 8.0 8.0 Obs sbp dbp bmi pr sex prot dose resid rerl lerl rl2 id 311 118 73 3.346 44 1 0 2 1 2 2 3 311 312 125 63 3.184 36 1 0 2 1 1 1 1 312 313 141 100 3.340 44 2 1 1 1 2 2 3 313 314 107 61 3.475 35 1 0 2 2 1 1 1 314 315 112 87 2.734 45 2 0 1 1 2 2 3 315 316 126 71 3.628 41 2 0 1 2 1 2 2 316 317 136 85 3.118 39 1 0 1 2 4 4 10 317 318 136 71 4.018 42 2 0 2 2 2 2 3 318 319 134 78 4.275 38 2 0 2 2 2 2 3 319 320 122 99 2.577 50 2 0 1 2 2 3 5 320 321 144 81 3.620 49 1 0 1 2 3 3 6 321 322 115 76 3.265 33 1 0 1 2 1 1 1 322 323 114 59 2.646 30 1 0 1 2 1 1 1 323 324 101 75 2.080 46 2 0 1 2 1 1 1 324 325 126 81 3.600 40 1 0 1 2 3 3 6 325 326 155 76 2.721 53 2 0 1 2 2 2 3 326 327 116 84 3.439 49 2 0 1 2 2 1 2 327 328 132 63 2.906 56 1 0 2 2 1 1 1 328 329 123 81 3.662 39 1 0 2 2 2 2 3 329 330 127 80 3.274 43 1 0 1 2 2 2 3 330 331 124 82 3.260 36 2 0 1 1 2 2 3 331 332 119 72 2.862 41 1 0 1 2 1 2 2 332 333 112 83 3.048 65 2 0 1 2 1 1 1 333 334 112 81 3.311 41 1 0 2 2 3 3 6 334 335 115 75 2.208 37 1 0 2 2 1 1 1 335 336 163 78 3.809 32 2 0 1 2 2 1 2 336 337 117 71 2.704 41 2 0 2 2 1 1 1 337 338 106 70 2.222 37 1 0 1 2 1 1 1 338 339 129 85 3.815 37 1 0 1 2 2 3 5 339 340 130 73 4.129 48 1 0 2 2 3 3 6 340 341 130 71 3.951 36 2 0 2 2 2 2 3 341 Obs rrl lrl rl rme lme rre lre riop liop age diab gh 342 1 1 40 0 0 -2.125 -2.500 19 19 13.9 4.9 14.5 343 3 3 31 0 0 -2.375 -0.625 18 14 27.9 18.4 13.8 344 1 1 40 0 0 0.000 0.000 16 19 8.8 5.4 9.6 345 3 3 31 0 0 -0.500 -0.500 18 20 11.2 10.0 15.0 346 2 2 35 0 0 -1.250 -1.250 18 14 7.1 11.8 11.1 347 2 3 32 0 0 -0.500 -0.625 17 16 3.8 14.0 12.5 348 1 1 40 0 0 2.750 -0.500 9 14 11.2 6.2 9.9 349 1 2 37 0 0 -1.625 -1.250 20 18 5.8 8.3 14.7 350 4 4 26 0 0 0.125 -0.250 19 18 10.6 13.7 13.6 351 4 4 26 0 0 -0.750 -0.250 11 12 26.8 10.4 10.2 352 2 2 35 0 0 0.125 0.125 19 17 7.1 11.8 12.7 353 2 1 37 0 0 0.375 0.000 18 18 15.7 4.3 15.1 354 1 1 40 0 0 -0.125 -0.750 13 17 16.6 10.8 10.3 355 3 2 32 0 0 0.000 -0.250 16 15 8.2 11.7 9.6 356 1 1 40 0 0 -1.500 -1.250 15 16 24.5 4.2 11.3 357 4 6 15 0 0 -0.625 -0.500 20 20 14.6 28.7 12.6 358 6 5 14 0 0 0.000 -0.250 17 17 7.9 21.7 16.3 359 3 1 34 0 0 -0.125 -0.125 16 16 12.0 7.7 14.8 360 3 2 32 0 0 -0.500 -0.250 14 12 14.4 11.7 12.6 361 1 2 37 0 0 2.750 2.125 13 13 29.8 29.7 11.9 362 3 3 31 0 0 -3.875 -4.375 14 18 6.9 17.7 12.7 363 1 1 40 0 0 0.125 0.000 14 16 24.4 3.7 11.8 364 1 2 37 0 0 0.250 0.250 11 14 27.0 6.7 13.2 365 2 2 35 0 0 1.375 1.250 14 16 6.8 21.4 12.4 366 3 2 32 0 0 -0.250 -0.750 12 13 25.1 6.2 12.1 367 4 6 15 0 0 0.000 0.000 19 16 12.7 18.7 15.3 368 1 2 37 0 0 -2.125 -0.875 21 21 19.1 5.4 13.2 369 3 4 27 0 0 -7.500 -8.625 21 18 21.0 30.7 10.8 370 2 2 35 0 0 -0.500 -0.500 21 19 18.3 10.7 12.9 371 3 2 32 0 0 -1.125 -1.375 16 16 22.0 11.9 13.1 372 6 6 13 0 0 0.125 -0.125 13 13 28.1 24.9 12.8 Obs sbp dbp bmi pr sex prot dose resid rerl lerl rl2 id 342 123 49 3.148 35 2 0 1 2 1 1 1 342 343 133 74 4.274 40 2 0 1 2 2 2 3 343 344 112 73 2.809 31 1 0 1 2 1 1 1 344 345 119 87 3.472 34 1 1 1 2 2 2 3 345 346 123 58 3.731 33 1 0 1 2 2 2 3 346 347 114 58 3.424 34 1 0 2 2 2 2 3 347 348 142 66 3.055 38 1 0 2 2 1 1 1 348 349 113 77 3.013 43 1 0 2 2 1 2 2 349 350 112 71 2.975 43 1 0 1 2 3 3 6 350 351 109 58 3.580 41 1 0 2 2 3 3 6 351 352 123 77 3.076 52 2 0 2 2 2 2 3 352 353 108 72 2.947 38 1 0 1 2 2 1 2 353 354 97 56 2.938 29 2 0 1 2 1 1 1 354 355 150 90 5.007 48 2 1 1 2 2 2 3 355 356 118 77 2.959 35 1 0 1 2 1 1 1 356 357 175 99 3.760 42 1 1 1 2 3 4 9 357 358 168 106 4.058 42 1 1 1 2 4 3 9 358 359 118 79 3.990 66 2 0 1 2 2 1 2 359 360 121 76 3.283 47 2 0 1 2 2 2 3 360 361 139 77 3.724 31 1 0 1 2 1 2 2 361 362 118 70 3.029 48 2 0 1 2 2 2 3 362 363 119 75 3.755 35 1 0 1 2 1 1 1 363 364 114 74 3.892 37 1 0 1 2 1 2 2 364 365 125 78 4.742 35 2 0 1 2 2 2 3 365 366 105 69 3.962 35 1 0 1 2 2 2 3 366 367 171 117 3.300 48 1 1 2 2 3 4 9 367 368 94 71 3.139 42 1 0 1 2 1 2 2 368 369 145 81 4.466 25 1 0 1 2 2 3 5 369 370 131 86 3.477 45 1 0 2 2 2 2 3 370 371 147 91 3.736 34 1 0 1 2 2 2 3 371 372 156 83 3.391 48 1 0 1 2 4 4 10 372 Obs rrl lrl rl rme lme rre lre riop liop age diab gh 373 6 6 13 0 0 -0.750 0.000 18 18 10.1 27.7 11.6 374 2 2 35 0 0 -0.750 -0.750 24 22 7.8 35.8 11.3 375 3 3 31 0 0 -0.125 0.125 17 16 26.3 29.7 12.3 376 2 2 35 0 0 0.000 0.000 16 14 11.3 5.0 15.0 377 4 2 28 0 0 0.250 -0.375 16 16 24.0 5.2 10.6 378 1 1 40 0 0 0.000 0.000 16 16 6.4 4.3 11.9 379 2 2 35 0 0 -4.250 -3.250 14 16 12.2 8.6 12.7 380 1 1 40 0 0 0.000 -0.250 17 16 8.3 1.3 12.0 381 6 6 13 1 1 0.250 -2.750 19 19 28.1 21.7 12.6 382 3 3 31 0 0 -0.875 -0.500 21 18 4.9 13.1 12.6 383 6 4 15 1 0 -0.625 -0.625 12 16 29.6 12.5 16.6 384 3 4 27 0 0 -4.375 -3.500 22 22 13.7 10.8 14.3 385 4 4 26 0 0 -1.000 -1.625 14 14 25.3 17.8 11.0 386 2 1 37 0 0 -0.250 0.000 16 14 25.2 5.8 11.0 387 1 1 40 0 0 -0.625 -1.000 16 16 17.5 2.6 10.7 388 1 1 40 0 0 -3.750 -4.000 19 19 7.7 8.0 11.6 389 1 1 40 0 0 -3.250 -3.375 12 13 15.9 6.4 10.8 390 1 1 40 0 0 -2.125 -2.750 15 14 20.7 1.9 8.1 391 3 2 32 0 0 -0.750 -1.125 14 14 10.6 41.8 11.1 392 2 2 35 0 0 -3.625 -2.000 18 20 18.8 8.9 8.9 393 1 1 40 0 0 -0.625 0.375 13 14 23.4 7.0 10.2 394 4 2 28 0 0 -0.125 0.000 17 17 26.8 7.4 15.3 395 1 1 40 0 0 -3.000 -3.000 14 14 14.6 5.9 16.5 396 6 6 13 1 1 -2.000 -2.250 14 18 10.9 15.0 12.7 397 3 4 27 0 0 0.000 -0.500 21 22 3.6 28.0 14.0 398 1 1 40 0 0 -1.500 -1.375 16 16 15.7 10.0 12.9 399 2 2 35 0 0 0.000 0.000 21 19 12.1 24.0 10.5 400 4 2 28 0 0 -2.875 -3.375 11 14 7.3 16.9 7.9 401 2 3 32 0 0 4.000 5.000 14 14 2.9 19.5 17.5 402 3 3 31 0 0 -1.000 -0.750 20 18 8.3 11.0 15.3 403 1 1 40 0 0 0.250 -0.125 13 14 11.7 3.9 15.2 Obs sbp dbp bmi pr sex prot dose resid rerl lerl rl2 id 373 129 89 3.351 49 2 0 1 2 4 4 10 373 374 153 81 3.058 45 1 1 1 2 2 2 3 374 375 140 70 3.066 38 1 0 2 2 2 2 3 375 376 130 82 3.669 35 1 0 2 2 2 2 3 376 377 120 70 3.571 41 1 0 1 2 3 2 5 377 378 108 82 2.754 50 1 0 1 2 1 1 1 378 379 117 77 3.246 35 2 0 1 2 2 2 3 379 380 110 76 2.677 43 1 0 1 2 1 1 1 380 381 185 89 5.509 42 1 1 1 2 4 4 10 381 382 116 92 3.677 57 2 1 2 2 2 2 3 382 383 116 77 3.926 34 1 0 2 2 4 3 9 383 384 148 99 3.694 47 2 0 1 2 2 3 5 384 385 132 82 3.424 39 1 0 1 2 3 3 6 385 386 118 68 3.065 40 1 0 1 2 2 1 2 386 387 116 78 3.094 36 1 0 2 2 1 1 1 387 388 140 61 3.852 34 1 0 1 2 1 1 1 388 389 147 98 4.132 57 1 1 1 2 1 1 1 389 390 118 65 3.536 33 1 0 2 2 1 1 1 390 391 135 64 3.640 41 1 0 2 2 2 2 3 391 392 119 83 3.725 39 1 0 1 2 2 2 3 392 393 127 85 5.666 51 1 0 1 2 1 1 1 393 394 106 81 2.794 41 2 0 1 2 3 2 5 394 395 115 63 3.304 45 2 0 1 2 1 1 1 395 396 151 96 3.987 47 1 1 1 2 4 4 10 396 397 142 79 3.309 46 2 0 1 2 2 3 5 397 398 111 77 2.896 32 2 0 1 2 1 1 1 398 399 130 84 4.318 49 2 0 1 2 2 2 3 399 400 108 58 3.247 37 2 0 2 2 3 2 5 400 401 108 83 4.021 38 2 0 1 2 2 2 3 401 402 118 79 3.642 49 2 0 2 2 2 2 3 402 403 108 81 3.484 32 1 0 2 2 1 1 1 403 Obs rrl lrl rl rme lme rre lre riop liop age diab gh 404 1 1 40 0 0 -2.750 -2.875 22 22 17.6 2.5 13.6 405 1 1 40 0 0 0.000 0.000 12 16 5.2 5.2 11.9 406 2 3 32 0 0 -2.875 -3.125 16 12 12.2 10.0 11.5 407 1 1 40 0 0 0.625 -0.125 18 18 28.3 4.6 17.7 408 6 6 13 0 0 -5.625 -5.125 17 19 24.7 10.0 11.3 409 3 3 31 0 0 -1.625 -1.500 20 20 13.0 8.1 15.0 410 2 2 35 0 0 0.000 0.125 11 12 18.8 30.0 11.3 411 1 1 40 0 0 0.000 0.000 9 9 28.6 3.9 8.1 412 2 2 35 0 0 0.000 0.000 14 20 3.4 15.0 12.1 413 1 1 40 0 0 -0.250 -2.125 17 18 11.6 2.2 10.4 414 3 2 32 0 0 -0.500 -0.500 23 20 10.6 4.5 12.8 415 2 2 35 0 0 0.125 0.125 21 18 20.1 6.2 13.4 416 4 4 26 0 0 -3.500 -3.875 20 20 11.9 15.8 11.1 417 4 4 26 0 0 0.375 0.000 11 10 17.2 9.7 15.0 418 3 3 31 0 0 -1.000 -0.750 16 19 29.1 24.1 13.3 419 3 4 27 0 0 -1.125 -1.750 14 15 26.8 28.1 12.4 420 3 2 32 0 0 0.250 -1.000 16 14 19.6 13.1 11.8 421 4 3 27 0 0 -9.750 -9.875 15 14 9.9 13.6 13.8 422 1 1 40 0 0 0.000 0.000 10 11 26.1 5.0 11.0 423 3 3 31 0 0 -0.250 -0.125 16 17 27.0 13.1 9.1 424 3 3 31 0 0 0.000 0.125 14 14 1.7 26.1 14.1 425 1 1 40 0 0 -0.125 -0.125 15 15 27.2 9.5 15.6 426 6 6 13 0 0 0.250 0.500 16 18 7.7 27.3 11.9 427 3 3 31 0 0 -0.625 -0.750 10 11 8.7 27.5 11.5 428 3 3 31 0 0 0.000 0.000 10 10 17.2 26.1 10.8 429 4 4 26 0 0 -1.000 -1.000 13 17 23.2 21.8 11.6 430 2 2 35 0 0 -1.125 -1.250 18 17 13.7 19.0 9.5 431 2 1 37 0 0 -0.500 -1.000 15 15 13.5 38.3 12.6 432 2 2 35 0 0 0.750 0.625 12 14 23.8 35.1 11.9 433 2 3 32 0 0 0.000 0.000 10 10 9.8 9.2 9.7 434 2 2 35 0 0 0.000 0.000 12 10 22.4 6.0 15.1 Obs sbp dbp bmi pr sex prot dose resid rerl lerl rl2 id 404 142 98 4.907 45 2 0 1 2 1 1 1 404 405 87 62 2.687 43 1 0 2 2 1 1 1 405 406 92 68 2.748 48 2 0 1 2 2 2 3 406 407 137 69 3.964 48 1 0 1 1 1 1 1 407 408 152 97 3.478 52 2 0 1 2 4 4 10 408 409 125 72 3.529 35 2 0 1 2 2 2 3 409 410 104 69 3.050 41 2 0 1 2 2 2 3 410 411 118 69 3.932 51 1 0 0 2 1 1 1 411 412 112 86 3.135 46 2 1 2 2 2 2 3 412 413 115 72 2.591 49 2 0 2 2 1 1 1 413 414 125 61 3.414 45 2 0 2 2 2 2 3 414 415 117 85 3.349 45 1 0 1 2 2 2 3 415 416 141 88 3.194 48 2 1 1 2 3 3 6 416 417 122 80 3.054 43 1 0 2 2 3 3 6 417 418 148 85 2.946 44 1 0 1 1 2 2 3 418 419 116 83 3.540 34 1 0 2 2 2 3 5 419 420 112 76 3.265 34 2 0 2 2 2 2 3 420 421 126 76 3.404 42 1 0 1 2 3 2 5 421 422 119 74 4.682 37 2 0 1 2 1 1 1 422 423 101 72 4.422 48 2 0 1 2 2 2 3 423 424 111 85 4.617 34 2 0 1 2 2 2 3 424 425 146 74 3.823 37 1 0 1 2 1 1 1 425 426 143 73 2.947 53 2 0 1 1 4 4 10 426 427 105 70 3.585 34 2 0 1 2 2 2 3 427 428 144 94 3.852 41 1 0 1 2 2 2 3 428 429 153 96 3.895 43 1 0 2 1 3 3 6 429 430 131 71 3.442 41 1 0 2 2 2 2 3 430 431 142 72 4.150 46 2 0 1 1 2 1 2 431 432 131 76 3.000 41 1 0 2 2 2 2 3 432 433 113 69 3.909 45 2 0 2 2 2 2 3 433 434 137 84 3.394 51 1 0 2 1 2 2 3 434 Obs rrl lrl rl rme lme rre lre riop liop age diab gh 435 1 2 37 0 0 -2.875 -1.875 14 13 7.6 9.1 17.3 436 2 2 35 0 0 -2.375 -2.875 16 19 19.1 10.2 12.4 437 4 4 26 0 0 -2.750 -3.250 20 24 9.4 8.8 21.3 438 1 1 40 0 0 -1.125 0.250 12 13 14.8 1.9 8.0 439 1 1 40 0 0 0.000 0.000 17 19 13.2 4.1 18.4 440 4 2 28 0 0 0.000 0.000 20 22 6.8 11.8 13.4 441 3 3 31 0 0 0.000 0.000 16 16 7.0 14.1 13.6 442 3 2 32 0 0 1.000 0.750 31 31 21.4 27.2 12.0 443 1 1 40 0 0 -1.125 -0.375 22 21 11.8 27.7 9.1 444 1 1 40 0 0 -0.250 0.125 14 13 10.6 9.1 13.0 445 3 3 31 0 0 -1.250 -0.875 12 13 26.6 13.1 11.0 446 1 1 40 0 0 -1.250 -1.125 14 14 14.0 2.0 11.1 447 1 1 40 0 0 0.125 0.000 13 13 13.0 2.9 15.8 448 4 4 26 0 0 0.250 0.000 15 17 14.9 8.9 13.1 449 2 1 37 0 0 0.375 0.375 16 13 26.0 24.2 10.3 450 3 6 16 0 0 -0.375 -0.250 10 11 24.7 25.4 13.3 451 4 4 26 0 0 0.125 0.000 22 20 19.9 14.0 15.8 452 1 1 40 0 0 -1.750 -2.000 19 18 11.0 7.0 10.0 453 2 2 35 0 0 -1.500 -1.125 11 11 16.6 8.0 9.6 454 7 7 5 1 1 0.000 0.000 12 12 12.1 12.7 12.4 455 2 1 37 0 0 -0.250 0.375 17 20 17.7 29.4 9.7 456 1 2 37 0 0 -6.000 -5.500 20 18 27.0 7.3 10.4 457 2 2 35 0 0 -1.500 -1.000 20 24 27.5 32.9 9.7 458 6 4 15 0 0 -1.250 -2.625 22 21 6.3 13.9 13.8 459 2 2 35 0 0 -0.500 -0.500 19 15 10.7 7.5 18.6 460 1 1 40 0 0 0.000 -0.125 17 17 14.8 2.0 8.2 461 6 4 15 0 0 -5.625 -5.375 21 19 11.7 15.9 10.8 462 3 4 27 0 0 0.000 0.000 13 16 8.0 9.5 8.2 463 4 4 26 0 0 14.375 12.125 14 14 28.1 28.9 12.6 464 1 1 40 0 0 -0.750 -0.500 24 25 13.2 3.9 15.3 465 6 4 15 0 1 -5.125 -5.125 21 22 8.9 21.8 9.3 Obs sbp dbp bmi pr sex prot dose resid rerl lerl rl2 id 435 118 58 2.602 40 2 0 1 2 1 2 2 435 436 128 80 3.268 41 1 0 1 2 2 2 3 436 437 109 77 2.814 40 1 0 1 2 3 3 6 437 438 116 65 2.775 36 1 0 2 2 1 1 1 438 439 120 76 3.303 39 2 0 2 2 1 1 1 439 440 116 82 3.912 41 1 0 1 2 3 2 5 440 441 114 74 3.153 47 1 0 1 2 2 2 3 441 442 116 66 3.328 30 2 0 1 2 2 2 3 442 443 118 76 3.620 37 1 0 1 1 1 1 1 443 444 141 89 3.470 44 1 0 2 1 1 1 1 444 445 121 79 4.022 41 2 0 1 2 2 2 3 445 446 113 78 2.895 34 2 0 1 2 1 1 1 446 447 121 71 2.991 39 2 0 2 2 1 1 1 447 448 112 68 3.743 36 2 0 1 2 3 3 6 448 449 144 84 4.341 46 2 0 1 2 2 1 2 449 450 119 56 4.247 30 1 1 2 2 2 4 8 450 451 124 88 3.909 41 2 0 1 2 3 3 6 451 452 113 69 4.323 47 2 0 1 2 1 1 1 452 453 116 68 3.354 33 2 0 2 2 2 2 3 453 454 142 99 3.043 44 1 1 1 2 4 4 10 454 455 120 69 2.804 42 2 0 1 1 2 1 2 455 456 118 88 3.838 37 2 0 2 2 1 2 2 456 457 166 88 3.479 55 1 1 2 2 2 2 3 457 458 117 83 3.591 43 2 1 1 2 4 3 9 458 459 109 80 3.033 53 2 0 1 2 2 2 3 459 460 114 60 3.104 36 1 0 1 2 1 1 1 460 461 113 90 3.728 36 2 1 1 2 4 3 9 461 462 104 70 3.224 41 1 0 2 2 2 3 5 462 463 164 77 4.920 44 2 1 2 2 3 3 6 463 464 137 97 3.394 51 1 0 1 2 1 1 1 464 465 120 72 2.961 40 2 0 2 2 4 3 9 465 Obs rrl lrl rl rme lme rre lre riop liop age diab gh 466 1 1 40 0 0 0.250 0.000 15 15 10.2 1.7 12.2 467 1 1 40 0 0 0.000 -0.375 11 12 9.9 6.8 14.9 468 2 2 35 0 0 0.000 0.125 19 17 12.2 9.3 13.6 469 1 1 40 0 0 -2.125 -2.500 14 11 28.6 12.6 7.3 470 1 1 40 0 0 -9.375 -8.875 16 18 10.4 6.0 12.1 471 1 1 40 0 0 0.250 0.250 12 12 7.7 4.5 14.4 472 1 1 40 0 0 -1.375 -1.750 13 11 18.2 1.9 14.0 473 1 1 40 0 0 -2.125 -2.000 13 15 11.0 8.4 13.9 474 3 3 31 0 0 -7.000 -6.500 19 21 16.0 16.4 11.1 475 3 2 32 0 0 0.000 0.000 14 14 12.6 16.4 8.6 476 1 1 40 0 0 1.000 0.000 10 12 28.3 3.4 11.5 477 2 2 35 0 0 -1.750 -1.750 12 11 9.7 11.4 12.7 478 4 4 26 0 0 -3.125 -3.000 13 13 16.2 23.4 9.7 479 2 2 35 0 0 -0.875 -1.250 11 12 2.2 10.4 14.8 480 1 1 40 0 0 -0.625 -1.750 13 13 3.2 7.8 10.9 481 1 1 40 0 0 -4.500 -5.125 15 15 8.0 6.3 9.7 482 1 1 40 0 0 -1.875 -0.875 15 15 13.7 2.8 12.9 483 1 1 40 0 0 -0.250 0.000 12 14 27.9 1.5 8.9 484 2 3 32 0 0 -4.000 -4.125 14 13 5.5 14.2 12.2 485 3 4 27 0 0 -7.750 -7.875 20 21 19.3 20.4 10.0 486 3 3 31 0 0 -4.000 -4.750 15 15 4.7 17.4 13.6 487 1 1 40 0 0 0.000 0.000 13 12 8.4 2.8 16.2 488 1 1 40 0 0 0.000 0.000 19 18 26.7 7.8 14.8 489 3 2 32 0 0 -0.375 -0.375 18 14 26.8 3.3 11.5 490 2 2 35 0 0 0.000 0.000 16 15 8.5 8.4 10.3 491 2 2 35 0 0 0.000 0.000 12 13 10.2 5.7 11.0 492 4 4 26 0 0 -2.875 -2.000 18 16 26.3 36.4 13.8 493 3 4 27 0 1 -0.500 -0.500 11 12 20.2 25.4 12.4 494 3 1 34 1 0 0.000 0.000 13 12 17.9 19.4 14.7 495 2 2 35 0 0 0.000 0.250 13 16 12.2 17.4 12.8 496 1 1 40 0 0 -2.875 -3.500 15 15 14.7 4.9 13.1 Obs sbp dbp bmi pr sex prot dose resid rerl lerl rl2 id 466 129 72 2.537 55 2 1 1 2 1 1 1 466 467 122 65 3.160 32 1 0 1 2 1 1 1 467 468 116 75 3.491 38 2 1 2 2 2 2 3 468 469 138 100 3.926 37 2 0 0 2 1 1 1 469 470 111 71 3.217 44 2 0 1 2 1 1 1 470 471 114 50 3.520 36 2 0 2 2 1 1 1 471 472 103 72 2.645 45 2 0 1 2 1 1 1 472 473 105 70 3.730 44 2 0 2 2 1 1 1 473 474 142 88 3.299 44 1 0 2 2 2 2 3 474 475 112 74 3.930 37 2 0 2 2 2 2 3 475 476 96 69 3.012 36 2 0 2 2 1 1 1 476 477 130 59 3.461 28 1 1 1 2 2 2 3 477 478 116 65 3.645 35 2 0 1 2 3 3 6 478 479 128 71 2.536 43 1 1 1 2 2 2 3 479 480 137 57 2.457 47 2 0 2 2 1 1 1 480 481 118 51 2.884 30 1 0 2 2 1 1 1 481 482 126 81 3.177 45 1 0 1 2 1 1 1 482 483 111 82 5.381 36 2 0 1 2 1 1 1 483 484 124 65 3.768 37 1 0 2 2 2 2 3 484 485 169 101 3.752 53 2 1 2 2 2 3 5 485 486 114 77 3.011 34 2 0 2 2 2 2 3 486 487 98 69 2.404 44 1 0 1 2 1 1 1 487 488 134 95 7.228 51 2 0 1 2 1 1 1 488 489 124 73 3.486 41 1 0 1 2 2 2 3 489 490 125 85 2.912 39 1 0 2 2 2 2 3 490 491 115 73 2.861 30 1 0 2 2 2 2 3 491 492 162 86 3.376 47 2 0 1 2 3 3 6 492 493 113 61 3.287 35 2 0 2 2 2 3 5 493 494 119 84 3.227 43 2 0 2 2 2 1 2 494 495 127 83 3.468 31 2 0 1 2 2 2 3 495 496 120 71 3.259 25 1 0 2 2 1 1 1 496 Obs rrl lrl rl rme lme rre lre riop liop age diab gh 497 1 1 40 0 0 -2.750 -2.500 13 11 21.2 1.9 11.4 498 1 1 40 0 0 -0.500 -0.250 18 18 10.4 5.1 14.2 499 1 1 40 0 0 -5.500 -5.250 16 15 12.3 6.2 17.2 500 4 4 26 0 0 -0.875 -1.750 15 15 20.1 24.2 12.1 501 4 3 27 0 0 -0.375 0.000 15 17 8.4 16.6 11.1 502 2 3 32 0 0 0.375 0.125 20 19 4.3 49.2 14.1 503 4 4 26 0 1 2.125 1.500 14 14 12.8 13.1 14.6 504 4 4 26 0 0 0.000 -0.500 15 17 22.6 12.8 12.6 505 4 4 26 1 1 0.000 0.625 20 19 2.1 55.2 17.1 506 2 1 37 0 0 -2.250 -2.375 13 13 27.2 5.7 9.2 507 1 1 40 0 0 0.000 0.000 18 18 18.3 5.6 8.8 508 4 3 27 0 0 0.125 0.250 12 13 29.5 10.4 10.8 509 1 1 40 0 0 -5.000 -4.500 13 14 10.9 12.1 13.2 510 4 4 26 0 0 -1.750 -1.375 16 16 6.9 12.1 10.4 511 1 1 40 0 0 -3.125 -3.500 10 10 25.9 9.2 14.4 512 4 4 26 0 0 -2.000 -2.875 11 11 12.7 9.1 13.5 513 6 6 13 0 0 -1.500 -1.125 17 16 12.0 13.2 15.3 514 4 4 26 0 0 -0.375 0.125 16 19 14.0 33.2 13.5 515 6 6 13 0 0 -0.875 -0.375 16 16 5.1 18.1 11.2 516 1 1 40 0 0 -2.000 -2.000 13 15 16.4 4.2 7.6 517 2 2 35 0 0 -1.000 -0.750 13 13 13.7 22.1 9.5 518 6 6 13 1 0 -0.500 -2.000 13 13 28.8 37.2 11.4 519 1 1 40 0 0 0.125 0.000 9 8 27.6 38.2 12.6 520 3 3 31 0 0 0.625 -0.250 15 13 3.8 14.2 13.5 521 3 4 27 0 0 -1.250 -1.375 16 17 20.3 11.1 13.5 522 7 7 5 1 1 -1.500 -1.625 23 20 12.6 15.7 17.1 523 3 4 27 0 0 0.250 -0.375 13 13 24.0 11.2 17.1 524 1 1 40 0 0 -5.375 -5.375 17 17 25.7 7.9 11.6 525 1 1 40 0 0 -4.000 -4.625 15 15 29.1 1.6 12.2 526 1 1 40 0 0 -0.625 -0.375 17 14 18.0 1.8 16.2 527 1 3 34 0 0 13.000 2.375 9 15 24.1 40.2 10.2 Obs sbp dbp bmi pr sex prot dose resid rerl lerl rl2 id 497 127 80 3.825 30 1 0 1 2 1 1 1 497 498 121 68 3.588 45 2 0 2 2 1 1 1 498 499 120 52 3.522 37 1 0 2 2 1 1 1 499 500 143 81 3.595 38 1 0 1 1 3 3 6 500 501 112 83 3.534 36 2 0 1 1 3 2 5 501 502 97 59 2.439 30 1 0 1 1 2 2 3 502 503 137 86 3.543 46 1 0 1 1 3 3 6 503 504 131 84 4.001 39 1 0 1 1 3 3 6 504 505 159 76 3.846 35 2 0 1 1 3 3 6 505 506 135 89 3.144 36 1 0 1 1 2 1 2 506 507 116 83 2.701 38 2 0 2 1 1 1 1 507 508 121 84 3.708 36 1 0 2 2 3 2 5 508 509 105 74 3.656 37 2 0 1 1 1 1 1 509 510 84 76 2.855 57 2 1 1 1 3 3 6 510 511 107 83 4.248 44 2 0 1 1 1 1 1 511 512 116 87 2.819 47 1 0 1 2 3 3 6 512 513 139 82 2.995 50 1 0 1 1 4 4 10 513 514 120 75 4.361 47 2 0 1 1 3 3 6 514 515 127 88 2.995 43 2 1 1 1 4 4 10 515 516 107 81 2.882 31 2 0 2 1 1 1 1 516 517 115 76 3.832 37 1 0 1 1 2 2 3 517 518 143 80 3.096 41 1 0 1 1 4 4 10 518 519 120 68 3.693 44 1 0 1 1 1 1 1 519 520 144 85 3.303 55 1 0 1 1 2 2 3 520 521 122 89 3.098 40 1 0 1 1 2 3 5 521 522 109 83 2.950 44 1 0 1 1 4 4 10 522 523 131 81 3.671 36 1 0 2 1 2 3 5 523 524 108 80 3.530 33 1 0 2 1 1 1 1 524 525 108 74 2.821 34 2 0 1 1 1 1 1 525 526 121 79 3.148 45 1 0 1 1 1 1 1 526 527 130 43 3.520 27 2 0 1 1 1 2 2 527 Obs rrl lrl rl rme lme rre lre riop liop age diab gh 528 3 1 34 0 0 -0.250 0.000 11 11 14.5 6.5 15.6 529 4 4 26 0 0 -2.125 -2.375 13 13 17.5 14.2 12.8 530 2 2 35 0 0 0.000 0.125 20 19 6.2 16.4 14.9 531 6 7 6 0 0 -1.500 -0.250 17 17 20.9 11.2 15.5 532 4 4 26 0 0 -1.500 -1.375 19 18 9.5 12.2 13.2 533 1 1 40 0 0 -0.250 0.125 12 15 5.7 1.8 10.2 534 1 1 40 0 0 -4.375 -3.625 21 18 7.0 5.6 10.1 535 1 1 40 0 0 -3.875 -5.000 19 20 15.7 2.1 14.8 536 2 2 35 0 0 -0.125 -0.250 23 24 12.6 6.1 17.9 537 1 1 40 0 0 -2.125 -3.625 15 16 7.4 5.6 15.9 538 1 1 40 0 0 0.000 0.000 11 13 16.0 2.4 13.4 539 4 4 26 1 1 -0.250 -0.500 17 18 4.8 16.1 15.9 540 1 1 40 0 0 -0.250 -0.125 11 10 4.9 3.9 13.3 541 1 1 40 0 0 1.000 0.875 14 14 1.1 5.6 15.7 542 1 1 40 0 0 0.000 0.000 11 11 3.2 1.8 11.2 543 1 1 40 0 0 -0.125 0.000 21 20 12.8 1.5 11.8 544 1 1 40 0 0 0.000 0.000 15 15 6.3 3.1 12.8 545 1 1 40 0 0 0.125 0.125 14 16 6.2 2.6 10.6 546 4 4 26 0 0 -0.250 -1.000 16 14 9.2 13.0 14.6 547 1 1 40 0 0 0.000 0.000 10 13 4.5 8.8 14.9 548 1 1 40 0 0 0.000 0.000 14 13 8.1 10.1 12.1 549 1 1 40 0 0 0.000 0.000 17 13 21.4 1.8 11.4 550 2 2 35 0 0 -1.875 -1.375 18 20 22.0 9.2 14.9 551 3 3 31 0 0 -0.250 -0.500 16 17 15.3 5.1 12.1 552 4 3 27 0 0 2.625 2.375 22 19 28.3 33.2 12.7 553 3 3 31 0 0 -3.375 -3.500 15 17 19.3 15.3 12.3 554 3 4 27 0 0 0.625 0.625 14 12 22.3 25.3 14.4 555 3 3 31 0 0 -2.375 -1.875 15 14 19.4 49.3 12.0 556 3 2 32 0 0 -1.375 -2.250 14 14 16.0 11.3 12.4 557 6 7 6 0 0 -0.125 0.375 19 16 26.1 14.3 10.7 558 4 4 26 0 0 -0.125 -0.125 21 24 16.3 15.0 11.5 Obs sbp dbp bmi pr sex prot dose resid rerl lerl rl2 id 528 123 81 3.402 41 1 0 1 1 2 1 2 528 529 116 78 3.244 52 1 0 1 1 3 3 6 529 530 106 79 3.520 45 2 0 1 2 2 2 3 530 531 141 71 2.983 41 1 0 2 1 4 4 10 531 532 141 86 2.664 42 1 1 1 1 3 3 6 532 533 88 66 2.198 34 2 0 1 1 1 1 1 533 534 112 78 2.803 45 2 0 2 2 1 1 1 534 535 113 75 2.951 42 1 0 1 1 1 1 1 535 536 122 87 3.125 38 1 0 1 1 2 2 3 536 537 110 76 2.393 39 2 0 2 1 1 1 1 537 538 112 65 2.955 31 1 0 1 1 1 1 1 538 539 116 77 3.365 42 1 0 1 1 3 3 6 539 540 97 68 2.570 43 1 0 1 1 1 1 1 540 541 105 60 2.313 51 2 0 1 1 1 1 1 541 542 88 52 2.049 43 1 0 1 1 1 1 1 542 543 116 76 2.661 41 2 0 1 1 1 1 1 543 544 106 69 2.641 30 1 0 1 1 1 1 1 544 545 78 53 2.663 54 1 0 1 2 1 1 1 545 546 111 74 3.062 43 1 0 1 1 3 3 6 546 547 104 67 2.875 34 1 0 1 1 1 1 1 547 548 118 73 3.595 28 1 0 1 1 1 1 1 548 549 95 76 2.849 40 2 0 2 1 1 1 1 549 550 113 85 3.207 53 2 0 1 1 2 2 3 550 551 103 66 4.167 41 2 0 1 1 2 2 3 551 552 139 76 2.474 33 2 0 1 1 3 2 5 552 553 116 79 3.200 41 2 0 1 1 2 2 3 553 554 153 77 2.630 46 2 0 1 1 2 3 5 554 555 166 71 3.075 37 2 0 1 1 2 2 3 555 556 113 80 3.640 39 1 1 1 1 2 2 3 556 557 174 109 3.257 42 1 1 1 1 4 4 10 557 558 119 81 3.415 34 1 0 1 1 3 3 6 558 Obs rrl lrl rl rme lme rre lre riop liop age diab gh 559 1.0 1 40 0 0 -1.000 -1.250 13 13 27.9 19.3 12.0 560 1.0 2 37 0 0 -2.000 -1.750 12 13 22.0 17.3 11.9 561 3.0 4 27 0 1 -0.750 -0.125 18 19 13.4 10.9 15.3 562 1.0 1 40 0 0 -0.750 -0.875 12 13 1.9 9.8 11.0 563 2.0 3 32 0 0 -4.000 -4.000 16 17 9.2 14.4 15.7 564 1.0 1 40 0 0 -2.250 -2.375 11 10 9.7 8.3 13.9 565 3.0 3 31 0 0 0.000 0.125 16 16 13.4 7.1 15.7 566 4.0 2 28 0 0 -0.125 -0.250 17 20 7.3 29.3 9.0 567 2.0 1 37 0 0 0.000 -0.750 9 10 11.9 14.3 8.9 568 1.0 2 37 0 0 0.000 0.000 12 11 18.6 3.0 11.4 569 4.0 4 26 0 0 -1.875 -2.000 15 12 8.3 11.9 12.3 570 4.0 2 28 0 0 1.500 1.250 21 20 26.0 17.6 12.2 571 2.0 3 32 0 0 -0.625 0.000 12 15 27.8 28.3 12.0 572 4.0 6 15 0 0 2.750 2.250 11 16 9.2 23.1 14.0 573 1.0 1 40 0 0 -0.125 -0.500 16 18 6.6 9.1 18.6 574 1.0 2 37 0 0 0.000 0.000 17 16 24.2 5.8 11.8 575 4.0 3 27 0 0 -7.000 -6.250 18 18 21.9 26.4 11.9 576 4.0 3 27 0 0 0.375 0.125 19 17 28.8 27.2 14.1 577 2.0 2 35 0 0 0.125 1.750 13 15 24.1 15.2 11.7 578 4.0 4 26 0 0 -2.500 -2.500 11 10 12.4 8.8 14.4 579 1.0 1 40 0 0 4.875 4.750 18 17 27.3 4.9 15.4 580 4.0 6 15 1 1 -0.250 0.125 17 18 8.6 32.0 10.4 581 4.0 4 26 1 0 -0.750 -0.375 17 13 12.7 13.1 14.2 582 1.0 1 40 0 0 -2.250 -2.500 12 9 27.5 8.2 11.3 583 2.0 2 35 0 0 -1.750 -0.375 17 16 23.9 6.5 17.2 584 3.0 4 27 0 0 -1.250 -1.375 11 10 13.7 22.8 6.4 585 1.0 1 40 0 0 0.125 0.000 19 16 15.1 2.5 11.2 586 1.0 1 40 0 0 -2.000 -2.125 17 19 16.7 5.1 10.6 587 4.0 3 27 0 0 -4.625 -4.250 15 14 24.6 21.8 12.0 588 1.0 1 40 0 0 0.000 -0.500 14 16 8.1 4.5 16.9 589 6.0 4 15 0 0 1.250 1.125 15 14 17.5 33.8 12.9 Obs sbp dbp bmi pr sex prot dose resid rerl lerl rl2 id 559 107 78 4.694 37 2 0 0 1 1 1 1 559 560 104 62 3.252 39 2 0 1 1 1 2 2 560 561 147 91 3.441 53 1 1 1 1 2 3 5 561 562 102 55 2.738 28 1 0 2 1 1 1 1 562 563 129 92 4.490 50 1 0 2 1 2 2 3 563 564 105 64 3.456 39 1 0 2 1 1 1 1 564 565 124 84 3.162 48 1 0 1 1 2 2 3 565 566 173 84 3.725 52 2 0 2 1 3 2 5 566 567 118 78 3.198 44 1 0 2 1 2 1 2 567 568 123 71 4.254 39 2 0 1 1 1 2 2 568 569 106 74 3.558 37 1 0 1 1 3 3 6 569 570 124 76 5.075 36 2 0 1 1 3 2 5 570 571 141 73 3.728 36 1 0 1 1 2 2 3 571 572 135 89 4.130 38 2 0 2 2 3 4 9 572 573 129 102 3.088 49 2 0 1 1 1 1 1 573 574 101 69 2.534 52 2 0 2 1 1 2 2 574 575 146 73 4.085 45 2 0 1 1 3 2 5 575 576 130 86 3.064 58 1 0 1 1 3 2 5 576 577 131 91 3.976 43 2 0 1 1 2 2 3 577 578 126 92 2.928 47 1 0 1 1 3 3 6 578 579 116 86 3.067 58 1 0 2 1 1 1 1 579 580 177 103 4.058 37 1 1 2 1 3 4 9 580 581 136 84 3.661 43 1 1 1 1 3 3 6 581 582 109 72 2.871 59 2 0 2 2 1 1 1 582 583 139 78 4.921 53 1 0 0 2 2 2 3 583 584 125 86 3.354 33 1 0 2 2 2 3 5 584 585 110 73 3.032 43 2 0 2 2 1 1 1 585 586 121 83 2.571 50 2 0 1 2 1 1 1 586 587 138 79 3.819 36 1 0 2 2 3 2 5 587 588 115 79 2.432 59 1 0 2 2 1 1 1 588 589 157 79 3.099 43 2 1 1 2 4 3 9 589 Obs rrl lrl rl rme lme rre lre riop liop age diab gh 590 2.0 2 35 0 0 -0.750 -1.500 21 18 2.5 36.9 11.3 591 1.0 1 40 0 0 -0.625 -0.625 17 18 23.7 1.8 9.0 592 1.0 1 40 0 0 0.500 0.500 13 13 5.9 3.0 13.6 593 2.0 4 28 0 0 -2.625 -2.500 14 15 25.0 7.3 9.7 594 1.0 2 37 0 0 -2.875 -3.125 18 18 11.8 5.0 16.9 595 1.5 1 39 0 0 0.250 -0.125 14 14 12.5 6.8 9.6 596 2.0 2 35 0 0 0.250 0.500 21 18 11.2 6.9 11.1 597 1.0 1 40 0 0 0.000 0.000 13 13 10.0 1.4 11.0 598 3.0 2 32 0 0 -1.125 -1.250 18 14 13.4 7.8 15.6 599 3.0 2 32 0 0 0.000 0.000 16 15 10.7 5.9 18.7 600 1.0 1 40 0 0 -3.250 -3.250 14 13 4.9 6.8 14.6 601 1.0 1 40 0 0 -4.250 -4.000 17 16 26.5 5.8 11.9 602 1.0 1 40 0 0 -4.000 -3.875 15 16 12.3 4.9 15.9 603 4.0 4 26 0 0 -1.625 -1.750 20 20 2.9 32.9 14.3 604 2.0 2 35 0 0 -1.125 -0.875 14 14 9.9 10.1 10.1 605 1.0 1 40 0 0 0.000 0.000 9 9 7.6 3.4 9.3 606 2.0 4 28 0 0 -4.500 -2.125 22 21 3.5 19.3 9.4 607 3.0 3 31 0 0 3.125 2.250 19 19 28.3 27.6 13.6 608 2.0 2 35 0 0 0.750 0.625 19 19 1.3 37.5 11.1 609 1.0 1 40 0 0 -5.000 -4.125 18 18 1.4 11.1 11.3 610 4.0 4 26 1 1 -2.250 -2.750 20 17 16.6 20.7 17.0 611 6.0 4 15 0 0 0.125 0.125 15 18 11.5 9.7 11.2 612 6.0 6 13 0 0 -1.125 -0.875 20 18 3.5 25.6 11.3 613 1.0 2 37 0 0 0.000 -0.375 17 17 12.1 10.7 11.2 614 1.0 2 37 0 0 -0.500 -0.500 20 18 23.6 5.3 11.0 615 6.0 6 13 0 0 -0.875 -0.750 18 15 21.6 5.9 11.0 616 1.5 1 39 0 0 -1.000 -1.000 8 9 21.4 5.8 13.5 617 1.0 1 40 0 0 0.125 0.000 18 17 5.9 10.1 12.6 618 6.0 4 15 0 1 0.000 0.000 11 12 7.6 14.1 14.8 619 1.0 1 40 0 0 0.500 0.500 13 14 10.7 3.3 16.0 620 2.0 3 32 0 0 -6.250 -6.875 9 9 27.8 8.9 10.2 Obs sbp dbp bmi pr sex prot dose resid rerl lerl rl2 id 590 120 65 3.059 47 2 0 1 2 2 2 3 590 591 108 59 3.233 37 1 0 2 2 1 1 1 591 592 87 67 3.036 49 2 0 2 2 1 1 1 592 593 112 72 2.709 42 2 1 2 2 2 3 5 593 594 122 94 3.736 53 2 1 2 2 1 2 2 594 595 125 87 3.944 34 2 0 2 2 2 1 2 595 596 88 73 2.416 49 1 0 1 2 2 2 3 596 597 107 69 2.704 39 1 0 1 1 1 1 1 597 598 107 83 2.969 39 2 1 2 2 2 2 3 598 599 131 85 3.082 58 1 0 1 2 2 2 3 599 600 110 66 2.169 48 2 0 2 2 1 1 1 600 601 128 90 4.663 43 2 0 1 2 1 1 1 601 602 127 71 2.672 32 1 0 2 2 1 1 1 602 603 139 77 3.220 45 1 0 1 2 3 3 6 603 604 114 74 3.737 37 2 0 2 2 2 2 3 604 605 102 73 3.553 45 1 0 2 1 1 1 1 605 606 88 72 3.047 40 2 0 2 1 2 3 5 606 607 207 112 2.732 48 1 1 1 1 2 2 3 607 608 153 70 3.188 43 2 0 1 1 2 2 3 608 609 119 65 2.684 32 1 0 2 2 1 1 1 609 610 113 90 2.567 48 2 0 1 2 3 3 6 610 611 100 72 2.977 41 1 0 1 2 4 3 9 611 612 138 74 4.439 39 1 0 1 2 4 4 10 612 613 125 78 4.063 43 1 0 2 2 1 2 2 613 614 104 77 3.433 45 2 1 1 2 1 2 2 614 615 107 76 3.267 33 2 1 1 2 4 4 10 615 616 118 84 3.682 37 1 0 2 2 2 1 2 616 617 114 69 2.692 38 2 0 2 1 1 1 1 617 618 95 68 3.300 46 2 0 1 2 4 3 9 618 619 121 79 2.980 47 2 0 1 1 1 1 1 619 620 105 69 3.331 39 1 0 1 1 2 2 3 620 Obs rrl lrl rl rme lme rre lre riop liop age diab gh 621 2 3 32 0 0 0.750 0.750 26 24 25.9 41.3 9.9 622 1 1 40 0 0 -2.000 -2.250 17 15 14.4 1.6 13.0 623 6 7 6 0 0 0.125 0.250 12 11 19.4 26.0 9.1 624 4 3 27 0 0 -3.500 -4.500 17 16 16.2 15.9 16.7 625 1 1 40 0 0 0.000 0.000 17 16 14.8 5.8 7.7 626 2 3 32 0 0 -1.625 -1.375 14 11 23.2 9.2 13.5 627 1 1 40 0 0 -6.375 -5.500 15 17 13.2 7.1 9.9 628 2 2 35 0 0 -1.000 -1.250 20 18 10.0 20.3 11.7 629 1 1 40 0 0 0.375 0.625 14 15 9.1 3.2 11.5 630 1 1 40 0 0 -1.500 -1.750 12 14 9.0 12.8 13.0 631 2 2 35 0 0 -4.500 -5.000 20 21 20.1 6.4 10.9 632 2 2 35 0 0 0.000 0.000 12 14 12.7 13.0 10.8 633 6 7 6 1 1 -1.000 -1.000 21 19 11.1 13.9 9.5 634 2 2 35 0 0 -1.625 -1.250 15 17 8.2 12.8 11.7 635 1 2 37 0 0 0.000 0.000 16 17 19.8 9.6 11.4 636 1 1 40 0 0 -0.625 -0.625 18 18 8.7 7.2 9.7 637 1 1 40 0 0 -1.250 -0.750 15 16 15.1 7.0 11.8 638 4 4 26 0 0 -1.375 -1.125 18 16 11.5 14.0 13.3 639 3 4 27 0 0 -1.250 -1.000 16 16 13.2 16.8 9.0 640 1 1 40 0 0 -3.500 -3.500 18 16 21.0 4.5 8.8 641 1 1 40 0 0 -3.500 -4.500 23 23 21.5 3.2 10.8 642 3 2 32 0 0 -1.500 -0.875 12 13 16.3 6.3 13.0 643 2 2 35 0 0 -2.125 -5.500 19 20 14.1 10.8 10.3 644 4 4 26 0 0 -5.125 -6.250 16 14 29.1 13.8 12.6 645 3 3 31 0 0 -3.500 -4.375 15 15 3.9 16.8 13.3 646 1 1 40 0 0 -3.250 -2.125 11 12 17.9 5.2 11.7 647 1 1 40 0 0 -3.125 -2.500 16 15 17.4 5.9 9.0 648 4 4 26 0 0 -1.000 -1.500 16 16 18.7 15.3 18.1 649 1 1 40 0 0 0.000 0.000 20 20 22.8 2.5 11.0 650 1 3 34 0 0 -7.375 -7.750 20 19 10.4 9.1 14.7 651 2 2 35 0 0 -1.250 -1.875 9 11 18.4 10.0 12.4 Obs sbp dbp bmi pr sex prot dose resid rerl lerl rl2 id 621 192 90 4.034 49 1 1 1 1 2 2 3 621 622 122 88 3.859 57 2 0 2 2 1 1 1 622 623 218 93 3.603 44 1 1 1 2 4 4 10 623 624 110 72 3.202 44 2 0 1 1 3 2 5 624 625 132 58 3.409 46 2 1 2 1 1 1 1 625 626 94 73 3.269 57 1 0 2 1 2 2 3 626 627 97 80 3.003 29 2 0 1 1 1 1 1 627 628 116 72 3.292 47 2 0 2 1 2 2 3 628 629 106 72 3.177 46 2 0 2 1 1 1 1 629 630 112 82 3.633 39 2 0 1 1 1 1 1 630 631 122 83 3.450 44 2 0 2 1 2 2 3 631 632 115 84 3.553 30 1 0 1 1 2 2 3 632 633 115 76 3.305 38 1 0 2 1 4 4 10 633 634 130 59 3.498 48 1 0 2 1 2 2 3 634 635 135 73 2.649 46 1 0 1 1 1 2 2 635 636 109 83 3.325 27 1 0 1 1 1 1 1 636 637 95 65 2.931 31 1 0 2 1 1 1 1 637 638 94 76 2.646 32 1 0 1 1 3 3 6 638 639 136 91 3.954 42 1 0 2 1 2 3 5 639 640 132 82 3.616 36 1 0 2 1 1 1 1 640 641 108 76 3.677 29 1 0 1 1 1 1 1 641 642 101 76 3.019 37 2 0 1 1 2 2 3 642 643 119 89 3.753 39 1 0 2 1 2 2 3 643 644 121 89 3.677 47 2 0 2 1 3 3 6 644 645 113 69 3.193 39 1 0 2 1 2 2 3 645 646 118 78 3.409 39 1 0 2 1 1 1 1 646 647 126 79 4.011 27 1 0 2 1 1 1 1 647 648 123 80 2.854 36 1 0 1 1 3 3 6 648 649 118 80 2.770 32 1 0 1 1 1 1 1 649 650 139 71 3.125 43 2 0 1 2 1 2 2 650 651 122 80 3.880 36 1 0 2 1 2 2 3 651 Obs rrl lrl rl rme lme rre lre riop liop age diab gh 652 4 4.0 26 0 0 -0.625 -0.250 20 20 6.7 17.7 10.2 653 1 1.0 40 0 0 0.000 0.250 12 12 6.1 1.5 15.8 654 3 1.0 34 0 0 -1.375 -1.250 7 8 7.6 10.1 12.3 655 2 3.0 32 0 1 -8.625 -8.125 19 20 19.3 8.1 8.6 656 6 5.0 14 1 1 0.250 0.000 17 16 17.2 22.5 11.7 657 1 1.0 40 0 0 -5.000 -4.375 14 15 17.5 11.9 9.1 658 1 1.0 40 0 0 -2.000 -1.000 20 20 8.7 1.9 10.5 659 1 1.0 40 0 0 -1.875 -1.500 16 15 5.3 5.1 11.3 660 1 1.0 40 0 0 0.875 0.000 11 11 4.8 13.1 11.6 661 4 2.0 28 0 0 -2.500 -2.750 19 20 4.8 16.3 15.0 662 2 1.0 37 0 0 -0.375 -0.125 16 17 12.5 39.2 9.9 663 1 1.0 40 0 0 -1.750 -2.125 13 10 16.3 2.8 11.9 664 4 2.0 28 0 0 -0.500 0.500 19 21 5.9 12.2 13.8 665 6 6.0 13 1 1 -8.750 -6.250 20 20 9.2 13.8 15.8 666 4 3.0 27 0 0 -1.750 -1.375 14 16 7.0 13.8 11.5 667 2 2.0 35 0 0 -2.125 -2.625 21 22 4.5 32.8 11.3 668 1 1.0 40 0 0 -1.375 -1.375 18 18 9.3 6.8 18.2 669 1 1.0 40 0 0 -2.500 -2.250 16 15 10.6 4.4 9.5 670 1 1.0 40 0 0 0.000 0.000 20 20 3.0 8.3 14.5 671 2 2.0 35 0 0 -2.750 -2.875 18 18 2.0 14.0 14.6 672 6 6.0 13 0 1 0.750 1.250 18 18 3.2 22.0 9.6 673 1 1.0 40 0 0 -0.500 -0.125 18 18 6.6 5.0 14.2 674 2 2.0 35 0 0 0.000 0.000 19 20 5.4 8.7 11.8 675 3 3.0 31 0 0 -1.750 -1.500 16 15 9.4 14.3 16.1 676 1 1.0 40 0 0 -0.250 -0.500 11 10 13.5 5.7 11.2 677 2 2.0 35 0 0 -5.875 -5.750 16 18 12.4 5.9 12.6 678 1 2.0 37 0 0 -1.750 -1.250 11 11 3.2 11.2 14.5 679 2 1.0 37 0 0 0.125 -0.375 15 15 5.9 4.7 9.8 680 2 2.0 35 0 0 -0.625 -0.250 15 13 10.8 9.4 9.8 681 3 4.0 27 0 0 -2.000 -2.000 15 14 2.9 21.4 11.9 682 1 2.0 37 0 0 -1.750 -1.750 20 18 8.2 7.8 15.8 Obs sbp dbp bmi pr sex prot dose resid rerl lerl rl2 id 652 113 76 3.015 32 1 0 2 1 3 3 6 652 653 103 51 2.263 43 1 0 1 1 1 1 1 653 654 118 45 3.394 31 1 0 1 2 2 1 2 654 655 119 80 3.127 37 2 0 1 2 2 2 3 655 656 102 70 3.551 38 1 0 1 2 4 3 9 656 657 133 96 3.318 38 2 1 2 2 1 1 1 657 658 109 55 2.118 40 2 0 1 2 1 1 1 658 659 100 61 2.365 34 2 0 2 2 1 1 1 659 660 139 88 2.980 38 1 0 2 2 1 1 1 660 661 136 85 3.110 50 2 0 1 2 3 2 5 661 662 116 69 3.367 46 2 0 2 2 2 1 2 662 663 123 62 2.724 38 2 0 2 2 1 1 1 663 664 107 79 2.866 45 2 0 2 2 3 2 5 664 665 122 84 2.578 51 2 0 1 2 4 4 10 665 666 135 100 3.291 38 1 1 2 2 3 2 5 666 667 123 71 3.158 35 2 0 1 2 2 2 3 667 668 115 65 3.191 43 2 0 1 2 1 1 1 668 669 121 83 3.389 37 2 1 0 2 1 1 1 669 670 87 69 3.029 41 1 1 2 2 1 1 1 670 671 133 95 3.051 54 1 0 2 2 2 2 3 671 672 137 92 3.919 46 1 1 1 2 4 4 10 672 673 103 70 2.555 37 1 0 2 2 1 1 1 673 674 115 59 2.869 41 1 0 2 2 2 2 3 674 675 129 98 3.414 46 2 0 2 2 2 2 3 675 676 117 73 3.072 46 1 0 2 2 1 1 1 676 677 105 67 3.561 39 1 1 2 2 2 2 3 677 678 114 88 3.374 55 2 0 2 2 1 2 2 678 679 83 46 2.551 36 1 0 2 2 2 1 2 679 680 108 75 3.629 29 1 0 2 2 2 2 3 680 681 119 70 3.181 52 2 0 1 2 2 3 5 681 682 112 82 3.245 42 2 0 2 2 1 2 2 682 Obs rrl lrl rl rme lme rre lre riop liop age diab gh 683 4 3.0 27 0 0 -1.250 0.250 20 16 12.3 8.5 13.4 684 1 1.0 40 0 0 0.000 0.000 13 14 9.0 4.6 10.9 685 1 2.0 37 0 0 -1.875 -1.375 22 21 9.7 9.8 10.6 686 1 1.0 40 0 0 0.000 0.000 12 12 1.7 7.5 7.1 687 1 1.0 40 0 0 -2.250 -2.625 18 17 4.2 13.7 9.7 688 1 1.5 39 0 0 0.750 0.125 18 17 9.8 6.4 17.3 689 1 2.0 37 0 0 -2.750 -2.500 17 17 6.0 10.9 14.9 690 2 3.0 32 0 0 0.375 0.000 12 10 14.6 7.0 15.8 691 1 2.0 37 0 0 -0.625 -1.000 10 12 0.4 18.8 8.5 692 1 2.0 37 0 0 0.000 0.000 12 17 8.3 7.2 9.2 693 2 2.0 35 0 0 -1.750 -2.500 12 11 14.1 16.9 12.2 694 1 1.0 40 0 0 0.125 0.000 15 13 15.3 3.0 15.0 695 2 2.0 35 0 0 -0.125 -0.125 16 15 12.9 5.7 16.8 696 1 1.0 40 0 0 -1.250 -1.250 15 15 15.2 9.0 10.0 697 3 3.0 31 0 0 0.000 -0.125 11 11 25.7 54.0 11.9 698 4 6.0 15 0 0 -0.125 -0.250 14 13 21.6 25.2 11.4 699 1 4.0 30 0 0 -5.625 -6.000 19 20 17.4 22.0 14.5 700 2 3.0 32 0 0 -2.750 -2.500 21 24 20.7 25.1 15.5 701 4 4.0 26 0 0 1.500 1.250 13 14 17.4 50.1 11.2 702 4 4.0 26 1 0 0.125 0.125 15 15 20.7 25.3 14.3 703 3 3.0 31 0 0 -7.750 -8.125 21 17 9.8 16.2 21.0 704 4 4.0 26 0 0 -0.500 0.125 16 16 28.8 13.1 14.7 705 4 4.0 26 0 0 -3.250 -3.000 13 15 13.9 18.1 15.8 706 3 6.0 16 0 1 -4.125 -4.000 13 14 8.4 34.1 11.9 707 3 3.0 31 0 0 -2.875 -3.375 13 11 9.8 10.7 13.4 708 2 1.0 37 0 0 -0.875 -1.250 10 12 29.5 7.3 14.8 709 5 4.0 21 1 0 -2.500 -1.750 11 11 6.2 27.9 11.0 710 1 1.0 40 0 0 0.000 0.000 13 15 16.8 1.9 7.7 711 4 4.0 26 0 0 -1.500 -1.375 16 12 12.1 21.2 7.3 712 1 1.0 40 0 0 -0.875 -1.125 12 11 13.0 3.7 6.8 713 2 3.0 32 0 0 -3.875 -3.875 13 15 22.1 5.8 12.0 Obs sbp dbp bmi pr sex prot dose resid rerl lerl rl2 id 683 131 88 3.765 47 1 0 2 2 3 2 5 683 684 127 78 2.867 34 1 0 2 2 1 1 1 684 685 114 86 4.278 39 2 0 2 2 1 2 2 685 686 95 68 2.285 30 1 0 2 2 1 1 1 686 687 106 60 3.216 33 1 0 2 2 1 1 1 687 688 111 79 3.085 38 1 0 2 2 1 2 2 688 689 108 65 2.973 38 2 0 2 2 1 2 2 689 690 103 78 3.146 46 2 0 1 2 2 2 3 690 691 104 79 3.329 37 2 0 2 2 1 2 2 691 692 116 71 3.776 38 1 0 2 2 1 2 2 692 693 108 72 4.599 42 2 0 1 2 2 2 3 693 694 117 78 2.989 42 2 0 1 2 1 1 1 694 695 136 67 3.704 41 1 0 1 2 2 2 3 695 696 132 77 3.285 38 1 0 2 2 1 1 1 696 697 146 60 3.840 28 1 1 1 2 2 2 3 697 698 123 59 2.923 38 2 0 2 2 3 4 9 698 699 134 80 3.314 46 2 0 1 2 1 3 4 699 700 158 77 3.431 46 2 0 1 2 2 2 3 700 701 116 61 3.225 45 2 0 2 2 3 3 6 701 702 127 81 3.080 53 1 1 1 2 3 3 6 702 703 114 85 3.158 49 2 1 1 2 2 2 3 703 704 115 90 5.602 31 1 0 1 2 3 3 6 704 705 122 76 2.790 41 1 0 2 2 3 3 6 705 706 152 74 3.308 52 2 0 1 2 2 4 8 706 707 128 66 3.788 34 1 0 2 2 2 2 3 707 708 102 73 3.456 39 1 0 2 2 2 1 2 708 709 112 73 3.257 38 1 0 1 2 3 3 6 709 710 133 91 3.024 38 1 0 1 2 1 1 1 710 711 103 72 4.365 50 2 0 1 2 3 3 6 711 712 110 78 3.072 42 2 0 2 2 1 1 1 712 713 108 71 2.698 44 2 0 1 2 2 2 3 713 Obs rrl lrl rl rme lme rre lre riop liop age diab gh 714 6 6.0 13 0 0 -3.250 -3.375 15 16 6.9 29.2 17.9 715 1 1.0 40 0 0 -8.000 -7.375 10 10 17.8 4.2 11.3 716 3 4.0 27 0 0 -0.500 -1.250 19 18 6.0 11.3 12.3 717 1 1.0 40 0 0 0.000 -0.375 12 11 12.3 3.2 16.3 718 4 6.0 15 0 0 -0.750 0.000 12 14 29.9 33.3 12.3 719 3 3.0 31 0 0 -2.750 -2.125 13 11 22.7 15.3 12.4 720 2 2.0 35 0 0 -7.875 -8.625 18 19 9.9 18.7 10.6 Obs sbp dbp bmi pr sex prot dose resid rerl lerl rl2 id 714 106 71 3.383 46 2 0 1 2 4 4 10 714 715 118 72 3.013 38 1 0 1 2 1 1 1 715 716 116 84 3.839 42 2 0 2 2 2 3 5 716 717 139 81 3.010 44 1 0 2 2 1 1 1 717 718 156 76 3.590 42 1 0 1 2 3 4 9 718 719 121 83 3.322 28 1 0 1 2 2 2 3 719 720 132 80 3.542 37 1 0 2 2 2 2 3 720
A.4 SAS Data file for Example in Section 10.2 on Zero-Inflated Count Data for Frequency of a Side Effect (variable sefreq), from Min and Agresti (2005).
options ls=75 ps=65;
data count;
input patient visit sefreq time trt2@;
1 1 0 21 1
1 2 0 5 1
1 3 0 8 1
1 4 0 16 1
1 5 0 19 1
1 6 0 36 1
3 1 0 23 0
3 2 0 24 0
3 3 0 7 0
3 4 0 14 0
3 5 0 28 0
3 6 0 28 0
4 1 1 23 0
4 2 4 10 0
4 3 1 8 0
4 4 0 13 0
4 5 0 29 0
4 6 0 27 0
5 1 0 27 1
5 2 0 7 1
5 3 0 11 1
5 4 1 15 1
5 5 1 27 1
5 6 0 31 1
6 1 0 27 0
6 2 0 7 0
6 3 0 7 0
6 4 0 14 0
6 5 0 36 0
6 6 0 22 0
7 1 0 27 0
7 2 0 5 0
7 3 0 7 0
7 4 0 14 0
7 5 0 36 0
7 6 0 27 0
8 1 0 28 0
8 2 0 10 0
8 3 0 11 0
8 4 0 14 0
8 5 0 28 0
8 6 0 31 0
21 1 0 31 0
21 2 0 8 0
21 3 0 7 0
21 4 0 14 0
21 5 0 28 0
21 6 0 27 0
22 1 0 31 1
22 2 0 8 1
22 3 0 7 1
22 4 0 14 1
22 5 0 28 1
22 6 0 27 1
23 1 0 38 1
23 2 0 7 1
23 3 0 7 1
23 4 0 15 1
23 5 0 28 1
23 6 0 28 1
24 1 0 39 0
24 2 1 7 0
24 3 0 10 0
24 4 2 11 0
24 5 0 28 0
24 6 0 35 0
25 1 0 25 1
25 2 2 7 1
25 3 1 7 1
25 4 2 15 1
25 5 0 28 1
25 6 3 27 1
500 1 0 28 0
500 2 0 8 0
500 3 0 8 0
500 4 0 12 0
500 5 0 28 0
500 6 0 34 0
501 1 0 28 1
501 2 1 8 1
501 3 1 8 1
501 4 0 12 1
501 5 2 28 1
501 6 4 34 1
504 1 0 33 0
504 2 0 8 0
504 3 0 9 0
504 4 0 11 0
504 5 0 30 0
504 6 0 28 0
505 1 0 33 1
505 2 0 7 1
505 3 0 10 1
505 4 0 11 1
505 5 0 30 1
505 6 0 33 1
506 1 0 30 1
506 2 0 7 1
506 3 0 7 1
506 4 2 19 1
506 5 5 23 1
506 6 2 27 1
507 1 0 28 0
507 2 0 7 0
507 3 0 9 0
507 4 0 15 0
507 5 0 27 0
507 6 0 32 0
509 1 0 33 0
509 2 0 9 0
509 3 0 8 0
509 4 0 14 0
509 5 0 25 0
509 6 0 27 0
1200 1 0 27 0
1200 2 0 5 0
1200 3 0 7 0
1200 4 0 14 0
1200 5 0 21 0
1200 6 0 28 0
1202 1 0 28 1
1202 2 0 7 1
1202 3 0 8 1
1202 4 0 15 1
1202 5 0 32 1
1202 6 0 21 1
1206 1 0 20 0
1206 2 0 7 0
1206 3 0 7 0
1206 4 0 14 0
1206 5 0 36 0
1206 6 0 21 0
1208 1 0 27 1
1208 2 0 7 1
1208 3 0 7 1
1208 4 0 13 1
1208 5 1 21 1
1208 6 0 14 1
1210 1 0 24 0
1210 2 0 7 0
1210 3 0 7 0
1210 4 0 14 0
1210 5 0 14 0
1210 6 0 21 0
1212 1 0 27 1
1212 2 1 7 1
1212 3 0 7 1
1212 4 0 14 1
1212 5 1 14 1
1212 6 0 21 1
1217 1 0 34 0
1217 2 0 7 0
1217 3 0 7 0
1217 4 0 14 0
1217 5 0 21 0
1217 6 0 28 0
1219 1 0 28 1
1219 2 1 7 1
1219 3 0 7 1
1219 4 0 11 1
1219 5 0 19 1
1219 6 0 29 1
1222 1 0 28 1
1222 2 0 8 1
1222 3 0 7 1
1222 4 0 18 1
1222 5 0 23 1
1222 6 0 26 1
1226 1 0 24 0
1226 2 0 7 0
1226 3 0 7 0
1226 4 0 14 0
1226 5 0 28 0
1226 6 0 32 0
1229 1 0 20 0
1229 2 2 7 0
1229 3 0 7 0
1229 4 0 14 0
1229 5 1 27 0
1229 6 0 37 0
1232 1 0 22 0
1232 2 0 7 0
1232 3 0 5 0
1232 4 0 15 0
1232 5 0 34 0
1232 6 0 21 0
1239 1 0 37 1
1239 2 0 6 1
1239 3 0 7 1
1239 4 0 14 1
1239 5 1 42 1
1239 6 0 28 1
1241 1 0 27 0
1241 2 0 8 0
1241 3 1 16 0
1241 4 0 14 0
1241 5 0 18 0
1241 6 0 31 0
1242 1 0 35 1
1242 2 0 7 1
1242 3 1 9 1
1242 4 0 21 1
1242 5 0 22 1
1242 6 0 40 1
1500 1 0 22 1
1500 2 0 7 1
1500 3 0 9 1
1500 4 0 6 1
1500 5 0 22 1
1500 6 0 40 1
1502 1 1 33 0
1502 2 0 8 0
1502 3 0 7 0
1502 4 0 15 0
1502 5 1 20 0
1502 6 0 21 0
1521 1 0 35 0
1521 2 0 7 0
1521 3 0 7 0
1521 4 1 14 0
1521 5 2 35 0
1521 6 1 35 0
1522 1 0 34 0
1522 2 0 9 0
1522 3 0 8 0
1522 4 0 12 0
1522 5 0 34 0
1522 6 0 20 0
1523 1 0 34 1
1523 2 0 7 1
1523 3 0 7 1
1523 4 1 7 1
1523 5 0 35 1
1523 6 0 35 1
1524 1 0 25 1
1524 2 0 7 1
1524 3 0 7 1
1524 4 0 7 1
1524 5 0 42 1
1524 6 0 14 1
1525 1 0 25 0
1525 2 0 7 0
1525 3 0 7 0
1525 4 0 14 0
1525 5 0 35 0
1525 6 0 35 0
1526 1 0 28 1
1526 2 0 7 1
1526 3 0 7 1
1526 4 0 7 1
1526 5 0 35 1
1526 6 0 35 1
3202 1 0 28 1
3202 2 0 7 1
3202 3 0 8 1
3202 4 0 13 1
3202 5 0 28 1
3202 6 1 26 1
3204 1 0 30 0
3204 2 0 8 0
3204 3 0 6 0
3204 4 0 14 0
3204 5 0 39 0
3204 6 0 15 0
3205 1 0 28 1
3205 2 0 8 1
3205 3 0 6 1
3205 4 3 14 1
3205 5 0 28 1
3205 6 0 32 1
3206 1 0 28 1
3206 2 0 8 1
3206 3 0 7 1
3206 4 0 18 1
3206 5 0 28 1
3206 6 0 29 1
3207 1 0 28 0
3207 2 0 7 0
3207 3 0 7 0
3207 4 0 14 0
3207 5 0 28 0
3207 6 0 42 0
3208 1 0 25 0
3208 2 0 7 0
3208 3 0 7 0
3208 4 1 16 0
3208 5 0 26 0
3208 6 0 38 0
3209 1 0 21 1
3209 2 0 7 1
3209 3 0 9 1
3209 4 0 14 1
3209 5 0 27 1
3209 6 0 27 1
3211 1 0 35 0
3211 2 0 7 0
3211 3 0 6 0
3211 4 0 13 0
3211 5 0 44 0
3211 6 0 21 0
3212 1 0 30 1
3212 2 0 7 1
3212 3 0 7 1
3212 4 0 13 1
3212 5 0 43 1
3212 6 0 21 1
3213 1 1 27 0
3213 2 0 6 0
3213 3 2 7 0
3213 4 0 14 0
3213 5 2 43 0
3213 6 1 21 0
3214 1 0 28 0
3214 2 0 7 0
3214 3 0 6 0
3214 4 0 16 0
3214 5 0 32 0
3214 6 0 30 0
3215 1 0 32 1
3215 2 0 7 1
3215 3 0 8 1
3215 4 0 16 1
3215 5 0 28 1
3215 6 0 28 1
3216 1 0 29 1
3216 2 0 8 1
3216 3 0 9 1
3216 4 0 13 1
3216 5 0 26 1
3216 6 0 28 1
3217 1 2 26 1
3217 2 1 8 1
3217 3 0 6 1
3217 4 1 14 1
3217 5 4 30 1
3217 6 2 30 1
4200 1 0 21 1
4200 2 1 7 1
4200 3 1 14 1
4200 4 1 7 1
4200 5 0 14 1
4200 6 0 21 1
4203 1 0 21 0
4203 2 0 7 0
4203 3 0 6 0
4203 4 0 15 0
4203 5 1 14 0
4203 6 0 21 0
4205 1 0 28 1
4205 2 0 7 1
4205 3 0 14 1
4205 4 0 13 1
4205 5 0 15 1
4205 6 0 21 1
4206 1 0 21 1
4206 2 0 7 1
4206 3 0 7 1
4206 4 0 14 1
4206 5 0 13 1
4206 6 0 29 1
4207 1 0 21 0
4207 2 0 7 0
4207 3 0 14 0
4207 4 0 13 0
4207 5 0 8 0
4207 6 0 28 0
4208 1 0 28 1
4208 2 0 7 1
4208 3 0 6 1
4208 4 0 15 1
4208 5 0 17 1
4208 6 0 25 1
4210 1 0 24 0
4210 2 0 8 0
4210 3 0 14 0
4210 4 0 4 0
4210 5 0 21 0
4210 6 0 36 0
4211 1 1 21 0
4211 2 0 7 0
4211 3 1 7 0
4211 4 0 7 0
4211 5 0 20 0
4211 6 0 33 0
4302 1 0 21 1
4302 2 2 7 1
4302 3 3 7 1
4302 4 6 15 1
4302 5 4 27 1
4302 6 3 21 1
4303 1 1 21 0
4303 2 0 7 0
4303 3 0 7 0
4303 4 0 15 0
4303 5 0 27 0
4303 6 0 21 0
4308 1 0 28 1
4308 2 0 7 1
4308 3 0 8 1
4308 4 0 10 1
4308 5 0 17 1
4308 6 0 28 1
4309 1 0 27 1
4309 2 1 7 1
4309 3 2 7 1
4309 4 1 7 1
4309 5 1 21 1
4309 6 2 28 1
4311 1 0 21 0
4311 2 0 6 0
4311 3 1 7 0
4311 4 1 14 0
4311 5 0 21 0
4311 6 0 28 0
4312 1 1 22 1
4312 2 0 5 1
4312 3 0 7 1
4312 4 0 14 1
4312 5 0 21 1
4312 6 1 28 1
4313 1 0 27 1
4313 2 0 7 1
4313 3 0 7 1
4313 4 0 7 1
4313 5 1 21 1
4313 6 0 28 1
4315 1 0 26 0
4315 2 1 7 0
4315 3 1 7 0
4315 4 0 7 0
4315 5 0 21 0
4315 6 0 28 0
4316 1 0 22 0
4316 2 1 7 0
4316 3 0 6 0
4316 4 0 14 0
4316 5 0 14 0
4316 6 0 21 0
4317 1 0 21 1
4317 2 0 8 1
4317 3 0 13 1
4317 4 0 7 1
4317 5 0 21 1
4317 6 0 21 1
6001 1 0 28 0
6001 2 0 7 0
6001 3 0 7 0
6001 4 0 14 0
6001 5 0 35 0
6001 6 0 21 0
6002 1 0 28 1
6002 2 2 7 1
6002 3 1 6 1
6002 4 3 8 1
6002 5 1 35 1
6002 6 1 28 1
6003 1 0 28 1
6003 2 0 7 1
6003 3 0 8 1
6003 4 0 18 1
6003 5 0 24 1
6003 6 2 28 1
6004 1 0 29 1
6004 2 1 7 1
6004 3 1 7 1
6004 4 0 18 1
6004 5 1 22 1
6004 6 0 35 1
6006 1 0 28 1
6006 2 0 10 1
6006 3 0 2 1
6006 4 1 23 1
6006 5 4 20 1
6006 6 4 28 1
6008 1 0 27 0
6008 2 0 7 0
6008 3 0 5 0
6008 4 0 9 0
6008 5 0 40 0
6008 6 0 23 0
6009 1 0 27 1
6009 2 0 7 1
6009 3 0 4 1
6009 4 2 17 1
6009 5 2 30 1
6009 6 3 30 1
6010 1 0 28 1
6010 2 0 7 1
6010 3 1 7 1
6010 4 0 14 1
6010 5 4 29 1
6010 6 5 26 1
6012 1 0 28 0
6012 2 0 7 0
6012 3 0 7 0
6012 4 0 20 0
6012 5 0 25 0
6012 6 0 32 0
6013 1 0 28 0
6013 2 0 7 0
6013 3 0 7 0
6013 4 0 15 0
6013 5 0 34 0
6013 6 0 27 0
6014 1 0 27 1
6014 2 0 7 1
6014 3 0 8 1
6014 4 0 14 1
6014 5 0 27 1
6014 6 0 32 1
6015 1 0 31 1
6015 2 0 7 1
6015 3 0 6 1
6015 4 1 15 1
6015 5 1 32 1
6015 6 2 29 1
6016 1 0 29 0
6016 2 0 8 0
6016 3 1 8 0
6016 4 0 17 0
6016 5 0 29 0
6016 6 0 28 0
6017 1 0 27 1
6017 2 0 7 1
6017 3 0 7 1
6017 4 0 22 1
6017 5 0 27 1
6017 6 0 26 1
6018 1 0 21 0
6018 2 0 7 0
6018 3 0 5 0
6018 4 0 19 0
6018 5 0 25 0
6018 6 0 31 0
6019 1 1 23 0
6019 2 0 5 0
6019 3 0 12 0
6019 4 0 14 0
6019 5 0 30 0
6019 6 0 37 0
6020 1 1 32 0
6020 2 0 10 0
6020 3 0 7 0
6020 4 0 11 0
6020 5 0 34 0
6020 6 0 16 0
6021 1 0 35 0
6021 2 0 5 0
6021 3 0 9 0
6021 4 0 21 0
6021 5 0 26 0
6021 6 0 34 0
6023 1 0 28 1
6023 2 0 9 1
6023 3 0 6 1
6023 4 0 12 1
6023 5 0 35 1
6023 6 0 22 1
6025 1 0 26 1
6025 2 0 9 1
6025 3 0 5 1
6025 4 0 21 1
6025 5 0 26 1
6025 6 0 28 1
6026 1 0 29 0
6026 2 2 6 0
6026 3 2 7 0
6026 4 2 23 0
6026 5 2 19 0
6026 6 0 35 0
6027 1 0 28 1
6027 2 0 7 1
6027 3 0 7 1
6027 4 0 7 1
6027 5 0 21 1
6027 6 0 35 1
6028 1 1 21 0
6028 2 0 9 0
6028 3 0 13 0
6028 4 1 13 0
6028 5 0 35 0
6028 6 1 17 0
6029 1 0 36 0
6029 2 0 7 0
6029 3 0 9 0
6029 4 0 13 0
6029 5 0 28 0
6029 6 0 37 0
6030 1 0 22 1
6030 2 0 8 1
6030 3 0 7 1
6030 4 0 14 1
6030 5 0 28 1
6030 6 6 28 1
6031 1 0 23 1
6031 2 0 7 1
6031 3 2 7 1
6031 4 0 14 1
6031 5 0 28 1
6031 6 0 28 1
6032 1 0 22 0
6032 2 0 13 0
6032 3 0 4 0
6032 4 0 12 0
6032 5 0 28 0
6032 6 0 28 0
6036 1 0 36 0
6036 2 0 7 0
6036 3 0 6 0
6036 4 0 15 0
6036 5 0 35 0
6036 6 1 22 0
6037 1 0 24 1
6037 2 0 7 1
6037 3 0 7 1
6037 4 0 18 1
6037 5 0 24 1
6037 6 1 25 1
6039 1 0 28 0
6039 2 0 7 0
6039 3 0 4 0
6039 4 0 14 0
6039 5 1 37 0
6039 6 2 29 0
7101 1 0 21 1
7101 2 0 13 1
7101 3 0 10 1
7101 4 0 14 1
7101 5 0 21 1
7101 6 0 31 1
7102 1 0 26 1
7102 2 1 12 1
7102 3 0 7 1
7102 4 2 8 1
7102 5 4 23 1
7102 6 1 35 1
7103 1 0 28 0
7103 2 1 10 0
7103 3 0 7 0
7103 4 0 11 0
7103 5 0 21 0
7103 6 0 35 0
7104 1 0 29 0
7104 2 0 6 0
7104 3 0 10 0
7104 4 0 14 0
7104 5 0 32 0
7104 6 0 14 0
7105 1 0 27 1
7105 2 0 5 1
7105 3 0 7 1
7105 4 0 9 1
7105 5 0 26 1
7105 6 0 35 1
7106 1 0 35 1
7106 2 0 3 1
7106 3 0 5 1
7106 4 0 13 1
7106 5 0 28 1
7106 6 0 23 1
7109 1 0 28 1
7109 2 1 7 1
7109 3 1 7 1
7109 4 2 21 1
7109 5 0 28 1
7109 6 0 28 1
7112 1 0 33 0
7112 2 0 7 0
7112 3 0 7 0
7112 4 0 14 0
7112 5 0 22 0
7112 6 2 35 0
7115 1 0 28 0
7115 2 0 7 0
7115 3 0 7 0
7115 4 0 14 0
7115 5 0 28 0
7115 6 0 28 0
7120 1 0 28 0
7120 2 0 7 0
7120 3 0 7 0
7120 4 0 15 0
7120 5 0 27 0
7120 6 0 28 0
7121 1 0 25 0
7121 2 0 7 0
7121 3 0 7 0
7121 4 0 13 0
7121 5 0 32 0
7121 6 0 28 0
7122 1 0 28 1
7122 2 0 7 1
7122 3 0 7 1
7122 4 0 14 1
7122 5 0 28 1
7122 6 0 28 1
7123 1 0 28 0
7123 2 0 8 0
7123 3 0 7 0
7123 4 0 14 0
7123 5 0 20 0
7123 6 0 22 0
7124 1 0 26 1
7124 2 0 7 1
7124 3 2 7 1
7124 4 0 14 1
7124 5 2 28 1
7124 6 0 28 1
/*Poisson GLMM*/
proc nlmixed data=count qpoints=50;
parms b0=-5 b1=1 b2=0.2 s=1.2 ;
eta=b0 + b1*trt2 + b2*log(time) + u;
m=exp(eta);
model sefreq ~ poisson(m);
random u~normal(0,s*s) subject= patient;
run;
/*Random effect type II negative binomial model*/
proc nlmixed data=count qpoints=50;
parms b0=-4 b1=1 b2=1 s=1.5 a=1;
eta=b0 + b1*trt2 + b2*log(time) + u;
m=exp(eta);
p=a*m/(1+a*m);
ll=lgamma(sefreq+1/a)-lgamma(sefreq+1)-lgamma(1/a);
ll=ll + 1/a*log(1-p) + sefreq*log(p);
model sefreq~general(ll);
random u~normal(0,s*s) subject= patient;
run;
/*Random effect zero-altered model*/
proc nlmixed data=count qpoints=50;
parms b0=-5 b1=1 b2=0.2 s1=1.256 b5=-1 ;
eta1 = b5+ (b0 + b1*trt2 + b2*log(time))+u1;
expeta1 = exp(eta1);
eta2 =b0 + b1*trt2 + b2*log(time) + u1;
expeta2 = exp(eta2);
if sefreq=0 then ll= -expeta1;
else ll=log(1-exp(-expeta1))+sefreq*eta2-log(exp(expeta2)-1)-lgamma(sefreq+1);
model sefreq ~ general(ll);
contrast 'b5=0' b5;
random u1 ~ normal(0,s1*s1) subject=patient;
run;
/*Random effect hurdle model (1) with log(time) in the first model*/
proc nlmixed data=count qpoints=50;
bounds -1<= s12 <=1;
parms b0=-2.874 b1=0.8945 b2=0.02137 s1=1.6469
b3=-2.844 b4=0.9626 b5=0.5398 s2=0.7060 s12=0.8481;
eta1 = b0 + b1*trt2 + b2*log(time) + u1;
expeta1 = exp(eta1);
eta2 = b3 + b4*trt2 + b5*log(time) + u2;
expeta2 = exp(eta2);
if sefreq=0 then ll= -log(1+expeta1);
else
ll=eta1-log(1+expeta1)+sefreq*eta2-log(exp(expeta2)-1)-lgamma(sefreq+1);
model sefreq ~ general(ll);
contrast 'b2=0' b2;
random u1 u2 ~ normal([0,0],[s1*s1,s12*s1*s2,s2*s2]) subject=patient;
run;
/*Random effect hurdle model (2) without log(time) in the first model*/
proc nlmixed data=count qpoints=50;
bounds -1<= s12 <=1;
parms b0=-2.874 b1=0.8945 s1=1.6469
b3=-2.844 b4=0.9626 b5=0.5398 s2=0.7060 s12=0.8481;
eta1 = b0 + b1*trt2 + u1;
expeta1 = exp(eta1);
eta2 = b3 + b4*trt2 + b5*log(time) + u2;
expeta2 = exp(eta2);
if sefreq=0 then ll= -log(1+expeta1);
else
ll=eta1-log(1+expeta1)+sefreq*eta2-log(exp(expeta2)-1)-lgamma(sefreq+1);
model sefreq ~ general(ll);
random u1 u2 ~ normal([0,0],[s1*s1,s12*s1*s2,s2*s2]) subject=patient;
run;
/*Random effect cumulative logit model with b2*/
proc nlmixed data=count qpoints=50;
parms i1=3.3 i2=4.62 i3=5.9 i4=6.37 i5=7.55
b1=0.9773 b2=0.153 s1=1.749;
eta1 = i1 - b1*trt2 - b2*log(time) - u;
eta2 = i2 - b1*trt2 - b2*log(time) - u;
eta3 = i3 - b1*trt2 - b2*log(time) - u;
eta4 = i4 - b1*trt2 - b2*log(time) - u;
eta5 = i5 - b1*trt2 - b2*log(time) - u;
if (sefreq=0) then z=1/(1+exp(-eta1));
else if (sefreq=1) then z=1/(1+exp(-eta2))-1/(1+exp(-eta1));
else if (sefreq=2) then z=1/(1+exp(-eta3))-1/(1+exp(-eta2));
else if (sefreq=3) then z=1/(1+exp(-eta4))-1/(1+exp(-eta3));
else if (sefreq=4) then z=1/(1+exp(-eta5))-1/(1+exp(-eta4));
else z=1-1/(1+exp(-eta5));
if (z>1e-8) then ll=log(z);
else ll=-1e100;
model sefreq ~ general(ll);
random u ~ normal(0,s1*s1) subject=patient;
run;
/*Random effect cumulative logit model without b2*/
proc nlmixed data=count qpoints=50;
parms i1=3.3 i2=4.62 i3=5.9 i4=6.37 i5=7.55
b1=0.9773 s1=1.749;
eta1 = i1 - b1*trt2 - u;
eta2 = i2 - b1*trt2 - u;
eta3 = i3 - b1*trt2 - u;
eta4 = i4 - b1*trt2 - u;
eta5 = i5 - b1*trt2 - u;
if (sefreq=0) then z=1/(1+exp(-eta1));
else if (sefreq=1) then z=1/(1+exp(-eta2))-1/(1+exp(-eta1));
else if (sefreq=2) then z=1/(1+exp(-eta3))-1/(1+exp(-eta2));
else if (sefreq=3) then z=1/(1+exp(-eta4))-1/(1+exp(-eta3));
else if (sefreq=4) then z=1/(1+exp(-eta5))-1/(1+exp(-eta4));
else z=1-1/(1+exp(-eta5));
if (z>1e-8) then ll=log(z);
else ll=-1e100;
model sefreq ~ general(ll);
random u ~ normal(0,s1*s1) subject=patient;
run;
A.5 SAS Data File for Example in Section 10.3 on Heterogeneity in Multi-Center Clinical Trials, from Hartzel, Liu, and Agresti (2001).
data ordinal;
do center = 1 to 8;
do trt = 1 to 0 by -1;
do resp = 3 to 1 by -1;
input count @@;
output;
end;
end;
end;
datalines;
13 7 6 1 1 10 2 5 10 2 2 1
11 23 7 2 8 2 7 11 8 0 3 2
15 3 5 1 1 5 13 5 5 4 0 1
7 4 13 1 1 11 15 9 2 3 2 2
run;
proc nlmixed data=ordinal qpoints=15;
** To maintain the threshold ordering define thresholds such that **;
** threshold 1 = 0 and threshold 2 = i2, where i2 > 0. **;
** Use starting value of 0 for sig_cb. **;
bounds i2>0; parms sig_cb=0;
eta1 = c - b*trt;
eta2 = i2 - c - b*trt;
if (resp=1) then z = 1/(1+exp(-eta1));
else if (resp=2) then z = 1/(1+exp(-eta2)) - 1/(1+exp(-eta1));
else z = 1 - 1/(1+exp(-eta2));
if (z > 1e-8) then ll = count*log(z); ** Check for small values of z **;
else ll=-1e100;
model resp ~ general(ll); ** Define general log-likelihood. **;
random c b ~ normal([gamma,beta],[sig_c*sig_c, sig_cb, sig_b*sig_b])
subject = center out = out1; ** OUT1 contains predicted center- **;
** specific cumulative log odds ratios **;
run;
A.6 SAS data file for using PROC NLMIXED to fit adjacent-categories logit random effects model to a 3x3x3x3 generalization of Table 8.7 on movie reviewers, for analyses conducted in Hartzel, Agresti, and Caffo (2001).
/*Input Data*/
data movie;
do lyons = 1 to 3;
do siskel = 1 to 3;
do medved = 1 to 3;
do ebert = 1 to 3;
input count @@;
if (count ne 0) then output;
end;
end;
end;
end;
datalines;
15 2 0 2 0 0 8 0 2
0 3 0 2 3 0 3 2 1
1 0 1 1 2 1 1 1 2
2 0 0 2 1 0 0 1 0
1 1 0 0 1 0 0 1 0
1 0 0 0 0 0 0 1 1
3 1 1 0 0 0 4 1 0
0 0 1 0 0 0 1 3 0
0 0 1 2 1 0 2 2 4
;
/*Create subject variable "movie" and code dummy variables*/
data movie;
set movie;
movie = _n_;
x1 = 1; x2 = 0; x3 = 0; resp = siskel; output;
x1 = 0; x2 = 1; x3 = 0; resp = ebert; output;
x1 = 0; x2 = 0; x3 = 1; resp = lyons; output;
x1 = 0; x2 = 0; x3 = 0; resp = medved; output;
drop lyons siskel medved ebert;
/*Recode responses as multivariate Bernoulli*/
data movie;
set movie;
if resp = 1 then do;
y1=1; y2=0; y3=0;
end;
else if resp = 2 then do;
y1=0; y2=1; y3=0;
end;
else if resp = 3 then do;
y1=0; y2=0; y3=1;
end;
else delete;
output;
drop resp;
run;
/*Run NLMIXED; add qpoints=? to change number of quadrature points*/
/*Adjacent Category Logit Model*/
proc nlmixed data=movie;
/*Code linear predictors*/
eta1 = alpha1 + x1 * beta1 + x2 * beta2 + x3 * beta3 + u1;
eta2 = alpha2 + x1 * beta1 + x2 * beta2 + x3 * beta3 + u2;
/*Adjacent category model*/
pi1 = exp(eta1+eta2) / (exp(eta1+eta2) + exp(eta2) + 1);
pi2 = exp(eta2) / (exp(eta1+eta2) + exp(eta2) + 1);
pi3 = 1 / (exp(eta1+eta2) + exp(eta2) + 1);
/*Define likelihood*/
z = (pi1**y1)*(pi2**y2)*(pi3**y3);
if (z > 1e-8) then ll = log(z);
else ll=-1e100;
model y1 ~ general(ll);
/*Specify random effect distribution*/
random u1 u2 ~ normal([0,0], [s1*s1, cov12 ,s2*s2]) subject = movie;
/*Replicate data points with identical likelihood values*/
replicate count;
run;
/*Adjacent Category Logit Model with Uncorrelated Random Effects*/
proc nlmixed data=movie qpoints = 50;
/*Code linear predictors*/
eta1 = alpha1 + x1 * beta1 + x2 * beta2 + x3 * beta3 + u1;
eta2 = alpha2 + x1 * beta1 + x2 * beta2 + x3 * beta3 + u2;
/*Adjacent category*/
pi1 = exp(eta1+eta2) / (exp(eta1+eta2) + exp(eta2) + 1);
pi2 = exp(eta2) / (exp(eta1+eta2) + exp(eta2) + 1);
pi3 = 1 / (exp(eta1+eta2) + exp(eta2) + 1);
z = (pi1**y1)*(pi2**y2)*(pi3**y3);
if (z > 1e-8) then ll = log(z);
else ll=-1e100;
model y1 ~ general(ll);
random u1 u2 ~ normal([0,0], [s1*s1, 0 ,s2*s2]) subject = movie;
replicate count;
run;
This page is maintained by Alan Agresti.
Last Updated: April 8, 2010.
(Our Department Home Page)