limma design matrix
3
0
Entering edit mode
Auer Michael ▴ 250
@auer-michael-953
Last seen 10.2 years ago
Hi I know this has been dealt with before, but I would really appreciate some help. I am anaylsing affy chips and I have three different sample types, lpl-high, lpl-low and brain. I construct a factor vector below TS <-c(rep(?lpl-high?,12),rep(?lpl-low?,10),rep(?brain?,17)) TS<- factor(TS, levels= c(?lpl-high?,?lpl-low?,?brain?)) I would simply want to make an analysis of variance, see wheter these three samples are differentially expressed (F value). I am not interested in any contrasts. How can I do that? Thanks a lot!!!!!
affy BRAIN affy BRAIN • 1.2k views
ADD COMMENT
0
Entering edit mode
Naomi Altman ★ 6.0k
@naomi-altman-380
Last seen 3.6 years ago
United States
Create the contrasts. Use limma and FStat on the results. Or, to do a classical ANOVA, use apply to do an aov to each gene and use anova to get the F-stat from the aov output. e.g. myanova=function(x,y){ #x=expression for 1 gene #y=TS anova(aov(x~-1+y)) } ANOVAlist=apply(exprsmat,1,myanova,y=TS) Each component of ANOVAlist is the ANOVA table for 1 row of exprsmat. At 05:42 PM 1/10/2005 +0100, Auer Michael wrote: >Hi > >I know this has been dealt with before, but I would really appreciate some >help. >I am anaylsing affy chips and I have three different sample types, >lpl-high, lpl-low and brain. >I construct a factor vector below >TS <-c(rep("lpl-high",12),rep("lpl-low",10),rep("brain",17)) >TS<- factor(TS, levels= c("lpl-high","lpl-low","brain")) > >I would simply want to make an analysis of variance, see wheter these >three samples are differentially expressed (F value). I am not interested >in any contrasts. How can I do that? > >Thanks a lot!!!!! > >_______________________________________________ >Bioconductor mailing list >Bioconductor@stat.math.ethz.ch >https://stat.ethz.ch/mailman/listinfo/bioconductor Naomi S. Altman 814-865-3791 (voice) Associate Professor Bioinformatics Consulting Center Dept. of Statistics 814-863-7114 (fax) Penn State University 814-865-1348 (Statistics) University Park, PA 16802-2111
ADD COMMENT
0
Entering edit mode
@gordon-smyth
Last seen 2 hours ago
WEHI, Melbourne, Australia
See limma User's Guide http://bioinf.wehi.edu.au/limma/usersguide.pdf middle of page 32. Gordon Mon Jan 10 17:42:20 CET 2005 >Hi > >I know this has been dealt with before, but I would really appreciate some >help. >I am anaylsing affy chips and I have three different sample types, >lpl-high, lpl-low and brain. >I construct a factor vector below >TS <-c(rep("lpl-high",12),rep("lpl-low",10),rep("brain",17)) >TS<- factor(TS, levels= c("lpl-high","lpl-low","brain")) > >I would simply want to make an analysis of variance, see wheter these >three samples are differentially expressed (F value). I am not interested >in any contrasts. How can I do that? > >Thanks a lot!!!!!
ADD COMMENT
0
Entering edit mode
Shi, Tao ▴ 720
@shi-tao-199
Last seen 9.4 years ago
Is there a way to extract the original (without moderating the residual mean squares) F stat or F.p.value before apply eBayes step? I know I can do this by writing a function and apply it to each gene, but lmFit is much faster.... ...Tao ========================================================== Message: 11 Date: Tue, 11 Jan 2005 17:46:40 +1100 From: Gordon Smyth Subject: [BioC] limma design matrix To: "Auer Michael"michael.auer@meduniwien.ac.at Cc: bioconductor@stat.math.ethz.ch Message-ID: <6.2.0.14.1.20050111174221.0215e5d0@imaphost.wehi.edu.au> Content-Type: text/plain; charset="us-ascii"; format=flowed See limma User's Guide http://bioinf.wehi.edu.au/limma/usersguide.pdf middle of page 32. Gordon Mon Jan 10 17:42:20 CET 2005 >Hi > >I know this has been dealt with before, but I would really appreciate some >help. >I am anaylsing affy chips and I have three different sample types, >lpl-high, lpl-low and brain. >I construct a factor vector below >TS <-c(rep("lpl-high",12),rep("lpl-low",10),rep("brain",17)) >TS<- factor(TS, levels= c("lpl-high","lpl-low","brain")) > >I would simply want to make an analysis of variance, see wheter these >three samples are differentially expressed (F value). I am not interested >in any contrasts. How can I do that? > >Thanks a lot!!!!! __________________________________________________ [[alternative HTML version deleted]]
ADD COMMENT

Login before adding your answer.

Traffic: 690 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6