New to ChipSeq getting 0 peaks using dba()
1
0
Entering edit mode
@fce3b503
Last seen 4 days ago
United States

I am very new to ChipSeq and am just so confused. I have my sequence data, the control and the reference from NiH and UCSC. I have successfully (I think) , mapped the reference, aligned them with bowtie and found peaks with MACS2. The output of MACS2 was a .r file, a .xls file, a .narrowPeak file and a .bed file. I now want to move onto finding differentials. Everything I read watch tells me I need a sample sheet. Is this the .xls file? Or do I need to create something else? Right now this is my R script:

library(ChIPQC)
library(DiffBind)
library(tidyverse)

db <- dba(sampleSheet =sampleSheet)


When I look at the created object the output is:

15741 Samples, 0 sites in matrix


Obviously I'm doing something wrong. Any help would be greatly appreciated.

ChipSeq ChIPSeqData chipseqDB • 153 views
0
Entering edit mode
Rory Stark ★ 4.6k
@rory-stark-5741
Last seen 5 days ago
CRUK, Cambridge, UK

The sample sheet is a separate .csv file that contains a row for each sample, with columns indicating the locations of the peak and bam files, as well as experimental metadata. The columns are document in the help page for ?dba.

You should think about working through the example in the vignette before you attempt your own analysis. You can see the example sample sheet as follows:

samplesheet <-  system.file("extra/tamoxifen.csv", package="DiffBind")
samples


For the peak files, you can use either the .xls files (PeakCaller="macs") or the .narrowPeak files (PeakCaller="narrow").

0
Entering edit mode

Thanks so much...I got it to work