Estimados integrantes de la lista. Disculpas por posteo cruzado. Estoy ajustando un modelo con lmer (lm4). La variable respuesta es un índice (ADI) que se midió en 3 áreas diferentes en 4 estaciones climáticas diferentes, así mis efectos fijos son area y estaciones climáticas. Cada área tiene 12 puntos de muestreo (Punto). En cada punto de muestreo el índice se midió repetidamente a diferentes horas, diferentes días y diferentes meses dentro de cada área y estación, pero no estoy interesado en la evolución del índice a través del tiempo. MI didea de parameterización son éstas: mod_adi_01 <- lmer(ADI ~ area + estacion + (1 | Punto/Mes/Dia/Hora), data = df_02, REML = FALSE) mod_adi_02 <- lmer(ADI ~ area + estacion + (1 | Punto) + (1 |Mes/Dia/Hora), data = df_02, REML = FALSE) Hay un de estas 2 alternativas correcta, o ninguna de las 2 lo son? Muchas gracias, Manuel -- *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]]
Javier Marcuzzi
2019-May-21 21:58 UTC
[R-es] Parameterización de modelo mixto (multilevel)
Estimado Manuel Spíndola Hay un problema, si no le interesa el tiempo, no lo ingrese en el modelo, si lo ingresa piense en que desea, tiempo como tiempo (H:M:S...) o que llueva a las 12 hs, a las 15 hs, como si las 15hs son totalmente independientes de las 12 hs. En muchas oportunidades mirando algunos ejemplos ayuda a escribir el modelo, por lo menos en mi caso luego de muchos años siempre me baso en algo escrito por otro, porque puede ser que yo piense en algo y descarte inconscientemente un modelo más adecuado. Javier Rubén Marcuzzi El dom., 19 may. 2019 a las 10:04, Manuel Spínola (<mspinola10 en gmail.com>) escribió:> Estimados integrantes de la lista. > > Disculpas por posteo cruzado. > > Estoy ajustando un modelo con lmer (lm4). > > La variable respuesta es un índice (ADI) que se midió en 3 áreas diferentes > en 4 estaciones climáticas diferentes, así mis efectos fijos son area y > estaciones climáticas. > > Cada área tiene 12 puntos de muestreo (Punto). > > En cada punto de muestreo el índice se midió repetidamente a diferentes > horas, diferentes días y diferentes meses dentro de cada área y estación, > pero no estoy interesado en la evolución del índice a través del tiempo. > > MI didea de parameterización son éstas: > > mod_adi_01 <- lmer(ADI ~ area + estacion + (1 | Punto/Mes/Dia/Hora), data > = df_02, REML = FALSE) > > mod_adi_02 <- lmer(ADI ~ area + estacion + (1 | Punto) + (1 > |Mes/Dia/Hora), data = df_02, REML = FALSE) > > Hay un de estas 2 alternativas correcta, o ninguna de las 2 lo son? > > Muchas gracias, > > Manuel > > -- > *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 >[[alternative HTML version deleted]]
Muchas gracias Javier. Voy a tener en consideración tu comentario. Manuel El mar., 21 may. 2019 a las 15:58, Javier Marcuzzi (< javier.ruben.marcuzzi en gmail.com>) escribió:> Estimado Manuel Spíndola > > Hay un problema, si no le interesa el tiempo, no lo ingrese en el modelo, > si lo ingresa piense en que desea, tiempo como tiempo (H:M:S...) o que > llueva a las 12 hs, a las 15 hs, como si las 15hs son totalmente > independientes de las 12 hs. En muchas oportunidades mirando algunos > ejemplos ayuda a escribir el modelo, por lo menos en mi caso luego de > muchos años siempre me baso en algo escrito por otro, porque puede ser que > yo piense en algo y descarte inconscientemente un modelo más adecuado. > > Javier Rubén Marcuzzi > > El dom., 19 may. 2019 a las 10:04, Manuel Spínola (<mspinola10 en gmail.com>) > escribió: > >> Estimados integrantes de la lista. >> >> Disculpas por posteo cruzado. >> >> Estoy ajustando un modelo con lmer (lm4). >> >> La variable respuesta es un índice (ADI) que se midió en 3 áreas >> diferentes >> en 4 estaciones climáticas diferentes, así mis efectos fijos son area y >> estaciones climáticas. >> >> Cada área tiene 12 puntos de muestreo (Punto). >> >> En cada punto de muestreo el índice se midió repetidamente a diferentes >> horas, diferentes días y diferentes meses dentro de cada área y estación, >> pero no estoy interesado en la evolución del índice a través del tiempo. >> >> MI didea de parameterización son éstas: >> >> mod_adi_01 <- lmer(ADI ~ area + estacion + (1 | Punto/Mes/Dia/Hora), data >> = df_02, REML = FALSE) >> >> mod_adi_02 <- lmer(ADI ~ area + estacion + (1 | Punto) + (1 >> |Mes/Dia/Hora), data = df_02, REML = FALSE) >> >> Hay un de estas 2 alternativas correcta, o ninguna de las 2 lo son? >> >> Muchas gracias, >> >> Manuel >> >> -- >> *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]]