BayesPeak

From SEQwiki
Jump to: navigation, search

Application data

Biological application domain(s) ChIP-seq, Simulation experiment
Principal bioinformatics method(s) Statistical calculation
Technology Illumina
Maintained? Yes
Input format(s) BED, GFF, BAM, or any other format supported by BioConductor
Output format(s) BED, GFF, BAM, or any other format supported by BioConductor
Software features Multicore
Programming language(s) R
Licence GPL

Summary: A Bayesian hidden Markov model to detect enriched locations in ChIP-seq data.

"Error: no local variable "counter" was set." is not a number.

Our proposed statistical algorithm, BayesPeak, uses a fully Bayesian hidden Markov model to detect enriched locations in the genome. The structure accommodates the natural features of the Solexa/Illumina sequencing data and allows for overdispersion in the abundance of reads in different regions. Moreover, a control sample can be incorporated in the analysis to account for experimental and sequence biases. Markov chain Monte Carlo algorithms are applied to estimate the posterior distributions of the model parameters, and posterior probabilities are used to detect the sites of interest. CONCLUSION: We have presented a flexible approach for identifying peaks from ChIP-seq reads, suitable for use on both transcription factor binding and histone modification data. Our method estimates probabilities of enrichment that can be used in downstream analysis.

Links


References

  1. . 2009. BMC Bioinformatics
  2. . 2011. Bioinformatics


To add a reference for BayesPeak, enter the PubMed ID in the field below and click 'Add'.

 


Search for "BayesPeak" in the SEQanswers forum / BioStar or:

Web Search Wiki Sites Scientific