[BioC] FW: normalization one color microarrays

Paz Tapia Ramirez verotapia at alumnos.utalca.cl
Tue Aug 30 00:05:23 CEST 2011












Step 14 is:
###################################################
### chunk number 14: filter.probes
###################################################
ddFILT=filter.probes(ddNORM,
                control=TRUE,
                wellaboveBG=TRUE,
                isfound=TRUE,
                wellaboveNEG=TRUE,
                sat=TRUE,
                PopnOL=TRUE,
                NonUnifOL=T,
                nas=TRUE,
                limWellAbove=75,
                limISF=75,
                limNEG=75,
                limSAT=75,
                limPopnOL=75,
                limNonUnifOL=75,
                limNAS=100,
                makePLOT=F,annotation.package="hgug4112a.db",flag.counts=T,targets)


I can run the step 13 (Correction background and Normalization).


sessionInfo() says:
R version 2.11.1 (2010-05-31) 
i686-pc-linux-gnu 

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

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

other attached packages:
 [1] vsn_3.16.0           annotate_1.26.1      hgug4112a.db_2.4.1  
 [4] org.Hs.eg.db_2.4.1   RSQLite_0.9-4        DBI_0.2-5           
 [7] AnnotationDbi_1.10.2 Biobase_2.8.0        marray_1.26.0       
[10] limma_3.4.5         

loaded via a namespace (and not attached):
[1] affy_1.26.1           affyio_1.16.0         grid_2.11.1          
[4] lattice_0.19-26       preprocessCore_1.10.0 xtable_1.5-6   

and specifically in step 14 says:

FILTERING BY ControlType FLAG 
------------------------------------------------------ 
 PROBES BEFORE FILTERING:  45220 
 PROBES AFTER ControlType FILTERING:  0 
 RAW DATA WITHOUT CONTROLS OUT : 0 
------------------------------------------------------ 
FILTERING BY IsWellAboveBG filterFLAG 

    FLAG FILTERING OPTIONS - FLAG OK = 1 - limWellAbove:  75 % 
    PROBES BEFORE FILTERING:  0 
    PROBES AFTER QC FILTERING:  0 
    IsNOTWellAboveBG OUT : 0 
------------------------------------------------------ 
FILTERING BY gIsFound filterFLAG 

    FLAG FILTERING OPTIONS - FLAG OK = 1 - limISF:  75 % 
    PROBES AFTER gIsFound FILTERING:  0 
    IsNOTFound OUT : 0 
------------------------------------------------------ 
Error en if (aux[jj] >= Limit[pos[jj]]) { : 
  valor ausente donde TRUE/FALSE es necesario

Thanks,
Paz


From: johannes.freudenberg at nih.gov
To: verotapia at alumnos.utalca.cl
CC: bioconductor at stat.math.ethz.ch
Date: Mon, 29 Aug 2011 13:12:42 -0400
Subject: RE: [BioC] normalization one color microarrays



Can you please be more specific?  I’m not really sure what step 14 is.  What exactly doesn’t work?   Please provide an example that doesn’t work for you and the kind of error messages it produces.  It’s also a good idea to provide the output of  > sessionInfo() when you post a question like this to the list. --Johannes  From: Paz Tapia Ramirez [mailto:verotapia at alumnos.utalca.cl] 
Sent: Monday, August 29, 2011 1:06 PM
To: Freudenberg, Johannes (NIH/NIEHS) [E]
Cc: bioconductor at stat.math.ethz.ch
Subject: RE: [BioC] normalization one color microarrays 
 when I work with Agi4x44PreProcess, from step 14 onwards it does not workFrom: johannes.freudenberg at nih.gov
To: verotapia at alumnos.utalca.cl
CC: bioconductor at stat.math.ethz.ch
Date: Mon, 29 Aug 2011 13:01:46 -0400
Subject: RE: [BioC] normalization one color microarraysHi Paz,  I’m not exactly sure what kind of microarrays you are trying to process but I’m guessing you are talking about Agilent 4x44 gene expression arrays.  If that is the case, I would recommend to use the Agi4x44Preprocess package.  You can download this package here:  http://www.bioconductor.org/packages/2.3/bioc/html/Agi4x44PreProcess.html  or by running the following in R > source("http://bioconductor.org/biocLite.R")> biocLite("Agi4x44PreProcess") This package comes with a nice manual which explains all the steps.  It also has plenty of examples you can try.  You can find the manual here: http://www.bioconductor.org/packages/2.3/bioc/vignettes/Agi4x44PreProcess/inst/doc/Agi4x44PreProcess.pdf I hope this helps for now.  Let me know if you more specific questions.   --Johannes  From: Paz Tapia Ramirez [mailto:verotapia at alumnos.utalca.cl] 
Sent: Monday, August 29, 2011 12:25 PM
To: Freudenberg, Johannes (NIH/NIEHS) [E]; bioconductor at stat.math.ethz.ch; bioconductor-bounces at r-project.org
Subject: RE: [BioC] normalization one color microarrays Hi,
  but if I make not summarization process, how do I build the set of expression?> From: johannes.freudenberg at nih.gov
> To: verotapia at alumnos.utalca.cl
> Date: Wed, 24 Aug 2011 16:17:07 -0400
> Subject: RE: [BioC] normalization one color microarrays
> 
> Summarization really only applies to Affymetrix arrays where you have several probes per probeset.
> 
> --Johannes
> 
> 
> -----Original Message-----
> From: Paz Tapia Ramirez [mailto:verotapia at alumnos.utalca.cl] 
> Sent: Wednesday, August 24, 2011 4:09 PM
> To: bioconductor at stat.math.ethz.ch; bioconductor-bounces at r-project.org
> Subject: [BioC] normalization one color microarrays
> 
> 
> 
> 
> 
> 
> 
> 
> 
> 
> Hello, I am starting to work with ONE color microarrays from Agilent, and I read that during the pre-processing is main tasks, such as background adjustment, normalization and summarization. Currently I have corrected the background and normalized samples, but do not understand the step of the summarization. Can you help me please?
> 
> Thks,
> Paz 
> [[alternative HTML version deleted]]
> 
> _______________________________________________
> 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 		 	   		   		 	   		  


More information about the Bioconductor mailing list