Description Usage Arguments Details Value See Also Examples
Helper function to create the input for the prosecution.
1 2 3 4 5 6 |
cspFile |
Path to the crime scene profile. |
refFile |
Path to the known profiles. |
ethnic |
Ethnicity, e.g. subgroup within the allele frequency database. |
nUnknowns |
Number of unknown contributors for which to perform calculation. |
adj |
Allele frequency adjustment parameter. |
fst |
Allele fraction adjustment F_{ST}. |
databaseFile |
Path to the allele database. If NULL, then defaults to
the NGMSelect database provided with |
linkageFile |
Path to recombination rate information. If NULL then defaults
to the linkage file provided with |
relationship |
Specified relationship between Q and X. Can take values of 0=unrelated, 1=parent/offspring, 2=siblings, 3=uncle/nephew, 4=half-uncle/half-nephew, 5=cousins, 6=grandparent/grandchild, 7=half-siblings. |
doDropin |
Whether or not to model drop-in. |
combineRare |
Whether or not to combine rare unobserved alleles into a single allele. |
rareThreshold |
If |
kit |
Parameter specifying which allele database supplied with likeLTD to use if
|
... |
Other parameters to determine how to perform calculations. |
It loads the CSP, known profiles, and allele database from file. It makes sure the relatedness factor is zero. It sets the reference individual to the first queried individual, unless it is not subject to dropout. In that case, the reference individual is the first individual with dropout.
A list of named input parameters, or hypothesis, suitable for the prosecution.
defence.hypothesis
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 | ## Not run:
# datapath to example files
datapath = file.path(system.file("extdata", package="likeLTD"),"hammer")
# File paths and case name for allele report
admin = pack.admin.input(
cspFile = file.path(datapath, 'hammer-CSP.csv'),
refFile = file.path(datapath, 'hammer-reference.csv'),
caseName = "hammer",
kit= "SGMplus"
)
# Enter arguments
args = list(
nUnknowns = 1,
doDropin = FALSE,
ethnic = "EA1",
adj = 1,
fst = 0.02,
relatedness = c(0,0)
)
# Create hypotheses
hypP = do.call(prosecution.hypothesis, append(admin,args))
hypD = do.call(defence.hypothesis, append(admin,args))
## End(Not run)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.