ASpli: Error in h(simpleError(msg, call))
0
0
Entering edit mode
Dmitri • 0
@e48040ed
Last seen 3.3 years ago

Hello ASpli maintainers/users!

I am currently having a problem running jCounts. Am using the latest devel version of ASpli as I was also getting the jCounts() Error in av[at] <- a[at] which was previously discussed here. Unsure what the new the error means though. Changing to " threshold = 1" does not help either.

Code I ran is below, have also attached head(gbcounts@junction.counts) if needed as well as session info.

Any guidance would be greatly appreciated.

-Dmitri

gtfFileName <- "Homo_sapiens.GRCh38.93.gtf"

genomeTxDb <- makeTxDbFromGFF(gtfFileName)
features <- binGenome( genomeTxDb )

gbcounts <- gbCounts(features=features, targets=targets,
minReadLength = 90, maxISize = 10000, libType="SE",strandMode=0)

asd <- jCounts(counts=gbcounts, features=features, minReadLength=90, libType="SE",strandMode=0)

Error in h(simpleError(msg, call)) : 
  error in evaluating the argument 'x' in selecting a method for function 'sort': row names supplied are of the wrong length
> head(gbcounts@junction.counts)
                          junction            gene strand multipleHit          symbol gene_coordinates
1.11211.12010                noHit           noHit      *           -               -                -
1.11671.12010                noHit           noHit      *           -               -                -
1.11844.12010                noHit           noHit      *           -               -                -
1.12227.12613 ENSG00000223972:J001 ENSG00000223972      +           - ENSG00000223972    1:11869-14409
1.12697.13221                noHit ENSG00000223972      +           - ENSG00000223972    1:11869-14409
1.12697.13403                noHit ENSG00000223972      +           - ENSG00000223972    1:11869-14409
                                                                                                           bin_spanned
1.11211.12010                                                                ENSG00000223972:E001;ENSG00000223972:E002
1.11671.12010                                                                ENSG00000223972:E001;ENSG00000223972:E002
1.11844.12010                                                                ENSG00000223972:E001;ENSG00000223972:E002
1.12227.12613                                                                ENSG00000223972:E004;ENSG00000223972:E005
1.12697.13221                      ENSG00000223972:E005;ENSG00000223972:E006;ENSG00000223972:E007;ENSG00000223972:E008
1.12697.13403 ENSG00000223972:E005;ENSG00000223972:E006;ENSG00000223972:E007;ENSG00000223972:E008;ENSG00000223972:E009
              j_within_bin DLD1 HCT15 SW480 CACO2 HCT116 HCC2998 HT29 IS2 SW948 SKCO1
1.11211.12010            -    0     0     6     0      0       0    0   0     0     0
1.11671.12010            -    0     0     0     0      0       0    1   0     0     0
1.11844.12010            -    0     0     2     1      0       0    0   0     0     0
1.12227.12613            -    0     0     2     0      0       0    0   0     0     0
1.12697.13221            -    0     0     4     0      1       3    1   0     0     0
1.12697.13403            -    0     0     1     0      0       0    0   0     0     0

Session info below

> sessionInfo()
R version 4.1.0 (2021-05-18)
Platform: x86_64-pc-linux-gnu (64-bit)
Running under: CentOS Linux 7 (Core)

Matrix products: default
BLAS:   /stornext/System/data/apps/R/R-4.1.0/lib64/R/lib/libRblas.so
LAPACK: /stornext/System/data/apps/R/R-4.1.0/lib64/R/lib/libRlapack.so

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

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

other attached packages:
 [1] GenomicFeatures_1.45.1 GenomicRanges_1.45.0   GenomeInfoDb_1.29.3    ASpli_2.3.1           
 [5] AnnotationDbi_1.55.1   IRanges_2.27.2         S4Vectors_0.31.1       Biobase_2.53.0        
 [9] BiocGenerics_0.39.2    edgeR_3.35.0           limma_3.49.4          

loaded via a namespace (and not attached):
  [1] colorspace_2.0-2            rjson_0.2.20                ellipsis_0.3.2             
  [4] htmlTable_2.2.1             biovizBase_1.41.0           XVector_0.33.0             
  [7] base64enc_0.1-3             dichromat_2.0-0             rstudioapi_0.13            
 [10] DT_0.18                     bit64_4.0.5                 fansi_0.5.0                
 [13] xml2_1.3.2                  splines_4.1.0               cachem_1.0.6               
 [16] knitr_1.33                  Formula_1.2-4               Rsamtools_2.9.1            
 [19] cluster_2.1.2               dbplyr_2.1.1                png_0.1-7                  
 [22] BiocManager_1.30.16         compiler_4.1.0              httr_1.4.2                 
 [25] backports_1.2.1             lazyeval_0.2.2              assertthat_0.2.1           
 [28] Matrix_1.3-4                fastmap_1.1.0               htmltools_0.5.1.1          
 [31] prettyunits_1.1.1           tools_4.1.0                 igraph_1.2.6               
 [34] gtable_0.3.0                glue_1.4.2                  GenomeInfoDbData_1.2.6     
 [37] dplyr_1.0.7                 rappdirs_0.3.3              Rcpp_1.0.7                 
 [40] vctrs_0.3.8                 Biostrings_2.61.2           rtracklayer_1.53.1         
 [43] xfun_0.25                   stringr_1.4.0               lifecycle_1.0.0            
 [46] ensembldb_2.17.4            restfulr_0.0.13             XML_3.99-0.7               
 [49] zlibbioc_1.39.0             MASS_7.3-54                 scales_1.1.1               
 [52] BiocStyle_2.21.3            BSgenome_1.61.0             VariantAnnotation_1.39.0   
 [55] ProtGenerics_1.25.1         hms_1.1.0                   MatrixGenerics_1.5.3       
 [58] SummarizedExperiment_1.23.1 AnnotationFilter_1.17.1     RColorBrewer_1.1-2         
 [61] yaml_2.2.1                  curl_4.3.2                  memoise_2.0.0              
 [64] gridExtra_2.3               ggplot2_3.3.5               UpSetR_1.4.0               
 [67] rpart_4.1-15                biomaRt_2.49.4              latticeExtra_0.6-29        
 [70] stringi_1.7.3               RSQLite_2.2.7               BiocIO_1.3.0               
 [73] checkmate_2.0.0             filelock_1.0.2              BiocParallel_1.27.3        
 [76] rlang_0.4.11                pkgconfig_2.0.3             matrixStats_0.60.0         
 [79] bitops_1.0-7                evaluate_0.14               lattice_0.20-44            
 [82] purrr_0.3.4                 htmlwidgets_1.5.3           GenomicAlignments_1.29.0   
 [85] bit_4.0.4                   tidyselect_1.1.1            plyr_1.8.6                 
 [88] magrittr_2.0.1              R6_2.5.1                    generics_0.1.0             
 [91] Hmisc_4.5-0                 DelayedArray_0.19.1         DBI_1.1.1                  
 [94] pillar_1.6.2                foreign_0.8-81              nnet_7.3-16                
 [97] survival_3.2-12             KEGGREST_1.33.0             RCurl_1.98-1.4             
[100] tibble_3.1.3                crayon_1.4.1                utf8_1.2.2                 
[103] BiocFileCache_2.1.1         rmarkdown_2.10              jpeg_0.1-9                 
[106] progress_1.2.2              locfit_1.5-9.4              grid_4.1.0                 
[109] data.table_1.14.0           blob_1.2.2                  digest_0.6.27              
[112] pbmcapply_1.5.0             tidyr_1.1.3                 munsell_0.5.0              
[115] Gviz_1.37.2
ASpli • 1.9k views
ADD COMMENT
0
Entering edit mode

Hi there. I have the same issue. Any update regarding this problem?

ADD REPLY

Login before adding your answer.

Traffic: 561 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6