RStudio Diffbind analysis shows no result
0
0
Entering edit mode
10 months ago
a_bis ▴ 40

Hi all, I'm trying to carry out a differential accessibility analysis on ATAC-seq data using DiffBind. I'm reading in a .csv file as my sample sheet and getting a decent correlation heatmap per the instructions in section 3.1 of the DiffBind manual. However, when I carry out dba.count(), I get "Error: No sites have activity greater than filter value." The same dataset that gave me thousands of differentially accessible sites when I ran a quick analysis on Galaxy, but I'm getting nothing here. Would anyone have an idea of what I could modify in these early stages of the analysis pipeline to address this issue? Thanks in advance.

galaxy diffbind rstudio • 495 views
ADD COMMENT
0
Entering edit mode

Let me kindly suggest to google "Error: No sites have activity greater than filter value." and you will find previous answers on that.

ADD REPLY
0
Entering edit mode

Thank you for your kind suggestion. I have indeed googled this in an effort to resolve the question, and found what appears to be a useful post on bioconductor (https://support.bioconductor.org/p/9137594/). It suggests running dba.peakset(samples, bRetrieve=TRUE) to look at the matrix after reading in the samples and before counting, but I seem to be getting Error in rep(T, nrow(pv$binding)) : invalid 'times' argument. Unfortunately I can't identify what this times argument refers to as I don't know where the rep() comes in play.

ADD REPLY

Login before adding your answer.

Traffic: 2570 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6