seeded.map: Construct the linkage map for a sequence of markers after...

In bschiffthaler/BatchMap: Software for the Creation of High Density Linkage Maps in Outcrossing Species

Description

Estimates the multipoint log-likelihood, linkage phases and recombination frequencies for a sequence of markers in a given order using seeded phases.

Usage

seeded.map(input.seq, tol = 1e-04, phase.cores = 1, seeds, verbosity = NULL)

Arguments

input.seq	an object of class sequence.
tol	tolerance for the C routine, i.e., the value used to evaluate convergence.
phase.cores	The number of parallel processes to use when estimating the phase of a marker. (Should be no more than 4)
seeds	A vector given the integer encoding of phases for the first N positions of the map
verbosity	A character vector that includes any or all of "batch", "order", "position", "time" and "phase" to output progress status information.

Details

Markers are mapped in the order defined in the object input.seq. The best combination of linkage phases is also estimated starting from the first position not in the given seeds.The multipoint likelihood is calculated according to Wu et al. (2002b)(Eqs. 7a to 11), assuming that the recombination fraction is the same in both parents. Hidden Markov chain codes adapted from Broman et al. (2008) were used.

Value

An object of class sequence, which is a list containing the following components:

seq.num	a vector containing the (ordered) indices of markers in the sequence, according to the input file.
seq.phases	a vector with the linkage phases between markers in the sequence, in corresponding positions. -1 means that there are no defined linkage phases.
seq.rf	a vector with the recombination frequencies between markers in the sequence. -1 means that there are no estimated recombination frequencies.
seq.like	log-likelihood of the corresponding linkage map.
data.name	name of the object of class outcross with the raw data.
twopt	name of the object of class rf.2pts with the 2-point analyses.

Author(s)

Adapted from Karl Broman (package 'qtl') by Gabriel R A Margarido, [e-mail] and Marcelo Mollinari, [e-mail]. Modified to use seeded phases by Bastian Schiffthaler [e-mail]

References

Broman, K. W., Wu, H., Churchill, G., Sen, S., Yandell, B. (2008) qtl: Tools for analyzing QTL experiments R package version 1.09-43

Jiang, C. and Zeng, Z.-B. (1997). Mapping quantitative trait loci with dominant and missing markers in various crosses from two inbred lines. Genetica 101: 47-58.

Lander, E. S., Green, P., Abrahamson, J., Barlow, A., Daly, M. J., Lincoln, S. E. and Newburg, L. (1987) MAPMAKER: An interactive computer package for constructing primary genetic linkage maps of experimental and natural populations. Genomics 1: 174-181.

Wu, R., Ma, C.-X., Painter, I. and Zeng, Z.-B. (2002a) Simultaneous maximum likelihood estimation of linkage and linkage phases in outcrossing species. Theoretical Population Biology 61: 349-363.

Wu, R., Ma, C.-X., Wu, S. S. and Zeng, Z.-B. (2002b). Linkage mapping of sex-specific differences. Genetical Research 79: 85-96

See Also

make.seq

Examples

  data(example.out)
  twopt <- rf.2pts(example.out)

  markers <- make.seq(twopt,c(30,12,3,14,2))
  seeded.map(markers, seeds = c(4,2))

bschiffthaler/BatchMap documentation built on Dec. 16, 2019, 2:22 a.m.