Skip to content

get_mapping_scores() IndexError #172

@SilviaZa

Description

@SilviaZa

Hi

I am comparing 2 species following the vignette.
Everything runs smoothely until I try to compute mapping scores with
D,MappingTable = get_mapping_scores(sm,keys,n_top = 0).
I consistently get this error:


IndexError Traceback (most recent call last)
Cell In[91], line 2
1 keys = {'nv':'ID.separate','sp':'cell_type'}
----> 2 D,MappingTable = get_mapping_scores(sm,keys,n_top = 0)

File c:\Users\z.SAMap_new_venv\Lib\site-packages\samap\analysis\scores.py:138, in get_mapping_scores(sm, keys, n_top)
135 label = "{}_mapping_scores".format(";".join([keys[sid] for sid in skeys]))
136 samap.adata.obs[label] = pd.Categorical(cl)
--> 138 CSIMth, clu = _compute_csim(samap, label, n_top=n_top, prepend=False)
140 A = pd.DataFrame(data=CSIMth, index=clu, columns=clu)
141 i = np.argsort(-A.values.max(0).flatten())

File c:\Users\z.SAMap_new_venv\Lib\site-packages\samap\analysis\scores.py:67, in _compute_csim(samap, key, X, prepend, n_top)
64 spyi = splabels[yi]
66 filt = spxi != spyi
---> 67 di = X.data[filt]
68 xi = xi[filt]
69 yi = yi[filt]

IndexError: boolean index did not match indexed array along axis 0; size of axis is 1348105 but size of corresponding boolean axis is 1347079


I would appreciate any help!

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions