SnapATAC2 icon indicating copy to clipboard operation
SnapATAC2 copied to clipboard

Unexpected behavior in AnnDataSet subset

Open mtvector opened this issue 8 months ago • 0 comments

Hi there,

When I run code like this:

data = snap.read_dataset("concatenated.h5ads",mode='r+')
data.obs['doubcall']=list(np.random.choice([True,False],data.shape[0]))
data.subset(obs_indices=pd.Series(data.obs_names)[~data.obs['doubcall']],out="concatenated_nodoublet")

It creates a directory concatenated_nodoublet, but this contains no files, while all of the h5ads within the AnnDataSet are copied in the locations where they are pointed to, with the suffix .h5ad.h5ad. I'm also wondering if this function is supposed write a new AnnDataSet file so that the datasets in this new directory can be easily accessed? Or does it inplace modify the original h5ads file?

Thanks! Matthew

mtvector avatar Jun 24 '24 04:06 mtvector