Thanks to visit codestin.com
Credit goes to github.com

Skip to content

Add fragment file path to a merged mudata object #133

@jolvhull

Description

@jolvhull

Dear muon-team,

I have converted my scATAC-seq SeuratObject to a mudata object but the information on the fragments files does not seem to be transferred. mdata.uns is empty.
My object is a merged object of multiple samples that each have a separate fragment file. I tried to run muon.atac.tl.locate_fragments by doing this for each subset in the object:
""""
samples = ['sample1', 'sample2', 'sample3']
fragment_dict = {
'sample1': os.path.join(work_dir, 'data/sample1_atac_fragments.tsv.gz),
sample2': os.path.join(work_dir, 'data/sample2_atac_fragments.tsv.gz),
sample3': os.path.join(work_dir, 'data/sample3_atac_fragments.tsv.gz)}

for x in samples:
muon.atac.tl.locate_fragments(mdata[mdata.obs.sample == x], fragments_dict[x])
""""

However, when I then try to continue with my QC, I still get the error "There is no fragments file located yet. Run muon.atac.tl.locate_fragments first."

When I subset my data, run muon.atac.tl.locate_fragments(sample1, fragments_dict['sample1'])
I am able to continue with QC (nucleosome signal)

How can I add the fragment file location for my complete dataset, where each subset, based on the value in an obs column, has its own fragment file?

Thank you!

Metadata

Metadata

Assignees

No one assigned

    Labels

    enhancementNew feature or request

    Type

    No type

    Projects

    No projects

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions