Estimada Lucía
Si no entiendo mal en su código anida semana dentro de año (=~1|year/week ),
como también (~plot*year) analiza las combinaciones entre parcelas y años.
Le comento que yoo tengo datos repetidos, el modelo también es en campos y
con animales, donde el efecto ambiental es importante...
Escribo muchos modelos en R, según todo lo que se me ocurra, algunas
variantes o alternativas posibles en R o estadísticamente medibles, no las
tengo en cuenta porque biológicamente me parecen sin sentido (la computadora
calcula cualquier cosa), en cambio observo entre los modelos distintos
algunas partes que me resultan interesantes, por ejemplo, al anidar si bien
la información que aporta R es superior, el AIC solo cambia en 1 respecto a
lo mismo sin anidar, los valores que a mi me interesan no se modifican
(sería la parcela en su caso), el tiempo de computo aumenta muchísimo, en
los números no me conviene anidar, aunque biológicamente es más lógico. Al
observar los resultados en algunas variables fijas, veo que tendría que
mejorar la parte estadística de esas variables, porque si bien en los
números es estadísticamente correcto, por biología se que eso no es del todo
lógico, en ese caso prefiero no mirar tanto el AIC y si dar respuesta en
términos biológicos a mis colegas (que como yo, ninguno es estadista).
La respuesta, pruebe, escriba, pruebe, tendrá dudas, pruebe otra
alternativa, ... al final se quedara con una respuesta para su caso en
particular, y muchas ideas para mejorar.
Javier Marcuzzi
-----Mensaje original-----
From: Lucía Torres
Sent: Friday, May 04, 2012 8:17 AM
To: r-help-es en r-project.org
Subject: [R-es] Medidas Repetidas
Hola a todos,
Estoy intentando realizar un análisis de medidas repetidas en R y me
gustaría pediros opinión porque estoy un poco atascada. Estoy utilizando
las funciones lme y lmer.
Tengo datos de censos de grillos (efi) realizados semanalmente (week) en
tres parcelas (plot) diferentes, cada una de las cuales tiene seis
transectos (transect) donde se realizan dichos censos. Los datos son de dos
años consecutivos (year)
mi modelo sería, en principio, algo así:
lmer.efi<-lmer(mean_efi_week ~plot*year+(1|transect),data = efi_sem, method
= "REML")
lme.efi<-lme(efi ~year* plot,random=~1|year/week | plot/transect,data
efi_sem, method = "REML")
Mi variable respuesta será el número de grillos, como variables fijas he
puesto parcela y año (solo hay dos años de muestreo), como variable
aleatoria he puesto el transecto (cada transecto se corresponde con una
única parcela). El problema es que no sé cómo definir el muestreo semanal
como medida repetida dentro de cada año, ¿podríais ayudarme?
Muchas Gracias de antemano por vuestra atención y colaboración,
Lucía
[[alternative HTML version deleted]]
_______________________________________________
R-help-es mailing list
R-help-es en r-project.org
https://stat.ethz.ch/mailman/listinfo/r-help-es