[R-es] Heterocedasticidad MLEM

[Javier Marcuzzi]

Estimado Jimmy Erney Reyes Velasco

No leí sus datos, justo ahora no tengo tiempo pero supe tener ese tipo de problemas.

No leí la última versión de lme4, pero para ciertas cosas nlme permite mejorar las el análisis cuándo hay curvas de producción que tienen que ver con la relación “vivo” en contraste a “tierra” y “clima”.

En lo personal me gustan las regresiones aleatorias, pero habría que leer sus datos y que está buscando.

Hay un paquete “agridat” si no recuerdo mal, lleno de ejemplos y creo que el autor también es de América del Sur, posiblemente ahí tiene alguien muy apropiado para su problema.

No puedo darle más ayuda en este momento.

Saludos.

Javier Marcuzzi

> El 26 jun 2026, a las 2:28 p. m., Jimmy Erney Reyes Velasco <jimmyreyesvelasco using gmail.com> escribió:
> 
> Buenos días, estimados, el siguiente correo lo envío con el objetivo de solicitarles un consejo, tengo un problema desde hace bastante y aún no lo he podido solucionar, es un poco complejo. Esta es la explicación:
> Hice un muestreo con una especie de planta en el periodo del 2022 hasta el 2024, medí la masa seca de todas las hojas de 32 plantas en 7 periodos, más especificamente, lo que produjo entre los días 0-193, 193-288, 288-389, 389-475, 475-556, 556-660, 660-731. Estas producciones se acumularon a través del tiempo, con el objetivo de estimar la producción de esta especie. El primer problema que surje es la autocorrelación de los residuos debido a que las mismas plantas fueron medidas a través del tiempo, esto lo intenté corregir con un modelo de efectos mixtos. El siguiente problema que es el que más me ha tenido tiempo en esto es la heterocedasticidad de todos los modelos que intento usar, intenté corregirlo con un modelo lineal de efectos mixtos robustos, sin embargo es difícil corregirlo por los siguientes problemas:
> Las plantas al parecer según mis resultados inhiben su crecimiento entre finales de noviembre y finales de febrero, ya que son los meses que marcan la temporada seca y el momento de mayor probabilidad de heladas. Al parecer es un mecanismo de adaptación. esto hace que se forme "microasíntotas" en el periodo de muestreo, mi hipótesis es que estos periodos de latencia de crecimiento es lo que hacen que ninguno de los modelos cumpla con el supuesto de homocedasticidad. Me interesa mucho la modelación de efectos mixtos, porque así puedo explicar que las tendencias propias de los individuos hace que la producción de masa seca foliar aporta a la variabilidad del modelo, a su vez también quiero demostrar que la temporada seca afecta la producción de masa seca foliar de la roseta de estas plantas, sin embargo ningún modelo hasta el momento se ajusta ni me permite probar mis hipótesis, el otro inconveniente es que no tengo datos propiamente dichos para ese tiempo del microclima, lo que hace difícil ingresar esas variables al modelo.
> Si alguno supiera una forma de solucionar el problema o forma de manejar el heterocedasticidad y la normalidad de los residuos de estos modelos, se los agradecería muchísimo. Esto permitirá culminar con la publicación de estos resultados que son bastante importantes.
> Este muestreo fue hecho en Colombia.
> Adjuntaré los datos por si alguno le gustaría echarle un vistazo.
> De antemano muchísimas gracias por su atención
> Saludos
> Jimmy
> <Producción de masa seca foliar.xlsx>_______________________________________________
> R-help-es mailing list
> R-help-es using r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-help-es