Error with dmrfinder bsseq after combining Bsseq obj
0
0
Entering edit mode
br • 0
@br-12537
Last seen 7.3 years ago

Dear all,

I have a problem occurring using the command dmrfinder
> dmrs0 <- dmrFinder(bsseqbismark_smooth_chrI_location3_5_tstat, cutoff = c(-4.6, 4.6))
Error in dmrFinder(bsseqbismark_smooth_chrI_location3_5_tstat, cutoff = c(-4.6, : 'stat' needs to be a column of 'bstat@stats'
I checked the
> colnames(bsseqbismark_smooth_chrI_location3_5_tstat@stats)
[1] "rawSds" "tstat.sd" "group2.means" "group1.means" "tstat"

this is the same as in your example cancer data set, except the
> colnames(BS.cancer.ex.tstat@stats)
[1] "rawSds" "tstat.sd" "group2.means" "group1.means" "tstat" [6] "tstat.corrected"

My object is
An object of type 'BSseqTstat' with 50684 methylation loci based on smoothed data: BSmooth (ns = 70, h = 1000, maxGap = 100000000) with parameters BSmooth.tstat (local.correct = FALSE, maxGap = 100000000) 'stats' slot is in-memory

The difference I find to your example data is that I have also:
> class(bsseqbismark_smooth_chrI_location3_5_tstat@stats)
[1] "DelayedMatrix" attr(,"package") [1] "DelayedArray"

and you have

> class(BS.cancer.ex.tstat@stats)
[1] "matrix"
Do you know what is happening? Can I provide you with more information? Should I prevent to get a Delayed Array? Does this come from how I subsetted my data to the type 3 and 5?

bsseqbismark_smooth_chrI_location3_5 <- (bsseqbismark_smooth_chrI[, which(bsseqbismark_smooth_chrI$Type %in% c(3,5))])

 

Thanks  a lot!! 

 

 

R version 3.4.0 (2017-04-21)
Platform: x86_64-pc-linux-gnu (64-bit)
Running under: Ubuntu 16.04.2 LTS

Matrix products: default
BLAS: /usr/lib/libblas/libblas.so.3.6.0
LAPACK: /usr/lib/lapack/liblapack.so.3.6.0

locale:
 [1] LC_CTYPE=en_US.UTF-8       LC_NUMERIC=C               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    LC_PAPER=en_US.UTF-8       LC_NAME=C                 
 [9] LC_ADDRESS=C               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] bsseqData_0.14.0           bsseq_1.12.1               SummarizedExperiment_1.6.3 DelayedArray_0.2.7        
 [5] matrixStats_0.52.2         Biobase_2.36.2             GenomicRanges_1.28.3       GenomeInfoDb_1.12.2       
 [9] IRanges_2.10.2             S4Vectors_0.14.3           BiocGenerics_0.22.0       

loaded via a namespace (and not attached):
 [1] Rcpp_0.12.11            XVector_0.16.0          zlibbioc_1.22.0         munsell_0.4.3           colorspace_1.3-2       
 [6] lattice_0.20-35         plyr_1.8.4              tools_3.4.0             grid_3.4.0              data.table_1.10.4      
[11] R.oo_1.21.0             gtools_3.5.0            permute_0.9-4           Matrix_1.2-10           GenomeInfoDbData_0.99.0
[16] R.utils_2.5.0           bitops_1.0-6            RCurl_1.95-4.8          limma_3.32.2            compiler_3.4.0         
[21] R.methodsS3_1.7.1       scales_0.4.1            locfit_1.5-9.1    ```

 

bsseq r package dmrfinder delayedarray combine • 1.5k views
ADD COMMENT
0
Entering edit mode

I found the answer myself, I have to change the statistics settings as I set local.correct = FALSE before. I did this choice as I read somewhere this might be better for RRBS data. 

dmrs0 <- dmrFinder(obj.tstat, cutoff = c(-4.6, 4.6),stat = "tstat")

Thanks anyways! :)

 

 

ADD REPLY

Login before adding your answer.

Traffic: 650 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