statMark | R Documentation |
Individual marker and interval statistics for an R/qtl cross object
statMark(cross, chr, stat.type = c("marker","interval"),
map.function = "kosambi")
cross |
An qtl |
chr |
Character vector of linkage group names used for subsetting the linkage map. |
stat.type |
Character string of either |
map.function |
Character string of either |
If "marker"
is chosen then a call to geno.table
from
qtl is used to return individual marker statistics for segregation distortion,
as well as allele and missing value proportions. For the current map
order the number of double crossovers at each marker are also returned.
If "interval"
is chosen then interval statistics are returned for
the current map order. These include the estimated recombination
fraction and LOD score between adjacent markers, calculated from
est.rf
in qtl. Also returned are the map interval distances and
converted map recombination fractions extracted from the "map"
component of each linkage group as well as the actual number of
recombinations between markers.
This function is used in profileMark
to plot any combination of
returned linkage map statistics on a single graphical display.
A list named by the stat.type
used in the call. Each element is a
data frame of statistics with columns named by the statistic.
Julian Taylor
Taylor, J., Butler, D. (2017) R Package ASMap: Efficient Genetic Linkage Map Construction and Diagnosis. Journal of Statistical Software, 79(6), 1–29.
profileMark
data(mapDH, package = "ASMap")
## produce all statistics
sm <- statMark(mapDH, stat.type = c("marker","interval"))
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.