define.peak: Defines the QTL with support interval and exclusionary window

In eqtl: Tools for analyzing eQTL experiments: A complementary to Karl Broman's 'qtl' package for genome-wide analysis

Description

Define QTL with LOD drop support interval by using the results of a single QTL genome scan scanone and using a genetical exclusionary window.

Usage

define.peak(scanone,lodcolumn=1,chr,th=2.3,si=1.5,graph=FALSE,window.size=20,
round,save.pict=FALSE, phe.name, ...)

Arguments

scanone	An object of class scanone. See qtl package manual for scanone function details.
lodcolumn	Indicates on which of the LOD score columns (phenotypes) the QTLs should be defined. This can be "all" which indicates all LOD score columns. This can also be a vector of integers indicating which of the columns should be used or a strings vector matching the names of the LOD columns, the phenotypes' name, to analyse. See qtl package manual for scanone function details.
chr	An optional vector indicating the chromosomes for which QTLs should be defined.
th	A single numeric value which sets the LOD score significance threshold. Only peaks with LOD score above this value will be analysed.
si	A single numeric value which sets the QTL's Support Interval. si is the value of the accepted drop of LOD score to estimate the likely region on which the QTL is localized.
graph	If TRUE, draws the LOD curve with LOD peaks and support interval for the detected QTLs.
window.size	The exclusionary window size: A single numeric value setting the minimum genetic distance between two distinct QTLs to be considered.
save.pict	If TRUE, save the LOD curves drawn with support interval as png files named like "trait name"\_"chromosome"\_"a number".png in the current folder.
round	An optional integer indicating the precision to be used for the LOD score. See round function for details.
phe.name	An optional character string specifying the name of the analysed trait. When performing scanone on a single trait, the lodcolumn is named 'lod' and as the analysed trait.
...	Passed to the functions plot and plot.scanone when they are called (if graph=TRUE). Passed the maximum size of the genomic region parameter: m=10 should set 2*10cM for the inferior and the superior SI bounds from the position of the peak

Details

This function is used to detect and report QTL regions from a one-QTL genome scan performed by the scanone function. A QTL is considered as a genomic region defined by a maximum LOD score peak value, its position and the position of its support interval (here called dQuoteSI). The SI is estimated by the accepted drop of LOD score from the maximum LOD value defining the QTL region (the LOD peak). The FDR falls as the QTL SI size increases with lower LOD scores away from the peak. Usually we use si=1.5 or si=2. A genetic exclusionary window sets the minimum distance between two distinct QTLs which we consider being able to detect and depends directly on the size of the population. Due to the shape of the LOD curve, the drop of LOD score cannot be reached in some cases. Therefore a maximum SI size is set at 20 cM by default. m=10 will set 2*10cM for the inferior and the superior SI bounds. graph=TRUE specify to draw the LOD curves and the LOD SI on different chart for each QTL on their chromosome. No graphical setup has been defined and therefore they will be drawn one above the other in the same R graphical window. To setup the graph page and print all the charts in same window, one may use the graphical parameter mfrow of the R function par() according specific needs before launching define.peak. You may not want to set graph=TRUE and lodcolumn="all" at the same time depending on the amount of data. The parameter save.pict is useful to save systematically all charts generated by define.peak. These graphs are already page setted by the usual graphical functions (like mfrow).

Value

Returns an object of class peak which is a list of components corresponding to traits. names(peak) contains the names of the traits. Each trait is itself a list with elements corresponding to chromosomes. For chromosomes on which no QTL have been detected, peak\$trait\$chromosome contains a NA value (where chromosome is the number identifying the chromosome). For those on which a QTL has been detected peak\$trait\$chromosome contains a data frame where rows are detected QTLs and columns are peak features (which describe QTLs). names(peak\$trait\$chromosome) contains the peak features:

lod	The peak's LOD score.
mname.peak	The maximum LOD peak's (pseudo-)marker name.
peak.cM	The maximum LOD peak's genetic position in centiMorgan (cM).
mname.inf	The (pseudo-)marker's name corresponding to the inferior si bound.
inf.cM	The genetic position of the inferior SI bound in centiMorgan (cM).
mname.sup	The (pseudo-)marker's name corresponding to the superior SI bound.
sup.cM	The genetic position of the superior SI bound in centiMorgan (cM).
si.quality	The subjective quality if the support interval. Due to the shape of the LOD curves and the methods used to define the LOD peaks, the subjective quality of the QTLs are various.

The subjective quality of the support interval

A QTL whose support interval can be reached and defined, has more weight than a QTL whose support interval cannot and has been defined by its maximum size (argument m). This quality information corresponds to symbols indicating, how each were defined by the bounds of the QTL support interval. The symbols on the right side gives the information for the superior SI bounds and so on for the left sided bounds. '+' indicates that the LOD-drop support interval has been reached. '<-' and '->' indicates that the LOD-drop SI hasn't been reached before the maximum SI size (defined by m argument) for the inferior and the superior bounds respectively. '|' indicates that the LOD-drop SI has been delimited by the beginning or the end of the LOD curve either for the inferior or superior bounds respectively. Therefore, the quality symbols '|->' indicates that the SI has been delimited on the left by the beginning of the LOD curve and on the right by the maximum SI size. Therefore, the drop of LOD score is not reached on either the left or right. '+|' indicates that the SI has been reached on the left but has been delimited on the right by the end of the LOD curve.

Symbols -

Signification

"++" -

The QTL is bounded by a LOD-drop with both both SI sides reached.

"<-->" -

The QTL is bounded by the m parameter. The SI is not reached.

"+|" -

The QTL is bounded by the end of the chromosome on the right and by a SI on the left.

"|+" -

The QTL is bounded by the beginning of the chromosome on the left and by a SI on the right.

"<-|" -

The QTL is bounded by the end of the chromosome on the right and m/2 on the left.

"|->" -

The QTL is bounded by the beginning of the chromosome on the left and m/2 on the right.

Author(s)

Hamid A. Khalili

References

Broman KW, Wu H, Sen S, Churchill GA (2003) R/qtl: QTL mapping in experimental crosses. Bioinformatics 19:889-890

See Also

scanone,read.cross

Examples

data(seed10);

out.em <- scanone( seed10, pheno.col=1:50, model='normal', method='hk');

################################################
# Detecting QTL with LOD drop support interval #
################################################

# Defining QTLs for all traits and saving the curves in png files.
out.peak <- define.peak(out.em, 'all',graph=TRUE,save.pict=TRUE,round=3);

# Defining QTLs for few traits and drawing the curves.
par(mfrow=c(1,5));
out.peak <- define.peak(out.em,lodcolumn=c(3,4,40,49),graph=TRUE,round=3);
par(mfrow=c(1,1));

# Defining QTLs for one trait and drawing the curves.
out.peak <- define.peak(out.em,lodcolumn='CATrck',graph=TRUE,round=3);

Example output

Loading required package: qtl
no. of traits: 50
trait: ID	lodcolumn: 1
trait: CATcmt	lodcolumn: 2
trait: CATmjm	lodcolumn: 3
trait: CATppt	lodcolumn: 4
trait: CATvvr	lodcolumn: 5
trait: CATxni	lodcolumn: 6
trait: CATfgb	lodcolumn: 7
trait: CATqoq	lodcolumn: 8
trait: CATsgu	lodcolumn: 9
trait: CATwwh	lodcolumn: 10
trait: CATpxm	lodcolumn: 11
trait: CATkks	lodcolumn: 12
trait: CATbaa	lodcolumn: 13
trait: CATham	lodcolumn: 14
trait: MITO312	lodcolumn: 15
trait: CATudg	lodcolumn: 16
trait: CATwgo	lodcolumn: 17
trait: CATdbn	lodcolumn: 18
trait: CATipy	lodcolumn: 19
trait: CATpyn	lodcolumn: 20
trait: CATfnu	lodcolumn: 21
trait: CATjoj	lodcolumn: 22
trait: CATpsl	lodcolumn: 23
trait: CATntj	lodcolumn: 24
trait: CATjwf	lodcolumn: 25
trait: CATrck	lodcolumn: 26
trait: CATane	lodcolumn: 27
trait: CAToti	lodcolumn: 28
trait: CATbzx	lodcolumn: 29
trait: CATrjo	lodcolumn: 30
trait: CATjxh	lodcolumn: 31
trait: CATlkr	lodcolumn: 32
trait: CATnrd	lodcolumn: 33
trait: CATtez	lodcolumn: 34
trait: CATjpz	lodcolumn: 35
trait: CATbko	lodcolumn: 36
trait: CATavz	lodcolumn: 37
trait: CATvkm	lodcolumn: 38
trait: CATuhl	lodcolumn: 39
trait: CATlbg	lodcolumn: 40
trait: CATvbb	lodcolumn: 41
trait: CATowo	lodcolumn: 42
trait: CATlzh	lodcolumn: 43
trait: CATepg	lodcolumn: 44
trait: CATbuu	lodcolumn: 45
trait: CATwrn	lodcolumn: 46
trait: CATtgk	lodcolumn: 47
trait: CATxmc	lodcolumn: 48
trait: CATdmq	lodcolumn: 49
trait: CATbzp	lodcolumn: 50
define.peak in process...
no. of traits: 4
trait: CATmjm	lodcolumn: 3
trait: CATppt	lodcolumn: 4
trait: CATlbg	lodcolumn: 40
trait: CATdmq	lodcolumn: 49
define.peak in process...
no. of traits: 1
trait: CATrck	lodcolumn: 26
define.peak in process...

eqtl documentation built on May 2, 2019, 5:42 a.m.