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]]
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
Muchas gracias Olivier, Voy a preparar el set de datos y lo envío, pero una rápida respuesta es que cada tratamiento tiene 4 parcelas y as u vez cada una de las 4 parcelas por tratamiento se subdivide en 4 subparcelas. Cada tratamiento tiene 16 observaciones en un momento determinado (por ejemplo la primera semana) y 16 observaciones en la segunda semana y así hasta 5 semanas. 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 <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]]
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 <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/> ------------ próxima parte ------------ Se ha borrado un adjunto en formato HTML... URL: <https://stat.ethz.ch/pipermail/r-help-es/attachments/20161101/e8829d66/attachment.html> ------------ próxima parte ------------ A non-text attachment was scrubbed... Name: mixed.csv Type: text/csv Size: 4791 bytes Desc: no disponible URL: <https://stat.ethz.ch/pipermail/r-help-es/attachments/20161101/e8829d66/attachment.bin>