[BioC] topTable coefficient

Atul Kakrana atulkakrana at outlook.com
Tue Jun 11 02:00:46 CEST 2013


Hi Ed,

The 'coef' argument let user's specify the contrast for which toptable 
should be generated. So, if you made multiple contrasts like:

contrast.matrix<-makeContrasts(Ten-Control,Eleven-Control,Twelve-Control, levels=design)## 
We made three contrasts here ten-control, eleven-control and twelve-control
fit2<- contrasts.fit(fit,contrast.matrix)
fit2 <- eBayes(fit2)

Now if you want the toptable for first contrast in your matrix 
'Ten-Control' specify coef=1, like this:
topTable(fit2, coef = 1, adjust = "fdr", number = 50,sort.by='logFC')

If you want toptable for 2nd contrast i.e. 'Eleven-Control' than specify 
coef=2, like this:
topTable(fit2, coef = 2, adjust = "fdr", number = 50,sort.by='logFC')

So, its just a way specifying your contrasts in numbers rather than 
writing whole contrast. The 'coef' count starts from 1 and not 0, 
therefore you get an error with 'coef=0'.

Best

Atul




On 10-Jun-13 6:00 PM, Ed Siefker wrote:
> What does the 'coef' argument to topTable() actually do?
> I know it says in the manual:
>
> column number or column name specifying which coefficient or
> contrast of the linear model is of interest
>
> But what does that mean mathematically?  Why does specifying
> the coefficient alter the results of the top table?  Either a probe
> is differentially expressed or it isn't, right?
>
> I'm having trouble conceptualizing why there are multiple coefficients
> at this stage at all.  Suppose my contrast is "MUvsWT=MU-WT".
> How do you get more than one result when you subtract WT from MU?
>
> If I run topTable(fit,coef=1) does that get me MU, WT, or MU-WT?
> If I run topTable(fit, coef=2) does that get me MU, WT, or MU-WT?
> If I run topTable(fit, coef=0) or topTable(fit, coef=2), I get an error.
> Which of those isn't a coefficient and why not?
>
> 	[[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