[BioC] Sampsizelibrary.R

[Hassane, Duane]

Mira,

>From aov() ... one way might be:

aov.out<-aov(outcome~factor)
aov.table<-summary(aov.out)
F.stat1<-as.vector(aov.table)[[1]][1,4])
# alternatively ... unlist() aov.table to give a named vector of values of interest
# F.stat1<-unlist(aov.table)['F value1']

df1<-as.vector(aov.table)[[1]][1,1])
# alternatively as above
# df1<-unlist(aov.table)['Df1']

Or from lm() + anova() ... another way might be:

lm.out<-lm(outcome~factor)
aov.out<-anova(lm.out)
F.stat1<-aov.out[1,4]
df1<-aov.out[1,1]

Hope this helps,

Duane Hassane


-----Original Message-----
From: bioconductor-bounces at stat.math.ethz.ch
[mailto:bioconductor-bounces at stat.math.ethz.ch]On Behalf Of
Cuperlovic-Culf, Miroslava
Sent: Tuesday, April 03, 2007 2:41 PM
To: bioconductor at stat.math.ethz.ch
Subject: [BioC] Sampsizelibrary.R


Dear All,

I am trying to get a sample size estimate from sampsizelibrary.R (Pounds, S, Cheng, C. Bioinformatics, 21, 4263, 2005). I am trying to use ANOVA F-statistics data for bg.theta function. Unfortunately I am having hard time getting F-statistics vectors as well as df1 and df2 out of aov. So the question is how can I get this information from aov or another anova application.

Thanks

Mira

 

_____________________________________________________________________________

 

Miroslava Èuperloviæ-Culf, Ph.D. 

National Research Council                                      Conseil national de recherches Canada

Research Officer                                                    Agente de Recherche

55 Crowley Farm Road, Suite 1100                        55, Chemin Crowley Farm, Suite 1100

Moncton, New Brunswick E1A 7R1                       Moncton (Nouveau-Brunswick) E1A 7R1

 

(506) 861-0952                      

Docufax (506) 851-3630

miroslava.cuperlovic-culf at nrc-cnrc.gc.ca <mailto:miroslava.cuperlovic-culf at nrc-cnrc.gc.ca> 

www.iit.nrc-cnrc.gc.ca <http://www.iit.nrc-cnrc.gc.ca/> 

 

_____________________________________________________________________________

 

 

 

 


	[[alternative HTML version deleted]]