RegulatoryGenomicsGroup / chicdiff

A differential caller for capture Hi-C data
4 stars 3 forks source link

Running IHWcorrection: error #10

Open odovgusha opened 4 years ago

odovgusha commented 4 years ago

Hi, I tried the test dataset from the tutorial and everything worked fine but when I try to analyze my own data I get this error.

*** Running IHWcorrection

Comparison against p-vals for out
Comparison against p-vals for outcontrol
Trained weights on the control sample
Saving 7 x 7 in image
Saving 7 x 7 in image
Learned distance dependency
applied to test data
Error: Aesthetics must be either length 1 or the same as the data (1): colour, x and y
Run `rlang::last_error()` to see where the error occurred.

sessionInfo() output

R version 3.6.2 (2019-12-12)
Platform: x86_64-pc-linux-gnu (64-bit)
Running under: Scientific Linux 7.3 (Nitrogen)

Matrix products: default
BLAS:   /usr/lib64/libblas.so.3.4.2
LAPACK: /usr/lib64/liblapack.so.3.4.2

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

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

other attached packages:
[1] Chicago_1.14.0    Chicdiff_0.5      data.table_1.12.8

loaded via a namespace (and not attached):
 [1] Biobase_2.46.0              IHW_1.14.0
 [3] bit64_0.9-7                 splines_3.6.2
 [5] Formula_1.2-3               stats4_3.6.2
 [7] latticeExtra_0.6-29         blob_1.2.1
 [9] GenomeInfoDbData_1.2.2      slam_0.1-47
[11] pillar_1.4.4                RSQLite_2.2.0
[13] backports_1.1.8             lattice_0.20-38
[15] glue_1.4.1                  digest_0.6.25
[17] GenomicRanges_1.38.0        RColorBrewer_1.1-2
[19] XVector_0.26.0              checkmate_2.0.0
[21] colorspace_1.4-1            cowplot_1.0.0
[23] htmltools_0.5.0             Matrix_1.2-18
[25] DESeq2_1.26.0               XML_3.99-0.3
[27] pkgconfig_2.0.3             genefilter_1.68.0
[29] zlibbioc_1.32.0             xtable_1.8-4
[31] scales_1.1.1                jpeg_0.1-8.1
[33] fdrtool_1.2.15              BiocParallel_1.20.1
[35] htmlTable_2.0.0             tibble_3.0.1
[37] annotate_1.64.0             farver_2.0.3
[39] IRanges_2.20.2              ggplot2_3.3.2
[41] ellipsis_0.3.1              SummarizedExperiment_1.16.1
[43] hexbin_1.28.1               nnet_7.3-12
[45] BiocGenerics_0.32.0         survival_3.1-8
[47] magrittr_1.5                crayon_1.3.4
[49] memoise_1.1.0               MASS_7.3-51.4
[51] foreign_0.8-72              tools_3.6.2
[53] lifecycle_0.2.0             matrixStats_0.56.0
[55] stringr_1.4.0               S4Vectors_0.24.4
[57] locfit_1.5-9.4              munsell_0.5.0
[59] cluster_2.1.0               DelayedArray_0.12.3
[61] AnnotationDbi_1.48.0        compiler_3.6.2
[63] GenomeInfoDb_1.22.1         rlang_0.4.6
[65] grid_3.6.2                  RCurl_1.98-1.2
[67] rstudioapi_0.11             Delaporte_7.0.3
[69] htmlwidgets_1.5.1           labeling_0.3
[71] bitops_1.0-6                base64enc_0.1-3
[73] FMStable_0.1-2              gtable_0.3.0
[75] DBI_1.1.0                   R6_2.4.1
[77] gridExtra_2.3               knitr_1.29
[79] bit_1.1-15.2                Hmisc_4.4-0
[81] lpsymphony_1.14.0           stringi_1.4.6
[83] parallel_3.6.2              Rcpp_1.0.4.6
[85] harmonicmeanp_3.0           geneplotter_1.64.0
[87] vctrs_0.3.1                 rpart_4.1-15
[89] acepack_1.4.1               png_0.1-7
[91] xfun_0.15

Thanks!

Shall-We-Dance commented 2 years ago

I've met the same error tonight. But basically, I think it may because some difference in design files (.baitmap & .rmap), I'm still trying to solve it.

Shall-We-Dance commented 2 years ago

Changing chromosome numbers from chr1,chr2,chr3... to1,2,3... will help.