[BioC] Extract top-ranked genes from model fit

James W. MacDonald jmacdon at med.umich.edu
Tue Oct 13 16:02:57 CEST 2009


Hi Chintanu,

Chintanu wrote:
> Hi,
> 
> I wish to get the UniGene annotation of the top ranked genes for two sets of
> Affymetrix files (Diseased vs. Reference; chip type: hgu133a). I am
> uncertain as to how to approach it.
> 
> A few of the simple background codes are:
> 
> dataset <- ReadAffy() ;  myRMA <- justRMA()
> 
> fit <- lmFit (myRMA, design)  # Analysis design is in object, “design”
> 
> contrasts_matrix <- makeContrasts(Diseased-Healthy, levels = design)
> 
>  fit2 <- contrasts.fit (fit, contrasts_matrix)
> 
> fit3 <- eBayes (fit2)
> 
> test_results <- topTable(fit3, number=10, adjust="BH")
> 
> # Then, I was uncertain regarding how to get the output with UniGene ID.
> Wonder if a bit of tweaking in the function, probes2table should somehow
> solve this !
> 
> # probes2table(****, featureNames(***), "hgu133a.db", anncols =
> aaf.handler()[c(1,2,7,9)], html = FALSE, text = TRUE, filename = "****")

Did you try it? This should output probe, symbol, gene name and UniGene 
ID, without any 'tweaking' required.

Best,

Jim


> 
> Many thanks,
> 
> Chintanu
> 
> 
> 
>> sessionInfo ()
> 
> R version 2.9.2 (2009-08-24)
> 
> i386-pc-mingw32
> 
> locale:
> 
> LC_COLLATE=English_New Zealand.1252;LC_CTYPE=English_New
> Zealand.1252;LC_MONETARY=English_New
> Zealand.1252;LC_NUMERIC=C;LC_TIME=English_New Zealand.1252
> 
> attached base packages:
> 
> [1] grDevices datasets  splines   graphics  stats     tcltk     utils
> methods   base
> 
> other attached packages:
> 
>  [1] annaffy_1.16.0       hgu133a.db_2.2.12    affycoretools_1.16.3
> KEGG.db_2.2.11       GO.db_2.2.11         RSQLite_0.7-2
> DBI_0.2-4
> 
> 
>  [8] AnnotationDbi_1.6.1  limma_2.18.3         hgu133acdf_2.4.0
> affy_1.22.1          Biobase_2.4.1        survival_2.35-7
> 
> loaded via a namespace (and not attached):
> 
>  [1] affyio_1.12.0        annotate_1.22.0      biomaRt_2.0.0
> Category_2.10.1      gcrma_2.16.0         genefilter_1.24.2
> GOstats_2.10.0
> 
>  [8] graph_1.22.2         GSEABase_1.6.1       preprocessCore_1.6.0
> RBGL_1.20.0          RCurl_1.2-0          tools_2.9.2
> XML_2.6-0
> 
> 
> [15] xtable_1.5-5
> 
> 	[[alternative HTML version deleted]]
> 
> 
> 
> ------------------------------------------------------------------------
> 
> _______________________________________________
> Bioconductor mailing list
> Bioconductor at stat.math.ethz.ch
> https://stat.ethz.ch/mailman/listinfo/bioconductor
> Search the archives: http://news.gmane.org/gmane.science.biology.informatics.conductor

-- 
James W. MacDonald, M.S.
Biostatistician
Douglas Lab
University of Michigan
Department of Human Genetics
5912 Buhl
1241 E. Catherine St.
Ann Arbor MI 48109-5618
734-615-7826



More information about the Bioconductor mailing list