[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