[R-es] anova medidas repetidas con lme

[Argel Gastélum Arellánez]

     Hola compañeros de la lista.

     Tengo el siguiente set de datos:

Repeticiones <- c(rep("RI", 14), rep("RII", 14), rep("RIII", 14))
Tiempo <- rep(c(0, 2, 4, 6, 8, 10, 12, 24, 36, 48, 60, 72, 96, 120), 3)
Concentracion_celular <- c(0.4862, 0.5375, 0.4309, 0.4390, 0.4603, 
0.4733, 0.3936, 0.9085, 0.5838, 0.5477, 0.6331, 0.8693, 1.0092, 0.6341, 
0.5350, 0.4951, 0.4120, 0.5032, 0.5968, 0.5734, 0.5258, 4.2884, 1.0763, 
0.6967, 1.6744, 0.6124, 1.1766, 0.6999, 0.5774, 0.6434, 0.4921, 0.5627, 
0.7045, 0.7065, 0.5667, 1.8556, 1.3378, 1.2013, 1.6612, 0.5042, 1.2924, 
1.0463)

DATOS <- data.frame(REPETICIONES = Repeticiones, TIEMPO = 
factor(Tiempo), CONCENTRACION_CELULAR = Concentracion_celular)

     Se trata de cultivos bacterianos, realizado con tres repeticiones, 
a los que se les monitoreó la concentración celular en el tiempo (en 
horas). Para determinar si existe diferencia significativa en la 
concentración celular en los diferentes tiempos necesito aplicar un 
anova de medidas repetidas, lo que estoy haciendo con la siguiente 
instrucción (tomado de 
http://blog.gribblelab.org/2009/03/09/repeated-measures-anova-using-r/):

# Modelo medidas repetidas:
library(nlme)

modelo.medidas.repetidas <- lme(CONCENTRACION_CELULAR ~ TIEMPO, random = 
~1|REPETICIONES/TIEMPO, na.action = na.exclude, data = DATOS)

¿Es adecuado aplicar el modelo de esta forma para realizar un anova de 
medidas repetidas? ¿Tiene sentido que aparezca la variable TIEMPO en el 
argumento "random" o sólo con dejar REPETICIONES es suficiente? Lo he 
aplicado de las dos forma y me da el mismo resultado para el valor de p.

     De antemano muchas gracias por la ayuda.

     Saludos.

--
     Argel.