EVOLUTION-MANAGER
Edit File: mqmscan.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: Genome scan with a multiple QTL model (MQM)</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 mqmscan {qtl}"><tr><td>mqmscan {qtl}</td><td style="text-align: right;">R Documentation</td></tr></table> <h2>Genome scan with a multiple QTL model (MQM)</h2> <h3>Description</h3> <p>Genome scan with a multiple QTL model. </p> <h3>Usage</h3> <pre> mqmscan(cross, cofactors=NULL, pheno.col = 1, model=c("additive","dominance"), forceML=FALSE, cofactor.significance=0.02, em.iter=1000, window.size=25.0, step.size=5.0, logtransform = FALSE, estimate.map = FALSE, plot=FALSE, verbose=FALSE, outputmarkers=TRUE, multicore=TRUE, batchsize=10, n.clusters=1, test.normality=FALSE,off.end=0 ) </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>cofactors</code></td> <td> <p> List of cofactors to be analysed as cofactors in backward elimination procedure when building the QTL model. See <code><a href="mqmsetcofactors.html">mqmsetcofactors</a></code> on how-to manually set cofactors for backward elimination. Or use <code><a href="mqmautocofactors.html">mqmautocofactors</a></code> for automatic selection of cofactors. Only three kind of (integer) values are allowed in the cofactor list. (0: no cofactor at this marker, 1: Use this marker as an additive cofactor, 2: Use this marker as an sexfactor (Dominant cofactor)) </p> </td></tr> <tr valign="top"><td><code>pheno.col</code></td> <td> <p> Column number in the phenotype matrix which should be used as the phenotype. This can be a vector of integers; One may also give a character strings matching the phenotype names. Finally, one may give a numeric vector of phenotypeIDs. This should consist of integers with 0 < value < no. phenotypes. </p> </td></tr> <tr valign="top"><td><code>model</code></td> <td> <p> When scanning for QTLs should haplotype dominance be considered in an F2 intercross. Using the dominance model we scan for additive effects but also allow an additional effect where AA+AB versus BB and AA versus AB+BB. This setting is ignored for BC and RIL populations </p> </td></tr> <tr valign="top"><td><code>forceML</code></td> <td> <p> Specify which statistical method to use to estimate variance components to use when QTL modeling and mapping. Default usage is the Restricted maximum likelihood approach (REML). With this option a user can disable REML and use maximum likelihood. </p> </td></tr> <tr valign="top"><td><code>cofactor.significance</code></td> <td> <p> Significance level at which a cofactor is considered significant. This is estimated using an analysis of deviance, and compared to the level specified by the user. The cofactors that dont reach this level of statistical significance are NOT used in the mapping stage. Value between 0 and 1 </p> </td></tr> <tr valign="top"><td><code>em.iter</code></td> <td> <p> Maximum number of iterations for the EM algorithm to converge </p> </td></tr> <tr valign="top"><td><code>window.size</code></td> <td> <p> Window size for mapping QTL locations, this parameter is used in the interval mapping stage. When calculating LOD scores at a genomic position all cofactors within window.size are dropped to estimate the (unbiased) effect of the location under interest. </p> </td></tr> <tr valign="top"><td><code>step.size</code></td> <td> <p> Step size used in interval mapping. A lower step.size parameter increases the number of output points, this creates a smoother QTL profile </p> </td></tr> <tr valign="top"><td><code>off.end</code></td> <td> <p>Distance (in cM) past the terminal markers on each chromosome to which the genotype simulations will be carried. </p> </td></tr> <tr valign="top"><td><code>logtransform</code></td> <td> <p> Indicate if the algorithm should do a log transformation on the trait data in the pheno.col </p> </td></tr> <tr valign="top"><td><code>estimate.map</code></td> <td> <p> Should Re-estimation of the marker locations on the genetic map occur before mapping QTLs. This method is deprecated rather use the <code><a href="est.map.html">est.map</a></code> function in R/qtl. This is because no map is returned into the crossobject. The old map remains in the cross object. </p> </td></tr> <tr valign="top"><td><code>plot</code></td> <td> <p> plot the results (default FALSE)</p> </td></tr> <tr valign="top"><td><code>verbose</code></td> <td> <p> verbose output</p> </td></tr> <tr valign="top"><td><code>outputmarkers</code></td> <td> <p>If TRUE (the default), the results include the marker locations as well as along a grid of pseudomarkers; if FALSE, the results include only the grid positions.</p> </td></tr> <tr valign="top"><td><code>multicore</code></td> <td> <p>Use multicore (if available)</p> </td></tr> <tr valign="top"><td><code>batchsize</code></td> <td> <p>Number of traits being analyzed as a batch.</p> </td></tr> <tr valign="top"><td><code>n.clusters</code></td> <td> <p>Number of child processes to split the job into.</p> </td></tr> <tr valign="top"><td><code>test.normality</code></td> <td> <p>If TRUE, test whether the phenotype follows a normal distribution via <code><a href="mqmtestnormal.html">mqmtestnormal</a></code>. </p> </td></tr> </table> <h3>Value</h3> <p>When scanning a single phenotype the function returns a <code><a href="scanone.html">scanone</a></code> object. </p> <p>The object contains a matrix of three columns for LOD scores, information content and LOD*information content with pseudo markers sorted in increasing order. For more information on the scanone object see: <code><a href="scanone.html">scanone</a></code> </p> <h3>Note</h3> <p>The resulting scanone object itself can be visualized using the standard R/qtl plotting routines (<code><a href="plot.scanone.html">plot.scanone</a></code>) or specialized function to show the mqm model (<code><a href="mqmplotsingletrait.html">mqmplot.singletrait</a></code>) and QTL profile. If cofactors were specified the QTL model used in scanning is also returned as a named attribute of the scanone object called mqmmodel. It can be extracted from the resulting scanone object by using the <code><a href="mqmgetmodel.html">mqmgetmodel</a></code> function or the <code><a href="../../base/html/attr.html">attr</a></code> function. </p> <p>Also note the <code>estimate.map</code> parameter does not return its re-estimated genetic map, altough it is used internally. When scanning multiple genotypes a <code>mqmmulti</code> object is created. This object is just a list composed of scanone objects. The results for a single trait can be obtained from the <code>mqmmulti</code> object, in scanone format. </p> <h3>Author(s)</h3> <p>Ritsert C Jansen; Danny Arends; Pjotr Prins; Karl W Broman <a href="mailto:broman@wisc.edu">broman@wisc.edu</a> </p> <h3>See Also</h3> <ul> <li><p> The MQM tutorial: <a href="https://rqtl.org/tutorials/MQM-tour.pdf">https://rqtl.org/tutorials/MQM-tour.pdf</a> </p> </li> <li> <p><code><a href="MQM.html">MQM</a></code> - MQM description and references </p> </li> <li> <p><code><a href="mqmscan.html">mqmscan</a></code> - Main MQM single trait analysis </p> </li> <li> <p><code><a href="mqmscanall.html">mqmscanall</a></code> - Parallellized traits analysis </p> </li> <li> <p><code><a href="mqmaugment.html">mqmaugment</a></code> - Augmentation routine for estimating missing data </p> </li> <li> <p><code><a href="mqmautocofactors.html">mqmautocofactors</a></code> - Set cofactors using marker density </p> </li> <li> <p><code><a href="mqmsetcofactors.html">mqmsetcofactors</a></code> - Set cofactors at fixed locations </p> </li> <li> <p><code><a href="mqmpermutation.html">mqmpermutation</a></code> - Estimate significance levels </p> </li> <li> <p><code><a href="scanone.html">scanone</a></code> - Single QTL scanning </p> </li></ul> <h3>Examples</h3> <pre> data(map10) # Genetic map modeled after mouse # simulate a cross (autosomes 1-10) qtl <- c(3,15,1,0) # QTL model: chr, pos'n, add've & dom effects cross <- sim.cross(map10[1:10],qtl,n=100,missing.prob=0.01) # MQM crossaug <- mqmaugment(cross) # Augmentation cat(crossaug$mqm$Nind,'real individuals retained in dataset', crossaug$mqm$Naug,'individuals augmented\n') result <- mqmscan(crossaug) # Scan # show LOD interval of the QTL on chr 3 lodint(result,chr=3) </pre> <hr /><div style="text-align: center;">[Package <em>qtl</em> version 1.66 <a href="00Index.html">Index</a>]</div> </body></html>