Description Usage Arguments Details Value Examples
Estimate map distances based on the estimated recombination fractions
1 2 3 4 5 6 7 | estimateMap(
mpcrossLG,
mapFunction = rfToHaldane,
maxOffset = 1,
maxMarkers = 2000,
verbose = FALSE
)
|
mpcrossLG |
An object of class |
mapFunction |
The map function to use to compute recombination fractions to centiMorgan distances. |
maxOffset |
The maximum separation between pairs of markers used for map construction, in terms of position within the ordering. Recombination fractions between pairs of markers, which are further apart than this, will not be used to estimate the map distances. |
maxMarkers |
The (approximate) number of markers for which distances are estimated simultaneously. |
verbose |
Should verbose output be produced? |
Once a marker order has been chosen, one possible way of estimating a genetic map is to convert the recombination fractions between adjacent markers into centiMorgan distances. This tends not to work well, because individual recombination fraction estimates can be highly variable, depending on the experimental design used, and the distribution of the marker alleles among the founders. It also wastes much of the information contained in the data; we can estimate recombination fractions between all pairs of markers, rather than just adjacent markers, and this information should be used in the estimation of map distances
This function uses non-linear least squares to estimate map distances as follows. Assume that there are n
markers on a chromosome, and for all pairs of markers there is an available estimate of the recombination fraction. For every pair of markers which differ by maxOffset
or less, in terms of their position within the ordering, the recombination fraction between these markers is turned into a centiMorgan distance. This centiMorgan distance is expressed as a sum of distances between adjacent markers, which is a simple equation. The set of all the equations generated in this way is represented as a matrix equation, and solved via non-linear least squares. As these non-linear least squares problems can become very large, input maxMarkers
allows the non-linear least squares problem to be broken into several smaller problems.
For example, assume that there are five markers, for which an order has been determined. The distance between markers i and j, as estimated by the recombination fractions, is d(i, j). The genetic distance between markers i and i + 1 in the final genetic map is a(i). So in this case, the parameters that are to be estimated are a(1), a(2), a(3) and a(4). If maxOffset
is 3, then the set of equations generated is
d(1, 3) = a(1) + a(2)
d(1, 4) = a(1) + a(2) + a(3)
d(2, 4) = a(2) + a(3)
d(3, 5) = a(3) + a(4)
d(2, 5) = a(2) + a(3) + a(4)
These constraints are represented as a matrix equation and solved for a(1), a(2), a(3) and a(4) using non-linear least squares. However, if maxOffset
is set to 2
, then the set of equations is
d(1, 3) = a(1) + a(2)
d(2, 4) = a(2) + a(3)
d(3, 5) = a(3) + a(4)
A map object, in the format specified by the qtl-package
package. This format is a list of chromosomes, with each entry being a named numeric vector of marker positions.
1 2 3 4 5 6 7 8 9 | data(simulatedFourParentData)
#Estimate recombination fractions
rf <- estimateRF(simulatedFourParentData)
#Assign all markers to one linkage group / chromosome
grouped <- formGroups(rf, groups = 1)
#Estimate map
estimatedMap <- estimateMap(grouped, maxOffset = 10)
#Create object that includes the map
mapped <- new("mpcrossMapped", grouped, map = estimatedMap)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.