[R] anova results?
dechao wang
dechwang at yahoo.co.uk
Thu Jan 24 17:08:29 CET 2002
Dear R users,
I did 16 experiments, with 16 responses (yield), and 4
factors(hidden, theres, lrate, cycl) with each having
2 levels as shown below. I want to do analysis of
variance to see which factor affect the response. The
coses are as follows:
hidden<-c(2,2,10,10,2,2,10,10,2,2,10,10,2,2,10,10)
theres<-c(0.0001, 0.1, 0.0001, 0.1, 0.0001, 0.1,
0.0001, 0.1, 0.0001, 0.1, 0.0001, 0.1, 0.0001, 0.1,
0.0001, 0.1)
lrate<-c(0.2,0.2,0.2,0.2,0.2,0.2,0.2,0.2,0.9,0.9,0.9,0.9,0.9,0.9,0.9,0.9)
cycl<-c(5000,5000,5000,5000,10000,10000,10000,10000,5000,5000,5000,5000,10000,10000,10000,10000)
yield<-c(0.26, 1.77, 0.29, 1.75, 0.26, 1.77, 0.31,
1.75, 0.29, 1.78, 0.35, 1.83, 0.29, 1.79, 0.35, 1.83)
npk <- data.frame(H=factor(hidden),
T=factor(theres),
L=factor(lrate),
C=factor(cycl), yield=yield)
( npk.aov <- aov(yield ~ H*T*L*C, npk) )
summary(npk.aov)
when I run it, I got the following results: but there
is no any F values and p-values. Can you please advise
me what I missed? Many thanks in advance, dechao
Df Sum Sq Mean Sq
H 1 0.0039 0.0039
T 1 8.8061 8.8061
L 1 0.0077 0.0077
C 1 0.0001 0.0001
H:T 1 0.0014 0.0014
H:L 1 0.0018 0.0018
H:C 1 6.25e-06 6.25e-06
T:L 1 0.0001 0.0001
T:C 1 6.25e-06 6.25e-06
L:C 1 6.25e-06 6.25e-06
H:T:L 1 0.0005 0.0005
H:T:C 1 0.0001 0.0001
H:L:C 1 0.0001 0.0001
T:L:C 1 0.0001 0.0001
H:T:L:C 1 6.25e-06 6.25e-06
__________________________________________________
Everything you'll ever need on one web page
from News and Sport to Email and Music Charts
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html
Send "info", "help", or "[un]subscribe"
(in the "body", not the subject !) To: r-help-request at stat.math.ethz.ch
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._
More information about the R-help
mailing list