The mispitools package provides tools for missing person identification combining genetic and non-genetic evidence within a Bayesian framework. It implements likelihood ratio (LR) calculations for various types of evidence commonly used in forensic identification.
In missing person identification, we compare two hypotheses:
The likelihood ratio (LR) quantifies the relative support for H1 vs H2:
\[LR = \frac{P(Evidence | H1)}{P(Evidence | H2)}\]
Multiple pieces of evidence can be combined by multiplying their LRs:
\[LR_{total} = LR_{genetic} \times LR_{sex} \times LR_{age} \times LR_{color}\]
Calculate the LR for sex evidence when the missing person is female:
library(mispitools)
# LR for sex evidence
# H1: MP is female, POI observed as female
# eps = probability of sex observation error
lr_sex(LR = TRUE, H = 1, eps = 0.05)
#> Sexo LRs
#> 1 F 1.9
#> 2 F 1.9
#> 3 F 1.9
#> 4 F 1.9
#> 5 F 1.9
#> 6 F 1.9
#> 7 F 1.9
#> 8 F 1.9
#> 9 F 1.9
#> 10 F 1.9
#> 11 F 1.9
#> 12 F 1.9
#> 13 F 1.9
#> 14 F 1.9
#> 15 F 1.9
#> 16 F 1.9
#> 17 F 1.9
#> 18 F 1.9
#> 19 F 1.9
#> 20 F 1.9
#> 21 F 1.9
#> 22 F 1.9
#> 23 F 1.9
#> 24 F 1.9
#> 25 F 1.9
#> 26 F 1.9
#> 27 F 1.9
#> 28 F 1.9
#> 29 F 1.9
#> 30 F 1.9
#> 31 F 1.9
#> 32 F 1.9
#> 33 F 1.9
#> 34 F 1.9
#> 35 F 1.9
#> 36 F 1.9
#> 37 F 1.9
#> 38 F 1.9
#> 39 M 0.1
#> 40 F 1.9
#> 41 F 1.9
#> 42 F 1.9
#> 43 F 1.9
#> 44 F 1.9
#> 45 F 1.9
#> 46 F 1.9
#> 47 F 1.9
#> 48 F 1.9
#> 49 F 1.9
#> 50 F 1.9
#> 51 F 1.9
#> 52 F 1.9
#> 53 F 1.9
#> 54 F 1.9
#> 55 F 1.9
#> 56 F 1.9
#> 57 F 1.9
#> 58 F 1.9
#> 59 F 1.9
#> 60 F 1.9
#> 61 F 1.9
#> 62 F 1.9
#> 63 F 1.9
#> 64 F 1.9
#> 65 F 1.9
#> 66 F 1.9
#> 67 F 1.9
#> 68 F 1.9
#> 69 F 1.9
#> 70 F 1.9
#> 71 F 1.9
#> 72 F 1.9
#> 73 F 1.9
#> 74 F 1.9
#> 75 F 1.9
#> 76 F 1.9
#> 77 F 1.9
#> 78 F 1.9
#> 79 F 1.9
#> 80 F 1.9
#> 81 F 1.9
#> 82 F 1.9
#> 83 F 1.9
#> 84 F 1.9
#> 85 F 1.9
#> 86 F 1.9
#> 87 F 1.9
#> 88 F 1.9
#> 89 F 1.9
#> 90 F 1.9
#> 91 F 1.9
#> 92 F 1.9
#> 93 F 1.9
#> 94 F 1.9
#> 95 F 1.9
#> 96 F 1.9
#> 97 F 1.9
#> 98 F 1.9
#> 99 F 1.9
#> 100 F 1.9
#> 101 F 1.9
#> 102 F 1.9
#> 103 F 1.9
#> 104 F 1.9
#> 105 F 1.9
#> 106 F 1.9
#> 107 F 1.9
#> 108 F 1.9
#> 109 F 1.9
#> 110 F 1.9
#> 111 F 1.9
#> 112 F 1.9
#> 113 M 0.1
#> 114 F 1.9
#> 115 F 1.9
#> 116 F 1.9
#> 117 F 1.9
#> 118 F 1.9
#> 119 F 1.9
#> 120 F 1.9
#> 121 F 1.9
#> 122 F 1.9
#> 123 M 0.1
#> 124 F 1.9
#> 125 F 1.9
#> 126 F 1.9
#> 127 F 1.9
#> 128 F 1.9
#> 129 F 1.9
#> 130 F 1.9
#> 131 M 0.1
#> 132 F 1.9
#> 133 F 1.9
#> 134 F 1.9
#> 135 F 1.9
#> 136 F 1.9
#> 137 F 1.9
#> 138 F 1.9
#> 139 F 1.9
#> 140 F 1.9
#> 141 F 1.9
#> 142 M 0.1
#> 143 F 1.9
#> 144 F 1.9
#> 145 F 1.9
#> 146 F 1.9
#> 147 F 1.9
#> 148 F 1.9
#> 149 M 0.1
#> 150 F 1.9
#> 151 F 1.9
#> 152 F 1.9
#> 153 F 1.9
#> 154 F 1.9
#> 155 F 1.9
#> 156 M 0.1
#> 157 F 1.9
#> 158 M 0.1
#> 159 F 1.9
#> 160 F 1.9
#> 161 F 1.9
#> 162 F 1.9
#> 163 F 1.9
#> 164 F 1.9
#> 165 F 1.9
#> 166 F 1.9
#> 167 F 1.9
#> 168 F 1.9
#> 169 F 1.9
#> 170 F 1.9
#> 171 F 1.9
#> 172 F 1.9
#> 173 F 1.9
#> 174 F 1.9
#> 175 F 1.9
#> 176 F 1.9
#> 177 F 1.9
#> 178 F 1.9
#> 179 F 1.9
#> 180 F 1.9
#> 181 F 1.9
#> 182 F 1.9
#> 183 F 1.9
#> 184 F 1.9
#> 185 M 0.1
#> 186 F 1.9
#> 187 F 1.9
#> 188 F 1.9
#> 189 F 1.9
#> 190 F 1.9
#> 191 F 1.9
#> 192 M 0.1
#> 193 F 1.9
#> 194 F 1.9
#> 195 F 1.9
#> 196 F 1.9
#> 197 F 1.9
#> 198 F 1.9
#> 199 M 0.1
#> 200 F 1.9
#> 201 F 1.9
#> 202 F 1.9
#> 203 F 1.9
#> 204 F 1.9
#> 205 F 1.9
#> 206 F 1.9
#> 207 F 1.9
#> 208 F 1.9
#> 209 F 1.9
#> 210 M 0.1
#> 211 F 1.9
#> 212 F 1.9
#> 213 F 1.9
#> 214 F 1.9
#> 215 F 1.9
#> 216 F 1.9
#> 217 F 1.9
#> 218 F 1.9
#> 219 F 1.9
#> 220 F 1.9
#> 221 F 1.9
#> 222 F 1.9
#> 223 F 1.9
#> 224 F 1.9
#> 225 F 1.9
#> 226 F 1.9
#> 227 F 1.9
#> 228 F 1.9
#> 229 F 1.9
#> 230 F 1.9
#> 231 F 1.9
#> 232 F 1.9
#> 233 F 1.9
#> 234 F 1.9
#> 235 F 1.9
#> 236 F 1.9
#> 237 F 1.9
#> 238 F 1.9
#> 239 F 1.9
#> 240 M 0.1
#> 241 F 1.9
#> 242 F 1.9
#> 243 F 1.9
#> 244 F 1.9
#> 245 M 0.1
#> 246 F 1.9
#> 247 F 1.9
#> 248 F 1.9
#> 249 M 0.1
#> 250 F 1.9
#> 251 F 1.9
#> 252 F 1.9
#> 253 F 1.9
#> 254 F 1.9
#> 255 F 1.9
#> 256 F 1.9
#> 257 F 1.9
#> 258 F 1.9
#> 259 F 1.9
#> 260 F 1.9
#> 261 F 1.9
#> 262 F 1.9
#> 263 F 1.9
#> 264 F 1.9
#> 265 F 1.9
#> 266 F 1.9
#> 267 F 1.9
#> 268 F 1.9
#> 269 F 1.9
#> 270 F 1.9
#> 271 F 1.9
#> 272 F 1.9
#> 273 F 1.9
#> 274 F 1.9
#> 275 F 1.9
#> 276 F 1.9
#> 277 F 1.9
#> 278 F 1.9
#> 279 F 1.9
#> 280 F 1.9
#> 281 F 1.9
#> 282 F 1.9
#> 283 F 1.9
#> 284 F 1.9
#> 285 F 1.9
#> 286 F 1.9
#> 287 F 1.9
#> 288 F 1.9
#> 289 F 1.9
#> 290 F 1.9
#> 291 F 1.9
#> 292 F 1.9
#> 293 M 0.1
#> 294 F 1.9
#> 295 F 1.9
#> 296 F 1.9
#> 297 F 1.9
#> 298 F 1.9
#> 299 F 1.9
#> 300 F 1.9
#> 301 F 1.9
#> 302 M 0.1
#> 303 F 1.9
#> 304 F 1.9
#> 305 F 1.9
#> 306 F 1.9
#> 307 F 1.9
#> 308 M 0.1
#> 309 F 1.9
#> 310 F 1.9
#> 311 F 1.9
#> 312 F 1.9
#> 313 F 1.9
#> 314 F 1.9
#> 315 F 1.9
#> 316 F 1.9
#> 317 F 1.9
#> 318 F 1.9
#> 319 F 1.9
#> 320 F 1.9
#> 321 F 1.9
#> 322 F 1.9
#> 323 F 1.9
#> 324 F 1.9
#> 325 F 1.9
#> 326 F 1.9
#> 327 F 1.9
#> 328 F 1.9
#> 329 F 1.9
#> 330 F 1.9
#> 331 F 1.9
#> 332 F 1.9
#> 333 M 0.1
#> 334 F 1.9
#> 335 F 1.9
#> 336 F 1.9
#> 337 F 1.9
#> 338 F 1.9
#> 339 F 1.9
#> 340 F 1.9
#> 341 M 0.1
#> 342 F 1.9
#> 343 F 1.9
#> 344 M 0.1
#> 345 F 1.9
#> 346 F 1.9
#> 347 F 1.9
#> 348 F 1.9
#> 349 F 1.9
#> 350 F 1.9
#> 351 F 1.9
#> 352 F 1.9
#> 353 F 1.9
#> 354 F 1.9
#> 355 M 0.1
#> 356 F 1.9
#> 357 F 1.9
#> 358 F 1.9
#> 359 F 1.9
#> 360 F 1.9
#> 361 F 1.9
#> 362 F 1.9
#> 363 F 1.9
#> 364 F 1.9
#> 365 F 1.9
#> 366 M 0.1
#> 367 F 1.9
#> 368 F 1.9
#> 369 F 1.9
#> 370 F 1.9
#> 371 M 0.1
#> 372 F 1.9
#> 373 F 1.9
#> 374 F 1.9
#> 375 F 1.9
#> 376 F 1.9
#> 377 F 1.9
#> 378 F 1.9
#> 379 F 1.9
#> 380 F 1.9
#> 381 F 1.9
#> 382 F 1.9
#> 383 F 1.9
#> 384 F 1.9
#> 385 F 1.9
#> 386 F 1.9
#> 387 F 1.9
#> 388 F 1.9
#> 389 M 0.1
#> 390 F 1.9
#> 391 F 1.9
#> 392 F 1.9
#> 393 F 1.9
#> 394 F 1.9
#> 395 F 1.9
#> 396 F 1.9
#> 397 F 1.9
#> 398 F 1.9
#> 399 F 1.9
#> 400 F 1.9
#> 401 F 1.9
#> 402 F 1.9
#> 403 F 1.9
#> 404 F 1.9
#> 405 F 1.9
#> 406 F 1.9
#> 407 F 1.9
#> 408 F 1.9
#> 409 F 1.9
#> 410 F 1.9
#> 411 F 1.9
#> 412 F 1.9
#> 413 M 0.1
#> 414 F 1.9
#> 415 F 1.9
#> 416 F 1.9
#> 417 F 1.9
#> 418 F 1.9
#> 419 F 1.9
#> 420 F 1.9
#> 421 F 1.9
#> 422 F 1.9
#> 423 F 1.9
#> 424 F 1.9
#> 425 F 1.9
#> 426 F 1.9
#> 427 F 1.9
#> 428 F 1.9
#> 429 F 1.9
#> 430 F 1.9
#> 431 F 1.9
#> 432 F 1.9
#> 433 F 1.9
#> 434 F 1.9
#> 435 F 1.9
#> 436 F 1.9
#> 437 F 1.9
#> 438 F 1.9
#> 439 F 1.9
#> 440 F 1.9
#> 441 F 1.9
#> 442 F 1.9
#> 443 F 1.9
#> 444 F 1.9
#> 445 M 0.1
#> 446 F 1.9
#> 447 F 1.9
#> 448 F 1.9
#> 449 F 1.9
#> 450 F 1.9
#> 451 F 1.9
#> 452 F 1.9
#> 453 M 0.1
#> 454 M 0.1
#> 455 F 1.9
#> 456 F 1.9
#> 457 F 1.9
#> 458 F 1.9
#> 459 F 1.9
#> 460 F 1.9
#> 461 F 1.9
#> 462 F 1.9
#> 463 F 1.9
#> 464 F 1.9
#> 465 F 1.9
#> 466 F 1.9
#> 467 F 1.9
#> 468 F 1.9
#> 469 F 1.9
#> 470 F 1.9
#> 471 F 1.9
#> 472 F 1.9
#> 473 F 1.9
#> 474 M 0.1
#> 475 F 1.9
#> 476 M 0.1
#> 477 F 1.9
#> 478 F 1.9
#> 479 F 1.9
#> 480 F 1.9
#> 481 F 1.9
#> 482 F 1.9
#> 483 F 1.9
#> 484 F 1.9
#> 485 F 1.9
#> 486 F 1.9
#> 487 F 1.9
#> 488 F 1.9
#> 489 F 1.9
#> 490 F 1.9
#> 491 F 1.9
#> 492 F 1.9
#> 493 F 1.9
#> 494 F 1.9
#> 495 F 1.9
#> 496 F 1.9
#> 497 F 1.9
#> 498 F 1.9
#> 499 F 1.9
#> 500 F 1.9
#> 501 F 1.9
#> 502 F 1.9
#> 503 F 1.9
#> 504 F 1.9
#> 505 F 1.9
#> 506 F 1.9
#> 507 F 1.9
#> 508 F 1.9
#> 509 F 1.9
#> 510 F 1.9
#> 511 F 1.9
#> 512 F 1.9
#> 513 F 1.9
#> 514 F 1.9
#> 515 F 1.9
#> 516 F 1.9
#> 517 F 1.9
#> 518 F 1.9
#> 519 F 1.9
#> 520 F 1.9
#> 521 F 1.9
#> 522 F 1.9
#> 523 F 1.9
#> 524 F 1.9
#> 525 F 1.9
#> 526 F 1.9
#> 527 F 1.9
#> 528 F 1.9
#> 529 F 1.9
#> 530 F 1.9
#> 531 F 1.9
#> 532 F 1.9
#> 533 F 1.9
#> 534 F 1.9
#> 535 F 1.9
#> 536 M 0.1
#> 537 F 1.9
#> 538 F 1.9
#> 539 F 1.9
#> 540 F 1.9
#> 541 F 1.9
#> 542 F 1.9
#> 543 M 0.1
#> 544 F 1.9
#> 545 F 1.9
#> 546 F 1.9
#> 547 F 1.9
#> 548 F 1.9
#> 549 M 0.1
#> 550 F 1.9
#> 551 F 1.9
#> 552 F 1.9
#> 553 F 1.9
#> 554 F 1.9
#> 555 F 1.9
#> 556 M 0.1
#> 557 F 1.9
#> 558 F 1.9
#> 559 F 1.9
#> 560 F 1.9
#> 561 F 1.9
#> 562 M 0.1
#> 563 F 1.9
#> 564 F 1.9
#> 565 F 1.9
#> 566 F 1.9
#> 567 F 1.9
#> 568 F 1.9
#> 569 F 1.9
#> 570 F 1.9
#> 571 F 1.9
#> 572 F 1.9
#> 573 F 1.9
#> 574 F 1.9
#> 575 F 1.9
#> 576 F 1.9
#> 577 F 1.9
#> 578 F 1.9
#> 579 F 1.9
#> 580 F 1.9
#> 581 F 1.9
#> 582 F 1.9
#> 583 F 1.9
#> 584 F 1.9
#> 585 F 1.9
#> 586 F 1.9
#> 587 F 1.9
#> 588 M 0.1
#> 589 M 0.1
#> 590 F 1.9
#> 591 M 0.1
#> 592 F 1.9
#> 593 F 1.9
#> 594 F 1.9
#> 595 F 1.9
#> 596 F 1.9
#> 597 F 1.9
#> 598 F 1.9
#> 599 F 1.9
#> 600 F 1.9
#> 601 F 1.9
#> 602 F 1.9
#> 603 F 1.9
#> 604 F 1.9
#> 605 F 1.9
#> 606 F 1.9
#> 607 F 1.9
#> 608 F 1.9
#> 609 F 1.9
#> 610 F 1.9
#> 611 F 1.9
#> 612 F 1.9
#> 613 F 1.9
#> 614 F 1.9
#> 615 F 1.9
#> 616 F 1.9
#> 617 F 1.9
#> 618 F 1.9
#> 619 F 1.9
#> 620 F 1.9
#> 621 F 1.9
#> 622 F 1.9
#> 623 F 1.9
#> 624 F 1.9
#> 625 F 1.9
#> 626 F 1.9
#> 627 F 1.9
#> 628 F 1.9
#> 629 F 1.9
#> 630 F 1.9
#> 631 F 1.9
#> 632 F 1.9
#> 633 F 1.9
#> 634 F 1.9
#> 635 F 1.9
#> 636 F 1.9
#> 637 F 1.9
#> 638 F 1.9
#> 639 F 1.9
#> 640 F 1.9
#> 641 F 1.9
#> 642 F 1.9
#> 643 F 1.9
#> 644 F 1.9
#> 645 F 1.9
#> 646 F 1.9
#> 647 F 1.9
#> 648 F 1.9
#> 649 F 1.9
#> 650 F 1.9
#> 651 F 1.9
#> 652 F 1.9
#> 653 F 1.9
#> 654 F 1.9
#> 655 F 1.9
#> 656 F 1.9
#> 657 F 1.9
#> 658 F 1.9
#> 659 F 1.9
#> 660 F 1.9
#> 661 F 1.9
#> 662 F 1.9
#> 663 F 1.9
#> 664 F 1.9
#> 665 F 1.9
#> 666 F 1.9
#> 667 F 1.9
#> 668 F 1.9
#> 669 F 1.9
#> 670 F 1.9
#> 671 F 1.9
#> 672 F 1.9
#> 673 F 1.9
#> 674 F 1.9
#> 675 F 1.9
#> 676 F 1.9
#> 677 F 1.9
#> 678 F 1.9
#> 679 F 1.9
#> 680 M 0.1
#> 681 F 1.9
#> 682 F 1.9
#> 683 F 1.9
#> 684 F 1.9
#> 685 F 1.9
#> 686 F 1.9
#> 687 F 1.9
#> 688 F 1.9
#> 689 F 1.9
#> 690 F 1.9
#> 691 F 1.9
#> 692 F 1.9
#> 693 F 1.9
#> 694 F 1.9
#> 695 F 1.9
#> 696 F 1.9
#> 697 F 1.9
#> 698 F 1.9
#> 699 F 1.9
#> 700 F 1.9
#> 701 F 1.9
#> 702 F 1.9
#> 703 F 1.9
#> 704 F 1.9
#> 705 F 1.9
#> 706 F 1.9
#> 707 F 1.9
#> 708 F 1.9
#> 709 F 1.9
#> 710 F 1.9
#> 711 F 1.9
#> 712 M 0.1
#> 713 M 0.1
#> 714 F 1.9
#> 715 F 1.9
#> 716 F 1.9
#> 717 F 1.9
#> 718 F 1.9
#> 719 M 0.1
#> 720 F 1.9
#> 721 F 1.9
#> 722 F 1.9
#> 723 M 0.1
#> 724 F 1.9
#> 725 F 1.9
#> 726 M 0.1
#> 727 F 1.9
#> 728 F 1.9
#> 729 F 1.9
#> 730 M 0.1
#> 731 F 1.9
#> 732 F 1.9
#> 733 F 1.9
#> 734 F 1.9
#> 735 F 1.9
#> 736 F 1.9
#> 737 F 1.9
#> 738 F 1.9
#> 739 M 0.1
#> 740 F 1.9
#> 741 F 1.9
#> 742 F 1.9
#> 743 F 1.9
#> 744 F 1.9
#> 745 F 1.9
#> 746 F 1.9
#> 747 F 1.9
#> 748 F 1.9
#> 749 F 1.9
#> 750 F 1.9
#> 751 F 1.9
#> 752 F 1.9
#> 753 F 1.9
#> 754 F 1.9
#> 755 F 1.9
#> 756 F 1.9
#> 757 F 1.9
#> 758 F 1.9
#> 759 F 1.9
#> 760 F 1.9
#> 761 F 1.9
#> 762 F 1.9
#> 763 F 1.9
#> 764 F 1.9
#> 765 F 1.9
#> 766 F 1.9
#> 767 F 1.9
#> 768 F 1.9
#> 769 F 1.9
#> 770 F 1.9
#> 771 F 1.9
#> 772 F 1.9
#> 773 F 1.9
#> 774 F 1.9
#> 775 F 1.9
#> 776 F 1.9
#> 777 F 1.9
#> 778 F 1.9
#> 779 F 1.9
#> 780 F 1.9
#> 781 F 1.9
#> 782 F 1.9
#> 783 F 1.9
#> 784 F 1.9
#> 785 F 1.9
#> 786 F 1.9
#> 787 F 1.9
#> 788 F 1.9
#> 789 F 1.9
#> 790 F 1.9
#> 791 F 1.9
#> 792 F 1.9
#> 793 F 1.9
#> 794 F 1.9
#> 795 F 1.9
#> 796 F 1.9
#> 797 F 1.9
#> 798 F 1.9
#> 799 F 1.9
#> 800 F 1.9
#> 801 F 1.9
#> 802 F 1.9
#> 803 F 1.9
#> 804 F 1.9
#> 805 F 1.9
#> 806 F 1.9
#> 807 F 1.9
#> 808 F 1.9
#> 809 F 1.9
#> 810 F 1.9
#> 811 F 1.9
#> 812 F 1.9
#> 813 F 1.9
#> 814 F 1.9
#> 815 F 1.9
#> 816 F 1.9
#> 817 F 1.9
#> 818 F 1.9
#> 819 F 1.9
#> 820 F 1.9
#> 821 F 1.9
#> 822 F 1.9
#> 823 F 1.9
#> 824 F 1.9
#> 825 F 1.9
#> 826 F 1.9
#> 827 F 1.9
#> 828 F 1.9
#> 829 F 1.9
#> 830 F 1.9
#> 831 F 1.9
#> 832 F 1.9
#> 833 F 1.9
#> 834 F 1.9
#> 835 F 1.9
#> 836 F 1.9
#> 837 F 1.9
#> 838 F 1.9
#> 839 F 1.9
#> 840 F 1.9
#> 841 F 1.9
#> 842 F 1.9
#> 843 F 1.9
#> 844 F 1.9
#> 845 F 1.9
#> 846 F 1.9
#> 847 F 1.9
#> 848 F 1.9
#> 849 F 1.9
#> 850 F 1.9
#> 851 F 1.9
#> 852 M 0.1
#> 853 F 1.9
#> 854 F 1.9
#> 855 F 1.9
#> 856 F 1.9
#> 857 F 1.9
#> 858 F 1.9
#> 859 F 1.9
#> 860 M 0.1
#> 861 F 1.9
#> 862 F 1.9
#> 863 F 1.9
#> 864 F 1.9
#> 865 F 1.9
#> 866 F 1.9
#> 867 F 1.9
#> 868 M 0.1
#> 869 F 1.9
#> 870 F 1.9
#> 871 F 1.9
#> 872 F 1.9
#> 873 F 1.9
#> 874 F 1.9
#> 875 F 1.9
#> 876 F 1.9
#> 877 F 1.9
#> 878 F 1.9
#> 879 F 1.9
#> 880 F 1.9
#> 881 F 1.9
#> 882 F 1.9
#> 883 F 1.9
#> 884 F 1.9
#> 885 F 1.9
#> 886 F 1.9
#> 887 F 1.9
#> 888 F 1.9
#> 889 F 1.9
#> 890 F 1.9
#> 891 F 1.9
#> 892 F 1.9
#> 893 F 1.9
#> 894 F 1.9
#> 895 F 1.9
#> 896 F 1.9
#> 897 F 1.9
#> 898 F 1.9
#> 899 M 0.1
#> 900 F 1.9
#> 901 F 1.9
#> 902 F 1.9
#> 903 F 1.9
#> 904 F 1.9
#> 905 F 1.9
#> 906 F 1.9
#> 907 F 1.9
#> 908 F 1.9
#> 909 F 1.9
#> 910 F 1.9
#> 911 F 1.9
#> 912 F 1.9
#> 913 F 1.9
#> 914 F 1.9
#> 915 F 1.9
#> 916 F 1.9
#> 917 F 1.9
#> 918 M 0.1
#> 919 M 0.1
#> 920 F 1.9
#> 921 F 1.9
#> 922 F 1.9
#> 923 F 1.9
#> 924 F 1.9
#> 925 F 1.9
#> 926 F 1.9
#> 927 F 1.9
#> 928 F 1.9
#> 929 F 1.9
#> 930 F 1.9
#> 931 F 1.9
#> 932 M 0.1
#> 933 F 1.9
#> 934 F 1.9
#> 935 F 1.9
#> 936 F 1.9
#> 937 F 1.9
#> 938 F 1.9
#> 939 F 1.9
#> 940 F 1.9
#> 941 F 1.9
#> 942 F 1.9
#> 943 F 1.9
#> 944 F 1.9
#> 945 F 1.9
#> 946 F 1.9
#> 947 F 1.9
#> 948 F 1.9
#> 949 F 1.9
#> 950 F 1.9
#> 951 F 1.9
#> 952 M 0.1
#> 953 F 1.9
#> 954 F 1.9
#> 955 F 1.9
#> 956 F 1.9
#> 957 F 1.9
#> 958 F 1.9
#> 959 F 1.9
#> 960 F 1.9
#> 961 F 1.9
#> 962 F 1.9
#> 963 F 1.9
#> 964 F 1.9
#> 965 F 1.9
#> 966 F 1.9
#> 967 F 1.9
#> 968 F 1.9
#> 969 F 1.9
#> 970 F 1.9
#> 971 M 0.1
#> 972 F 1.9
#> 973 F 1.9
#> 974 F 1.9
#> 975 F 1.9
#> 976 F 1.9
#> 977 F 1.9
#> 978 F 1.9
#> 979 M 0.1
#> 980 F 1.9
#> 981 F 1.9
#> 982 F 1.9
#> 983 F 1.9
#> 984 F 1.9
#> 985 M 0.1
#> 986 F 1.9
#> 987 F 1.9
#> 988 M 0.1
#> 989 F 1.9
#> 990 F 1.9
#> 991 F 1.9
#> 992 F 1.9
#> 993 M 0.1
#> 994 F 1.9
#> 995 F 1.9
#> 996 F 1.9
#> 997 F 1.9
#> 998 F 1.9
#> 999 M 0.1
#> 1000 F 1.9Calculate the LR for age evidence:
# LR for age evidence
# MP age = 25, tolerance range = 5 years
# POI observed age falls within range
lr_age(LR = TRUE, H = 1, MPa = 25, MPr = 5, epa = 0.05)
#> group Age UHRmin UHRmax LRa
#> 1 T1 30 26.90 33.10 7.60000000
#> 2 T1 22 19.46 24.54 7.60000000
#> 3 T1 25 22.25 27.75 7.60000000
#> 4 T1 21 18.53 23.47 7.60000000
#> 5 T1 30 26.90 33.10 7.60000000
#> 6 T1 21 18.53 23.47 7.60000000
#> 7 T1 23 20.39 25.61 7.60000000
#> 8 T1 21 18.53 23.47 7.60000000
#> 9 T1 27 24.11 29.89 7.60000000
#> 10 T1 28 25.04 30.96 7.60000000
#> 11 T1 26 23.18 28.82 7.60000000
#> 12 T1 20 17.60 22.40 7.60000000
#> 13 T1 25 22.25 27.75 7.60000000
#> 14 T1 28 25.04 30.96 7.60000000
#> 15 T1 23 20.39 25.61 7.60000000
#> 16 T1 22 19.46 24.54 7.60000000
#> 17 T1 26 23.18 28.82 7.60000000
#> 18 T1 24 21.32 26.68 7.60000000
#> 19 T1 27 24.11 29.89 7.60000000
#> 20 T1 25 22.25 27.75 7.60000000
#> 21 T1 23 20.39 25.61 7.60000000
#> 22 T1 24 21.32 26.68 7.60000000
#> 23 T1 26 23.18 28.82 7.60000000
#> 24 T1 24 21.32 26.68 7.60000000
#> 25 T1 26 23.18 28.82 7.60000000
#> 26 T1 21 18.53 23.47 7.60000000
#> 27 T1 20 17.60 22.40 7.60000000
#> 28 T1 27 24.11 29.89 7.60000000
#> 29 T1 24 21.32 26.68 7.60000000
#> 30 T1 29 25.97 32.03 7.60000000
#> 31 T1 20 17.60 22.40 7.60000000
#> 32 T1 25 22.25 27.75 7.60000000
#> 33 T1 24 21.32 26.68 7.60000000
#> 34 T1 25 22.25 27.75 7.60000000
#> 35 T1 20 17.60 22.40 7.60000000
#> 36 T1 28 25.04 30.96 7.60000000
#> 37 T1 27 24.11 29.89 7.60000000
#> 38 T1 22 19.46 24.54 7.60000000
#> 39 T0 75 68.75 81.25 0.05714286
#> 40 T1 28 25.04 30.96 7.60000000
#> 41 T1 29 25.97 32.03 7.60000000
#> 42 T1 30 26.90 33.10 7.60000000
#> 43 T1 23 20.39 25.61 7.60000000
#> 44 T1 28 25.04 30.96 7.60000000
#> 45 T1 30 26.90 33.10 7.60000000
#> 46 T1 23 20.39 25.61 7.60000000
#> 47 T1 20 17.60 22.40 7.60000000
#> 48 T1 20 17.60 22.40 7.60000000
#> 49 T1 21 18.53 23.47 7.60000000
#> 50 T1 28 25.04 30.96 7.60000000
#> 51 T1 27 24.11 29.89 7.60000000
#> 52 T1 30 26.90 33.10 7.60000000
#> 53 T1 21 18.53 23.47 7.60000000
#> 54 T1 22 19.46 24.54 7.60000000
#> 55 T1 28 25.04 30.96 7.60000000
#> 56 T1 20 17.60 22.40 7.60000000
#> 57 T1 27 24.11 29.89 7.60000000
#> 58 T1 26 23.18 28.82 7.60000000
#> 59 T1 30 26.90 33.10 7.60000000
#> 60 T1 30 26.90 33.10 7.60000000
#> 61 T1 20 17.60 22.40 7.60000000
#> 62 T1 29 25.97 32.03 7.60000000
#> 63 T1 27 24.11 29.89 7.60000000
#> 64 T1 22 19.46 24.54 7.60000000
#> 65 T1 30 26.90 33.10 7.60000000
#> 66 T1 26 23.18 28.82 7.60000000
#> 67 T1 28 25.04 30.96 7.60000000
#> 68 T1 21 18.53 23.47 7.60000000
#> 69 T1 28 25.04 30.96 7.60000000
#> 70 T1 20 17.60 22.40 7.60000000
#> 71 T1 21 18.53 23.47 7.60000000
#> 72 T1 27 24.11 29.89 7.60000000
#> 73 T1 28 25.04 30.96 7.60000000
#> 74 T1 27 24.11 29.89 7.60000000
#> 75 T1 30 26.90 33.10 7.60000000
#> 76 T1 23 20.39 25.61 7.60000000
#> 77 T1 24 21.32 26.68 7.60000000
#> 78 T1 30 26.90 33.10 7.60000000
#> 79 T1 20 17.60 22.40 7.60000000
#> 80 T1 22 19.46 24.54 7.60000000
#> 81 T1 21 18.53 23.47 7.60000000
#> 82 T1 22 19.46 24.54 7.60000000
#> 83 T1 20 17.60 22.40 7.60000000
#> 84 T1 28 25.04 30.96 7.60000000
#> 85 T1 25 22.25 27.75 7.60000000
#> 86 T1 26 23.18 28.82 7.60000000
#> 87 T1 22 19.46 24.54 7.60000000
#> 88 T1 28 25.04 30.96 7.60000000
#> 89 T1 20 17.60 22.40 7.60000000
#> 90 T1 24 21.32 26.68 7.60000000
#> 91 T1 22 19.46 24.54 7.60000000
#> 92 T1 25 22.25 27.75 7.60000000
#> 93 T1 26 23.18 28.82 7.60000000
#> 94 T1 21 18.53 23.47 7.60000000
#> 95 T1 29 25.97 32.03 7.60000000
#> 96 T1 28 25.04 30.96 7.60000000
#> 97 T1 27 24.11 29.89 7.60000000
#> 98 T1 25 22.25 27.75 7.60000000
#> 99 T1 27 24.11 29.89 7.60000000
#> 100 T1 20 17.60 22.40 7.60000000
#> 101 T1 30 26.90 33.10 7.60000000
#> 102 T1 26 23.18 28.82 7.60000000
#> 103 T1 24 21.32 26.68 7.60000000
#> 104 T1 23 20.39 25.61 7.60000000
#> 105 T1 20 17.60 22.40 7.60000000
#> 106 T1 25 22.25 27.75 7.60000000
#> 107 T1 22 19.46 24.54 7.60000000
#> 108 T1 21 18.53 23.47 7.60000000
#> 109 T1 22 19.46 24.54 7.60000000
#> 110 T1 30 26.90 33.10 7.60000000
#> 111 T1 24 21.32 26.68 7.60000000
#> 112 T1 21 18.53 23.47 7.60000000
#> 113 T0 59 53.87 64.13 0.05714286
#> 114 T1 25 22.25 27.75 7.60000000
#> 115 T1 27 24.11 29.89 7.60000000
#> 116 T1 21 18.53 23.47 7.60000000
#> 117 T1 25 22.25 27.75 7.60000000
#> 118 T1 26 23.18 28.82 7.60000000
#> 119 T1 30 26.90 33.10 7.60000000
#> 120 T1 22 19.46 24.54 7.60000000
#> 121 T1 24 21.32 26.68 7.60000000
#> 122 T1 27 24.11 29.89 7.60000000
#> 123 T0 8 6.44 9.56 0.05714286
#> 124 T1 21 18.53 23.47 7.60000000
#> 125 T1 22 19.46 24.54 7.60000000
#> 126 T1 27 24.11 29.89 7.60000000
#> 127 T1 26 23.18 28.82 7.60000000
#> 128 T1 25 22.25 27.75 7.60000000
#> 129 T1 26 23.18 28.82 7.60000000
#> 130 T1 25 22.25 27.75 7.60000000
#> 131 T0 6 4.58 7.42 0.05714286
#> 132 T1 30 26.90 33.10 7.60000000
#> 133 T1 24 21.32 26.68 7.60000000
#> 134 T1 30 26.90 33.10 7.60000000
#> 135 T1 24 21.32 26.68 7.60000000
#> 136 T1 27 24.11 29.89 7.60000000
#> 137 T1 25 22.25 27.75 7.60000000
#> 138 T1 25 22.25 27.75 7.60000000
#> 139 T1 29 25.97 32.03 7.60000000
#> 140 T1 24 21.32 26.68 7.60000000
#> 141 T1 25 22.25 27.75 7.60000000
#> 142 T0 51 46.43 55.57 0.05714286
#> 143 T1 27 24.11 29.89 7.60000000
#> 144 T1 22 19.46 24.54 7.60000000
#> 145 T1 23 20.39 25.61 7.60000000
#> 146 T1 29 25.97 32.03 7.60000000
#> 147 T1 20 17.60 22.40 7.60000000
#> 148 T1 29 25.97 32.03 7.60000000
#> 149 T0 67 61.31 72.69 0.05714286
#> 150 T1 26 23.18 28.82 7.60000000
#> 151 T1 29 25.97 32.03 7.60000000
#> 152 T1 27 24.11 29.89 7.60000000
#> 153 T1 25 22.25 27.75 7.60000000
#> 154 T1 24 21.32 26.68 7.60000000
#> 155 T1 23 20.39 25.61 7.60000000
#> 156 T0 51 46.43 55.57 0.05714286
#> 157 T1 30 26.90 33.10 7.60000000
#> 158 T0 38 34.34 41.66 0.05714286
#> 159 T1 23 20.39 25.61 7.60000000
#> 160 T1 30 26.90 33.10 7.60000000
#> 161 T1 20 17.60 22.40 7.60000000
#> 162 T1 29 25.97 32.03 7.60000000
#> 163 T1 24 21.32 26.68 7.60000000
#> 164 T1 30 26.90 33.10 7.60000000
#> 165 T1 28 25.04 30.96 7.60000000
#> 166 T1 23 20.39 25.61 7.60000000
#> 167 T1 23 20.39 25.61 7.60000000
#> 168 T1 25 22.25 27.75 7.60000000
#> 169 T1 24 21.32 26.68 7.60000000
#> 170 T1 28 25.04 30.96 7.60000000
#> 171 T1 26 23.18 28.82 7.60000000
#> 172 T1 26 23.18 28.82 7.60000000
#> 173 T1 29 25.97 32.03 7.60000000
#> 174 T1 29 25.97 32.03 7.60000000
#> 175 T1 23 20.39 25.61 7.60000000
#> 176 T1 21 18.53 23.47 7.60000000
#> 177 T1 20 17.60 22.40 7.60000000
#> 178 T1 20 17.60 22.40 7.60000000
#> 179 T1 24 21.32 26.68 7.60000000
#> 180 T1 24 21.32 26.68 7.60000000
#> 181 T1 24 21.32 26.68 7.60000000
#> 182 T1 21 18.53 23.47 7.60000000
#> 183 T1 30 26.90 33.10 7.60000000
#> 184 T1 21 18.53 23.47 7.60000000
#> 185 T0 8 6.44 9.56 0.05714286
#> 186 T1 29 25.97 32.03 7.60000000
#> 187 T1 22 19.46 24.54 7.60000000
#> 188 T1 23 20.39 25.61 7.60000000
#> 189 T1 28 25.04 30.96 7.60000000
#> 190 T1 29 25.97 32.03 7.60000000
#> 191 T1 26 23.18 28.82 7.60000000
#> 192 T0 52 47.36 56.64 0.05714286
#> 193 T1 27 24.11 29.89 7.60000000
#> 194 T1 22 19.46 24.54 7.60000000
#> 195 T1 21 18.53 23.47 7.60000000
#> 196 T1 28 25.04 30.96 7.60000000
#> 197 T1 25 22.25 27.75 7.60000000
#> 198 T1 21 18.53 23.47 7.60000000
#> 199 T0 77 70.61 83.39 0.05714286
#> 200 T1 26 23.18 28.82 7.60000000
#> 201 T1 28 25.04 30.96 7.60000000
#> 202 T1 27 24.11 29.89 7.60000000
#> 203 T1 25 22.25 27.75 7.60000000
#> 204 T1 22 19.46 24.54 7.60000000
#> 205 T1 22 19.46 24.54 7.60000000
#> 206 T1 21 18.53 23.47 7.60000000
#> 207 T1 26 23.18 28.82 7.60000000
#> 208 T1 23 20.39 25.61 7.60000000
#> 209 T1 21 18.53 23.47 7.60000000
#> 210 T0 59 53.87 64.13 0.05714286
#> 211 T1 22 19.46 24.54 7.60000000
#> 212 T1 25 22.25 27.75 7.60000000
#> 213 T1 21 18.53 23.47 7.60000000
#> 214 T1 20 17.60 22.40 7.60000000
#> 215 T1 27 24.11 29.89 7.60000000
#> 216 T1 26 23.18 28.82 7.60000000
#> 217 T1 27 24.11 29.89 7.60000000
#> 218 T1 30 26.90 33.10 7.60000000
#> 219 T1 27 24.11 29.89 7.60000000
#> 220 T1 20 17.60 22.40 7.60000000
#> 221 T1 27 24.11 29.89 7.60000000
#> 222 T1 27 24.11 29.89 7.60000000
#> 223 T1 24 21.32 26.68 7.60000000
#> 224 T1 20 17.60 22.40 7.60000000
#> 225 T1 27 24.11 29.89 7.60000000
#> 226 T1 27 24.11 29.89 7.60000000
#> 227 T1 21 18.53 23.47 7.60000000
#> 228 T1 30 26.90 33.10 7.60000000
#> 229 T1 24 21.32 26.68 7.60000000
#> 230 T1 30 26.90 33.10 7.60000000
#> 231 T1 21 18.53 23.47 7.60000000
#> 232 T1 21 18.53 23.47 7.60000000
#> 233 T1 29 25.97 32.03 7.60000000
#> 234 T1 30 26.90 33.10 7.60000000
#> 235 T1 21 18.53 23.47 7.60000000
#> 236 T1 26 23.18 28.82 7.60000000
#> 237 T1 20 17.60 22.40 7.60000000
#> 238 T1 24 21.32 26.68 7.60000000
#> 239 T1 27 24.11 29.89 7.60000000
#> 240 T0 45 40.85 49.15 0.05714286
#> 241 T1 23 20.39 25.61 7.60000000
#> 242 T1 27 24.11 29.89 7.60000000
#> 243 T1 27 24.11 29.89 7.60000000
#> 244 T1 20 17.60 22.40 7.60000000
#> 245 T0 55 50.15 59.85 0.05714286
#> 246 T1 20 17.60 22.40 7.60000000
#> 247 T1 22 19.46 24.54 7.60000000
#> 248 T1 29 25.97 32.03 7.60000000
#> 249 T0 37 33.41 40.59 0.05714286
#> 250 T1 21 18.53 23.47 7.60000000
#> 251 T1 21 18.53 23.47 7.60000000
#> 252 T1 23 20.39 25.61 7.60000000
#> 253 T1 22 19.46 24.54 7.60000000
#> 254 T1 26 23.18 28.82 7.60000000
#> 255 T1 24 21.32 26.68 7.60000000
#> 256 T1 25 22.25 27.75 7.60000000
#> 257 T1 30 26.90 33.10 7.60000000
#> 258 T1 21 18.53 23.47 7.60000000
#> 259 T1 21 18.53 23.47 7.60000000
#> 260 T1 26 23.18 28.82 7.60000000
#> 261 T1 28 25.04 30.96 7.60000000
#> 262 T1 20 17.60 22.40 7.60000000
#> 263 T1 25 22.25 27.75 7.60000000
#> 264 T1 20 17.60 22.40 7.60000000
#> 265 T1 29 25.97 32.03 7.60000000
#> 266 T1 25 22.25 27.75 7.60000000
#> 267 T1 23 20.39 25.61 7.60000000
#> 268 T1 26 23.18 28.82 7.60000000
#> 269 T1 27 24.11 29.89 7.60000000
#> 270 T1 27 24.11 29.89 7.60000000
#> 271 T1 24 21.32 26.68 7.60000000
#> 272 T1 28 25.04 30.96 7.60000000
#> 273 T1 28 25.04 30.96 7.60000000
#> 274 T1 24 21.32 26.68 7.60000000
#> 275 T1 26 23.18 28.82 7.60000000
#> 276 T1 28 25.04 30.96 7.60000000
#> 277 T1 23 20.39 25.61 7.60000000
#> 278 T1 28 25.04 30.96 7.60000000
#> 279 T1 30 26.90 33.10 7.60000000
#> 280 T1 22 19.46 24.54 7.60000000
#> 281 T1 24 21.32 26.68 7.60000000
#> 282 T1 23 20.39 25.61 7.60000000
#> 283 T1 20 17.60 22.40 7.60000000
#> 284 T1 30 26.90 33.10 7.60000000
#> 285 T1 22 19.46 24.54 7.60000000
#> 286 T1 23 20.39 25.61 7.60000000
#> 287 T1 29 25.97 32.03 7.60000000
#> 288 T1 28 25.04 30.96 7.60000000
#> 289 T1 28 25.04 30.96 7.60000000
#> 290 T1 24 21.32 26.68 7.60000000
#> 291 T1 22 19.46 24.54 7.60000000
#> 292 T1 24 21.32 26.68 7.60000000
#> 293 T0 65 59.45 70.55 0.05714286
#> 294 T1 23 20.39 25.61 7.60000000
#> 295 T1 24 21.32 26.68 7.60000000
#> 296 T1 24 21.32 26.68 7.60000000
#> 297 T1 22 19.46 24.54 7.60000000
#> 298 T1 22 19.46 24.54 7.60000000
#> 299 T1 23 20.39 25.61 7.60000000
#> 300 T1 29 25.97 32.03 7.60000000
#> 301 T1 21 18.53 23.47 7.60000000
#> 302 T0 3 1.79 4.21 0.05714286
#> 303 T1 30 26.90 33.10 7.60000000
#> 304 T1 24 21.32 26.68 7.60000000
#> 305 T1 27 24.11 29.89 7.60000000
#> 306 T1 25 22.25 27.75 7.60000000
#> 307 T1 22 19.46 24.54 7.60000000
#> 308 T0 72 65.96 78.04 0.05714286
#> 309 T1 22 19.46 24.54 7.60000000
#> 310 T1 28 25.04 30.96 7.60000000
#> 311 T1 20 17.60 22.40 7.60000000
#> 312 T1 23 20.39 25.61 7.60000000
#> 313 T1 25 22.25 27.75 7.60000000
#> 314 T1 20 17.60 22.40 7.60000000
#> 315 T1 30 26.90 33.10 7.60000000
#> 316 T1 28 25.04 30.96 7.60000000
#> 317 T1 24 21.32 26.68 7.60000000
#> 318 T1 22 19.46 24.54 7.60000000
#> 319 T1 23 20.39 25.61 7.60000000
#> 320 T1 25 22.25 27.75 7.60000000
#> 321 T1 24 21.32 26.68 7.60000000
#> 322 T1 21 18.53 23.47 7.60000000
#> 323 T1 30 26.90 33.10 7.60000000
#> 324 T1 21 18.53 23.47 7.60000000
#> 325 T1 22 19.46 24.54 7.60000000
#> 326 T1 27 24.11 29.89 7.60000000
#> 327 T1 24 21.32 26.68 7.60000000
#> 328 T1 29 25.97 32.03 7.60000000
#> 329 T1 22 19.46 24.54 7.60000000
#> 330 T1 26 23.18 28.82 7.60000000
#> 331 T1 28 25.04 30.96 7.60000000
#> 332 T1 22 19.46 24.54 7.60000000
#> 333 T0 7 5.51 8.49 0.05714286
#> 334 T1 25 22.25 27.75 7.60000000
#> 335 T1 22 19.46 24.54 7.60000000
#> 336 T1 21 18.53 23.47 7.60000000
#> 337 T1 24 21.32 26.68 7.60000000
#> 338 T1 24 21.32 26.68 7.60000000
#> 339 T1 22 19.46 24.54 7.60000000
#> 340 T1 22 19.46 24.54 7.60000000
#> 341 T0 59 53.87 64.13 0.05714286
#> 342 T1 25 22.25 27.75 7.60000000
#> 343 T1 25 22.25 27.75 7.60000000
#> 344 T0 77 70.61 83.39 0.05714286
#> 345 T1 22 19.46 24.54 7.60000000
#> 346 T1 29 25.97 32.03 7.60000000
#> 347 T1 26 23.18 28.82 7.60000000
#> 348 T1 26 23.18 28.82 7.60000000
#> 349 T1 28 25.04 30.96 7.60000000
#> 350 T1 24 21.32 26.68 7.60000000
#> 351 T1 26 23.18 28.82 7.60000000
#> 352 T1 24 21.32 26.68 7.60000000
#> 353 T1 23 20.39 25.61 7.60000000
#> 354 T1 27 24.11 29.89 7.60000000
#> 355 T0 46 41.78 50.22 0.05714286
#> 356 T1 27 24.11 29.89 7.60000000
#> 357 T1 27 24.11 29.89 7.60000000
#> 358 T1 22 19.46 24.54 7.60000000
#> 359 T1 25 22.25 27.75 7.60000000
#> 360 T1 23 20.39 25.61 7.60000000
#> 361 T1 29 25.97 32.03 7.60000000
#> 362 T1 26 23.18 28.82 7.60000000
#> 363 T1 30 26.90 33.10 7.60000000
#> 364 T1 25 22.25 27.75 7.60000000
#> 365 T1 30 26.90 33.10 7.60000000
#> 366 T0 10 8.30 11.70 0.05714286
#> 367 T1 27 24.11 29.89 7.60000000
#> 368 T1 25 22.25 27.75 7.60000000
#> 369 T1 25 22.25 27.75 7.60000000
#> 370 T1 27 24.11 29.89 7.60000000
#> 371 T0 44 39.92 48.08 0.05714286
#> 372 T1 28 25.04 30.96 7.60000000
#> 373 T1 23 20.39 25.61 7.60000000
#> 374 T1 23 20.39 25.61 7.60000000
#> 375 T1 28 25.04 30.96 7.60000000
#> 376 T1 21 18.53 23.47 7.60000000
#> 377 T1 30 26.90 33.10 7.60000000
#> 378 T1 28 25.04 30.96 7.60000000
#> 379 T1 27 24.11 29.89 7.60000000
#> 380 T1 27 24.11 29.89 7.60000000
#> 381 T1 30 26.90 33.10 7.60000000
#> 382 T1 26 23.18 28.82 7.60000000
#> 383 T1 28 25.04 30.96 7.60000000
#> 384 T1 26 23.18 28.82 7.60000000
#> 385 T1 20 17.60 22.40 7.60000000
#> 386 T1 21 18.53 23.47 7.60000000
#> 387 T1 22 19.46 24.54 7.60000000
#> 388 T1 20 17.60 22.40 7.60000000
#> 389 T0 68 62.24 73.76 0.05714286
#> 390 T1 21 18.53 23.47 7.60000000
#> 391 T1 27 24.11 29.89 7.60000000
#> 392 T1 27 24.11 29.89 7.60000000
#> 393 T1 26 23.18 28.82 7.60000000
#> 394 T1 26 23.18 28.82 7.60000000
#> 395 T1 30 26.90 33.10 7.60000000
#> 396 T1 24 21.32 26.68 7.60000000
#> 397 T1 27 24.11 29.89 7.60000000
#> 398 T1 29 25.97 32.03 7.60000000
#> 399 T1 27 24.11 29.89 7.60000000
#> 400 T1 23 20.39 25.61 7.60000000
#> 401 T1 22 19.46 24.54 7.60000000
#> 402 T1 20 17.60 22.40 7.60000000
#> 403 T1 24 21.32 26.68 7.60000000
#> 404 T1 20 17.60 22.40 7.60000000
#> 405 T1 27 24.11 29.89 7.60000000
#> 406 T1 25 22.25 27.75 7.60000000
#> 407 T1 23 20.39 25.61 7.60000000
#> 408 T1 28 25.04 30.96 7.60000000
#> 409 T1 23 20.39 25.61 7.60000000
#> 410 T1 25 22.25 27.75 7.60000000
#> 411 T1 28 25.04 30.96 7.60000000
#> 412 T1 21 18.53 23.47 7.60000000
#> 413 T0 33 29.69 36.31 0.05714286
#> 414 T1 26 23.18 28.82 7.60000000
#> 415 T1 27 24.11 29.89 7.60000000
#> 416 T1 26 23.18 28.82 7.60000000
#> 417 T1 26 23.18 28.82 7.60000000
#> 418 T1 29 25.97 32.03 7.60000000
#> 419 T1 30 26.90 33.10 7.60000000
#> 420 T1 25 22.25 27.75 7.60000000
#> 421 T1 27 24.11 29.89 7.60000000
#> 422 T1 30 26.90 33.10 7.60000000
#> 423 T1 30 26.90 33.10 7.60000000
#> 424 T1 25 22.25 27.75 7.60000000
#> 425 T1 25 22.25 27.75 7.60000000
#> 426 T1 21 18.53 23.47 7.60000000
#> 427 T1 30 26.90 33.10 7.60000000
#> 428 T1 30 26.90 33.10 7.60000000
#> 429 T1 20 17.60 22.40 7.60000000
#> 430 T1 24 21.32 26.68 7.60000000
#> 431 T1 22 19.46 24.54 7.60000000
#> 432 T1 21 18.53 23.47 7.60000000
#> 433 T1 26 23.18 28.82 7.60000000
#> 434 T1 30 26.90 33.10 7.60000000
#> 435 T1 25 22.25 27.75 7.60000000
#> 436 T1 24 21.32 26.68 7.60000000
#> 437 T1 23 20.39 25.61 7.60000000
#> 438 T1 23 20.39 25.61 7.60000000
#> 439 T1 22 19.46 24.54 7.60000000
#> 440 T1 20 17.60 22.40 7.60000000
#> 441 T1 22 19.46 24.54 7.60000000
#> 442 T1 26 23.18 28.82 7.60000000
#> 443 T1 29 25.97 32.03 7.60000000
#> 444 T1 23 20.39 25.61 7.60000000
#> 445 T0 63 57.59 68.41 0.05714286
#> 446 T1 23 20.39 25.61 7.60000000
#> 447 T1 23 20.39 25.61 7.60000000
#> 448 T1 23 20.39 25.61 7.60000000
#> 449 T1 30 26.90 33.10 7.60000000
#> 450 T1 30 26.90 33.10 7.60000000
#> 451 T1 30 26.90 33.10 7.60000000
#> 452 T1 30 26.90 33.10 7.60000000
#> 453 T0 17 14.81 19.19 0.05714286
#> 454 T0 67 61.31 72.69 0.05714286
#> 455 T1 28 25.04 30.96 7.60000000
#> 456 T1 28 25.04 30.96 7.60000000
#> 457 T1 23 20.39 25.61 7.60000000
#> 458 T1 23 20.39 25.61 7.60000000
#> 459 T1 24 21.32 26.68 7.60000000
#> 460 T1 30 26.90 33.10 7.60000000
#> 461 T1 30 26.90 33.10 7.60000000
#> 462 T1 20 17.60 22.40 7.60000000
#> 463 T1 24 21.32 26.68 7.60000000
#> 464 T1 25 22.25 27.75 7.60000000
#> 465 T1 24 21.32 26.68 7.60000000
#> 466 T1 30 26.90 33.10 7.60000000
#> 467 T1 28 25.04 30.96 7.60000000
#> 468 T1 20 17.60 22.40 7.60000000
#> 469 T1 29 25.97 32.03 7.60000000
#> 470 T1 28 25.04 30.96 7.60000000
#> 471 T1 27 24.11 29.89 7.60000000
#> 472 T1 26 23.18 28.82 7.60000000
#> 473 T1 20 17.60 22.40 7.60000000
#> 474 T0 57 52.01 61.99 0.05714286
#> 475 T1 26 23.18 28.82 7.60000000
#> 476 T0 42 38.06 45.94 0.05714286
#> 477 T1 29 25.97 32.03 7.60000000
#> 478 T1 20 17.60 22.40 7.60000000
#> 479 T1 28 25.04 30.96 7.60000000
#> 480 T1 26 23.18 28.82 7.60000000
#> 481 T1 22 19.46 24.54 7.60000000
#> 482 T1 30 26.90 33.10 7.60000000
#> 483 T1 29 25.97 32.03 7.60000000
#> 484 T1 30 26.90 33.10 7.60000000
#> 485 T1 25 22.25 27.75 7.60000000
#> 486 T1 22 19.46 24.54 7.60000000
#> 487 T1 25 22.25 27.75 7.60000000
#> 488 T1 25 22.25 27.75 7.60000000
#> 489 T1 30 26.90 33.10 7.60000000
#> 490 T1 20 17.60 22.40 7.60000000
#> 491 T1 28 25.04 30.96 7.60000000
#> 492 T1 29 25.97 32.03 7.60000000
#> 493 T1 20 17.60 22.40 7.60000000
#> 494 T1 20 17.60 22.40 7.60000000
#> 495 T1 22 19.46 24.54 7.60000000
#> 496 T1 26 23.18 28.82 7.60000000
#> 497 T1 29 25.97 32.03 7.60000000
#> 498 T1 21 18.53 23.47 7.60000000
#> 499 T1 22 19.46 24.54 7.60000000
#> 500 T1 27 24.11 29.89 7.60000000
#> 501 T1 25 22.25 27.75 7.60000000
#> 502 T1 30 26.90 33.10 7.60000000
#> 503 T1 29 25.97 32.03 7.60000000
#> 504 T1 30 26.90 33.10 7.60000000
#> 505 T1 28 25.04 30.96 7.60000000
#> 506 T1 21 18.53 23.47 7.60000000
#> 507 T1 27 24.11 29.89 7.60000000
#> 508 T1 20 17.60 22.40 7.60000000
#> 509 T1 26 23.18 28.82 7.60000000
#> 510 T1 27 24.11 29.89 7.60000000
#> 511 T1 24 21.32 26.68 7.60000000
#> 512 T1 23 20.39 25.61 7.60000000
#> 513 T1 23 20.39 25.61 7.60000000
#> 514 T1 22 19.46 24.54 7.60000000
#> 515 T1 29 25.97 32.03 7.60000000
#> 516 T1 23 20.39 25.61 7.60000000
#> 517 T1 24 21.32 26.68 7.60000000
#> 518 T1 25 22.25 27.75 7.60000000
#> 519 T1 29 25.97 32.03 7.60000000
#> 520 T1 27 24.11 29.89 7.60000000
#> 521 T1 27 24.11 29.89 7.60000000
#> 522 T1 21 18.53 23.47 7.60000000
#> 523 T1 24 21.32 26.68 7.60000000
#> 524 T1 22 19.46 24.54 7.60000000
#> 525 T1 30 26.90 33.10 7.60000000
#> 526 T1 21 18.53 23.47 7.60000000
#> 527 T1 23 20.39 25.61 7.60000000
#> 528 T1 27 24.11 29.89 7.60000000
#> 529 T1 23 20.39 25.61 7.60000000
#> 530 T1 24 21.32 26.68 7.60000000
#> 531 T1 30 26.90 33.10 7.60000000
#> 532 T1 20 17.60 22.40 7.60000000
#> 533 T1 20 17.60 22.40 7.60000000
#> 534 T1 30 26.90 33.10 7.60000000
#> 535 T1 23 20.39 25.61 7.60000000
#> 536 T0 9 7.37 10.63 0.05714286
#> 537 T1 28 25.04 30.96 7.60000000
#> 538 T1 22 19.46 24.54 7.60000000
#> 539 T1 21 18.53 23.47 7.60000000
#> 540 T1 25 22.25 27.75 7.60000000
#> 541 T1 27 24.11 29.89 7.60000000
#> 542 T1 27 24.11 29.89 7.60000000
#> 543 T0 10 8.30 11.70 0.05714286
#> 544 T1 30 26.90 33.10 7.60000000
#> 545 T1 24 21.32 26.68 7.60000000
#> 546 T1 29 25.97 32.03 7.60000000
#> 547 T1 26 23.18 28.82 7.60000000
#> 548 T1 26 23.18 28.82 7.60000000
#> 549 T0 36 32.48 39.52 0.05714286
#> 550 T1 24 21.32 26.68 7.60000000
#> 551 T1 28 25.04 30.96 7.60000000
#> 552 T1 23 20.39 25.61 7.60000000
#> 553 T1 27 24.11 29.89 7.60000000
#> 554 T1 21 18.53 23.47 7.60000000
#> 555 T1 20 17.60 22.40 7.60000000
#> 556 T0 6 4.58 7.42 0.05714286
#> 557 T1 21 18.53 23.47 7.60000000
#> 558 T1 22 19.46 24.54 7.60000000
#> 559 T1 24 21.32 26.68 7.60000000
#> 560 T1 26 23.18 28.82 7.60000000
#> 561 T1 25 22.25 27.75 7.60000000
#> 562 T0 80 73.40 86.60 0.05714286
#> 563 T1 26 23.18 28.82 7.60000000
#> 564 T1 26 23.18 28.82 7.60000000
#> 565 T1 26 23.18 28.82 7.60000000
#> 566 T1 26 23.18 28.82 7.60000000
#> 567 T1 24 21.32 26.68 7.60000000
#> 568 T1 29 25.97 32.03 7.60000000
#> 569 T1 24 21.32 26.68 7.60000000
#> 570 T1 23 20.39 25.61 7.60000000
#> 571 T1 30 26.90 33.10 7.60000000
#> 572 T1 21 18.53 23.47 7.60000000
#> 573 T1 26 23.18 28.82 7.60000000
#> 574 T1 25 22.25 27.75 7.60000000
#> 575 T1 29 25.97 32.03 7.60000000
#> 576 T1 28 25.04 30.96 7.60000000
#> 577 T1 27 24.11 29.89 7.60000000
#> 578 T1 24 21.32 26.68 7.60000000
#> 579 T1 23 20.39 25.61 7.60000000
#> 580 T1 25 22.25 27.75 7.60000000
#> 581 T1 22 19.46 24.54 7.60000000
#> 582 T1 21 18.53 23.47 7.60000000
#> 583 T1 25 22.25 27.75 7.60000000
#> 584 T1 29 25.97 32.03 7.60000000
#> 585 T1 20 17.60 22.40 7.60000000
#> 586 T1 30 26.90 33.10 7.60000000
#> 587 T1 29 25.97 32.03 7.60000000
#> 588 T0 19 16.67 21.33 0.05714286
#> 589 T0 40 36.20 43.80 0.05714286
#> 590 T1 25 22.25 27.75 7.60000000
#> 591 T0 11 9.23 12.77 0.05714286
#> 592 T1 21 18.53 23.47 7.60000000
#> 593 T1 23 20.39 25.61 7.60000000
#> 594 T1 30 26.90 33.10 7.60000000
#> 595 T1 21 18.53 23.47 7.60000000
#> 596 T1 30 26.90 33.10 7.60000000
#> 597 T1 29 25.97 32.03 7.60000000
#> 598 T1 20 17.60 22.40 7.60000000
#> 599 T1 30 26.90 33.10 7.60000000
#> 600 T1 29 25.97 32.03 7.60000000
#> 601 T1 25 22.25 27.75 7.60000000
#> 602 T1 27 24.11 29.89 7.60000000
#> 603 T1 29 25.97 32.03 7.60000000
#> 604 T1 25 22.25 27.75 7.60000000
#> 605 T1 23 20.39 25.61 7.60000000
#> 606 T1 28 25.04 30.96 7.60000000
#> 607 T1 26 23.18 28.82 7.60000000
#> 608 T1 29 25.97 32.03 7.60000000
#> 609 T1 29 25.97 32.03 7.60000000
#> 610 T1 23 20.39 25.61 7.60000000
#> 611 T1 24 21.32 26.68 7.60000000
#> 612 T1 21 18.53 23.47 7.60000000
#> 613 T1 23 20.39 25.61 7.60000000
#> 614 T1 30 26.90 33.10 7.60000000
#> 615 T1 24 21.32 26.68 7.60000000
#> 616 T1 25 22.25 27.75 7.60000000
#> 617 T1 25 22.25 27.75 7.60000000
#> 618 T1 29 25.97 32.03 7.60000000
#> 619 T1 25 22.25 27.75 7.60000000
#> 620 T1 24 21.32 26.68 7.60000000
#> 621 T1 20 17.60 22.40 7.60000000
#> 622 T1 30 26.90 33.10 7.60000000
#> 623 T1 30 26.90 33.10 7.60000000
#> 624 T1 24 21.32 26.68 7.60000000
#> 625 T1 29 25.97 32.03 7.60000000
#> 626 T1 23 20.39 25.61 7.60000000
#> 627 T1 24 21.32 26.68 7.60000000
#> 628 T1 26 23.18 28.82 7.60000000
#> 629 T1 22 19.46 24.54 7.60000000
#> 630 T1 24 21.32 26.68 7.60000000
#> 631 T1 20 17.60 22.40 7.60000000
#> 632 T1 29 25.97 32.03 7.60000000
#> 633 T1 24 21.32 26.68 7.60000000
#> 634 T1 26 23.18 28.82 7.60000000
#> 635 T1 28 25.04 30.96 7.60000000
#> 636 T1 23 20.39 25.61 7.60000000
#> 637 T1 26 23.18 28.82 7.60000000
#> 638 T1 29 25.97 32.03 7.60000000
#> 639 T1 21 18.53 23.47 7.60000000
#> 640 T1 30 26.90 33.10 7.60000000
#> 641 T1 27 24.11 29.89 7.60000000
#> 642 T1 22 19.46 24.54 7.60000000
#> 643 T1 29 25.97 32.03 7.60000000
#> 644 T1 24 21.32 26.68 7.60000000
#> 645 T1 21 18.53 23.47 7.60000000
#> 646 T1 30 26.90 33.10 7.60000000
#> 647 T1 21 18.53 23.47 7.60000000
#> 648 T1 22 19.46 24.54 7.60000000
#> 649 T1 21 18.53 23.47 7.60000000
#> 650 T1 25 22.25 27.75 7.60000000
#> 651 T1 28 25.04 30.96 7.60000000
#> 652 T1 28 25.04 30.96 7.60000000
#> 653 T1 21 18.53 23.47 7.60000000
#> 654 T1 28 25.04 30.96 7.60000000
#> 655 T1 26 23.18 28.82 7.60000000
#> 656 T1 23 20.39 25.61 7.60000000
#> 657 T1 20 17.60 22.40 7.60000000
#> 658 T1 28 25.04 30.96 7.60000000
#> 659 T1 29 25.97 32.03 7.60000000
#> 660 T1 28 25.04 30.96 7.60000000
#> 661 T1 25 22.25 27.75 7.60000000
#> 662 T1 25 22.25 27.75 7.60000000
#> 663 T1 24 21.32 26.68 7.60000000
#> 664 T1 23 20.39 25.61 7.60000000
#> 665 T1 29 25.97 32.03 7.60000000
#> 666 T1 23 20.39 25.61 7.60000000
#> 667 T1 30 26.90 33.10 7.60000000
#> 668 T1 22 19.46 24.54 7.60000000
#> 669 T1 30 26.90 33.10 7.60000000
#> 670 T1 24 21.32 26.68 7.60000000
#> 671 T1 24 21.32 26.68 7.60000000
#> 672 T1 24 21.32 26.68 7.60000000
#> 673 T1 29 25.97 32.03 7.60000000
#> 674 T1 23 20.39 25.61 7.60000000
#> 675 T1 30 26.90 33.10 7.60000000
#> 676 T1 23 20.39 25.61 7.60000000
#> 677 T1 20 17.60 22.40 7.60000000
#> 678 T1 22 19.46 24.54 7.60000000
#> 679 T1 28 25.04 30.96 7.60000000
#> 680 T0 61 55.73 66.27 0.05714286
#> 681 T1 24 21.32 26.68 7.60000000
#> 682 T1 24 21.32 26.68 7.60000000
#> 683 T1 27 24.11 29.89 7.60000000
#> 684 T1 22 19.46 24.54 7.60000000
#> 685 T1 26 23.18 28.82 7.60000000
#> 686 T1 24 21.32 26.68 7.60000000
#> 687 T1 27 24.11 29.89 7.60000000
#> 688 T1 20 17.60 22.40 7.60000000
#> 689 T1 21 18.53 23.47 7.60000000
#> 690 T1 30 26.90 33.10 7.60000000
#> 691 T1 23 20.39 25.61 7.60000000
#> 692 T1 26 23.18 28.82 7.60000000
#> 693 T1 20 17.60 22.40 7.60000000
#> 694 T1 20 17.60 22.40 7.60000000
#> 695 T1 20 17.60 22.40 7.60000000
#> 696 T1 24 21.32 26.68 7.60000000
#> 697 T1 24 21.32 26.68 7.60000000
#> 698 T1 25 22.25 27.75 7.60000000
#> 699 T1 24 21.32 26.68 7.60000000
#> 700 T1 25 22.25 27.75 7.60000000
#> 701 T1 27 24.11 29.89 7.60000000
#> 702 T1 22 19.46 24.54 7.60000000
#> 703 T1 21 18.53 23.47 7.60000000
#> 704 T1 21 18.53 23.47 7.60000000
#> 705 T1 24 21.32 26.68 7.60000000
#> 706 T1 21 18.53 23.47 7.60000000
#> 707 T1 22 19.46 24.54 7.60000000
#> 708 T1 30 26.90 33.10 7.60000000
#> 709 T1 25 22.25 27.75 7.60000000
#> 710 T1 24 21.32 26.68 7.60000000
#> 711 T1 23 20.39 25.61 7.60000000
#> 712 T0 51 46.43 55.57 0.05714286
#> 713 T0 39 35.27 42.73 0.05714286
#> 714 T1 21 18.53 23.47 7.60000000
#> 715 T1 20 17.60 22.40 7.60000000
#> 716 T1 26 23.18 28.82 7.60000000
#> 717 T1 22 19.46 24.54 7.60000000
#> 718 T1 28 25.04 30.96 7.60000000
#> 719 T0 75 68.75 81.25 0.05714286
#> 720 T1 26 23.18 28.82 7.60000000
#> 721 T1 24 21.32 26.68 7.60000000
#> 722 T1 24 21.32 26.68 7.60000000
#> 723 T0 51 46.43 55.57 0.05714286
#> 724 T1 26 23.18 28.82 7.60000000
#> 725 T1 26 23.18 28.82 7.60000000
#> 726 T0 14 12.02 15.98 0.05714286
#> 727 T1 28 25.04 30.96 7.60000000
#> 728 T1 28 25.04 30.96 7.60000000
#> 729 T1 23 20.39 25.61 7.60000000
#> 730 T0 50 45.50 54.50 0.05714286
#> 731 T1 27 24.11 29.89 7.60000000
#> 732 T1 29 25.97 32.03 7.60000000
#> 733 T1 24 21.32 26.68 7.60000000
#> 734 T1 28 25.04 30.96 7.60000000
#> 735 T1 23 20.39 25.61 7.60000000
#> 736 T1 26 23.18 28.82 7.60000000
#> 737 T1 30 26.90 33.10 7.60000000
#> 738 T1 20 17.60 22.40 7.60000000
#> 739 T0 73 66.89 79.11 0.05714286
#> 740 T1 21 18.53 23.47 7.60000000
#> 741 T1 21 18.53 23.47 7.60000000
#> 742 T1 27 24.11 29.89 7.60000000
#> 743 T1 23 20.39 25.61 7.60000000
#> 744 T1 21 18.53 23.47 7.60000000
#> 745 T1 21 18.53 23.47 7.60000000
#> 746 T1 28 25.04 30.96 7.60000000
#> 747 T1 26 23.18 28.82 7.60000000
#> 748 T1 30 26.90 33.10 7.60000000
#> 749 T1 26 23.18 28.82 7.60000000
#> 750 T1 28 25.04 30.96 7.60000000
#> 751 T1 20 17.60 22.40 7.60000000
#> 752 T1 24 21.32 26.68 7.60000000
#> 753 T1 28 25.04 30.96 7.60000000
#> 754 T1 30 26.90 33.10 7.60000000
#> 755 T1 30 26.90 33.10 7.60000000
#> 756 T1 23 20.39 25.61 7.60000000
#> 757 T1 29 25.97 32.03 7.60000000
#> 758 T1 21 18.53 23.47 7.60000000
#> 759 T1 25 22.25 27.75 7.60000000
#> 760 T1 23 20.39 25.61 7.60000000
#> 761 T1 27 24.11 29.89 7.60000000
#> 762 T1 27 24.11 29.89 7.60000000
#> 763 T1 25 22.25 27.75 7.60000000
#> 764 T1 28 25.04 30.96 7.60000000
#> 765 T1 23 20.39 25.61 7.60000000
#> 766 T1 24 21.32 26.68 7.60000000
#> 767 T1 22 19.46 24.54 7.60000000
#> 768 T1 29 25.97 32.03 7.60000000
#> 769 T1 24 21.32 26.68 7.60000000
#> 770 T1 23 20.39 25.61 7.60000000
#> 771 T1 25 22.25 27.75 7.60000000
#> 772 T1 27 24.11 29.89 7.60000000
#> 773 T1 23 20.39 25.61 7.60000000
#> 774 T1 24 21.32 26.68 7.60000000
#> 775 T1 27 24.11 29.89 7.60000000
#> 776 T1 20 17.60 22.40 7.60000000
#> 777 T1 23 20.39 25.61 7.60000000
#> 778 T1 22 19.46 24.54 7.60000000
#> 779 T1 26 23.18 28.82 7.60000000
#> 780 T1 24 21.32 26.68 7.60000000
#> 781 T1 23 20.39 25.61 7.60000000
#> 782 T1 26 23.18 28.82 7.60000000
#> 783 T1 25 22.25 27.75 7.60000000
#> 784 T1 25 22.25 27.75 7.60000000
#> 785 T1 26 23.18 28.82 7.60000000
#> 786 T1 29 25.97 32.03 7.60000000
#> 787 T1 23 20.39 25.61 7.60000000
#> 788 T1 29 25.97 32.03 7.60000000
#> 789 T1 30 26.90 33.10 7.60000000
#> 790 T1 24 21.32 26.68 7.60000000
#> 791 T1 29 25.97 32.03 7.60000000
#> 792 T1 25 22.25 27.75 7.60000000
#> 793 T1 24 21.32 26.68 7.60000000
#> 794 T1 22 19.46 24.54 7.60000000
#> 795 T1 28 25.04 30.96 7.60000000
#> 796 T1 21 18.53 23.47 7.60000000
#> 797 T1 29 25.97 32.03 7.60000000
#> 798 T1 26 23.18 28.82 7.60000000
#> 799 T1 26 23.18 28.82 7.60000000
#> 800 T1 28 25.04 30.96 7.60000000
#> 801 T1 23 20.39 25.61 7.60000000
#> 802 T1 21 18.53 23.47 7.60000000
#> 803 T1 25 22.25 27.75 7.60000000
#> 804 T1 29 25.97 32.03 7.60000000
#> 805 T1 30 26.90 33.10 7.60000000
#> 806 T1 27 24.11 29.89 7.60000000
#> 807 T1 27 24.11 29.89 7.60000000
#> 808 T1 27 24.11 29.89 7.60000000
#> 809 T1 28 25.04 30.96 7.60000000
#> 810 T1 21 18.53 23.47 7.60000000
#> 811 T1 29 25.97 32.03 7.60000000
#> 812 T1 27 24.11 29.89 7.60000000
#> 813 T1 30 26.90 33.10 7.60000000
#> 814 T1 30 26.90 33.10 7.60000000
#> 815 T1 26 23.18 28.82 7.60000000
#> 816 T1 22 19.46 24.54 7.60000000
#> 817 T1 21 18.53 23.47 7.60000000
#> 818 T1 27 24.11 29.89 7.60000000
#> 819 T1 24 21.32 26.68 7.60000000
#> 820 T1 28 25.04 30.96 7.60000000
#> 821 T1 29 25.97 32.03 7.60000000
#> 822 T1 22 19.46 24.54 7.60000000
#> 823 T1 21 18.53 23.47 7.60000000
#> 824 T1 23 20.39 25.61 7.60000000
#> 825 T1 25 22.25 27.75 7.60000000
#> 826 T1 23 20.39 25.61 7.60000000
#> 827 T1 27 24.11 29.89 7.60000000
#> 828 T1 21 18.53 23.47 7.60000000
#> 829 T1 20 17.60 22.40 7.60000000
#> 830 T1 24 21.32 26.68 7.60000000
#> 831 T1 20 17.60 22.40 7.60000000
#> 832 T1 29 25.97 32.03 7.60000000
#> 833 T1 23 20.39 25.61 7.60000000
#> 834 T1 30 26.90 33.10 7.60000000
#> 835 T1 27 24.11 29.89 7.60000000
#> 836 T1 21 18.53 23.47 7.60000000
#> 837 T1 28 25.04 30.96 7.60000000
#> 838 T1 24 21.32 26.68 7.60000000
#> 839 T1 21 18.53 23.47 7.60000000
#> 840 T1 25 22.25 27.75 7.60000000
#> 841 T1 20 17.60 22.40 7.60000000
#> 842 T1 29 25.97 32.03 7.60000000
#> 843 T1 29 25.97 32.03 7.60000000
#> 844 T1 27 24.11 29.89 7.60000000
#> 845 T1 26 23.18 28.82 7.60000000
#> 846 T1 20 17.60 22.40 7.60000000
#> 847 T1 27 24.11 29.89 7.60000000
#> 848 T1 21 18.53 23.47 7.60000000
#> 849 T1 24 21.32 26.68 7.60000000
#> 850 T1 22 19.46 24.54 7.60000000
#> 851 T1 24 21.32 26.68 7.60000000
#> 852 T0 67 61.31 72.69 0.05714286
#> 853 T1 27 24.11 29.89 7.60000000
#> 854 T1 26 23.18 28.82 7.60000000
#> 855 T1 22 19.46 24.54 7.60000000
#> 856 T1 28 25.04 30.96 7.60000000
#> 857 T1 26 23.18 28.82 7.60000000
#> 858 T1 24 21.32 26.68 7.60000000
#> 859 T1 22 19.46 24.54 7.60000000
#> 860 T0 40 36.20 43.80 0.05714286
#> 861 T1 25 22.25 27.75 7.60000000
#> 862 T1 25 22.25 27.75 7.60000000
#> 863 T1 26 23.18 28.82 7.60000000
#> 864 T1 29 25.97 32.03 7.60000000
#> 865 T1 21 18.53 23.47 7.60000000
#> 866 T1 30 26.90 33.10 7.60000000
#> 867 T1 24 21.32 26.68 7.60000000
#> 868 T0 12 10.16 13.84 0.05714286
#> 869 T1 30 26.90 33.10 7.60000000
#> 870 T1 20 17.60 22.40 7.60000000
#> 871 T1 20 17.60 22.40 7.60000000
#> 872 T1 21 18.53 23.47 7.60000000
#> 873 T1 26 23.18 28.82 7.60000000
#> 874 T1 22 19.46 24.54 7.60000000
#> 875 T1 29 25.97 32.03 7.60000000
#> 876 T1 27 24.11 29.89 7.60000000
#> 877 T1 26 23.18 28.82 7.60000000
#> 878 T1 28 25.04 30.96 7.60000000
#> 879 T1 25 22.25 27.75 7.60000000
#> 880 T1 26 23.18 28.82 7.60000000
#> 881 T1 30 26.90 33.10 7.60000000
#> 882 T1 24 21.32 26.68 7.60000000
#> 883 T1 26 23.18 28.82 7.60000000
#> 884 T1 25 22.25 27.75 7.60000000
#> 885 T1 24 21.32 26.68 7.60000000
#> 886 T1 29 25.97 32.03 7.60000000
#> 887 T1 28 25.04 30.96 7.60000000
#> 888 T1 26 23.18 28.82 7.60000000
#> 889 T1 22 19.46 24.54 7.60000000
#> 890 T1 24 21.32 26.68 7.60000000
#> 891 T1 20 17.60 22.40 7.60000000
#> 892 T1 27 24.11 29.89 7.60000000
#> 893 T1 27 24.11 29.89 7.60000000
#> 894 T1 30 26.90 33.10 7.60000000
#> 895 T1 29 25.97 32.03 7.60000000
#> 896 T1 22 19.46 24.54 7.60000000
#> 897 T1 29 25.97 32.03 7.60000000
#> 898 T1 30 26.90 33.10 7.60000000
#> 899 T0 49 44.57 53.43 0.05714286
#> 900 T1 25 22.25 27.75 7.60000000
#> 901 T1 25 22.25 27.75 7.60000000
#> 902 T1 24 21.32 26.68 7.60000000
#> 903 T1 30 26.90 33.10 7.60000000
#> 904 T1 29 25.97 32.03 7.60000000
#> 905 T1 28 25.04 30.96 7.60000000
#> 906 T1 24 21.32 26.68 7.60000000
#> 907 T1 29 25.97 32.03 7.60000000
#> 908 T1 30 26.90 33.10 7.60000000
#> 909 T1 26 23.18 28.82 7.60000000
#> 910 T1 25 22.25 27.75 7.60000000
#> 911 T1 25 22.25 27.75 7.60000000
#> 912 T1 21 18.53 23.47 7.60000000
#> 913 T1 28 25.04 30.96 7.60000000
#> 914 T1 26 23.18 28.82 7.60000000
#> 915 T1 24 21.32 26.68 7.60000000
#> 916 T1 28 25.04 30.96 7.60000000
#> 917 T1 25 22.25 27.75 7.60000000
#> 918 T0 59 53.87 64.13 0.05714286
#> 919 T0 73 66.89 79.11 0.05714286
#> 920 T1 21 18.53 23.47 7.60000000
#> 921 T1 20 17.60 22.40 7.60000000
#> 922 T1 20 17.60 22.40 7.60000000
#> 923 T1 26 23.18 28.82 7.60000000
#> 924 T1 24 21.32 26.68 7.60000000
#> 925 T1 30 26.90 33.10 7.60000000
#> 926 T1 24 21.32 26.68 7.60000000
#> 927 T1 28 25.04 30.96 7.60000000
#> 928 T1 24 21.32 26.68 7.60000000
#> 929 T1 20 17.60 22.40 7.60000000
#> 930 T1 25 22.25 27.75 7.60000000
#> 931 T1 21 18.53 23.47 7.60000000
#> 932 T0 54 49.22 58.78 0.05714286
#> 933 T1 28 25.04 30.96 7.60000000
#> 934 T1 23 20.39 25.61 7.60000000
#> 935 T1 20 17.60 22.40 7.60000000
#> 936 T1 25 22.25 27.75 7.60000000
#> 937 T1 23 20.39 25.61 7.60000000
#> 938 T1 28 25.04 30.96 7.60000000
#> 939 T1 28 25.04 30.96 7.60000000
#> 940 T1 22 19.46 24.54 7.60000000
#> 941 T1 25 22.25 27.75 7.60000000
#> 942 T1 25 22.25 27.75 7.60000000
#> 943 T1 24 21.32 26.68 7.60000000
#> 944 T1 23 20.39 25.61 7.60000000
#> 945 T1 21 18.53 23.47 7.60000000
#> 946 T1 22 19.46 24.54 7.60000000
#> 947 T1 25 22.25 27.75 7.60000000
#> 948 T1 22 19.46 24.54 7.60000000
#> 949 T1 23 20.39 25.61 7.60000000
#> 950 T1 26 23.18 28.82 7.60000000
#> 951 T1 23 20.39 25.61 7.60000000
#> 952 T0 45 40.85 49.15 0.05714286
#> 953 T1 29 25.97 32.03 7.60000000
#> 954 T1 24 21.32 26.68 7.60000000
#> 955 T1 23 20.39 25.61 7.60000000
#> 956 T1 24 21.32 26.68 7.60000000
#> 957 T1 20 17.60 22.40 7.60000000
#> 958 T1 30 26.90 33.10 7.60000000
#> 959 T1 26 23.18 28.82 7.60000000
#> 960 T1 20 17.60 22.40 7.60000000
#> 961 T1 22 19.46 24.54 7.60000000
#> 962 T1 28 25.04 30.96 7.60000000
#> 963 T1 22 19.46 24.54 7.60000000
#> 964 T1 20 17.60 22.40 7.60000000
#> 965 T1 26 23.18 28.82 7.60000000
#> 966 T1 27 24.11 29.89 7.60000000
#> 967 T1 21 18.53 23.47 7.60000000
#> 968 T1 25 22.25 27.75 7.60000000
#> 969 T1 22 19.46 24.54 7.60000000
#> 970 T1 23 20.39 25.61 7.60000000
#> 971 T0 75 68.75 81.25 0.05714286
#> 972 T1 20 17.60 22.40 7.60000000
#> 973 T1 25 22.25 27.75 7.60000000
#> 974 T1 22 19.46 24.54 7.60000000
#> 975 T1 23 20.39 25.61 7.60000000
#> 976 T1 25 22.25 27.75 7.60000000
#> 977 T1 21 18.53 23.47 7.60000000
#> 978 T1 26 23.18 28.82 7.60000000
#> 979 T0 40 36.20 43.80 0.05714286
#> 980 T1 24 21.32 26.68 7.60000000
#> 981 T1 30 26.90 33.10 7.60000000
#> 982 T1 23 20.39 25.61 7.60000000
#> 983 T1 24 21.32 26.68 7.60000000
#> 984 T1 23 20.39 25.61 7.60000000
#> 985 T0 2 0.86 3.14 0.05714286
#> 986 T1 21 18.53 23.47 7.60000000
#> 987 T1 22 19.46 24.54 7.60000000
#> 988 T0 59 53.87 64.13 0.05714286
#> 989 T1 24 21.32 26.68 7.60000000
#> 990 T1 27 24.11 29.89 7.60000000
#> 991 T1 29 25.97 32.03 7.60000000
#> 992 T1 30 26.90 33.10 7.60000000
#> 993 T0 33 29.69 36.31 0.05714286
#> 994 T1 20 17.60 22.40 7.60000000
#> 995 T1 20 17.60 22.40 7.60000000
#> 996 T1 27 24.11 29.89 7.60000000
#> 997 T1 22 19.46 24.54 7.60000000
#> 998 T1 23 20.39 25.61 7.60000000
#> 999 T0 66 60.38 71.62 0.05714286
#> 1000 T1 27 24.11 29.89 7.60000000The package uses Conditional Probability Tables (CPTs) to model evidence under each hypothesis:
# CPT under H2 (population hypothesis)
cpt_h2 <- cpt_population(
propS = c(0.5, 0.5), # 50% female, 50% male
MPa = 30, # MP age
MPr = 5, # Age range
propC = c(0.3, 0.25, 0.2, 0.15, 0.1) # Hair color proportions
)
# CPT under H1 (MP hypothesis)
cpt_h1 <- cpt_missing_person(
MPs = 1, # Female
MPc = 2, # Hair color 2
eps = 0.05, # Sex error
epa = 0.05, # Age error
epc = error_matrix_hair() # Hair color error matrix
)
# View dimensions
dim(cpt_h1)
#> [1] 4 5| Function | Description |
|---|---|
lr_sex() |
LR for biological sex evidence |
lr_age() |
LR for age evidence |
lr_hair_color() |
LR for hair color evidence |
lr_birthdate() |
LR for birthdate evidence (Dirichlet model) |
lr_pigmentation() |
LR for multiple pigmentation traits |
| Function | Description |
|---|---|
sim_lr_genetic() |
Simulate genetic LRs from pedigrees |
sim_lr_prelim() |
Simulate non-genetic LRs |
sim_posterior() |
Combine prior and LRs for posterior odds |
sim_poi_genetic() |
Generate random DNA profiles |
sim_poi_prelim() |
Generate random non-genetic data |
| Function | Description |
|---|---|
decision_threshold() |
Find optimal LR threshold |
threshold_rates() |
Calculate TPR/FPR at different thresholds |
lr_combine() |
Combine genetic and non-genetic LRs |
| Function | Description |
|---|---|
plot_lr_distribution() |
Plot LR distributions under H1/H2 |
plot_decision_curve() |
ROC-like decision curves |
plot_cpt() |
Heatmap of CPT values |
The package includes allele frequency databases for genetic LR calculations:
# Available databases
data(Argentina)
data(Europe)
data(USA)
data(Asia)
data(Austria)
data(BosniaHerz)
data(China)
data(Japan)
# View structure
dim(Argentina)
#> [1] 93 25
names(Argentina)[1:10]
#> [1] "Allele" "D8S1179" "D21S11" "D7S820" "CSF1PO" "D3S1358" "THO1"
#> [8] "D13S317" "D16S539" "D2S1338"See the “Complete Workflow” vignette for a full example combining genetic and non-genetic evidence in a missing person case.
Marsico FL, Vigeland MD, Egeland T, Herrera Pinero F (2021). “Making decisions in missing person identification cases with low statistical power.” Forensic Science International: Genetics, 52, 102519. https://doi.org/10.1016/j.fsigen.2021.102519
Marsico FL, et al. (2023). “Likelihood ratios for non-genetic evidence in missing person cases.” Forensic Science International: Genetics, 66, 102891. https://doi.org/10.1016/j.fsigen.2023.102891