bumphunter-methods | R Documentation |
bumphunter
in Package minfiEstimate regions for which a genomic profile deviates from its baseline value. Originally implemented to detect differentially methylated genomic regions between two populations, but can be applied to any CpG-level coefficient of interest.
## S4 method for signature 'GenomicRatioSet'
bumphunter(object, design, cluster=NULL,
coef=2, cutoff=NULL, pickCutoff=FALSE, pickCutoffQ=0.99,
maxGap=500, nullMethod=c("permutation","bootstrap"),
smooth=FALSE, smoothFunction=locfitByCluster,
useWeights=FALSE, B=ncol(permutations), permutations=NULL,
verbose=TRUE, type = c("Beta","M"), ...)
object |
An object of class GenomicRatioSet. |
design |
Design matrix with rows representing samples and columns representing covariates. Regression is applied to each row of mat. |
cluster |
The clusters of locations that are to be analyzed
together. In the case of microarrays, the clusters are many times
supplied by the manufacturer. If not available the function
|
coef |
An integer denoting the column of the design matrix containing the covariate of interest. The hunt for bumps will be only be done for the estimate of this coefficient. |
cutoff |
A numeric value. Values of the estimate of the genomic profile above the cutoff or below the negative of the cutoff will be used as candidate regions. It is possible to give two separate values (upper and lower bounds). If one value is given, the lower bound is minus the value. |
pickCutoff |
Should bumphunter attempt to pick a cutoff using the permutation distribution? |
pickCutoffQ |
The quantile used for picking the cutoff using the permutation distribution. |
maxGap |
If cluster is not provided this maximum location gap
will be used to define cluster via the |
nullMethod |
Method used to generate null candidate regions, must be one of ‘boots trap’ or ‘permutation’ (defaults to ‘permutation’). However, if covariates in addition to the outcome of interest are included in the design matrix (ncol(design)>2), the ‘permutation’ approach is not recommended. See vignette and original paper for more information. |
smooth |
A logical value. If TRUE the estimated profile will be smoothed with the
smoother defined by |
smoothFunction |
A function to be used for smoothing the estimate of the genomic
profile. Two functions are provided by the package: |
useWeights |
A logical value. If |
B |
An integer denoting the number of resamples to use when computing
null distributions. This defaults to 0. If |
permutations |
is a matrix with columns providing indexes to be used to
scramble the data and create a null distribution when
|
verbose |
logical value. If |
type |
Should bumphunting be performed on M-values ("M") or Beta values ("Beta")? |
... |
further arguments to be passed to the smoother functions. |
See help file for bumphunter
method in the
bumphunter
package for for details.
An object of class bumps
with the following components:
tab |
The table with candidate regions and annotation for these. |
coef |
The single loci coefficients. |
fitted |
The estimated genomic profile used to determine the regions. |
pvaluesMarginal |
marginal p-value for each genomic location. |
null |
The null distribution. |
algorithm |
details on the algorithm. |
Rafael A. Irizarry, Martin J. Aryee and Kasper D. Hansen
AE Jaffe, P Murakami, H Lee, JT Leek, MD Fallin, AP Feinberg, and RA Irizarry. Bump hunting to identify differentially methylated regions in epigenetic epidemiology studies. International Journal of Epidemiology (2012) 41(1):200-209. doi:10.1093/ije/dyr238
bumphunter
if(require(minfiData)) {
gmSet <- preprocessQuantile(MsetEx)
design <- model.matrix(~ gmSet$status)
bumps <- bumphunter(gmSet, design = design, B = 0,
type = "Beta", cutoff = 0.25)
}
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.