Change location of 'fragments' in multiome Seurat file
1
0
Entering edit mode
3.1 years ago
V ▴ 380

Hello,

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:
  '/Users/m/Data/Multiome_1220/atac_fragments.tsv.gz'
  '/Users/m/Data/Multiome_1220/atac_fragments.tsv.gz.tbi'

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 • 3.1k views
ADD COMMENT
4
Entering edit mode
2.6 years ago
V ▴ 380

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
ADD COMMENT

Login before adding your answer.

Traffic: 1897 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

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

Powered by the version 2.3.6