Open mfiers opened 4 years ago
Weird, this sure looks like a bug! Can you give use code to reproduce it with the builtin datasets?
previously filed on anndata: theislab/anndata#234
Not sure if still required (the code in the anndata repo is close to this), but here is a minimal example:
import scanpy as sc
adata = sc.datasets.pbmc3k()
h5adfile = 'pbmc3k.h5ad'
adata.write(h5adfile)
a1 = sc.read_h5ad(h5adfile)
a2 = sc.read_h5ad(h5adfile, backed='r+')
sc.tl.score_genes(a1, ['KIR3DL2-1', 'AL590523.1', 'CT476828.1'])
# OK
sc.tl.score_genes(a2, ['KIR3DL2-1', 'AL590523.1', 'CT476828.1'])
# ERROR
thanks..
Has anyone solved it?
Hi - excellent software, thanks! -
but I do have a problem. If i load a disk backed dataset, I cannot run
sc.tl.score_genes
.Given these two sets:
and
this works perfectly:
but, this:
yields the following error:
thanks Mark