Skip to content

Question about global permutation test #31

@sealx017

Description

@sealx017

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!

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type
    No fields configured for issues without a type.

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions