Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Error in adding RNA assay #178

Open
eblchen opened this issue Mar 15, 2024 · 1 comment
Open

Error in adding RNA assay #178

eblchen opened this issue Mar 15, 2024 · 1 comment

Comments

@eblchen
Copy link

eblchen commented Mar 15, 2024

Having difficulty converting h5Seurat to Anndata.
Had similar issues to #27 in the SaveH5Seurat step(), resolved by FindVariableFeatures() suggestion.
Now having issues with the Convert() step:

Convert("data/ut_seurat.h5Seurat",` dest = "h5ad", assay = "RNA", overwrite = T)
Validating h5Seurat file
Adding data from RNA as X
Error in assay.group$obj_copy_to(dst_loc = dfile, dst_name = "X", src_name = x.data) : 
  HDF5-API Errors:
    error #000: ../../src/H5Ocopy.c in H5Ocopy(): line 233: unable to copy object
        class: HDF5
        major: Object header
        minor: Unable to copy object

    error #001: ../../src/H5Ocopy.c in H5O__copy(): line 286: source object not found
        class: HDF5
        major: Symbol table
        minor: Object not found

    error #002: ../../src/H5Gloc.c in H5G_loc_find(): line 422: can't find object
        class: HDF5
        major: Symbol table
        minor: Object not found

    error #003: ../../src/H5Gtraverse.c in H5G_traverse(): line 851: internal path traversal failed
        class: HDF5
        major: Symbol table
        minor: Object not found

    error #004: ../../src/H5Gtraverse.c in H5G__traverse_real(): line 627: traversal operator failed
        class: HDF5
        major: Symbol table
        minor: Callback failed

    error #005: ../../src/H5Gloc.c in H5G__loc_find_cb(): line 378: object 'data' doesn't exist
        class: HD

Same error as quoted by Rohit-Satyam in #166. Have tried downgrading to Seurat and SeuratObject v4, also didn't solve the issue.

N.B. I can add the SCT assay, the error is only when I try to add the RNA assay and I need the uncorrected counts in the AnnData object.

sessionInfo()
R version 4.3.2 (2023-10-31)
Platform: x86_64-redhat-linux-gnu (64-bit)
Running under: Rocky Linux 8.9 (Green Obsidian)

Matrix products: default
BLAS/LAPACK: /usr/lib64/libopenblaso-r0.3.15.so;  LAPACK version 3.9.0

locale:
[1] C

time zone: UTC
tzcode source: system (glibc)

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

other attached packages:
 [1] SeuratDisk_0.0.0.9021       hues_0.2.0                  ggbeeswarm_0.7.2            scattermore_1.2             patchwork_1.2.0.9000       
 [6] slingshot_2.10.0            princurve_2.1.6             TSCAN_1.38.0                TrajectoryUtils_1.10.1      scater_1.30.1              
[11] ggplot2_3.5.0               scuttle_1.10.2              SingleCellExperiment_1.22.0 SummarizedExperiment_1.32.0 Biobase_2.62.0             
[16] GenomicRanges_1.54.1        GenomeInfoDb_1.38.6         IRanges_2.36.0              S4Vectors_0.40.2            BiocGenerics_0.48.1        
[21] MatrixGenerics_1.14.0       matrixStats_1.1.0           Seurat_5.0.1                SeuratObject_5.0.0          sp_2.1-1                   

loaded via a namespace (and not attached):
  [1] RcppAnnoy_0.0.21          splines_4.3.2             later_1.3.1               bitops_1.0-7              tibble_3.2.1             
  [6] polyclip_1.10-6           fastDummies_1.7.3         lifecycle_1.0.4           hdf5r_1.3.9               processx_3.8.2           
 [11] globals_0.16.2            lattice_0.21-9            pals_1.8                  MASS_7.3-60               magrittr_2.0.3           
 [16] plotly_4.10.4             remotes_2.4.2.1           httpuv_1.6.12             sctransform_0.4.1         spam_2.10-0              
 [21] pkgbuild_1.4.3            spatstat.sparse_3.0-3     reticulate_1.34.0         cowplot_1.1.1             mapproj_1.2.11           
 [26] pbapply_1.7-2             RColorBrewer_1.1-3        maps_3.4.1                abind_1.4-5               zlibbioc_1.48.0          
 [31] Rtsne_0.16                purrr_1.0.2               RCurl_1.98-1.14           GenomeInfoDbData_1.2.11   ggrepel_0.9.4            
 [36] irlba_2.3.5.1             listenv_0.9.0             spatstat.utils_3.0-4      goftest_1.2-3             RSpectra_0.16-1          
 [41] spatstat.random_3.2-1     fitdistrplus_1.1-11       parallelly_1.36.0         DelayedMatrixStats_1.22.6 leiden_0.4.3             
 [46] codetools_0.2-19          DelayedArray_0.28.0       tidyselect_1.2.0          ScaledMatrix_1.8.1        viridis_0.6.5            
 [51] spatstat.explore_3.2-5    jsonlite_1.8.7            BiocNeighbors_1.18.0      ellipsis_0.3.2            progressr_0.14.0         
 [56] ggridges_0.5.4            survival_3.5-8            tools_4.3.2               ica_1.0-3                 Rcpp_1.0.12              
 [61] glue_1.7.0                gridExtra_2.3             SparseArray_1.2.4         mgcv_1.9-0                dplyr_1.1.3              
 [66] withr_3.0.0               combinat_0.0-8            fastmap_1.1.1             fansi_1.0.6               callr_3.7.3              
 [71] caTools_1.18.2            digest_0.6.33             rsvd_1.0.5                R6_2.5.1                  mime_0.12                
 [76] colorspace_2.1-0          gtools_3.9.4              tensor_1.5                dichromat_2.0-0.1         spatstat.data_3.0-3      
 [81] utf8_1.2.4                tidyr_1.3.0               generics_0.1.3            data.table_1.14.8         httr_1.4.7               
 [86] htmlwidgets_1.6.2         S4Arrays_1.2.0            uwot_0.1.16               pkgconfig_2.0.3           gtable_0.3.4             
 [91] lmtest_0.9-40             XVector_0.42.0            htmltools_0.5.6.1         dotCall64_1.1-0           scales_1.3.0             
 [96] png_0.1-8                 rstudioapi_0.15.0         reshape2_1.4.4            curl_5.1.0                nlme_3.1-163             
[101] zoo_1.8-12                stringr_1.5.0             KernSmooth_2.23-22        parallel_4.3.2            miniUI_0.1.1.1           
[106] vipor_0.4.5               desc_1.4.3                pillar_1.9.0              grid_4.3.2                fastICA_1.2-3            
[111] vctrs_0.6.5               RANN_2.6.1                gplots_3.1.3              promises_1.2.1            BiocSingular_1.16.0      
[116] beachmat_2.16.0           xtable_1.8-4              cluster_2.1.4             beeswarm_0.4.0            cli_3.6.2                
[121] compiler_4.3.2            rlang_1.1.3               crayon_1.5.2              future.apply_1.11.0       mclust_6.0.0             
[126] ps_1.7.5                  plyr_1.8.9                stringi_1.7.12            viridisLite_0.4.2         deldir_1.0-9             
[131] BiocParallel_1.36.0       munsell_0.5.0             lazyeval_0.2.2            spatstat.geom_3.2-7       Matrix_1.6-1.1           
[136] RcppHNSW_0.5.0            bit64_4.0.5               sparseMatrixStats_1.12.2  future_1.33.0             shiny_1.7.5.1            
[141] ROCR_1.0-11               igraph_2.0.2              bit_4.0.5   
@yz3n18
Copy link

yz3n18 commented May 27, 2024

I had the same issue.
Following the suggestion from Alexis-Varin #27 solved my issue.

adata[["RNA"]] <- as(object = adata[["RNA"]], Class = "Assay")
SaveH5Seurat(adata,overwrite = T,verbose = T, filename ='./adata.h5Seurat')
Convert('./adata.h5Seurat',dest='h5ad') 

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants