This is a test version of Biostars. For the public version, visit https://www.biostars.org.
How to start with bamsignals

Hi everybody! I'm trying to discover new ChIP-seq tools, I'm trying bamsignals but in this tutorial fro me is not clear what is bampath and genes (and how to perform them with my data) does anyone have experience in this package and could give me his opinion or advice?

Thanks in advance!

chip-seq bam bamsignal bioconductor r
• 2,034 views • link
updated by GenoMax GenoMax • written by Lila M Lila M

The info is right there in the link you posted

bampath <- system.file("extdata", "randomBam.bam", package="bamsignals")
genes <-  get(load(system.file("extdata", "randomAnnot.Rdata", package="bamsignals")))

Are you not familiar with how R works?

• 0 views • link
written by GenoMax GenoMax

Yes, I know that , but this info is to get the "toy data" and I din't know what I need for my personal data, please, if you understand the code, could you explain what is extdata, randomAnnot.Rdata and how can I get them in my own?

Thanks

• 0 views • link
written by Lila M Lila M

You are reading paths for randomBam.bam/randomAnnot.Rdata from a directory extdata (for package bamsignals). bampath will contain the absolute path to the bam file on your disk. extdata generally contains the raw data included with a R-package.

• 0 views • link
written by GenoMax GenoMax

is still not clear for me how load the BAM files and processes them, I've also tried Rsamtools but not success at all, help please!

• 0 views • link
written by Lila M Lila M

0 answers

No answers yet.

Log in to answer this question.

More posts like this
  • Eigencor plot after limma::removeBatchEffect
    written by f2ff7a46 • f2ff7a46
    I'm trying to create an eigencor plot of clinical variables against principle components using the PCAtools guidelines. [https://www.bioconductor.org/packages/devel/bioc/vignettes/PCAtools/inst/doc/PCAtools.html#correlate-the-principal-components-back-to-the-clinical-data][1] I want to first eliminate batch effects …
  • GSVA from DESeq2
    written by herny.bahus • herny.bahus
    Hi, I'm analyzing RNA-seq data and i followed DESeq2 tutorial (http://www.bioconductor.org/packages/release/bioc/vignettes/DESeq2/inst/doc/DESeq2.html#why-un-normalized-counts) starting from the raw counts for gene expression, I want now to perform GSVA …
  • metaSeq: Meta-analysis of RNA-seq count data
    written by Anand • Anand
    Hi everyone. Has anyone here used the [metaSeq package][1] to carry out a meta-analysis of raw counts data from RNA-seq runs from different studies? I …
  • how to select peaks at not transcribed genes
    written by Lila M Lila M
    Hi there, Maybe this is a silly question. I have several peaks called with macss (n = 12808) and I want to create two groups …
  • Does it make sense to use methylation data for GESA
    written by CY CY
    I came across [ChAMP][1] which analyzes methylation data using GSEA. I know that DNA methylation and corresponding mRNA expression are interconnected in some degree. However, …
  • Is there a name for these charts?
    written by russhh russhh
    There's some plots on the fgsea [vignette][1] that I'd like to refer to by name. It's the chart returned by the plotGseaTable function that I'm …
  • What is the meaning of Bumphunter output?
    written by boyu93 • boyu93
    I tried using ChAMP to analyze differential Methylation Regions (DMR).[ChAMP][1] The method I used was based on Bumphunter. But its output made me confused. Especially …
  • How to annotate TTS or TES
    written by Lila M Lila M
    Hi everybody, I'm performing some gene profiles from my ChiP seq experiment. To do that, I've used [homer][1] or [ChipSeeker][2] for the annotation. I am …
  • How do I find out regions of methylation enrichment using MEDIPS bioconductor package?
    written by lakhujanivijay lakhujanivijay
    I have a single MEDIP-seq sample and I am only interested in finding out regions of methylation. How do I go about it? Since, I …
  • DNAcopy plateaus plot
    written by 2nelly 2nelly
    Hi all, Could someone briefly explain me why a plateau plot in DNAcopy R package can be used to define thresholds in CNA analysis? [DNAcopy …