CHIP Seq Analysis by DiffBind Package
Hi I am working on ChipSeq data of histone marks. I am using Diffbind package. I have 8 samples in 2 groups each. I am trying to see differential peak enrichment in these two groups. I have used this to read one file at a time. How can I club 4 samples so that I can do a group comparison?
bedJ <- as.data.frame(read.table("CEMT_178.H3K27ac-J.bed",header = FALSE, sep="\t",stringsAsFactors=FALSE, quote=""))
Thanks Shrinka
• 1,235 views
•
link
updated
by
Rory Stark
•
written
by
shrinka.genetics •
•
written
by
shrinka.genetics •
1 answer
The easiest way to get these data into DiffBind is to to use a sample sheet with entries for the aligned reads, peaks, and other metadata for each sample.
• 0 views
•
link
written
by
Rory Stark
Log in to answer this question.
More posts like this
-
.BigWig file has more peaks than reads in .Bam file?
written by Aryan •
-
Diffbind output file default contrasts
written by c.eskiw •
-
chromoMap Error: .rowNamesDF<-`(x, value = value) : invalid 'row.names' length
written by kamanovae •
-
Diffbind adding my own consensus peaks
written by lara.bedeyan •
-
CHIP Seq Analysis by DiffBind Package
written by shrinka.genetics •
-
my file do not have duplicates, but it still shows duplicate are not allow
written by mikysyc2016
-
strange results from differential methylation analysis using limma
written by RC •
-
reading a bed file in R
written by zizigolu
-
Problem with finding overlapping chips peaks with ChIPpeakAnno
written by Debbie •
-
Geoquery And Gse Error
written by Avoks •
What is unclear after reading the manual, especially section 3? https://www.bioconductor.org/packages/release/bioc/vignettes/DiffBind/inst/doc/DiffBind.pdf