Entering edit mode
2.1 years ago
bioinformatics.girl
▴
20
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.
You define the object
frag.104
with a 4; but you reference objectfrag.101
with a 1; is this intended?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"