GreenleafLab / ArchR

ArchR : Analysis of Regulatory Chromatin in R (www.ArchRProject.com)
MIT License
388 stars 140 forks source link

Error with addCombinedDims function: Error in f(init, x[[i]]) : number of rows of matrices must match (see arg 2) #1779

Closed millerkaplanlab closed 1 year ago

millerkaplanlab commented 1 year ago

Error with addCombinedDims function: Error in f(init, x[[i]]) : number of rows of matrices must match (see arg 2) #1779

Attach your log file

ArchR-addIterativeLSI-6debc678dfc2f-Date-2022-12-20_Time-13-06-44.log

Describe the bug Error with addCombinedDims function: Error in f(init, x[[i]]) : number of rows of matrices must match (see arg 2)

To Reproduce

Expected behavior Expected dimensions to be added to project. They are of equal size.

Screenshots

Session Info

sessionInfo() R version 4.2.1 (2022-06-23) Platform: x86_64-pc-linux-gnu (64-bit) Running under: Ubuntu 20.04.5 LTS

Matrix products: default BLAS: /usr/lib/x86_64-linux-gnu/blas/libblas.so.3.9.0 LAPACK: /usr/lib/x86_64-linux-gnu/lapack/liblapack.so.3.9.0

Random number generation: RNG: L'Ecuyer-CMRG Normal: Inversion Sample: Rejection

locale: [1] LC_CTYPE=en_CA.UTF-8 LC_NUMERIC=C LC_TIME=en_CA.UTF-8 LC_COLLATE=en_CA.UTF-8 LC_MONETARY=en_CA.UTF-8 LC_MESSAGES=en_CA.UTF-8 LC_PAPER=en_CA.UTF-8
[8] LC_NAME=C LC_ADDRESS=C LC_TELEPHONE=C LC_MEASUREMENT=en_CA.UTF-8 LC_IDENTIFICATION=C

attached base packages: [1] parallel grid stats4 stats graphics grDevices utils datasets methods base

other attached packages: [1] nabor_0.5.0 Rsamtools_2.12.0 BSgenome.Mmusculus.UCSC.mm10_1.4.3 BSgenome_1.64.0 rtracklayer_1.56.1
[6] Biostrings_2.64.1 XVector_0.36.0 rhdf5_2.40.0 stringr_1.5.0 plyr_1.8.8
[11] magrittr_2.0.3 gtable_0.3.1 gtools_3.9.4 gridExtra_2.3 ArchR_1.0.2
[16] cluster_2.1.4 presto_1.0.0 data.table_1.14.6 Rcpp_1.0.9 colorspace_2.0-3
[21] scClustViz_1.3.11 shiny_1.7.3 SeuratObject_4.1.3 Seurat_4.3.0 scran_1.24.1
[26] scater_1.24.0 ggplot2_3.4.0 scuttle_1.6.3 DropletUtils_1.16.0 SingleCellExperiment_1.18.1
[31] SummarizedExperiment_1.26.1 MatrixGenerics_1.8.1 matrixStats_0.63.0 Matrix_1.5-1 org.Dr.eg.db_3.15.0
[36] biovizBase_1.44.0 Signac_1.9.0 ensembldb_2.20.2 AnnotationFilter_1.20.0 GenomicFeatures_1.48.4
[41] AnnotationDbi_1.58.0 Biobase_2.56.0 GenomicRanges_1.48.0 GenomeInfoDb_1.32.4 IRanges_2.32.0
[46] S4Vectors_0.36.1 AnnotationHub_3.4.0 BiocFileCache_2.4.0 dbplyr_2.2.1 BiocGenerics_0.44.0

loaded via a namespace (and not attached): [1] rsvd_1.0.5 Hmisc_4.7-2 ica_1.0-3 RcppRoll_0.3.0 ps_1.7.2 lmtest_0.9-40
[7] crayon_1.5.2 MASS_7.3-58 rhdf5filters_1.8.0 nlme_3.1-159 backports_1.4.1 rlang_1.0.6
[13] ROCR_1.0-11 fontawesome_0.4.0 irlba_2.3.5.1 callr_3.7.3 limma_3.54.0 filelock_1.0.2
[19] BiocParallel_1.30.4 rjson_0.2.21 bit64_4.0.5 harmony_0.1.1 glue_1.6.2 sctransform_0.3.5
[25] processx_3.8.0 vipor_0.4.5 spatstat.sparse_3.0-0 spatstat.geom_3.0-3 tidyselect_1.2.0 usethis_2.1.6
[31] fitdistrplus_1.1-8 XML_3.99-0.13 tidyr_1.2.1 zoo_1.8-11 GenomicAlignments_1.32.1 xtable_1.8-4
[37] cli_3.4.1 zlibbioc_1.42.0 rstudioapi_0.14 miniUI_0.1.1.1 sp_1.5-1 rhandsontable_0.3.8
[43] bslib_0.4.1 rpart_4.1.16 fastmatch_1.1-3 BiocSingular_1.12.0 xfun_0.35 pkgbuild_1.4.0
[49] KEGGREST_1.36.3 tibble_3.1.8 interactiveDisplayBase_1.34.0 ggrepel_0.9.2 listenv_0.8.0 png_0.1-8
[55] future_1.29.0 withr_2.5.0 bitops_1.0-7 dqrng_0.3.0 pillar_1.8.1 cachem_1.0.6
[61] fs_1.5.2 DelayedMatrixStats_1.18.2 vctrs_0.5.1 ellipsis_0.3.2 generics_0.1.3 devtools_2.4.5
[67] tools_4.2.1 foreign_0.8-82 beeswarm_0.4.0 munsell_0.5.0 DelayedArray_0.22.0 fastmap_1.1.0
[73] compiler_4.2.1 pkgload_1.3.2 abind_1.4-5 httpuv_1.6.6 sessioninfo_1.2.2 plotly_4.10.1
[79] GenomeInfoDbData_1.2.8 edgeR_3.38.4 lattice_0.20-45 deldir_1.0-6 utf8_1.2.2 later_1.3.0
[85] dplyr_1.0.10 jsonlite_1.8.4 scales_1.2.1 ScaledMatrix_1.4.1 pbapply_1.6-0 sparseMatrixStats_1.8.0
[91] lazyeval_0.2.2 promises_1.2.0.1 latticeExtra_0.6-30 R.utils_2.12.2 goftest_1.2-3 spatstat.utils_3.0-1
[97] reticulate_1.26 checkmate_2.1.0 cowplot_1.1.1 textshaping_0.3.6 statmod_1.4.37 Rtsne_0.16
[103] dichromat_2.0-0.1 uwot_0.1.14 igraph_1.3.5 HDF5Array_1.24.2 survival_3.4-0 rsconnect_0.8.28
[109] yaml_2.3.6 systemfonts_1.0.4 htmltools_0.5.4 memoise_2.0.1 VariantAnnotation_1.42.1 profvis_0.3.7
[115] BiocIO_1.6.0 locfit_1.5-9.6 viridisLite_0.4.1 digest_0.6.31 assertthat_0.2.1 mime_0.12
[121] commonmark_1.8.1 rappdirs_0.3.3 RSQLite_2.2.19 future.apply_1.10.0 remotes_2.4.2 urlchecker_1.0.1
[127] blob_1.2.3 R.oo_1.25.0 ragg_1.2.4 shinythemes_1.2.0 splines_4.2.1 Formula_1.2-4
[133] labeling_0.4.2 Rhdf5lib_1.18.2 Cairo_1.6-0 ProtGenerics_1.28.0 RCurl_1.98-1.9 hms_1.1.2
[139] base64enc_0.1-3 BiocManager_1.30.19 ggbeeswarm_0.6.0 nnet_7.3-17 sass_0.4.4 RANN_2.6.1
[145] fansi_1.0.3 parallelly_1.32.1 R6_2.5.1 ggridges_0.5.4 lifecycle_1.0.3 bluster_1.6.0
[151] curl_4.3.3 leiden_0.4.3 jquerylib_0.1.4 RcppAnnoy_0.0.20 RColorBrewer_1.1-3 spatstat.explore_3.0-5
[157] htmlwidgets_1.5.4 beachmat_2.12.0 polyclip_1.10-4 biomaRt_2.52.0 markdown_1.4 purrr_0.3.5
[163] globals_0.16.2 htmlTable_2.4.1 patchwork_1.1.2 spatstat.random_3.0-1 progressr_0.12.0 codetools_0.2-18
[169] metapod_1.4.0 prettyunits_1.1.1 R.methodsS3_1.8.2 DBI_1.1.3 tensor_1.5 httr_1.4.4
[175] KernSmooth_2.23-20 stringi_1.7.8 progress_1.2.2 reshape2_1.4.4 farver_2.1.1 viridis_0.6.2
[181] xml2_1.3.3 BiocNeighbors_1.14.0 restfulr_0.0.15 interp_1.1-3 scattermore_0.8 BiocVersion_3.15.2
[187] bit_4.0.5 jpeg_0.1-10 spatstat.data_3.0-0 pkgconfig_2.0.3 knitr_1.41

Additional context

rcorces commented 1 year ago

Hi @millerkaplanlab! Thanks for using ArchR! Please make sure that your post belongs in the Issues section. Only bugs and error reports belong in the Issues section. Usage questions and feature requests should be posted in the Discussions section, not in Issues.
Before we help you, you must respond to the following questions unless your original post already contained this information: 1. If you've encountered an error, have you already searched previous Issues to make sure that this hasn't already been solved? 2. Can you recapitulate your error using the tutorial code and dataset? If so, provide a reproducible example. 3. Did you post your log file? If not, add it now. 4. Remove any screenshots that contain text and instead copy and paste the text using markdown's codeblock syntax (three consecutive backticks). You can do this by editing your original post.

millerkaplanlab commented 1 year ago

1. If you've encountered an error, have you already searched previous Issues to make sure that this hasn't already been solved? Yes 2. Can you recapitulate your error using the tutorial code and dataset? If so, provide a reproducible example. No 3. Did you post your log file? If not, add it now. Yes 4. Remove any screenshots that contain text and instead copy and paste the text using markdown's codeblock syntax (three consecutive backticks). You can do this by editing your original post.

rcorces commented 1 year ago

The log file that you uploaded in your original post is for addIterativeLSI rather than addCombinedDims

Since this isnt a problem with the tutorial data, its likely some incompatibility of your datasets or something you've done prior. In order to help, I'll need the correct log file and the code used from start to finish.

millerkaplanlab commented 1 year ago

I don't know exactly how I fixed the issue; however, it seems to be working now!

ZayaAyush commented 1 year ago

UMAP-scATAC-scRNA.pdf

I don't know exactly how I fixed the issue; however, it seems to be working now!

I have a same issue! Could you let me know how to fix this issue! projMulti3 <- addCombinedDims(projMulti3, reducedDims = c("LSI_ATAC", "LSI_RNA"), name = "LSI_Combined")

Error in f(init, x[[i]]) : number of rows of matrices must match (see arg 2)

I can see individual UMAP of LSI_ATAC and LSI_RNA. Please help me! thank you