Question: Conditional Quantile Normalization (cqn) error
gravatar for Eric.Martin
4 months ago by
Eric.Martin0 wrote:

I generated a count list using featureCounts of my ATAC-seq data and fed it into R. I got the GC bias and lenght for each peak queried, and made a list for the sequencing depth. When I run the following command:

ataccountscqn <- cqn(ataccounts0715, x = ataccovar3$gccontent, lengths = ataccovar3$length, sizeFactors = atac.sizefactor, verbose = TRUE)

I get the following error:

Error in qr.default(t(const)) : NA/NaN/Inf in foreign function call (arg 1)

I checked for NAs in the datasets, there aren't any. I checked the types of data in each file to make sure there are numbers and not characters for the values, and a checked the type() of each dataset to make sure it aligns with the sample data. The sample dataset works fine and produces an output, so i'm at a loss as to what (I most likely did wrong) is going on.

sessionInfo() R version 3.6.0 (2019-04-26) Platform: x86_64-apple-darwin15.6.0 (64-bit) Running under: macOS Sierra 10.12.6

Matrix products: default BLAS: /System/Library/Frameworks/Accelerate.framework/Versions/A/Frameworks/vecLib.framework/Versions/A/libBLAS.dylib LAPACK: /Library/Frameworks/R.framework/Versions/3.6/Resources/lib/libRlapack.dylib

Random number generation: RNG: Mersenne-Twister Normal: Inversion Sample: Rounding

locale: [1] enUS.UTF-8/enUS.UTF-8/enUS.UTF-8/C/enUS.UTF-8/en_US.UTF-8

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

other attached packages: [1] dplyr0.8.3 scales1.0.0 SummarizedExperiment1.14.0 [4] DelayedArray0.10.0 BiocParallel1.18.0 matrixStats0.54.0
[7] Biobase2.44.0 GenomicRanges1.36.0 GenomeInfoDb1.20.0
[10] IRanges
2.18.1 S4Vectors0.22.0 BiocGenerics0.30.0
[13] Matrix1.2-17 readr1.3.1 readxl1.3.1
[16] cqn
1.30.0 quantreg5.42 SparseM1.77
[19] preprocessCore1.46.0 nor1mix1.3-0 mclust_5.4.5

loaded via a namespace (and not attached): [1] bitops1.0-6 DirichletMultinomial1.26.0 TFBSTools1.22.0
[4] bit64
0.9-7 httr1.4.0 backports1.1.4
[7] tools3.6.0 utf81.1.4 R62.4.0
[10] DT
0.7 seqLogo1.50.0 DBI1.0.0
[13] lazyeval0.2.2 colorspace1.4-1 tidyselect0.2.5
[16] bit
1.1-14 compiler3.6.0 cli1.1.0
[19] plotly4.9.0 rtracklayer1.44.0 caTools1.17.1.2
[22] stringr
1.4.0 digest0.6.20 Rsamtools2.0.0
[25] R.utils2.9.0 XVector0.24.0 pkgconfig2.0.2
[28] htmltools
0.3.6 BSgenome1.52.0 htmlwidgets1.3
[31] rlang0.4.0 rstudioapi0.10 RSQLite2.1.1
[34] VGAM
1.1-1 shiny1.3.2 jsonlite1.6
[37] gtools3.8.1 R.oo1.22.0 RCurl1.95-4.12
[40] magrittr
1.5 GO.db3.8.2 GenomeInfoDbData1.2.1
[43] fansi0.4.0 Rcpp1.0.1 munsell0.5.0
[46] R.methodsS3
1.7.1 stringi1.4.3 zlibbioc1.30.0
[49] plyr1.8.4 grid3.6.0 blob1.2.0
[52] promises
1.0.1 crayon1.3.4 miniUI0.1.1.1
[55] CNEr1.20.0 lattice0.20-38 Biostrings2.52.0
[58] annotate
1.62.0 hms0.5.0 KEGGREST1.24.0
[61] zeallot0.1.0 pillar1.4.2 reshape21.4.3
[64] TFMPvalue
0.0.8 XML3.98-1.20 glue1.3.1
[67] BiocManager1.30.4 data.table1.12.2 vctrs0.2.0
[70] png
0.1-7 chromVAR1.6.0 httpuv1.5.1
[73] cellranger1.1.0 MatrixModels0.4-1 tidyr0.8.3
[76] gtable
0.3.0 poweRlaw0.70.2 purrr0.3.2
[79] assertthat0.2.1 ggplot23.2.0 mime0.7
[82] xtable
1.8-4 later0.8.0 viridisLite0.3.0
[85] tibble2.1.3 GenomicAlignments1.20.1 AnnotationDbi1.46.0
[88] memoise

ADD COMMENTlink modified 4 months ago by Kasper Daniel Hansen6.4k • written 4 months ago by Eric.Martin0
Answer: Conditional Quantile Normalization (cqn) error
gravatar for Kasper Daniel Hansen
4 months ago by
United States
Kasper Daniel Hansen6.4k wrote:

If you send me the input, I can take a slow look.

Best, Kasper (Sent from my phone.)

ADD COMMENTlink written 4 months ago by Kasper Daniel Hansen6.4k
Please log in to add an answer.


Use of this site constitutes acceptance of our User Agreement and Privacy Policy.
Powered by Biostar version 16.09
Traffic: 192 users visited in the last hour