EVOLUTION-MANAGER
Edit File: addqtl.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: Scan for an additional QTL in a multiple-QTL model</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 addqtl {qtl}"><tr><td>addqtl {qtl}</td><td style="text-align: right;">R Documentation</td></tr></table> <h2>Scan for an additional QTL in a multiple-QTL model</h2> <h3>Description</h3> <p>Scan for an additional QTL in the context of a multiple QTL model. </p> <h3>Usage</h3> <pre> addqtl(cross, chr, pheno.col=1, qtl, covar=NULL, formula, method=c("imp","hk"), model=c("normal", "binary"), incl.markers=TRUE, verbose=FALSE, tol=1e-4, maxit=1000, forceXcovar=FALSE, require.fullrank=FALSE) </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>chr</code></td> <td> <p>Optional vector indicating the chromosomes to be scanned. If missing, all chromosomes are scanned. Refer to chromosomes by name. Refer to chromosomes with a preceding <code>-</code> to have all chromosomes but those considered. A logical (TRUE/FALSE) vector may also be used.</p> </td></tr> <tr valign="top"><td><code>pheno.col</code></td> <td> <p>Column number in the phenotype matrix to be used as the phenotype. One may also give a character string matching a phenotype name. Finally, one may give a numeric vector of phenotypes, in which case it must have the length equal to the number of individuals in the cross, and there must be either non-integers or values < 1 or > no. phenotypes; this last case may be useful for studying transformations.</p> </td></tr> <tr valign="top"><td><code>qtl</code></td> <td> <p>An object of class <code>qtl</code>, as output from <code><a href="makeqtl.html">makeqtl</a></code>.</p> </td></tr> <tr valign="top"><td><code>covar</code></td> <td> <p>A matrix or data.frame of covariates. These must be strictly numeric.</p> </td></tr> <tr valign="top"><td><code>formula</code></td> <td> <p>An object of class <code><a href="../../stats/html/formula.html">formula</a></code> indicating the model to be fitted. (It can also be the character string representation of a formula.) QTLs are referred to as <code>Q1</code>, <code>Q2</code>, etc. Covariates are referred to by their names in the data frame <code>covar</code>. If the new QTL is not included in the formula, its main effect is added.</p> </td></tr> <tr valign="top"><td><code>method</code></td> <td> <p>Indicates whether to use multiple imputation or Haley-Knott regression.</p> </td></tr> <tr valign="top"><td><code>model</code></td> <td> <p>The phenotype model: the usual model or a model for binary traits</p> </td></tr> <tr valign="top"><td><code>incl.markers</code></td> <td> <p>If FALSE, do calculations only at points on an evenly spaced grid. If <code><a href="calc.genoprob.html">calc.genoprob</a></code> or <code><a href="sim.geno.html">sim.geno</a></code> were run with <code>stepwidth="variable"</code> or <code>stepwidth="max"</code>, we force <code>incl.markers=TRUE</code>.</p> </td></tr> <tr valign="top"><td><code>verbose</code></td> <td> <p>If TRUE, display information about the progress of calculations. If <code>verbose</code> is an integer > 1, further messages from <code><a href="scanqtl.html">scanqtl</a></code> are also displayed.</p> </td></tr> <tr valign="top"><td><code>tol</code></td> <td> <p>Tolerance for convergence for the binary trait model.</p> </td></tr> <tr valign="top"><td><code>maxit</code></td> <td> <p>Maximum number of iterations for fitting the binary trait model.</p> </td></tr> <tr valign="top"><td><code>forceXcovar</code></td> <td> <p>If TRUE, force inclusion of X-chr-related covariates (like sex and cross direction).</p> </td></tr> <tr valign="top"><td><code>require.fullrank</code></td> <td> <p>If TRUE, give LOD=0 when covariate matrix in the linear regression is not of full rank.</p> </td></tr> </table> <h3>Details</h3> <p>The formula is used to specified the model to be fit. In the formula, use <code>Q1</code>, <code>Q2</code>, etc., or <code>q1</code>, <code>q2</code>, etc., to represent the QTLs, and the column names in the covariate data frame to represent the covariates. </p> <p>We enforce a hierarchical structure on the model formula: if a QTL or covariate is in involved in an interaction, its main effect must also be included. </p> <p>If one wishes to scan for QTL that interact with another QTL, include it in the formula (with an index of one more than the number of QTL in the input <code>qtl</code> object). </p> <h3>Value</h3> <p>An object of class <code>scanone</code>, as produced by the <code><a href="scanone.html">scanone</a></code> function. LOD scores are relative to the base model (with any terms that include the new QTL omitted). </p> <h3>Author(s)</h3> <p>Karl W Broman, <a href="mailto:broman@wisc.edu">broman@wisc.edu</a></p> <h3>References</h3> <p>Haley, C. S. and Knott, S. A. (1992) A simple regression method for mapping quantitative trait loci in line crosses using flanking markers. <em>Heredity</em> <b>69</b>, 315–324. </p> <p>Sen, Ś. and Churchill, G. A. (2001) A statistical framework for quantitative trait mapping. <em>Genetics</em> <b>159</b>, 371–387. </p> <h3>See Also</h3> <p><code><a href="scanone.html">scanone</a></code>, <code><a href="fitqtl.html">fitqtl</a></code>, <code><a href="scanqtl.html">scanqtl</a></code>, <code><a href="refineqtl.html">refineqtl</a></code>, <code><a href="makeqtl.html">makeqtl</a></code>, <code><a href="addtoqtl.html">addtoqtl</a></code>, <code><a href="addpair.html">addpair</a></code>, <code><a href="addint.html">addint</a></code> </p> <h3>Examples</h3> <pre> data(fake.f2) # take out several QTLs and make QTL object qc <- c(1, 8, 13) qp <- c(26, 56, 28) fake.f2 <- subset(fake.f2, chr=c(1,2,3,8,13)) fake.f2 <- calc.genoprob(fake.f2, step=2, err=0.001) qtl <- makeqtl(fake.f2, qc, qp, what="prob") # scan for an additional QTL out1 <- addqtl(fake.f2, qtl=qtl, formula=y~Q1+Q2+Q3, method="hk") max(out1) # scan for an additional QTL that interacts with the locus on chr 1 out2 <- addqtl(fake.f2, qtl=qtl, formula=y~Q1*Q4+Q2+Q3, method="hk") max(out2) # plot interaction LOD scores plot(out2-out1) </pre> <hr /><div style="text-align: center;">[Package <em>qtl</em> version 1.66 <a href="00Index.html">Index</a>]</div> </body></html>