cran / simDNAmixtures

:exclamation: This is a read-only mirror of the CRAN R package repository. simDNAmixtures — Simulate Forensic DNA Mixtures
0 stars 1 forks source link

Two markers not simulated? #1

Open thoree opened 1 year ago

thoree commented 1 year ago

Nice and useful library! Apparently Penta D and Penta E are not simulated:

library(simDNAmixtures)
# Essentially first example in sample_mixtures documentation 
# (changed to n = 1 and added seed = 1729)
# remotes::install_github("https://github.com/cran/simDNAmixtures/", force = T)

# freqs <- read_allele_freqs(system.file("extdata","FBI_extended_Cauc.csv",
#                            package = "simDNAmixtures")) #problem with reprex
freqs <- read_allele_freqs("C:/Rlibrary/simDNAmixtures/extdata/FBI_extended_Cauc.csv")
data(gf)

sampling_parameters <- list(min_mu = 50., max_mu = 5e3,
                           min_cv = 0.05, max_cv = 0.35,
                           degradation_shape1 = 10, degradation_shape2 = 1)

mixtures <- sample_mixtures(n = 1, contributors = c("U1", "U2"), freqs = freqs,
                            sampling_parameters = sampling_parameters,
                            model_settings = gf$gamma_settings_no_stutter,
                            sample_model = sample_gamma_model, seed = 1729)

# Markers in in database of frequencies and in mixture
markersInDB  <- names(freqs)
markersInMixture <- unique(mixtures$smash$Marker)
setdiff(markersInDB, markersInMixture)
#> [1] "Penta E" "Penta D"

Created on 2023-09-18 with reprex v2.0.2

gaborcsardi commented 1 year ago

Hi, this is a read only mirror of CRAN, please see the package authors in the DESCRIPTION file. Look for Maintainer, BugReports and URL. Thanks!