[R-es] Modelo mixto en R
Olivier Nuñez
onunez en unex.es
Mie Nov 2 12:06:08 CET 2016
Bueno,
dando por hecho lo anterior, tu modelo me parece corecto para contestar a la pregunta que planteas.
Sin embargo, es posible que experimentes problemas de convergencia con estos datos.
Te aconsejo tipificar la variable "colecta":
datos$colecta=scale(datos$colecta)
Un saludo. Olivier
De: "Manuel Spínola" <mspinola10 en gmail.com>
Para: "Olivier Nuñez" <onunez en unex.es>
CC: "R" <r-help-es en r-project.org>
Enviados: Martes, 1 de Noviembre 2016 13:53:01
Asunto: Re: [R-es] Modelo mixto en R
Hola Olivier,
Adjunto los datos.
Manuel
El 31 de octubre de 2016, 3:29, Olivier Nuñez < onunez en unex.es > escribió:
Manuel,
no estoy seguro de entender el diseño.
Cada parcela recibe los 4 tratamientos (un tratamiento por subparcela de una misma parcela)?
Las "ocasiones" (supongo que instantes) son los mismos para cada parcela*tratamiento?
Si mandas un ejemplo de base de datos igual ayudaría a entender el diseño....
Un saludo. Olivier
----- Mensaje original -----
De: "Manuel Spínola" < mspinola10 en gmail.com >
Para: "R" < r-help-es en r-project.org >
Enviados: Sábado, 29 de Octubre 2016 14:56:02
Asunto: [R-es] Modelo mixto en R
Hola quisiera saber como parameterizar el siguiente diseño con lme4.
Tengo 4 tratamientos, con 4 parcelas cada uno, divididas a su vez en 4
subparcelas, y se han tomado 5 medidas repetidas (5 ocasiones) y se quiere
ver el cambio en la variable respuesta producido por los tratamientos a
través de las 5 ocasiones.
Mi intento con lmer es el siguiente:
mod <- lmer(variable respuesta ~ ocasión * tratamiento + (ocasión |
parcela/subparcela)
Es correcta esta forma de parameterización del diseño planteado?
Saludos,
Manuel Spínola
--
*Manuel Spínola, Ph.D.*
Instituto Internacional en Conservación y Manejo de Vida Silvestre
Universidad Nacional
Apartado 1350-3000
Heredia
COSTA RICA
mspinola en una.cr < mspinola en una.ac.cr >
mspinola10 en gmail.com
Teléfono: (506) 8706 - 4662
Personal website: Lobito de río < https://sites.google.com/site/lobitoderio/ >
Institutional website: ICOMVIS < http://www.icomvis.una.ac.cr/ >
[[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
--
Manuel Spínola, Ph.D.
Instituto Internacional en Conservación y Manejo de Vida Silvestre
Universidad Nacional
Apartado 1350-3000
Heredia
COSTA RICA
mspinola en una.cr
mspinola10 en gmail.com
Teléfono: (506) 8706 - 4662
Personal website: Lobito de río
Institutional website: ICOMVIS
[[alternative HTML version deleted]]
Más información sobre la lista de distribución R-help-es