I am new to DiffBind. I have tried working with your test data and I have seen the tutorial online.
I am facing some issues with DiffBind.
- I am working on Rat data. when I enter my data, 'no genome detected' option comes up. I have seen that dba.blacklist does not have rat genome in background. How can I provide rat genome?
- I am not getting any peaks using DEseq2 and EdgeR in my data.
- Can we compare 4 different conditions in DiffBind? Is there an example available with you with more than 'Resistance' and 'Restrained' conditions?
- For MACS2, do we need to provide specifically 'broad' peak or 'narrow' peak data as input for DiffBind?
Code should be placed in three backticks as shown below
```dba.analyze
test_data.analyzed <- dba.analyze(test_data.counted)
Output: Applying Blacklist/Greylists... No genome detected. Normalize DESeq2 with defaults... Analyzing... gene-wise dispersion estimates mean-dispersion relationship final dispersion estimates
I have also tried:
test_data.analyzed <- dba.analyze(test_data.counted, method=DBA_ALL_METHODS, bBlacklist = FALSE, bGreylist = FALSE)
Output: Normalize edgeR with defaults... Normalize DESeq2 with defaults... Analyzing... gene-wise dispersion estimates mean-dispersion relationship final dispersion estimates
For both the commands, I do not get any differential peaks.
dba.show(test_data.analyzed, bContrasts = T)
Factor Group Samples Group2 Samples2 DB.edgeR DB.DESeq2 1 Condition TREATMENT 2 CONTROL 2 0 0
Kindly help.