Why can't I quantify the peaks (scATAC)?
0
0
Entering edit mode
2.1 years ago

I generated /outs files following Cellranger 10X and for the fragment file did:

>  create fragment objects
frags.104 <- CreateFragmentObject(
  path = "fragments_101.tsv.gz",
  cells = rownames(md.101))


> pbmc101.counts <- FeatureMatrix(
  fragments = frags.101,
   features = combined.peaks,
cells = rownames(md.101)

**Error in h(simpleError(msg, call)) : 
error in evaluating the argument 'y' in selecting a method for function 'intersect': [internal] bgzf_seek() failed**

What am I doing wrong here?? I've been trying for days and sincerely do not want to re-generate a new set of files from Cellranger due to time.

Thanks.

Signac ATAC snATAC scATAC • 876 views
ADD COMMENT
0
Entering edit mode

You define the object frag.104 with a 4; but you reference object frag.101 with a 1; is this intended?

ADD REPLY
0
Entering edit mode

Yes. It was just mistyped but the error is "Error in h(simpleError(msg, call)) : error in evaluating the argument 'y' in selecting a method for function 'intersect': [internal] bgzf_seek() failed"

ADD REPLY

Login before adding your answer.

Traffic: 1088 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