How does diffbind extract the original coordinates of each difference peak
After using dba. analyze() to calculate the difference, use dba. report() to extract the difference peak, but the difference peak coordinates have been merged. How to output the real peak coordinates corresponding to each difference peak in different sample species?
> library(DiffBind)
> dbObj_S3_S2 <-dba(sampleSheet="D5_D3_Linux.csv")
> dbObj_S3_S2 <- dba.count(dbObj_S3_S2 , bUseSummarizeOverlaps=TRUE)
> dbObj_S3_S2 <- dba.normalize(dbObj_S3_S2)
> dbObj_S3_S2 <- dba.contrast(dbObj_S3_S2,reorderMeta=list(Condition="S2"),minMembers = 2)
> dbObj_S3_S2 <- dba.analyze(dbObj_S3_S2, method=DBA_ALL_METHODS)
> comp1.deseq2_S3_S2 <- dba.report(dbObj_S3_S2, method=DBA_DESEQ2,contrast=1,th=1)
• 11 views
Read more here: Source link