[R-es] Heterocedasticidad MLEM
Jimmy Erney Reyes Velasco
j|mmyreye@ve|@@co @end|ng |rom gm@||@com
Vie Jun 26 19:28:56 CEST 2026
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
------------ próxima parte ------------
Se ha borrado un adjunto en formato HTML...
URL: <https://stat.ethz.ch/pipermail/r-help-es/attachments/20260626/32f30268/attachment-0001.html>
------------ próxima parte ------------
A non-text attachment was scrubbed...
Name: Producción de masa seca foliar.xlsx
Type: application/vnd.openxmlformats-officedocument.spreadsheetml.sheet
Size: 23050 bytes
Desc: no disponible
URL: <https://stat.ethz.ch/pipermail/r-help-es/attachments/20260626/32f30268/attachment-0001.xlsx>
Más información sobre la lista de distribución R-help-es