Hi all,
 
I have a list of gene symbols generated from the differential expression analysis below. How do I convert these symbols to emsembl id? Thanks 
 
 
Regards,
Michelle
 
 
 
 
R version 2.14.1 (2011-12-22)
Platform: x86_64-pc-mingw32/x64 (64-bit)


&gt; library(affy)
&gt;  library(limma)
&gt; pd=read.AnnotatedDataFrame("phenodata.txt",header=TRUE,sep="",row.names=1)
&gt;  a=ReadAffy(filenames=rownames(pData(pd)),phenoData=pd,verbose=TRUE)
1 reading Control-1.cel ...instantiating an AffyBatch (intensity a 1004004x6 matrix)...done.
Reading in : Control-1.cel
Reading in : Control-2.cel
Reading in : Dicer-1.cel
Reading in : Dicer-2.cel
Reading in : Drosha-1.cel
Reading in : Drosha-2.cel
&gt; x=rma(a)
Loading required package: AnnotationDbi
Background correcting
Normalizing
Calculating Expression
Warning message:
package ‘AnnotationDbi’ was built under R version 2.14.2 
&gt; c=paste(pd$treatment,pd$n,sep="")
&gt;  f=factor(c)
&gt;  design=model.matrix(~0+f)
&gt; colnames(design)=levels(f)
&gt; fit=lmFit(x,design)
&gt; library(mouse4302.db)
Loading required package: org.Mm.eg.db
Loading required package: DBI

Warning messages:
1: package ‘RSQLite’ was built under R version 2.14.2 
2: package ‘DBI’ was built under R version 2.14.2 
&gt; library(annotate)
Warning message:
package ‘annotate’ was built under R version 2.14.2 
&gt;  fit$genes$Symbol &lt;- getSYMBOL(fit$genes$ID,"mouse4302.db")
&gt; contrast.matrix=makeContrasts(E1="present-absent.Dicer",E2="present-absent.Drosha",E3="absent.Drosha-absent.Dicer",levels=design)
&gt;  fit2=contrasts.fit(fit,contrast.matrix)
&gt; 
&gt; fit2=eBayes(fit2)
&gt; 
&gt; results1 &lt;-topTable (fit2, coef=1, p.value=0.0001,number=nrow(fit2))
&gt; write.table(results1, file="control-Dicer5.txt")
&gt; results2 &lt;-topTable (fit2, coef=2, p.value=0.0001,number=nrow(fit2))
&gt; write.table(results2, file="control-Drosha5.txt")
&gt; results3 &lt;-topTable (fit2, coef=3, p.value=0.0001,number=nrow(fit2))
&gt; results=decideTests(fit2)
&gt; summary(results2)
&gt; b=venncounts(results2)
&gt; print(b)
&gt; vennDiagram(results)
&gt; a=vennDiagram(results,include=c("up","down"),counts.col=c("red","green")) 



	[[alternative HTML version deleted]]

