Promedio bayesiano de modelos para la predicción del consumo máximo de oxígeno en deportistas con datos de no-ejercicio

Históricamente, los modelos de no-ejercicio para predecir el consumo máximo de oxígeno (VO2max) han sido construidos mediante regresión lineal frecuentista, usando técnicas estándar de selección de modelos. Sin embargo, existe incertidumbre acerca de la estructura estadística en el proceso de selección del modelo. En este estudio se propuso construir un modelo de no-ejercicio para predecir el VO2max en deportistas orientados al rendimiento, considerando la incertidumbre de modelo a través del Promedio Bayesiano de Modelos (BMA). Un objetivo adicional fue comparar la performance predictiva del BMA con las de los modelos derivados de varias técnicas frecuentistas usuales de selección de variables. Con tal fin, se implementó un submuestreo aleatorio estratificado repetido. Los datos incluyeron observaciones de la variable respuesta (en L·min-1), así como registros de Género, Deporte, Edad, Peso, Talla e Índice de masa corporal (BMI) (Edad = 22.1 ± 4.9 años, media ± SD; n = 272). Se propuso una clasificación de deportes con el objetivo de incluirla dentro del proceso de construcción del modelo: Combate, Juego, Resistencia 1 y Resistencia 2. El enfoque BMA se implementó en base a dos métodos: Occam's window y Composición de Modelo mediante el método de Monte Carlo con Cadenas de Markov (MC²). Se observaron discrepancias en la selección de variables entre los procedimientos frecuentistas. Ambos métodos de BMA produjeron resultados muy similares. Los modelos que incluyeron Género y las variables dummies para Resistencia 1 y Resistencia 2 acumularon virtualmente toda la probabilidad de modelo a posteriori. El Peso fue el predictor continuo con la más alta probabilidad de inclusión a posteriori (menor a 0.8). Las combinaciones de variables que involucraron predictores con un alto nivel de multicolinealidad fueron desacreditadas. Los modelos con sustancial contribución para el BMA presentaron un ajuste apreciable (R² ajustado menor a 0.8). Entre los modelos seleccionados por estrategias frecuentistas, el obtenido mediante el método de regresión por pasos (Stepwise regression method) con alfa igual a 0.05 fue el más respaldado por los datos, en términos de probabilidad de modelo a posteriori. En concordancia con la literatura, el BMA tuvo mejor performance predictiva de los datos fuera de la muestra que los modelos seleccionados por técnicas frecuentistas, medida por la cobertura del intervalo de predicción de 90 por ciento. La clasificación de deportes reveló resultados consistentes.

Flujos de carbono y nitrógeno bajo distintos sistemas de labranza en la región semiárida pampeana

Los contenidos de carbono y nitrógeno orgánicos del suelo están relacionados con su capacidad productiva. Esta relación es esperable que sea más estrecha en suelos de zonas áridas y semiáridas. La cantidad de carbono que puede secuestrar el suelo es, a la vez, de interés por la posibilidad de usar a los suelos como sumideros de carbono y mitigar el efecto invernadero. Esta tesis tiene como objetivo determinar los efectos del sistema de labranza sobre la dinámica del carbono y el nitrógeno en un suelo de una región semiárida. Se realizó un experimento de campo de larga duración en el que se contrastó un manejo tradicional, que empleaba un arado de disco en la labor principal, con un manejo tradicional, que empleaba un arado de disco en la labor principal, con un manejo de siembra directa. Los contenidos de humedad del suelo fueron superiores bajo siembra directa y mayores también los rendimientos de los cultivos que integraban la rotación siendo, consecuentemente, los aportes de carbono en forma de residuos al suelo también mayores. El flujo anual de carbono emitido a la atmósfera como C - CO2 fue similar entre sistemas de labranza, de lo que resultó un balance de carbono neutro bajo siembra directa y negativo bajo labranza con arado. La producción de C - CO2 aumentó con la cantidad de residuos presentes en el suelo y la temperatura y decrecía al incrementarse la cantidad de nitratos. La descomposición de los residuos vegetales y la liberación al suelo de su nitrógeno fue más rápida en el suelo labrado, lo que llevó a mayores niveles de nitratos a la siembra de algunos cultivos en el tratamiento con labranza en casi todos los cultivos integrantes de la rotación. Bajo siembra directa el suelo no perdió carbono durante el experimento mientras que bajo arado de disco el stock de carbono orgánico del suelo decreció. En regiones semiáridas la adopción de siembra directa determina balances de carbono menos negativos que bajo manejos con labranza, por incrementar las entradas de carbono al suelo y no afectar las salidas por respiración heterótrofa.

Estimación de heredabilidades y correlaciones genéticas para caracteres funcionales mediante análisis de supervivencia en Holstein Colombiano.

Las evaluaciones genéticas para caracteres funcionales requieren metodologías como el análisis de supervivencia (SA), que contemplen registros pertenecientes a hembras que no han presentado el evento de interes (descarte o prenez) al momento de la evaluacion. Así, los objetivos de esta tesis fueron: 1) analizar los factores que afectan el riesgo de descarte en vida productiva (VP), y de concepción en días abiertos (DA) e intervalo primer-último servicio (IP), 2) estimar parámetros genéticos de dispersión para los caracteres evaluados. Se emplearon 44652 registros de lactancia de 15706 vacas Holstein Colombiano, ajustando un modelo animal frágil Weibull para VP, con el número de partos (NP), producción de leche (PL) y edad al primer parto (EPP) como efectos fijos. En el caso de DA e IPU se emplearon 14789 servicios de 6205 vacas, ajustando un modelo de datos agrupados con efectos fijos de NP, EPP y número de servicios (SC). En todos los casos se incluyeron efectos aleatorios de animal y de hato. Tanto PL como NP tuvieron gran influencia en la VP, al igual que NP en DA e IPU. Se estimaron valores de h2 de 0,104, 0,086 y 0,1013 para VP, DA, e IPU, respectivamente, asi como correlaciones genéticas simples (rs) entre el nivel de PL con VP y DA de ƒ{0,03 y 0,47, respectivamente. Teniendo en cuenta tanto los valores de cría expresados como riesgo relativo, asi como la magnitud de rs, un aumento en PL conllevaría a disminución en el riesgo de descarte y aumento en el riesgo de concepción, dando lugar a aumentos en la VP y reducción en DA. Las h2s estimadas sugieren la existencia de variabilidad genetica para VP, DA e IPU en Holstein Colombiano, sirviendo de sustento para la implementación de evaluaciones genética de la raza, aprovechando las ventajas de metodologías como SA.

Modelo de ecuaciones estructurales para la evaluación de la calidad de una vacuna bovina

Los modelos de ecuaciones estructurales fueron introducidos en el contexto de la genética y la econometría y han sido extensamente empleados en las ciencias sociales. Sin embargo han sido poco empleados en el área de salud animal. En este trabajo se propone la construcción de un modelo de ecuaciones estructurales para la validación de un sistema de evaluación de la calidad de una vacuna contra la Rinotraqueítis infecciosa bovina. La dificultad para encontrar rodeos de bovinos seronegativos libres del virus BoHV-1 plantea la necesidad de desarrollar pruebas estandarizadas y armonizadas en animales de laboratorio que permitan definir la consistencia de lote a lote de la potencia de la vacuna, de un modo rápido y menos costoso para asegurar la presencia de productos confiables en el mercado. En este sentido, es de interés verificar la concordancia de la respuesta serológica en el animal de laboratorio con la protección a la enfermedad y a la infección en bovinos. Este sistema involucra variables latentes y relaciones complejas entre variables. El modelo de ecuaciones estructurales demostró ser una técnica útil para la validación de un modelo que emplea mediciones serológicas. en animales de laboratorio como alternativa a las pruebas clásicas de desafío en bovinos para la evaluación de la efectividad de las vacunas empleadas en el control y prevención de la enfermedad causada por un herpesvirus. En este trabajo se pudo comprobar que las variables que representaban los conceptos de interés - protección a la infección, protección a la enfermedad y calidad determinada por el modelo animal en el laboratorio (medición alternativa de la efectividad de la vacuna) explicaron adecuadamente la variabilidad de las variables observables. Se destacaron las variables área bajo las curvas de excreción y área bajo la curva de severidad para representar las variables protección a la infección y protección a la enfermedad, respectivamente, y las variables serológicas medidas por ELISA en bovinos y cobayos representando la calidad por laboratorio. Por otra parte, el modelo permitió la incorporación de correlación de errores de las variables serológicas, que hubiese sido imposible con otras técnicas multivariadas. Asimismo, la técnica de análisis permitió explorar por partes las relaciones entre variables para detectar la fuente de falta de ajuste. Se pudo comprobar el buen ajuste de las relaciones de las variables latentes protección a la infección y calidad por laboratorio, obteniéndose una alta correlación entre ambos conceptos. Teniendo en cuenta la alta fiabilidad de las variables serológicas medidas por ELISA tanto en bovinos como en cobayos, se puede inferir que la medición de estas variables de laboratorio puede ser útil en la predicción de la protección de la vacuna. El efecto de la protección a la infección sobre la protección a la enfermedad fue significativo pero débil. La parte estructural del modelo presentó un ajuste adecuado, validando así la teoría inicial de interrelaciones entre las variables latentes. Se plantea la posibilidad de extender los modelos de ecuaciones estructurales para la evaluación de vacunas que involucran otros virus que afectan el ganado bovino.

Recomendaciones de fertilización en maíz según la forma algebraica del modelo predictivo empleado

p.155-164

Flujos de carbono y nitrógeno bajo distintos sistemas de labranza en la región semiárida pampeana

Los contenidos de carbono y nitrógeno orgánicos del suelo están relacionados con su capacidad productiva. Esta relación es esperable que sea más estrecha en suelos de zonas áridas y semiáridas. La cantidad de carbono que puede secuestrar el suelo es, a la vez, de interés por la posibilidad de usar a los suelos como sumideros de carbono y mitigar el efecto invernadero. Esta tesis tiene como objetivo determinar los efectos del sistema de labranza sobre la dinámica del carbono y el nitrógeno en un suelo de una región semiárida. Se realizó un experimento de campo de larga duración en el que se contrastó un manejo tradicional, que empleaba un arado de disco en la labor principal, con un manejo tradicional, que empleaba un arado de disco en la labor principal, con un manejo de siembra directa. Los contenidos de humedad del suelo fueron superiores bajo siembra directa y mayores también los rendimientos de los cultivos que integraban la rotación siendo, consecuentemente, los aportes de carbono en forma de residuos al suelo también mayores. El flujo anual de carbono emitido a la atmósfera como C - CO2 fue similar entre sistemas de labranza, de lo que resultó un balance de carbono neutro bajo siembra directa y negativo bajo labranza con arado. La producción de C - CO2 aumentó con la cantidad de residuos presentes en el suelo y la temperatura y decrecía al incrementarse la cantidad de nitratos. La descomposición de los residuos vegetales y la liberación al suelo de su nitrógeno fue más rápida en el suelo labrado, lo que llevó a mayores niveles de nitratos a la siembra de algunos cultivos en el tratamiento con labranza en casi todos los cultivos integrantes de la rotación. Bajo siembra directa el suelo no perdió carbono durante el experimento mientras que bajo arado de disco el stock de carbono orgánico del suelo decreció. En regiones semiáridas la adopción de siembra directa determina balances de carbono menos negativos que bajo manejos con labranza, por incrementar las entradas de carbono al suelo y no afectar las salidas por respiración heterótrofa.

Estimación de heredabilidades y correlaciones genéticas para caracteres funcionales mediante análisis de supervivencia en Holstein Colombiano.

Las evaluaciones genéticas para caracteres funcionales requieren metodologías como el análisis de supervivencia (SA), que contemplen registros pertenecientes a hembras que no han presentado el evento de interes (descarte o prenez)al momento de la evaluacion. Así, los objetivos de esta tesis fueron: 1)analizar los factores que afectan el riesgo de descarte en vida productiva (VP), y de concepción en días abiertos (DA)e intervalo primer-último servicio (IP), 2)estimar parámetros genéticos de dispersión para los caracteres evaluados. Se emplearon 44652 registros de lactancia de 15706 vacas Holstein Colombiano, ajustando un modelo animal frágil Weibull para VP, con el número de partos (NP), producción de leche (PL)y edad al primer parto (EPP)como efectos fijos. En el caso de DA e IPU se emplearon 14789 servicios de 6205 vacas, ajustando un modelo de datos agrupados con efectos fijos de NP, EPP y número de servicios (SC). En todos los casos se incluyeron efectos aleatorios de animal y de hato. Tanto PL como NP tuvieron gran influencia en la VP, al igual que NP en DA e IPU. Se estimaron valores de h2 de 0,104, 0,086 y 0,1013 para VP, DA, e IPU, respectivamente, asi como correlaciones genéticas simples (rs)entre el nivel de PL con VP y DA de ƒ{0,03 y 0,47, respectivamente. Teniendo en cuenta tanto los valores de cría expresados como riesgo relativo, asi como la magnitud de rs, un aumento en PL conllevaría a disminución en el riesgo de descarte y aumento en el riesgo de concepción, dando lugar a aumentos en la VP y reducción en DA. Las h2s estimadas sugieren la existencia de variabilidad genetica para VP, DA e IPU en Holstein Colombiano, sirviendo de sustento para la implementación de evaluaciones genética de la raza, aprovechando las ventajas de metodologías como SA.

Un modelo de rendimiento para sorgo granífero (Sorghum Bicolor (L.) Moench.)

p.113-123

Inferencia bayesiana sobre los parámetros de dispersión genéticos y ambientales en modelos animales con efectos maternos

Los modelos 'modelos animales con efectos maternos' (MAM) son modelos lineales mixtos que se utilizan para ajustar registros de caracteres bajo la influencia de efectos maternos. Uno de los desafíos más importantes en el marco de los MAM es la estimación de los parámetros de dispersión o 'componentes de (co) varianza' (CVC). En esta tesis se introducen desde una perspectiva bayesiana contribuciones teóricas y metodológicas con relación a la estimación de CVC para MAM sujetos a estructuras de covarianza novedosas. En primer lugar, se describe una implementación del análisis bayesiano jerárquico vía el algoritmo del muestreo de Gibbs. Luego, se considera una especificación conjugada diferente para la distribución a priori de la matriz de covarianza genética, basada en la distribución Wishart invertida generalizada, y se presenta una estrategia para determinar los correspondientes hiperparámetros. Esta estrategia fue comparada contra otras especificaciones a priori mediante un estudio de simulación estocástica, y produjo estimaciones precisas de los parámetros genéticos, con menores errores estándares y mejor tasa de convergencia. En segundo lugar, se presenta una formulación alternativa del MAM que incluye un parámetro de correlación ambiental entre pares de observaciones madre-progenie, y se desarrolla un procedimiento de estimación basado en un algoritmo de muestreo por grilla. El procedimiento fue programado y ejecutado exitosamente, y se obtuvo la primera estimación del parámetro de correlación con datos de campo para peso al destete en bovinos de carne. Por último, se considera el problema de la estimación de CVC en una población multirracial, donde en general es necesario especificar una estructura de covarianza heterogénea para los valores de cría. En particular, se demuestra que el modelo basado en la descomposición de la matriz de covarianza genética es equivalente al que deriva de la teoría genética cuantitativa. Además, se extiende el modelo para incluir efectos maternos y se describe la implementación de un análisis bayesiano jerárquico con el objetivo de estimar los CVC. El procedimiento fue implementado con éxito en datos experimentales de peso al destete y se obtuvieron por primera vez estimaciones para el conjunto completo de CVC.

Modelos animales para caracteres con efectos maternos, madres no identificadas e información incompleta

La evaluación genética para caracteres de crecimiento pre - destete requiere ajustar modelos animales con efectos maternos (MAM). Tanto la estimación paramétrica de la variabilidad como la evaluación genética mediante MAM son realizadas empleando datos de campo, muchos de los cuales no poseen información completa para todas las variables explicativas maternas. Es común no contar con la identificación de madres (biológicas y/o receptoras), de abuelas maternas y, consecuentemente, de la edad de la madre (EM). Este problema es bien marcado en razas compuestas como Brangus y Braford que tienen políticas para registrar animales de pedigrí "abierto". Además, no existe un consenso sobre cuál es el mejor modelo de predicción, y existen interrogantes sobre la magnitud de los componentes de (co) varianza genético-aditivos y ambientales del modelo de evaluación. La primera investigación de esta tesis consistió en la estimación, mediante métodos bayesianos de los parámetros de dispersión en MAMs con distintas estructuras de (co) varianza, para datos de peso al destete de animales Angus de pedigrí. El análisis se caracterizó por la originalidad en los muestreos de las distribuciones marginales posteriores de las covarianzas genéticas aditivas y de la correlación entre los efectos ambientales maternos permanentes de una vaca y sus hijas también madres. Con el objeto de especificar correctamente la fracción aditiva de las (co) varianzas cuando se desconocen las madres y/o abuelas maternas de los animales con datos, en otro capítulo se desarrollaron MAMs equivalentes que no requieren alargar los vectores de los valores de cría con madres o abuelas fantasmas. Finalmente, se desarrolló un modelo mixto que atenúa el sesgo por error de medición clásico en el efecto EM, e introduce splines penalizadas y una estructura de (co) variación autoregresiva de orden 1 para suavizar las covarianzas residuales Este modelo es apropiado para ajustar datos de animales nacidos por transplante embrionario con madres receptoras desconocidas

Sistemas de evaluación de tierras y elaboración de modelos de aptitud de uso agrícola, para distintos escenarios climáticos, en un sector de la subregión pampa arenosa (Provincia de Buenos Aires, Argentina)

La Pampa Arenosa ha sido escenario de cambios en el uso del territorio que respondieron principalmente al aumento de las precipitaciones a partir de la década del 70. La evaluación de las tierras es una etapa crítica en la planificación del uso sustentable. Por este motivo, se analizaron distintos sistemas de evaluación de tierras y se desarrollaron modelos expertos que consideren los factores ambientales heredados y las variaciones climáticas, para el sector de dunas longitudinales de la Pampa Arenosa en la Provincia de Buenos Aires, a escala 1: 50.000, considerando a los partidos de Nueve de Julio, Carlos Casares, Pehuajó y Trenque Lauquen. Las tierras fueron clasificadas por Capacidad de Uso, Indice de Productividad (IP) y se generaron sistemas expertos, utilizando el programa ALES, para los tipos de utilización de las tierras (TUTs) : maíz, soja y trigo. La homogeneidad de las series climáticas de precipitaciones se determinó mediante el test de Rachas. La aplicación del test de Pettitt permitió identificar la existencia de un cambio abrupto en las precipitaciones y el este de Mann Kendall mostró una tendencia creciente en relación a la precipitación anual. Las tierras con severas (clase III) y muy severas limitaciones (clase IV), fueron las más frecuentes ocupando el 42,6 por ciento y 29,8 por ciento respectivamente del área. Se comprobó que el IP de las tierras se incrementó con el aumento de las precipitaciones, alcanzando su máxima expresión climática en el período posterior al cambio abrupto. Las tierras de moderada capacidad productiva con valores de IP entre 65-51 ocuparon la mayor superficie de área de estudio. Los modelos expertos según los TUTs presentaron una aptitud de uso de las tierras variable, condicionada por la capacidad de retención hídrica de los suelos. Los modelos expertos fueron sensibles a las variaciones climáticas y el cambio abrupto en las precipitaciones.

Modelos explicativos y predictivos de la productividad, disponibilidad y calidad de forraje basados en sensores remotos

La productividad primaria neta aérea (PPNA) y la biomasa y la calidad forrajeras son clave en los sistemas pastoriles. Dada su variación espacial y temporal, contar con descripciones de dicha variación y con herramientas de monitorización facilita y perfecciona su manejo. Los sensores remotos permiten estimar frecuente y detalladamente la PPNA pero no la biomasa y la calidad. El objetivo de esta tesis es mejorar el conocimiento sobre la PPNA y las posibilidades de monitorización de los recursos forrajeros en general y los de la Pampa Deprimida en particular. Primero, se presentan patrones de variación de la PPNA de los principales recursos forrajeros y modelos de prospección de la PPNA primaveral para distintas unidades de vegetación. Segundo, se presentan modelos de estimación de la biomasa y calidad forrajera basados en sensores remotos y generados mediante ensayos en condiciones controladas y a campo. La variación espacial de la PPNA fue explicada por la precipitación a escala regional, los suelos a escala subregional y los suelos y el pastoreo a escala local. Entre recursos, las pasturas de loma tuvieron mayor PPNA que los pastizales de bajo; dicha diferencia aumentó con la precipitación media. Entre escalas, la variación espacial fue mayor a escala regional y la temporal fue mayor a escala local. La PPNA primaveral se relacionó con variables ambientales de meses previos diferentes según los pastizales estuvieran dominados por especies invernales o estivales, y según las pasturas fueran de suelos profundos o someros. Las estimaciones de la biomasa y calidad forrajera fueron satisfactorias en general, aunque imprecisas bajo ciertas condiciones experimentales como sequía, fertilización, y presencia de biomasa senescente, que requirieron modelos específicos. La estimación de la biomasa total a campo mediante un modelo basado en estimaciones independientes de la biomasa fotosintética y la no fotosintética resultó precisa a lo largo de un año. Los resultados aportan información acerca de los controles ambientales de la PPNA en pastizales subhúmedos y representan un avance en la estimación de la biomasa y calidad forrajera mediante sensores remotos

Sistema de observaciones fenológicas y fenométricas para la calibración de modelos fisiológicos de rendimiento

p.99-108

Sistemas de evaluación de tierras y elaboración de modelos de aptitud de uso agrícola, para distintos escenarios climáticos, en un sector de la subregión pampa arenosa (Provincia de Buenos Aires, Argentina)

La Pampa Arenosa ha sido escenario de cambios en el uso del territorio que respondieron principalmente al aumento de las precipitaciones a partir de la década del 70. La evaluación de las tierras es una etapa crítica en la planificación del uso sustentable. Por este motivo, se analizaron distintos sistemas de evaluación de tierras y se desarrollaron modelos expertos que consideren los factores ambientales heredados y las variaciones climáticas, para el sector de dunas longitudinales de la Pampa Arenosa en la Provincia de Buenos Aires, a escala 1: 50.000, considerando a los partidos de Nueve de Julio, Carlos Casares, Pehuajó y Trenque Lauquen. Las tierras fueron clasificadas por Capacidad de Uso, Indice de Productividad (IP)y se generaron sistemas expertos, utilizando el programa ALES, para los tipos de utilización de las tierras (TUTs): maíz, soja y trigo. La homogeneidad de las series climáticas de precipitaciones se determinó mediante el test de Rachas. La aplicación del test de Pettitt permitió identificar la existencia de un cambio abrupto en las precipitaciones y el este de Mann Kendall mostró una tendencia creciente en relación a la precipitación anual. Las tierras con severas (clase III)y muy severas limitaciones (clase IV), fueron las más frecuentes ocupando el 42,6 por ciento y 29,8 por ciento respectivamente del área. Se comprobó que el IP de las tierras se incrementó con el aumento de las precipitaciones, alcanzando su máxima expresión climática en el período posterior al cambio abrupto. Las tierras de moderada capacidad productiva con valores de IP entre 65-51 ocuparon la mayor superficie de área de estudio. Los modelos expertos según los TUTs presentaron una aptitud de uso de las tierras variable, condicionada por la capacidad de retención hídrica de los suelos. Los modelos expertos fueron sensibles a las variaciones climáticas y el cambio abrupto en las precipitaciones.

Eficiencia in vitro de antagonistas de Rhizoctonia solani y Sclerotium rolfsii y análisis camparativo de distintos modelos de crecimiento

p.59-65