cole-trapnell-lab / cicero-release

https://cole-trapnell-lab.github.io/cicero-release/
MIT License
55 stars 15 forks source link

A question about generate_ccans() results #61

Closed chent5 closed 4 years ago

chent5 commented 4 years ago

Thanks a lot for developing this package.

I checked the pairwise co-accessibility of a CCAN detected by generate_ccans() function. I found that there are also some negative values of co-accessibility. Here are my code:

image

I am a little confused... CCANs are clusters of sites that tend to be co-accessible. How to explain these negative values for the peaks in the same CCAN? Or is there anything wrong with my code or my understanding?

Thanks so much for your help in advance!

hpliner commented 4 years ago

Hello, CCANs are clusters of coaccessible sites as determined by finding modules within a graph where edges have positive coaccessibility. This doesn't exclude the possibility that there are pairs of sites within that module that have negative coaccessibility values. In the above code you are sorting by coaccessibility, I suspect there are many more positively coaccessible pairs further down the list!