[R] random effects in analysis of variance

[Monika Ferster]

I have a data.frame paint (below-mentioned). In SPlus I used the command 

raov( MOISTURE ~ BATCH / PROBE )

(you could do raov( MOISTURE ~ BATCH +  PROBE%in%BATCH) as well)

so that the factors are taken as random. In R this function raov doesn't
exist. How can I calculate the same? Maybe with lme, but how? 

str(paint)
----------
`data.frame':	60 obs. of  5 variables:
 $ BATCH   : Factor w/ 15 levels "1","2","3","4",..: 1 1 1 1 2 2 2 2 3 3 ...
 $ SAMPLE  : num  1 1 2 2 3 3 4 4 5 5 ...
 $ REP     : num  1 2 1 2 1 2 1 2 1 2 ...
 $ MOISTURE: num  40 39 30 30 26 28 25 26 29 28 ...
 $ PROBE   : Factor w/ 2 levels "1","2": 1 1 2 2 1 1 2 2 1 1 ...

dput(paint)
-----------
structure(list(BATCH = structure(c(1, 1, 1, 1, 2, 2, 2, 2, 3, 
3, 3, 3, 4, 4, 4, 4, 5, 5, 5, 5, 6, 6, 6, 6, 7, 7, 7, 7, 8, 8, 
8, 8, 9, 9, 9, 9, 10, 10, 10, 10, 11, 11, 11, 11, 12, 12, 12, 
12, 13, 13, 13, 13, 14, 14, 14, 14, 15, 15, 15, 15), .Label = c("1", 
"2", "3", "4", "5", "6", "7", "8", "9", "10", "11", "12", "13", 
"14", "15"), class = "factor"), SAMPLE = c(1, 1, 2, 2, 3, 3, 
4, 4, 5, 5, 6, 6, 7, 7, 8, 8, 9, 9, 10, 10, 11, 11, 12, 12, 13, 
13, 14, 14, 15, 15, 16, 16, 17, 17, 18, 18, 19, 19, 20, 20, 21, 
21, 22, 22, 23, 23, 24, 24, 25, 25, 26, 26, 27, 27, 28, 28, 29, 
29, 30, 30), REP = c(1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 
1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 
2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 1, 2, 
1, 2, 1, 2), MOISTURE = c(40, 39, 30, 30, 26, 28, 25, 26, 29, 
28, 14, 15, 30, 31, 24, 24, 19, 20, 17, 17, 33, 32, 26, 24, 23, 
24, 32, 33, 34, 34, 29, 29, 27, 27, 31, 31, 13, 16, 27, 24, 25, 
23, 25, 27, 29, 29, 31, 32, 19, 20, 29, 30, 23, 24, 25, 25, 39, 
37, 26, 28), PROBE = structure(c(1, 1, 2, 2, 1, 1, 2, 2, 1, 1, 
2, 2, 1, 1, 2, 2, 1, 1, 2, 2, 1, 1, 2, 2, 1, 1, 2, 2, 1, 1, 2, 
2, 1, 1, 2, 2, 1, 1, 2, 2, 1, 1, 2, 2, 1, 1, 2, 2, 1, 1, 2, 2, 
1, 1, 2, 2, 1, 1, 2, 2), .Label = c("1", "2"), class = "factor")), .Names = c("BATCH", 
"SAMPLE", "REP", "MOISTURE", "PROBE"), row.names = c("1", "2", 
"3", "4", "5", "6", "7", "8", "9", "10", "11", "12", "13", "14", 
"15", "16", "17", "18", "19", "20", "21", "22", "23", "24", "25", 
"26", "27", "28", "29", "30", "31", "32", "33", "34", "35", "36", 
"37", "38", "39", "40", "41", "42", "43", "44", "45", "46", "47", 
"48", "49", "50", "51", "52", "53", "54", "55", "56", "57", "58", 
"59", "60"), class = "data.frame")
-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
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
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._