[Bioc-devel] problem w/ "exprs<-" method and LumiBatch objects
dear maintainers of Biobase and lumi, this is a question raised from a different thread a few days ago here: https://stat.ethz.ch/pipermail/bioc-devel/2014-January/005129.html the matrix of genes by samples expression values in an 'ExpressionSet' object can be accessed and replaced with the "exprs<-" method and this changes automatically the dimension of the 'ExpressionSet' object as illustrated in the following example: library(Biobase) data(sample.ExpressionSet) dim(sample.ExpressionSet) Features Samples 500 26 exprs(sample.ExpressionSet) <- matrix(rnorm(2600), nrow=100) dim(sample.ExpressionSet) Features Samples 100 26 however, this does not seem to be happening with 'LumiBatch' objects which are an extension of 'ExpressionSet' objects: library(lumi) data(example.lumi) dim(example.lumi) Features Samples 8000 4 exprs(example.lumi) <- matrix(rnorm(400), nrow=100) dim(example.lumi) ## NOTE NOW THAT DIMENSION DOES NOT CHANGE!! Features Samples 8000 4 a couple of diagnostics are that the expression matrix itself it does seem to have been replaced: dim(exprs(example.lumi)) [1] 100 4 and the "exprs<-" method does not seem to be explicitely defined for 'LumiBatch' objects: showMethods("exprs") Function: exprs (package Biobase) object="AffyBatch" object="ExpressionSet" object="MethyLumiSet" object="SnpSet" so i guess something is not working at the interaction of Biobase::exprs<- and LumiBatch-class, whether this should be fixed at Biobase or lumi, have no idea. thanks! robert. ps: sessionInfo() > sessionInfo() R version 3.0.2 (2013-09-25) Platform: x86_64-unknown-linux-gnu (64-bit) locale: [1] LC_CTYPE=en_US.UTF8 LC_NUMERIC=C LC_TIME=en_US.UTF8 [4] LC_COLLATE=en_US.UTF8 LC_MONETARY=en_US.UTF8 LC_MESSAGES=en_US.UTF8 [7] LC_PAPER=en_US.UTF8 LC_NAME=C LC_ADDRESS=C [10] LC_TELEPHONE=C LC_MEASUREMENT=en_US.UTF8 LC_IDENTIFICATION=C attached base packages: [1] parallel stats graphics grDevices utils datasets methods base other attached packages: [1] lumi_2.14.1 Biobase_2.22.0 BiocGenerics_0.8.0 vimcom_0.9-92 [5] setwidth_1.0-3 colorout_1.0-1 loaded via a namespace (and not attached): [1] affy_1.40.0 affyio_1.30.0 annotate_1.40.0 [4] AnnotationDbi_1.24.0 base64_1.1 beanplot_1.1 [7] BiocInstaller_1.12.0 biomaRt_2.18.0 Biostrings_2.30.1 [10] bitops_1.0-6 BSgenome_1.30.0 bumphunter_1.2.0 [13] codetools_0.2-8 colorspace_1.2-4 DBI_0.2-7 [16] digest_0.6.4 doRNG_1.5.5 foreach_1.4.1 [19] genefilter_1.44.0 GenomicFeatures_1.14.2 GenomicRanges_1.14.4 [22] grid_3.0.2 illuminaio_0.4.0 IRanges_1.20.6 [25] iterators_1.0.6 itertools_0.1-1 KernSmooth_2.23-10 [28] lattice_0.20-24 limma_3.18.7 locfit_1.5-9.1 [31] MASS_7.3-29 Matrix_1.1-1.1 matrixStats_0.8.14 [34] mclust_4.2 methylumi_2.8.0 mgcv_1.7-27 [37] minfi_1.8.9 multtest_2.18.0 nleqslv_2.1 [40] nlme_3.1-113 nor1mix_1.1-4 pkgmaker_0.17.4 [43] preprocessCore_1.24.0 RColorBrewer_1.0-5 RCurl_1.95-4.1 [46] registry_0.2 reshape_0.8.4 R.methodsS3_1.6.1 [49] rngtools_1.2.3 Rsamtools_1.14.2 RSQLite_0.11.4 [52] rtracklayer_1.22.0 siggenes_1.36.0 splines_3.0.2 [55] stats4_3.0.2 stringr_0.6.2 survival_2.37-4 [58] tools_3.0.2 XML_3.98-1.1 xtable_1.7-1 [61] XVector_0.2.0 zlibbioc_1.8.0