Introduction to mispitools

Franco L. Marsico

2026-01-08

Overview

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.

Key Features

Installation

# Install from CRAN (when available)
install.packages("mispitools")

# Or install the development version from GitHub
# devtools::install_github("MarsicoFL/mispitools")

The Bayesian Framework

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}\]

Quick Start

Non-Genetic Evidence

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.9

Calculate 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.60000000

Conditional Probability Tables

The 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

Visualizing CPTs

# Visualize both CPTs and LR heatmap
plot_cpt(cpt_h2, cpt_h1)

Available Functions

LR Calculation Functions

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

Simulation Functions

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

Decision Analysis

Function Description
decision_threshold() Find optimal LR threshold
threshold_rates() Calculate TPR/FPR at different thresholds
lr_combine() Combine genetic and non-genetic LRs

Visualization

Function Description
plot_lr_distribution() Plot LR distributions under H1/H2
plot_decision_curve() ROC-like decision curves
plot_cpt() Heatmap of CPT values

Interactive Apps

Launch interactive Shiny applications:

# Basic CPT explorer
app_mispitools()

# Advanced LR comparison with ROC analysis
app_lr_comparison()

Population Frequency Databases

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"

Next Steps

See the “Complete Workflow” vignette for a full example combining genetic and non-genetic evidence in a missing person case.

References

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