Error when running Integrating scRNA-seq and scATAC-seq data vignette #1877
Unanswered
realmungmung
asked this question in
Q&A
Replies: 0 comments
Sign up for free
to join this conversation on GitHub.
Already have an account?
Sign in to comment
Uh oh!
There was an error while loading. Please reload this page.
-
Hi all,
I'm running this Integrating scRNA-seq and scATAC-seq data vignette and I was able to run through the tutorial until this command:
gene.activities <- GeneActivity(pbmc.atac, features = VariableFeatures(pbmc.rna))
where I got a huge error:
gene.activities <- GeneActivity(pbmc.atac, features = VariableFeatures(pbmc.rna))
Extracting gene coordinates
Extracting reads overlapping genomic regions
| | 0 % ~calculating [E::bgzf_read_block] Failed to read BGZF header at offset 1416828390
[E::hts_itr_next] Failed to seek to offset 92869838159767: Illegal seek
[E::hts_itr_next] Failed to seek to offset 92902743745887: Illegal seek
[E::hts_itr_next] Failed to seek to offset 92916919541141: Illegal seek
[E::hts_itr_next] Failed to seek to offset 93381738192310: Illegal seek
[E::hts_itr_next] Failed to seek to offset 93413213899239: Illegal seek
[E::hts_itr_next] Failed to seek to offset 93516938894158: Illegal seek
[E::hts_itr_next] Failed to seek to offset 93609721460164: Illegal seek
[E::hts_itr_next] Failed to seek to offset 93629593997030: Illegal seek
[E::hts_itr_next] Failed to seek to offset 93633988935576: Illegal seek
[E::hts_itr_next] Failed to seek to offset 93643503715391: Illegal seek
......
my session info is:
sessionInfo()
R version 4.4.2 (2024-10-31)
Platform: aarch64-apple-darwin20
Running under: macOS Sonoma 14.6.1
Matrix products: default
BLAS: /System/Library/Frameworks/Accelerate.framework/Versions/A/Frameworks/vecLib.framework/Versions/A/libBLAS.dylib
LAPACK: /Library/Frameworks/R.framework/Versions/4.4-arm64/Resources/lib/libRlapack.dylib; LAPACK version 3.12.0
locale:
[1] en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8
time zone: America/Los_Angeles
tzcode source: internal
attached base packages:
[1] stats4 stats graphics grDevices utils datasets methods base
other attached packages:
[1] cowplot_1.1.3 ggplot2_3.5.1 EnsDb.Hsapiens.v86_2.99.0 ensembldb_2.30.0 AnnotationFilter_1.30.0
[6] GenomicFeatures_1.58.0 AnnotationDbi_1.68.0 Biobase_2.66.0 GenomicRanges_1.58.0 GenomeInfoDb_1.42.1
[11] IRanges_2.40.1 S4Vectors_0.44.0 BiocGenerics_0.52.0 Signac_1.14.0 Seurat_5.1.0
[16] SeuratObject_5.0.2 sp_2.1-4 pbmcMultiome.SeuratData_0.1.4 SeuratData_0.2.2.9001
loaded via a namespace (and not attached):
[1] RcppAnnoy_0.0.22 splines_4.4.2 later_1.4.1 BiocIO_1.16.0 bitops_1.0-9 tibble_3.2.1
[7] polyclip_1.10-7 rpart_4.1.23 XML_3.99-0.18 fastDummies_1.7.4 lifecycle_1.0.4 globals_0.16.3
[13] lattice_0.22-6 MASS_7.3-63 backports_1.5.0 magrittr_2.0.3 rmarkdown_2.29 Hmisc_5.2-1
[19] plotly_4.10.4 yaml_2.3.10 httpuv_1.6.15 sctransform_0.4.1 spam_2.11-0 spatstat.sparse_3.1-0
[25] reticulate_1.40.0 pbapply_1.7-2 DBI_1.2.3 RColorBrewer_1.1-3 abind_1.4-8 zlibbioc_1.52.0
[31] Rtsne_0.17 purrr_1.0.2 biovizBase_1.54.0 RCurl_1.98-1.16 nnet_7.3-20 VariantAnnotation_1.52.0
[37] rappdirs_0.3.3 GenomeInfoDbData_1.2.13 ggrepel_0.9.6 irlba_2.3.5.1 listenv_0.9.1 spatstat.utils_3.1-1
[43] goftest_1.2-3 RSpectra_0.16-2 spatstat.random_3.3-2 fitdistrplus_1.2-1 parallelly_1.41.0 leiden_0.4.3.1
[49] codetools_0.2-20 DelayedArray_0.32.0 RcppRoll_0.3.1 tidyselect_1.2.1 UCSC.utils_1.2.0 farver_2.1.2
[55] base64enc_0.1-3 matrixStats_1.4.1 spatstat.explore_3.3-3 GenomicAlignments_1.42.0 jsonlite_1.8.9 Formula_1.2-5
[61] progressr_0.15.1 ggridges_0.5.6 survival_3.8-3 tools_4.4.2 ica_1.0-3 Rcpp_1.0.13-1
[67] glue_1.8.0 gridExtra_2.3 SparseArray_1.6.0 xfun_0.49 MatrixGenerics_1.18.0 dplyr_1.1.4
[73] withr_3.0.2 fastmap_1.2.0 digest_0.6.37 R6_2.5.1 mime_0.12 colorspace_2.1-1
[79] scattermore_1.2 tensor_1.5 dichromat_2.0-0.1 spatstat.data_3.1-4 RSQLite_2.3.9 tidyr_1.3.1
[85] generics_0.1.3 data.table_1.16.4 rtracklayer_1.66.0 httr_1.4.7 htmlwidgets_1.6.4 S4Arrays_1.6.0
[91] uwot_0.2.2 pkgconfig_2.0.3 gtable_0.3.6 blob_1.2.4 lmtest_0.9-40 XVector_0.46.0
[97] htmltools_0.5.8.1 dotCall64_1.2 ProtGenerics_1.38.0 scales_1.3.0 png_0.1-8 spatstat.univar_3.1-1
[103] knitr_1.49 rstudioapi_0.17.1 reshape2_1.4.4 rjson_0.2.23 checkmate_2.3.2 nlme_3.1-166
[109] curl_6.0.1 zoo_1.8-12 cachem_1.1.0 stringr_1.5.1 KernSmooth_2.23-26 parallel_4.4.2
[115] miniUI_0.1.1.1 foreign_0.8-87 restfulr_0.0.15 pillar_1.10.0 grid_4.4.2 vctrs_0.6.5
[121] RANN_2.6.2 promises_1.3.2 xtable_1.8-4 cluster_2.1.8 htmlTable_2.4.3 evaluate_1.0.1
[127] cli_3.6.3 compiler_4.4.2 Rsamtools_2.22.0 rlang_1.1.4 crayon_1.5.3 future.apply_1.11.3
[133] labeling_0.4.3 plyr_1.8.9 stringi_1.8.4 viridisLite_0.4.2 deldir_2.0-4 BiocParallel_1.40.0
[139] munsell_0.5.1 Biostrings_2.74.1 lazyeval_0.2.2 spatstat.geom_3.3-4 Matrix_1.7-1 BSgenome_1.74.0
[145] RcppHNSW_0.6.0 patchwork_1.3.0 bit64_4.5.2 future_1.34.0 KEGGREST_1.46.0 shiny_1.10.0
[151] SummarizedExperiment_1.36.0 ROCR_1.0-11 igraph_2.1.2 memoise_2.0.1 fastmatch_1.1-6 bit_4.5.0.1
Please let me know if you have any suggestions thank you so much!
Beta Was this translation helpful? Give feedback.
All reactions