magnusdv / forrel

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

`markers` problem in `missingPersonIP` #25

Closed thoree closed 4 years ago

thoree commented 4 years ago

Great functions: missingPersonIP and missingPersonEP. So far only one problem detected:

# Four siblings; the fourth is missing
library(forrel, quietly = T)
x = nuclearPed(4)

# Remaining sibs typed with 5 triallelic markers
x = markerSim(x, N = 5, ids = 3:5, alleles = 1:3, seed = 123, verbose = FALSE)

# Compute exclusion power statistics
missingPersonIP(x, missing = 6, markers = 1,
                nsim = 5, threshold = c(10, 100))
#> Simulating 5 profiles...done
#> Computing likelihood ratios...
#> Error in `rownames<-`(`*tmp*`, value = markers): attempt to set 'rownames' on an object with no dimensions

Created on 2019-10-26 by the reprex package (v0.3.0)