Open stevevanhooser opened 4 years ago
@dlwise : I added at_colocalization_summary. Can you check it to see if it provides the outputs you are seeking? The data I got only had 2-way colocalization.
It produces (for me):
[out,msg] = at_colocalization_summary(the_atdir, colname);
out
out =
struct with fields:
totalA: 89
totalB: 54
NumAPaired_ab: 57
NumBPaired_ab: 45
NumANotPaired_ab: 32
NumBNotPaired_ab: 9
NumAPaired_ba: 57
NumBPaired_ba: 45
NumANotPaired_ba: 32
NumBNotPaired_ba: 9
msg =
'AoB:Number of As paired with at least one B: 57; AoB:Number of Bs paired with at least one A: 45; AoB:Number of As not paired with at least one B: 32; AoB:Number of Bs not paired with at least one A: 9; BoA:Number of As paired with at least one B: 57; BoA:Number of Bs paired with at least one A: 45; BoA:Number of As not paired with at least one B: 32; BoA:Number of Bs not paired with at least one A: 9; '
Should report #paired and #unpaired, especially useful for 3-d colocalization