[R] random and group arguments in lme

[Eusebio Arenal Gutiérrez]

I want to fit the following model:

> aov.madera <- aov(tiempo~producto/panel, data = madera)
> anova(aov.madera)
Analysis of Variance Table

Response: tiempo
               Df Sum Sq Mean Sq F value    Pr(>F)    
producto        2 93.631  46.815 63.0511 4.304e-07 ***
producto:panel  9 43.533   4.837  6.5144  0.001881 ** 
Residuals      12  8.910   0.742                      
---
Signif. codes:  0  `***'  0.001  `**'  0.01  `*'  0.05  `.'  0.1  ` '  1 

As panel is random I would need
["producto", "Mean Sq"]/["producto:panel", "Mean Sq"] in the cell
["producto", "F value"]
and its p-value, that is

> 46.815/4.837
[1] 9.67852
> 2*pf(46.815/4.837, 2, 9, lower.tail=FALSE)
[1] 0.01143266

How can I fit this model with lme function (of nlme)? What arguments have I
to put in the random and group arguments?

Where can I find references about the nlme library? a tutorial?

Thaks in advance,

Eusebio


       _\|/_      
       (o o)                                         
*---oOO-(_)-OOo-----------------------------------,,,-------*
*                      use at eio.uva.es            (o o)      *
*---------------------------------------------oOO-(_)-OOo---*


-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-
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
_._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._._