[BioC] Lumi problem

Davis, Wade davisjwa at health.missouri.edu
Tue Mar 8 21:00:57 CET 2011


Mali,
One start is to use a lumi*IDMapping package instead, as shown below.

## load the data
matt_exp <- lumiR(fileName,lib.mapping="lumiHumanIDMapping")

I've also had intermittent problems with the function addControlData2lumi, and have used the alternate route below to get the info into the controlData slot. I have NOT looked at the source code of the addControlData2lumi function for a few years, but this should, in theory, be doing the same thing...or at least getting the same end result. But in my experience, this gives me less trouble.

################################
#Add QC Info
################################

controlFile <- paste(filedir,'QC.txt',sep="")
controlData<-getControlData(controlFile, type = "data.frame")

matt_exp at controlData<-controlData

Good luck,
Wade Davis


-----Original Message-----
From: mali salmon [mailto:shalmom1 at gmail.com] 
Sent: Tuesday, March 08, 2011 1:30 AM
To: bioconductor at r-project.org
Subject: Re: [BioC] Lumi problem

I forgot to add sessionInfo()

> sessionInfo()
R version 2.12.1 (2010-12-16)
Platform: i386-redhat-linux-gnu (32-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] grid      stats     graphics  grDevices utils     datasets  methods
[8] base

other attached packages:
 [1] genefilter_1.32.0      gplots_2.8.0           caTools_1.11
 [4] bitops_1.0-4.1         gdata_2.8.1            gtools_2.6.2
 [7] statmod_1.4.8          lumiMouseAll.db_1.12.0 org.Mm.eg.db_2.4.6
[10] lumiHumanAll.db_1.12.0 org.Hs.eg.db_2.4.6     RSQLite_0.9-4
[13] DBI_0.2-5              annotate_1.28.0        AnnotationDbi_1.12.0
[16] limma_3.6.9            lumi_2.2.1             Biobase_2.10.0

loaded via a namespace (and not attached):
 [1] affy_1.28.0           affyio_1.18.0         hdrcde_2.15
 [4] KernSmooth_2.23-4     lattice_0.19-13       MASS_7.3-9
 [7] Matrix_0.999375-46    methylumi_1.6.1       mgcv_1.7-2
[10] nlme_3.1-97           preprocessCore_1.12.0 splines_2.12.1
[13] survival_2.36-2       tcltk_2.12.1          tools_2.12.1
[16] xtable_1.5-6




On Tue, Mar 8, 2011 at 9:23 AM, mali salmon <shalmom1 at gmail.com> wrote:

> Hi All
> I am analysing illumina array data using "Lumi", and I found a problem with
> adding control data to Lumi.
> I have tried:
>
> >addControlData2lumi("illuminaControls.txt",x.lumi)
> and got the following message:
> "Annotation columns are not available in the data".
>
> Then I tried to create the LumiBatch by specifing the mouse annotation
> library
> >x.lumi<-lumiR("illuminaAll.txt",lib="lumiMouseAll.db")
> >addControlData2lumi("illuminaControls.txt",x.lumi)
>
> but it didn't help, when I try to do background correction, I still get the
> error:
> "There is no control probe information in the LumiBatch object!
>  No background adjustment will be performed."
>
> Any ideas?
> Thanks
> Mali
>
>
>


-- 
Mali Salmon-Divon PhD
EMBL European Bioinformatics Institute
Wellcome Trust Genome Campus
Cambridge CB10 1SD
tel +44 (0)1223 494444
fax +44 (0)1223 494468

	[[alternative HTML version deleted]]



More information about the Bioconductor mailing list