[BioC] question about topTable in limma
Sean Davis
sdavis2 at mail.nih.gov
Tue Apr 26 23:08:55 CEST 2011
On Mon, Apr 25, 2011 at 11:53 PM, samane fazeli <samanefazeli at gmail.com> wrote:
> Hello
> Dear all
> I have a question about topTable in limma package. Can i specify
> parameters of interest in
> topTable, i.e. it gives me specific parameters such as adj.P-value
> and logFC, not ID and names
> columns because both of them are caused to disorder output file of
> topTable. I would like to
> process our data in excel
>
> limma(library)
> targets=readTargets("ApoAITargets.txt")
> RG=read.maimages(targets$FileName, source="spot")
> MA=normalizeWithinArrays(RG, method="loess")
> controlRef=rep(1,16)
> KOcontrol=c(rep(0,8),rep(1,8))
> design=cbind(controlRef,KOcontrol)
> fit=lmFit(MA,design)
> fit=eBayes(fit)
> option(digits=3)
> result=topTable(fit, coef=2, number=6384, adjust="BH")
> as.matrix(result)
> write.csv(result,file="printtiploess.csv", quote=FALSE, col.names=TRUE,
> sep=",")
Hi. See the genelist argument to topTable().
Sean
More information about the Bioconductor
mailing list