Match bloodmeals and humans based on log10LR static threshold

Usage

match_static_thresh(log10_lrs, thresh)

Arguments

log10_lrs: Output from calc_log10_lrs() or from bistro with return_lrs = TRUE (lrs: the second element in the list)
thresh: log10LR threshold for matching. All bloodmeal-human pairs with a log10LR ≥ thresh will be considered a match.

Value

Dataframe including columns similar to the bistro() output.

Examples

bistro_output <- bistro(bloodmeal_profiles, human_profiles,
  pop_allele_freqs = pop_allele_freqs,
  kit = "ESX17", peak_thresh = 200, return_lrs = TRUE
)
#> 1/17 markers in kit but not in pop_allele_freqs: AMEL
#> Formatting bloodmeal profiles
#> Removing 6 peaks under the threshold of 200 RFU.
#> For 1/4 bloodmeal ids, all peaks are below the threshold
#> Formatting human profiles
#> Markers being used: D10S1248, D12S391, D16S539, D18S51, D19S433, D1S1656, D21S11, D22S1045, D2S1338, D2S441, D3S1358, D8S1179, FGA, SE33, TH01, VWA
#> Calculating log10LRs
#> # bloodmeal ids: 3
#> # human ids: 3
#> Bloodmeal id 1/3
#> Human id 1/3
#> Human id 2/3
#> Human id 3/3
#> Bloodmeal id 2/3
#> Human id 1/3
#> Human id 2/3
#> Human id 3/3
#> Bloodmeal id 3/3
#> Human id 1/3
#> Human id 2/3
#> Human id 3/3
#> Identifying matches
match_static_thresh(bistro_output$lrs, 10)
#> # A tibble: 4 × 6
#>   bloodmeal_id locus_count est_noc match human_id log10_lr
#>   <chr>              <int>   <dbl> <chr> <chr>       <dbl>
#> 1 evid1                 16       2 yes   P1           21.8
#> 2 evid1                 16       2 yes   P2           10.3
#> 3 evid2                  1       2 no    NA           NA  
#> 4 evid3                  8       1 no    NA           NA