r3fang / SnapATAC

Analysis Pipeline for Single Cell ATAC-seq
GNU General Public License v3.0
307 stars 126 forks source link

Error in `createSnap` #144

Open galib36 opened 4 years ago

galib36 commented 4 years ago

Hello, When I try to create a snap object with the following command

x.sp = createSnap(
    file=c("RJ_ATACseq_from_cellranger.snap"),
    sample=c("Pancreas_190724"),
    num.cores=1
    )

it gives me the following error

Error in data.frame(barcode, TN, UM, PP, UQ, CM) : arguments imply differing number of rows: 0, 8241
9.
stop(gettextf("arguments imply differing number of rows: %s", paste(unique(nrows), collapse = ", ")), domain = NA)
8.
data.frame(barcode, TN, UM, PP, UQ, CM)
7.
readMetaData.default(file)
6.
readMetaData(file)
5.
createSnapSingle(file = fileList[[i]], sample = sampleList[[i]])
4.
FUN(X[[i]], ...)
3.
lapply(as.list(seq(fileList)), function(i) { createSnapSingle(file = fileList[[i]], sample = sampleList[[i]]) })
2.
createSnap.default(file = c("RJ_ATACseq_from_cellranger.snap"), sample = c("Pancreas_190724"), num.cores = 1)
1.
createSnap(file = c("RJ_ATACseq_from_cellranger.snap"), sample = c("Pancreas_190724"), num.cores = 1)

In addition, in the console it is showing the following error

Epoch: reading the barcode session ...
[1] "Warning @readSnap: 'BD/name' not found in  RJ_ATACseq_from_cellranger.snap"
An open HDF5 file handle exists. If the file has changed on disk meanwhile, the function may not work properly. Run 'h5closeAll()' to close all open HDF5 object handles.An open HDF5 file handle exists. If the file has changed on disk meanwhile, the function may not work properly. Run 'h5closeAll()' to close all open HDF5 object handles.

I ran the h5closeAll() command but the problem continues. I looked into #78 and updated rhdf5 to rhdf5_2.28.1 but the problem continues. Can you please let me know how to solve the issue.

xiangrong7 commented 4 years ago

Hello, When I try to create a snap object with the following command

x.sp = createSnap(
  file=c("RJ_ATACseq_from_cellranger.snap"),
  sample=c("Pancreas_190724"),
  num.cores=1
  )

it gives me the following error

Error in data.frame(barcode, TN, UM, PP, UQ, CM) : arguments imply differing number of rows: 0, 8241
9.
stop(gettextf("arguments imply differing number of rows: %s", paste(unique(nrows), collapse = ", ")), domain = NA)
8.
data.frame(barcode, TN, UM, PP, UQ, CM)
7.
readMetaData.default(file)
6.
readMetaData(file)
5.
createSnapSingle(file = fileList[[i]], sample = sampleList[[i]])
4.
FUN(X[[i]], ...)
3.
lapply(as.list(seq(fileList)), function(i) { createSnapSingle(file = fileList[[i]], sample = sampleList[[i]]) })
2.
createSnap.default(file = c("RJ_ATACseq_from_cellranger.snap"), sample = c("Pancreas_190724"), num.cores = 1)
1.
createSnap(file = c("RJ_ATACseq_from_cellranger.snap"), sample = c("Pancreas_190724"), num.cores = 1)

In addition, in the console it is showing the following error

Epoch: reading the barcode session ...
[1] "Warning @readSnap: 'BD/name' not found in  RJ_ATACseq_from_cellranger.snap"
An open HDF5 file handle exists. If the file has changed on disk meanwhile, the function may not work properly. Run 'h5closeAll()' to close all open HDF5 object handles.An open HDF5 file handle exists. If the file has changed on disk meanwhile, the function may not work properly. Run 'h5closeAll()' to close all open HDF5 object handles.

I ran the h5closeAll() command but the problem continues. I looked into #78 and updated rhdf5 to rhdf5_2.28.1 but the problem continues. Can you please let me know how to solve the issue.

Me too~~~

xiangrong7 commented 4 years ago

updated rhdf5 to rhdf5_2.30.0 run ok~~

galib36 commented 4 years ago

I updated it to rhdf5_2.30.1 but the problem continues.

galib36 commented 4 years ago

I tried with rhd5_2.30.0 but still no luck, the problem continues. Any solution?