Description Details Author(s) References See Also Examples
The 'les' package estimates Loci of Enhanced Significance (LES) in tiling microarray data. These are regions of regulation such as found in differential transcription, CHiP-chip, or DNA modification analysis. The package provides a universal framework suitable for identifying differential effects in tiling microarray data sets, and is independent of the underlying statistics at the level of single probes.
The 'les' package provides a universal framework for detecting differential effects in tiling microarray experiments.
It is universal in the sense that one is free to choose any statistical test to estimate the effect of differential effect for each probe on the tiling microarray. Provided with the p-values for each probe and the corresponding positions of the probes, 'les' uses a sliding window approach to estimate the fraction of regulated probes in the local surrounding of each probe. The approach is related to computing a spatially resolved and weighted false discovery rate, and yields a interpretable statistical feature Lambda.
Resulting regions can be scored according to their overall effect. Methods for high-level plotting and export of the results to other software and genome browsers are provided.
The 'les' package is published under the GPL-3 license.
Julian Gehring, Clemens Kreutz, Jens Timmer
Maintainer: Julian Gehring <julian.gehring@fdm.uni-freiburg.de>
in preparation
This package is based on:
Kilian Bartholome, Clemens Kreutz, and Jens Timmer: Estimation of gene induction enables a relevance-based ranking of gene sets, Journal of Computational Biology: A Journal of Computational Molecular Cell Biology 16, no. 7 (July 2009): 959-967. http://www.liebertonline.com/doi/abs/10.1089/cmb.2008.0226
Class:
Les
Methods and functions:
Les
estimate
threshold
regions
ci
chi2
export
plot
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 | data(spikeInStat)
x <- Les(pos, pval)
x <- estimate(x, 200)
x <- threshold(x)
x <- regions(x)
subset <- pos >= 5232300 & pos <= 5233200
x <- ci(x, subset, conf=0.90, nBoot=50)
## plot data
plot(x, region=TRUE)
plot(x, region=TRUE, error="ci")
## Not run:
## export data of chromosome 1
export(x, file="les_out.bed", chr=1)
export(x, file="les_out.wig", format="wig", chr=1)
## End(Not run)
|
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.