Change location of 'fragments' in multiome Seurat file
Entering edit mode
2.6 years ago
V ▴ 360


I have a multiome dataset of RNA +ATAC seq that was given to me as an analysed .RData file by a collaborator. This is already integrated and fully analysed.

When I attempt to visualise something in the multiome Seurat object using something like the below :

CoveragePlot(test, region = 'CD8A', features = 'CD8A', assay = 'ATAC', expression.assay = 'SCT', peaks = FALSE)

I get an error saying :

Error: TabixFile: file(s) do not exist:

I am not 'm' so it would seem that Rstudio is looking for these files in the original directory on my collaborators machine rather on mine in the pwc. Does anyone know how I can edit the stored path to the fragments.tsv.gz file and replace it with my own? I've tried the 'SetFragments' argument that use to be active in Signac but it is no longer working. Any help would be much appreciated. Thank you

RNA-Seq R Seurat Signac ATAC-Seq • 2.4k views
Entering edit mode
2.1 years ago
V ▴ 360

To add an answer to this because I figured it out in the end :

  1. Reset the current fragments object in the Seurat
    1. Make a new path
  2. Add to your object
Fragments(HF1clean@assays$ATAC) <- NULL

 fragments <- CreateFragmentObject(path = "FullMultiome/atac_fragments.tsv.gz", cells = colnames(HF1), validate.fragments = TRUE)

Fragments(HF1clean@assays$ATAC) <- fragments

Login before adding your answer.

Traffic: 2577 users visited in the last hour
Help About
Access RSS

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6