EVOLUTION-MANAGER
Edit File: est.rf.html
<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Strict//EN" "http://www.w3.org/TR/xhtml1/DTD/xhtml1-strict.dtd"><html xmlns="http://www.w3.org/1999/xhtml"><head><title>R: Estimate pairwise recombination fractions</title> <meta http-equiv="Content-Type" content="text/html; charset=utf-8" /> <link rel="stylesheet" type="text/css" href="R.css" /> </head><body> <table width="100%" summary="page for est.rf {qtl}"><tr><td>est.rf {qtl}</td><td style="text-align: right;">R Documentation</td></tr></table> <h2>Estimate pairwise recombination fractions</h2> <h3>Description</h3> <p>Estimate the sex-averaged recombination fraction between all pairs of genetic markers. </p> <h3>Usage</h3> <pre> est.rf(cross, maxit=10000, tol=1e-6) </pre> <h3>Arguments</h3> <table summary="R argblock"> <tr valign="top"><td><code>cross</code></td> <td> <p>An object of class <code>cross</code>. See <code><a href="read.cross.html">read.cross</a></code> for details.</p> </td></tr> <tr valign="top"><td><code>maxit</code></td> <td> <p>Maximum number of iterations for the EM algorithm (not used with backcrosses).</p> </td></tr> <tr valign="top"><td><code>tol</code></td> <td> <p>Tolerance for determining convergence (not used with backcrosses).</p> </td></tr> </table> <h3>Details</h3> <p>For a backcross, one can simply count recombination events. For an intercross or 4-way cross, a version of the EM algorithm must be used to estimate recombination fractions. (Since, for example, in an intercross individual that is heterozygous at two loci, it is not known whether there were 0 or 2 recombination events.) Note that, for the 4-way cross, we estimate sex-averaged recombination fractions. </p> <h3>Value</h3> <p>The input <code>cross</code> object is returned with a component, <code>rf</code>, added. This is a matrix of size (tot.mar x tot.mar). The diagonal contains the number of typed meioses per marker, the lower triangle contains the estimated recombination fractions, and the upper triangle contains the LOD scores (testing rf = 0.5). </p> <h3>Author(s)</h3> <p>Karl W Broman, <a href="mailto:broman@wisc.edu">broman@wisc.edu</a> </p> <h3>See Also</h3> <p><code><a href="plot.rf.html">plotRF</a></code>, <code><a href="pull.rf.html">pull.rf</a></code>, <code><a href="plot.rfmatrix.html">plot.rfmatrix</a></code>, <code><a href="est.map.html">est.map</a></code>, <code><a href="badorder.html">badorder</a></code>, <code><a href="checkAlleles.html">checkAlleles</a></code> </p> <h3>Examples</h3> <pre> data(badorder) badorder <- est.rf(badorder) plotRF(badorder) </pre> <hr /><div style="text-align: center;">[Package <em>qtl</em> version 1.66 <a href="00Index.html">Index</a>]</div> </body></html>