Question: bsseq - dmrFinder error
0
gravatar for parker
2.3 years ago by
parker0
Switzerland
parker0 wrote:

Hello,

I have previously used the bsseq package to get a list of dmrs successfully. However, now I have encountered an error when I run the dmrFinder - please see below:

Any help is much appreciated.

Thank you!

Hannah

> data.tstat

An object of type 'BSseqTstat' with
  3142216 methylation loci
based on smoothed data:
  BSmooth (ns = 20, h = 500, maxGap = 100000000)
with parameters
  BSmooth.tstat (local.correct = TRUE, maxGap = 100000000)

> dmrs_quantil <- dmrFinder(data.tstat, qcutoff=c(0.005,0.995))
Error in (1 - h) * qs[i] : non-numeric argument to binary operator

 

R version 3.3.1 (2016-06-21)
Platform: x86_64-pc-linux-gnu (64-bit)
Running under: Ubuntu 14.04.5 LTS

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

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

other attached packages:
[1] bsseq_1.10.0               limma_3.30.13             
[3] SummarizedExperiment_1.4.0 Biobase_2.34.0            
[5] GenomicRanges_1.26.2       GenomeInfoDb_1.10.3       
[7] IRanges_2.8.2              S4Vectors_0.12.1          
[9] BiocGenerics_0.20.0       

loaded via a namespace (and not attached):
 [1] Rcpp_0.12.11       XVector_0.14.1     zlibbioc_1.20.0    munsell_0.4.3     
 [5] colorspace_1.3-2   lattice_0.20-35    plyr_1.8.4         grid_3.3.1        
 [9] data.table_1.10.4  R.oo_1.21.0        gtools_3.5.0       matrixStats_0.52.2
[13] permute_0.9-4      Matrix_1.2-8       R.utils_2.5.0      bitops_1.0-6      
[17] RCurl_1.95-4.8     R.methodsS3_1.7.1  scales_0.4.1       locfit_1.5-9.1 

ADD COMMENTlink modified 2.3 years ago • written 2.3 years ago by parker0
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: 339 users visited in the last hour