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@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]]
