abagen icon indicating copy to clipboard operation
abagen copied to clipboard

Qustion about the hemisphere filter method in drop_mismatch_samples function (samples_.py L223)

Open YCHuang0610 opened this issue 1 year ago • 0 comments

https://github.com/rmarkello/abagen/blob/6f5542522247b9a7fb00ded198c34f4dc25957f3/abagen/samples_.py#L220C1-L224C46

Line 220-224 in samples_.py

annot = annotation.assign(hemisphere=hemisphere, structure=structure) \
                  .query('(hemisphere == "L" & mni_x < 0) '
                         '| (hemisphere == "R" & mni_x > 0) '
                         '| (hemisphere == "B" & mni_x == 0)',
                         engine='python')

Some of the samples in AHBA were not annotated hemisphere, such as 4700, Ve-I-II; 4701, Ve-III in donor 9861. In samples_.py#L223, the samples that are not annotated hemispheres and have a non-zero mni_x will be removed. This probably causes that all the samples without annotated hemispheres will be filtered out at this step, because no sample has a mni_x of 0. I am wondering is this resonable?

image

image

YCHuang0610 avatar Aug 08 '23 09:08 YCHuang0610