[BioC] Help on experiment design

Mark Cowley m.cowley at garvan.org.au
Fri Aug 20 05:50:57 CEST 2010


Hi Marcelo,
take a look at the annotation packages here: http://bioconductor.org/packages/release/AnnotationData.html

you'll need to do the mapping from probe to gene before running heatmap.2, since heatmap.2 won't do that conversion for you

cheers,
Mark

On 20/08/2010, at 12:09 AM, Marcelo Brandão wrote:

> Thanks Mark.
> 
> Now, the data makes sense.
> I draw the results on a heatmap (heatmap.2 from gplots), do you have
> any idea on how to change the probes names to geneID? I am using
> Nimblegen data and have an NDF file with all annotation.
> 
> Cheers and thanks again
> 
> Marcelo
> 
> Em 19 de agosto de 2010 01:29, Mark Cowley <m.cowley at garvan.org.au> escreveu:
>> hi Marcelo,
>> 
>> If you use that design matrix, then you'll need to also fit contrasts
>>> designEsp <- cbind(controle=c(1,1,1,0,0,0),experimento=c(0,0,0,1,1,1))
>>> fitEsp <- lmFit(exprs(summarized),designEsp)
>> contr.matrix <- c(-1,1)
>> fitEsp <- contrasts.fit(fitEsp, contr.matrix)
>>> fitEsp <- eBayes(fitEsp)
>>> topTable(fitEsp,coef=1, adjust = "fdr", n = 10)
>> 
>> or use this design
>>> designEsp2 <- cbind(intercept=c(1,1,1,1,1,1),expvscon=c(0,0,0,1,1,1))
>>> fitEsp2 <- lmFit(exprs(summarized),designEsp2)
>>> fitEsp2 <- eBayes(fitEsp2)
>>> topTable(fitEsp2,coef=2, adjust = "fdr", n = 10)
>> 
>> cheers,
>> Mark
>> On 19/08/2010, at 10:59 AM, Marcelo Brandão wrote:
>> 
>>> Hello all!
>>> I am analyzing a set with 6 microarrays, 3 controls and 3
>>> experimentals. I am interested in look after differentially expressed
>>> genes. My doubt, among a lot of else, is if my experiment is designed
>>> correctly. I am currently using the following comands:
>>> 
>>> designEsp <- cbind(controle=c(1,1,1,0,0,0),experimento=c(0,0,0,1,1,1))
>>> fitEsp <- lmFit(exprs(summarized),designEsp)
>>> fitEsp <- eBayes(fitEsp)
>>> topTable(fitEsp,coef=2, adjust = "fdr", n = 10)
>>> 
>>> is it the best way to infer different expression? I am using nimblegen
>>> microarrays.
>>> 
>>> Thanks in advance.
>>> 
>>> Marcelo
>>> 
>>> --
>>> Marcelo Mendes Brandão
>>> Postdoc fellow
>>> Laboratório de Biologia Molecular de Plantas - ESALQ/USP
>>> Website: http://bioinfo.esalq.usp.br
>>> AtPIN: http://bioinfo.esalq.usp.br/atpin
>>> SKYPE: mmbrand
>>> Tel: (+55) 19 3429 4442
>>> 
>>> _______________________________________________
>>> 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
>> 
>> 
> 
> 
> 
> -- 
> Marcelo Mendes Brandão
> Postdoc fellow
> Laboratório de Biologia Molecular de Plantas - ESALQ/USP
> Website: http://bioinfo.esalq.usp.br
> AtPIN: http://bioinfo.esalq.usp.br/atpin
> SKYPE: mmbrand
> Tel: (+55) 19 3429 4442



More information about the Bioconductor mailing list