[BioC] combineBeadLevelLists() in Bioconductor 2.8 beadarray
Alex Gutteridge
alexg at ruggedtextile.com
Tue Dec 7 12:32:47 CET 2010
Hi Mark,
Sorry, maybe I'm being dim, but how do you call combine() correctly?
> dir1 <- "/5168671015"
> iBatch1 <- readIllumina(dir = paste(".", dir1, sep = ""),
+ useImages = FALSE,
+ metricsFile = "Metrics.txt",
+ illuminaAnnotation="Humanv4")
Processing section 5168671015_A_1
Processing section 5168671015_B_1
Processing section 5168671015_C_1
Processing section 5168671015_D_1
Processing section 5168671015_E_1
Processing section 5168671015_F_1
Processing section 5168671015_G_1
Processing section 5168671015_H_1
Processing section 5168671015_I_1
Processing section 5168671015_J_1
Processing section 5168671015_K_1
Processing section 5168671015_L_1
> dir2 <- "/5168671022"
> iBatch2 <- readIllumina(dir = paste(".", dir2, sep = ""),
+ useImages = FALSE,
+ metricsFile = "Metrics.txt",
+ illuminaAnnotation="Humanv4")
Processing section 5168671022_A_1
Processing section 5168671022_B_1
Processing section 5168671022_C_1
Processing section 5168671022_D_1
Processing section 5168671022_E_1
Processing section 5168671022_F_1
Processing section 5168671022_G_1
Processing section 5168671022_H_1
Processing section 5168671022_I_1
Processing section 5168671022_J_1
Processing section 5168671022_K_1
Processing section 5168671022_L_1
> iBatch = combine(iBatch1,iBatch2)
Error in combine(iBatch1, iBatch2) : unable to combine these objects
> class(iBatch1)
[1] "beadLevelData"
attr(,"package")
[1] "beadarray"
> class(iBatch2)
[1] "beadLevelData"
attr(,"package")
[1] "beadarray"
On Tue, 7 Dec 2010 10:46:21 +0000, Mark Dunning <mark.dunning at gmail.com>
wrote:
> Hi Alex,
>
> The functionality still exists in beadarray, but can be called using
> the generic 'combine' rather than combineBeadLevelLists.
>
> Best wishes,
>
> Mark
>
> On Tue, Dec 7, 2010 at 10:20 AM, Alex Gutteridge
> <alexg at ruggedtextile.com> wrote:
>> Hi,
>>
>> Older versions of beadarray have a combineBeadLevelList() for combining
>> bead level data lists read from two different Illumina chips into one
>> list.
>> In the latest beadarray this seems to have been removed, but I can't
find
>> anything in the documentation describing the correct way to do this
with
>> the new beadLevelData class. Am I missing something obvious or is this
no
>> longer supported?
>>
>>> sessionInfo()
>> R version 2.12.0 (2010-10-15)
>> Platform: x86_64-unknown-linux-gnu (64-bit)
>>
>> locale:
>> [1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C
>> [3] LC_TIME=en_US.UTF-8 LC_COLLATE=en_US.UTF-8
>> [5] LC_MONETARY=C LC_MESSAGES=en_US.UTF-8
>> [7] LC_PAPER=en_US.UTF-8 LC_NAME=C
>> [9] LC_ADDRESS=C LC_TELEPHONE=C
>> [11] LC_MEASUREMENT=en_US.UTF-8 LC_IDENTIFICATION=C
>>
>> attached base packages:
>> [1] splines stats graphics grDevices utils datasets methods
>> [8] base
>>
>> other attached packages:
>> [1] reshape_0.8.3 plyr_1.2.1 nlme_3.1-97
>> [4] lumiHumanAll.db_1.12.0 org.Hs.eg.db_2.4.6 RSQLite_0.9-3
>> [7] DBI_0.2-5 AnnotationDbi_1.12.0 lumi_2.2.0
>> [10] limma_3.6.6 latticeExtra_0.6-14 RColorBrewer_1.0-2
>> [13] lattice_0.19-13 contrast_0.13 Design_2.3-0
>> [16] Hmisc_3.8-3 survival_2.36-1 Cairo_1.4-5
>> [19] beadarray_2.0.2 Biobase_2.10.0
>>
>> loaded via a namespace (and not attached):
>> [1] affy_1.28.0 affyio_1.18.0 annotate_1.28.0
>> [4] cluster_1.13.2 grid_2.12.0 hdrcde_2.14
>> [7] KernSmooth_2.23-4 MASS_7.3-8 Matrix_0.999375-44
>> [10] methylumi_1.4.0 mgcv_1.7-1 preprocessCore_1.12.0
>> [13] tcltk_2.12.0 tools_2.12.0 xtable_1.5-6
>>
>> --
>> Alex Gutteridge
>>
>> _______________________________________________
>> Bioconductor mailing list
>> Bioconductor at r-project.org
>> https://stat.ethz.ch/mailman/listinfo/bioconductor
>> Search the archives:
>> http://news.gmane.org/gmane.science.biology.informatics.conductor
>>
--
Alex Gutteridge
More information about the Bioconductor
mailing list