js2264 / HiContacts

Analyzing Hi-C data in R with HiCExperiment objects
https://js2264.github.io/OHCA/
Other
10 stars 1 forks source link

unable to import pairs file #8

Closed ratheraarif closed 2 weeks ago

ratheraarif commented 8 months ago

When I try to import pairs file it throws the following error"

import(PairsFile(paste0(data_dir, 'DDSC.sorted.pairs')), format = 'pairs') Error in GenomeInfoDb:::getDanglingSeqlevels(x, new2old = new2old, pruning.mode = pruning.mode, :
The following seqlevels are to be dropped but are currently in use (i.e. have ranges on them): chr1, chr2, chr3, chr4, chr5, chr6, chr7, chr8, chr9, chr10, chr11, chr12, chr13, chr14, chr15, chr16, chrM. Please use the 'pruning.mode' argument to control how to prune 'x', that is, how to remove the ranges in 'x' that are on these sequences. For example, do something like: seqlevels(x, pruning.mode="coarse") <- new_seqlevels or: keepSeqlevels(x, new_seqlevels, pruning.mode="coarse") See ?seqinfo for a description of the pruning modes.


R version 4.3.2 (2023-10-31)
Platform: aarch64-apple-darwin20 (64-bit)
Running under: macOS Sonoma 14.3

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.3-arm64/Resources/lib/libRlapack.dylib;  LAPACK version 3.11.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: Asia/Kolkata
tzcode source: internal

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

other attached packages:
 [1] purrr_1.0.2                 diffHic_1.34.0              InteractionSet_1.30.0      
 [4] SummarizedExperiment_1.32.0 Biobase_2.62.0              MatrixGenerics_1.14.0      
 [7] matrixStats_1.2.0           GenomicRanges_1.54.1        GenomeInfoDb_1.38.5        
[10] IRanges_2.36.0              S4Vectors_0.40.2            dplyr_1.1.4                
[13] Cairo_1.6-2                 DNAZooData_1.2.0            fourDNData_1.2.0           
[16] HiContactsData_1.4.0        ExperimentHub_2.10.0        AnnotationHub_3.10.0       
[19] BiocFileCache_2.10.1        dbplyr_2.4.0                BiocGenerics_0.48.1        
[22] HiContacts_1.4.0            HiCool_1.2.0                ggplot2_3.4.4              
[25] HiCExperiment_1.3.1        

loaded via a namespace (and not attached):
  [1] strawr_0.0.91                 rstudioapi_0.15.0             jsonlite_1.8.8               
  [4] magrittr_2.0.3                ggbeeswarm_0.7.2              farver_2.1.1                 
  [7] rmarkdown_2.25                BiocIO_1.12.0                 zlibbioc_1.48.0              
 [10] vctrs_0.6.5                   memoise_2.0.1                 Rsamtools_2.18.0             
 [13] RCurl_1.98-1.14               htmltools_0.5.7               S4Arrays_1.2.0               
 [16] curl_5.2.0                    Rhdf5lib_1.24.1               SparseArray_1.2.3            
 [19] rhdf5_2.46.1                  htmlwidgets_1.6.4             basilisk_1.14.2              
 [22] plotly_4.10.4                 cachem_1.0.8                  GenomicAlignments_1.38.2     
 [25] mime_0.12                     lifecycle_1.0.4               pkgconfig_2.0.3              
 [28] Matrix_1.6-5                  R6_2.5.1                      fastmap_1.1.1                
 [31] GenomeInfoDbData_1.2.11       shiny_1.8.0                   digest_0.6.34                
 [34] colorspace_2.1-0              AnnotationDbi_1.64.1          plyinteractions_1.1.0        
 [37] RSpectra_0.16-1               pkgload_1.3.4                 RSQLite_2.3.5                
 [40] Rhtslib_2.4.1                 labeling_0.4.3                filelock_1.0.3               
 [43] fansi_1.0.6                   mgcv_1.9-1                    httr_1.4.7                   
 [46] abind_1.4-5                   compiler_4.3.2                bit64_4.0.5                  
 [49] withr_3.0.0                   csaw_1.36.1                   BiocParallel_1.36.0          
 [52] DBI_1.2.0                     rappdirs_0.3.3                DelayedArray_0.28.0          
 [55] sessioninfo_1.2.2             rjson_0.2.21                  tools_4.3.2                  
 [58] vipor_0.4.7                   beeswarm_0.4.0                interactiveDisplayBase_1.40.0
 [61] httpuv_1.6.14                 glue_1.7.0                    restfulr_0.0.15              
 [64] nlme_3.1-164                  promises_1.2.1                rhdf5filters_1.14.1          
 [67] grid_4.3.2                    generics_0.1.3                BSgenome_1.70.1              
 [70] gtable_0.3.4                  tzdb_0.4.0                    rmdformats_1.0.4             
 [73] tidyr_1.3.1                   data.table_1.14.10            hms_1.1.3                    
 [76] metapod_1.10.1                utf8_1.2.4                    XVector_0.42.0               
 [79] BiocVersion_3.18.1            pillar_1.9.0                  stringr_1.5.1                
 [82] limma_3.58.1                  vroom_1.6.5                   later_1.3.2                  
 [85] splines_4.3.2                 lattice_0.22-5                rtracklayer_1.62.0           
 [88] bit_4.0.5                     tidyselect_1.2.0              locfit_1.5-9.8               
 [91] Biostrings_2.70.1             knitr_1.45                    bookdown_0.37                
 [94] edgeR_4.0.12                  xfun_0.41                     statmod_1.5.0                
 [97] stringi_1.8.3                 lazyeval_0.2.2                yaml_2.3.8                   
[100] evaluate_0.23                 codetools_0.2-19              tibble_3.2.1                 
[103] BiocManager_1.30.22           cli_3.6.2                     xtable_1.8-4                 
[106] reticulate_1.34.0             munsell_0.5.0                 Rcpp_1.0.12                  
[109] dir.expiry_1.10.0             png_0.1-8                     XML_3.99-0.16.1              
[112] ggrastr_1.0.2                 parallel_4.3.2                ellipsis_0.3.2               
[115] readr_2.1.5                   blob_1.2.4                    basilisk.utils_1.14.1        
[118] plyranges_1.22.0              bitops_1.0-7                  viridisLite_0.4.2            
[121] scales_1.3.0                  crayon_1.5.2                  rlang_1.1.3                  
[124] KEGGREST_1.42.0               cowplot_1.1.3                
> 
js2264 commented 4 months ago

Hi @ratheraarif, have you tried importing the data with plyinteractions? See more details here: https://bioconductor.org/packages/release/bioc/vignettes/plyinteractions/inst/doc/plyinteractions.html#2_Importing_genomic_interactions_in_R