Question: Conditional Quantile Normalization (cqn) error
0
gravatar for Eric.Martin
4 weeks ago by
UCSC
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
1.1.0

ADD COMMENTlink modified 4 weeks ago by Kasper Daniel Hansen6.4k • written 4 weeks ago by Eric.Martin0
Answer: Conditional Quantile Normalization (cqn) error
0
gravatar for Kasper Daniel Hansen
4 weeks 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 weeks ago by Kasper Daniel Hansen6.4k
Please log in to add an answer.

Help
Access

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