Closed TingTingShao closed 6 months ago
This is an expected behavior. datas.to_adata()
will not copy those in datas.adatas.obsm
. In order to copy fragment_paired
, you need to do this before calling to_adata()
:
datas.obsm['fragment_paired'] = datas.adatas.obsm['fragment_paired']
Thanks!
Hi,
I used simply run this to convert anndataset into anndata
dataset=datas.to_adata()
However, when I call peaks with snap.tl.macs3, error shows:
neither 'fragment_single' nor 'fragment_paired' is present in the '.obsm'
Looking forward to your reply! thanks tingting