The hardware and bandwidth for this mirror is donated by dogado GmbH, the Webhosting and Full Service-Cloud Provider. Check out our Wordpress Tutorial.
If you wish to report a bug, or if you are interested in having us mirror your free-software or open-source project, please feel free to contact us at mirror[@]dogado.de.

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

These binaries (installable software) and packages are in development.
They may not be fully stable and should be used with caution. We make no claims about them.
Health stats visible at Monitor.