-
Notifications
You must be signed in to change notification settings - Fork 4
New issue
Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.
By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.
Already on GitHub? Sign in to your account
How to use replicates? #11
Comments
Hi Ivan, in case each condition has multiple Hi-C replicates, DiffDomain takes the combined Hi-C contact matrix from the replicates as the input which is a common practice to generate a large number of Hi-C interactions. |
Hi Tian, thank you for your answer. The problem is that I was explicitly asked to take into account replicates in order to have statistics about the different tads, like in RNA-seq, therefore merging is not an option. I was wondering if you are planning something in this sense, creating tad by tad a sort of distribution and then testing it case vs control. I know it's difficult for intervals, such as for atac-seq, but for example diffbind implements a similar approach. |
Hey Ivan. Current version of DiffDomain cannot take the replicates as the input. We are sorry about this. Indeed, it is challenging for comparing TADs using replicates Hi-C data. We will take this as a future work. Please stay tuned. Meanwhile, I would recommend trying DiffDomain on combined Hi-C data. From our extensive analysis, using only combined Hi-C data achieves reasonable results (see our bioRxiv paper here) |
Hi, I have replicates for the two conditions to test. How I can I use them? I don't see option to use replicates.
The text was updated successfully, but these errors were encountered: