GM and PM not found by snapAtac

[rLannes]

Hi,

I have this strange behaviour, I started by using bam of 10X data.

I added PM and GM using snaptools

snaptools snap-add-pmat --snap-file T20276_atac.header.snap --peak-file ../peaks.bed
error: cell x peak matrix already exists, delete it using snap-del first

 snaptools snap-add-gmat --snap-file T20276_atac.header.snap --gene-file ../genes_chr.bed 
error: cell x gene matrix already exists, delete it first using snap-del first

But in R I encounter the following Error:

snap = createSnap( file="T20276_atac.header.snap", sample="T20276" );

snap = createSnap(
+   file="Projet/atac_snap/snap_atac/T20276_atac.header.snap",
+   sample="T20276"
+ );
Epoch: reading the barcode session ...
> x.sp = addBmatToSnap(snap)
Epoch: reading cell-bin count matrix session ...
> x.sp = addPmatToSnap(x.sp)
Epoch: reading cell-peak count matrix session ...
Error in value[[3L]](cond) : 
  Warning @readSnap: 'PM/idx' not found in  T20276_atac.header.snap
> x.sp = addGmatToSnap(x.sp)
Epoch: reading cell-gene count matrix session ...
Error in value[[3L]](cond) : 
  Warning @addGmat: 'GM/idx' not found in  T20276_atac.header.snap

Any help appreciated!

[Yydtqqqg]

I encountered the same issue:

Error in value[[3L]](cond) : 
  Warning @addGmat: 'GM/idx' not found in ...... .snap
Calls: addGmatToSnap ... tryCatch -> tryCatchList -> tryCatchOne -> <Anonymous>
停止执行

Are you able to solve the problem?