EVOLUTION-MANAGER
Edit File: checkAlleles.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: Identify markers with switched alleles</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 checkAlleles {qtl}"><tr><td>checkAlleles {qtl}</td><td style="text-align: right;">R Documentation</td></tr></table> <h2>Identify markers with switched alleles</h2> <h3>Description</h3> <p>Identify markers whose alleles might have been switched by comparing the LOD score for linkage to all other autosomal markers with the original data to that when the alleles have been switched. </p> <h3>Usage</h3> <pre>checkAlleles(cross, threshold=3, verbose)</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>threshold</code></td> <td> <p>Only an increase in maximum 2-point LOD of at least this amount will lead to a marker being flagged.</p> </td></tr> <tr valign="top"><td><code>verbose</code></td> <td> <p>If TRUE and there are no markers above the threshold, print a message.</p> </td></tr> </table> <h3>Details</h3> <p>For each marker, we compare the maximum LOD score for the cases where the estimated recombination fraction > 0.5 to those where r.f. < 0.5. The function <code><a href="est.rf.html">est.rf</a></code> must first be run. </p> <p><b>Note</b>: Markers that are tightly linked to a marker whose alleles are switched are likely to also be flagged by this method. The real problem markers are likely those with the biggest difference in LOD scores. </p> <h3>Value</h3> <p>A data frame containing the flagged markers, having four columns: the marker name, chromosome ID, numeric index within chromosome, and the difference between the maximum two-point LOD score with the alleles switched to that from the original data. </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="est.rf.html">est.rf</a></code>, <code><a href="geno.crosstab.html">geno.crosstab</a></code>, <code><a href="switchAlleles.html">switchAlleles</a></code> </p> <h3>Examples</h3> <pre> data(fake.f2) # switch homozygotes at marker D5M391 fake.f2 <- switchAlleles(fake.f2, "D5M391") fake.f2 <- est.rf(fake.f2) checkAlleles(fake.f2) </pre> <hr /><div style="text-align: center;">[Package <em>qtl</em> version 1.46-2 <a href="00Index.html">Index</a>]</div> </body></html>