SpatialDM icon indicating copy to clipboard operation
SpatialDM copied to clipboard

Question about global permutation test

Open sealx017 opened this issue 1 year ago • 0 comments

Hello, I noticed that the package compares the actual global I with the permuted global I's to compute the permutation p-value:

    adata.uns['global_stat']['perm']['global_p'] = 1 - (adata.uns['global_I'] \
                         > adata.uns['global_stat']['perm']['global_perm'].T).sum(axis=0) / n_perm

My question is why do you compare only one direction (">"), i.e., greater than the mean under NULL? Global I can be negative, i.e., less than the mean under NULL indicating negative spatial association.

Please let me know if I misunderstood the package and the paper!

sealx017 avatar Aug 22 '24 15:08 sealx017