magnusdv / forrel

Forensic pedigree analysis and relatedness inference
GNU General Public License v2.0
10 stars 0 forks source link

LRpower output #40

Closed thoree closed 4 years ago

thoree commented 4 years ago

I believe that there is a typo in the last line of the output below and that it should be P(LR ≥ 0)

library(forrel)
#> Loading required package: pedtools
PO = nuclearPed(fa = "A", ch = "B")
UN = list(singleton("A"), singleton("B"))
LRpower(PO, UN, UN, ids = c("A", "B"), nsim = 1, threshold = 0, 
        alleles = 1:100, seed = 123, verbose = F)
#> LR power analysis using individuals A, B.
#> Mean total LR: 0 
#> Mean total log10(LR): -Inf 
#> Estimated inclusion power: 
#>   P(LR > 0) = 1

Created on 2020-08-26 by the reprex package (v0.3.0)