¿Es ANOVA de medidas repetidas un método correcto para mis datos?

Question

Mi experimento es un diseño de bloques completamente al azar. El objetivo es determinar si una variable es diferente entre las especies$a$,$b$,$c$. La medición se realizó 2 veces (junio, julio) en cada año (2011, 2012). Me preguntaba si ANOVA de medidas repetidas es un método correcto de usar. Si es así, ¿podría ayudarme a redactar la sintaxis para el análisis ANOVA y post-hoc?

Mis datos son como los siguientes:

 dat <- read.table(text = "species block   year    time    variable
a   1   2011    June    1
a   2   2011    June    2
a   3   2011    June    3
b   1   2011    June    4
b   2   2011    June    5
b   3   2011    June    6
c   1   2011    June    7
c   2   2011    June    8
c   3   2011    June    9
a   1   2011    July    10
a   2   2011    July    11
a   3   2011    July    12
b   1   2011    July    13
b   2   2011    July    14
b   3   2011    July    15
c   1   2011    July    16
c   2   2011    July    17
c   3   2011    July    18
a   1   2012    June    19
a   2   2012    June    20
a   3   2012    June    21
b   1   2012    June    22
b   2   2012    June    23
b   3   2012    June    24
c   1   2012    June    25
c   2   2012    June    26
c   3   2012    June    27
a   1   2012    July    28
a   2   2012    July    29
a   3   2012    July    30
b   1   2012    July    31
b   2   2012    July    32
b   3   2012    July    33
c   1   2012    July    34
c   2   2012    July    35
c   3   2012    July    36", header=TRUE)
 

Answer 1

La respuesta corta es sí, es una de medidas repetidas. Hay muchas maneras de hacerlo, en R. Ver este post para ver una lista de algunas de las posibilidades. Me gustaría utilizar el paquete nlme, suponiendo que los datos están distribuidos normalmente.

dat$block <- factor(dat$block)
dat$time = factor(paste(dat$time,dat$year,sep="."))
    dat$variable = rnorm(36)
library(nlme)
summary(lme(variable~species,data=dat,random=~1|block/time))

He combinado el año y la hora en un solo factor "tiempo". Si usted no hace esto solo tiene 2 niveles de año y dos de tiempo, y las cosas se ponen un poco absurdo - difícil obtener estimaciones de la variación con sólo 2 repeticiones.

Answer 2

El post atiretoo es bueno, y el enlace proporcionado es excelente. La función aov de base r estima el tipo I anova, y el SPSS por defecto es el tipo II. Es importante tener muy claro qué tipo de ano desea asegurar para obtener las estimaciones correctas.