Resumen: El principal objetivo del presente estudio es estimar los factores de seguridad y la probabilidad de falla planar en los taludes de la Comuna Universitaria de la ciudad de Manizales (Colombia). Se realiza un procesamiento estadístico de los datos recolectados durante 14 años de estudios geotécnicos. Se aplica la ecuación de talud infinito en una simulación de Monte Carlo para analizar escenarios típicos de entrada como cohesión, ángulo de fricción, nivel de aguas freáticas e inclinación de estratos con presencia de ceniza y llenos antrópicos. En la simulación se generan 5000 iteraciones para cada variable, obteniendo como resultado la representación de una herramienta para detección de condiciones que generan inestabilidad en los taludes de la comuna. Se concluye que para factores de seguridad aceptables calculados por métodos determinísticos, se pueden tener niveles de amenaza altos en función de probabilidad de falla por lo que los análisis probabilísticos se consideran conservadores y útiles en la detección temprana de desastres.

Abstract: The main objective of this study is to estimate the safety factors and the probability of planar failure in the slopes of the Comuna Universitaria of the city of Manizales (Colombia). Geotechnical data gathered over a period of 14 years is processed statistically. The infinite slope equation is applied through Monte Carlo simulations to assess typical input scenarios including cohesion, friction angle, groundwater level, and inclination of layers with the presence of ash and anthropogenic fills. The simulation involved 5000 iterations for each variable, resulting in a tool for detecting conditions causing instability in the slopes of the Comuna Universitaria. It is concluded that for acceptable safety factors calculated deterministically, it is possible to obtain high risk levels with relation to failure probability function, making probabilistic analyses conservative and useful for early warning detection.

Resumen La regulación de Basilea sobre grandes exposiciones está en proceso de implementación y, aun que es una regulación que reduce el límite del tamaño para las exposiciones de gran tamaño en los bancos y está orientada a un mejor control de los parámetros de riesgo que ocasiona que los portafolios se inclinen más a tener pérdidas grandes por concentración de crédito, es importante conocer la porción de riesgo de crédito y capital que lleva implícita. En este trabajo se determinan, mediante una simulación Montecarlo a través de un modelo de riesgo de crédito, los add-ons implícitos (o añadidos implícitos) debidos a la concentración de crédito, y se comparan con los requisitos actuales de capital. Al mismo tiempo, se analiza el marco completo de Basilea para entender cómo se aborda el riesgo de crédito de concentración, esto es, de forma integral, más que mediante un suplemento particular de capital.

Abstract The Basel large exposures standard, addressed to banks, is in the process of being implemented, and although the new rule is reducing the limit for the large credit exposures in banks, and is geared to better control the portfolio risk parameters that make a portfolio more prone to losses due to credit concentration, it is important to know the share of risk and capital implicit to this new rule. In this work, the implicit add-ons for credit risk concentration are determined through a Monte Carlo credit risk model, and the results are compared with current capital requirements. The author also analyzed the complete Basel framework to understand how the concentration risk is addressed in an integrated approach rather than with a specific capital supplement.

Resumen: En el presente artículo se presenta el estudio de la estimación de demanda de una estación de carga de vehículos eléctricos basado en el uso de simulación de Montecarlo. La modelación del sistema eléctrico se lo realiza a través del software de PowerFactory, por otra parte, para el desarrollo de las simulaciones de Montecarlo y el procesamiento de la información, se lo realiza a través del uso del software de Python. El análisis en general está enfocado en determinar el impacto de la integración de vehículos eléctricos en la red eléctrica, de forma que se generan escenarios que emulen la operación de la electrolinera, teniendo en cuenta la demanda de vehículos eléctricos dentro de sus puntos de carga, de forma que se generan escenarios aleatorios considerando el número de vehículos conectados y el porcentaje de carga de los vehículos. Otro aspecto que se considera es la integración de sistemas de generación fotovoltaica, de forma que se realiza el análisis de su impacto en la operación de la red eléctrica teniendo en cuenta la operación aleatoria de la electrolinera.

Abstract: This paper presents the study of the demand estimation of an electric vehicle charging station based on the use of Monte Carlo simulation. The modeling of the electrical system is done through the PowerFactory software, on the other hand, for the development of the Monte Carlo simulations and the processing of the information, it is done through the use of Python software. The analysis in general is focused on determining the impact of the integration of electric vehicles in the electric grid, so that scenarios that emulate the operation of the electric station are generated, taking into account the demand of electric vehicles within their charging points, so that random scenarios are generated considering the number of connected vehicles and the percentage of vehicle charging. Another aspect that is considered is the integration of photovoltaic generation systems, so that an analysis of their impact on the operation of the electric grid is performed, taking into account the random operation of the electric station.

Resumen La carga no controlada de vehículos eléctricos plantea un gran desafío para los operadores de redes de distribución y los planificadores de sistemas de potencia. En lugar de focalizarse en el control de esta carga no controlada, se propone un modelo que utiliza variables de análisis de contingencias para calcular la capacidad de potencia necesaria en el sistema de potencia. Se emplea la variable de potencia no servida para evaluar la cantidad de potencia de carga no cubierta en cada barra del sistema, seguido del cálculo de la capacidad de potencia adicional requerida, utilizando un sistema fotovoltaico y de almacenamiento y otra alternativa de generación constante en el sistema de potencia IEEE de 14 barras con información sobre algunos vehículos eléctricos y la carga diaria en el sistema de potencia de Perú. Los resultados obtenidos en el sistema de potencia con generación distribuida muestran que no hay presencia de potencia no servida, corroborando el éxito de la metodología utilizada. Este modelo brinda herramientas tanto a los operadores de redes de distribución como a los planificadores de sistemas de potencia, reduciendo el impacto en el sistema de potencia de los vehículos eléctricos y aportando una metodología aplicable a otros sistemas de distribución eléctrica con cargas no controladas.

Abstract The uncontrolled charging of electric vehicles poses a great challenge for distribution network operators and power system planners. Instead of focusing on controlling this uncontrolled load, a model that uses contingency analysis variables to calculate the power capacity needed in the power system is proposed. The unserved power variable is used to evaluate the amount of uncovered load power at each bus of the system, followed by the calculation of the additional power capacity required using a photovoltaic and storage system and another constant generation alternative in the 14-bus IEEE power system with information on some electric vehicles and daily load in the power system of Peru. The results obtained in the power system with distributed generation, the absence of unserved power, corroborate the success of the methodology used. This model provides tools to both distribution network operators and power system planners, reducing the impact on the power system of electric vehicles and providing a methodology applicable to other electric distribution systems with uncontrolled loads.

Abstract The purpose of this study is to generate the simulation and validation of a cobalt 60 beam of a radiotherapy equipment under reference conditions, using evaluation parameters corresponding to commercial planning systems. For this, we used dose profile and percentage depth dose (PDD) data for a field size of 10 cm x 10 cm at a surface source distance (DFS) of 80 cm from the water phantom. The program used to simulate the radiation beam was Geant4, including a collimation system and the phantom. When comparing the simulation results with the measured data, a difference of 3.50 % was obtained in the region defined between the maximum dose and 50 % in the PDD curve, as well as a maximum deviation of 0.96 mm in buildup region in the PDD. In addition, there is a deviation of 1.42 mm in the penumbra region for the dose profile. For the zone outside the central axis in the dose profile and in the zone of low doses, they were 2.41 % and 2.56 %, respectively. Based on the results obtained, this paper validates the simulated conditions.

Resumen El propósito de este estudio es generar la simulación y validación de un haz de cobalto 60 de un equipo de radioterapia en condiciones de referencia, mediante el uso de parámetros de evaluación correspondientes a sistemas de planificación comerciales. Para esto, se utilizaron datos de perfiles de dosis y de porcentaje de dosis en profundidad (PDD) para un tamaño de campo de 10 cm x 10 cm a una distancia fuente superficie (DFS) de 80 cm del maniquí de agua. El programa utilizado para la simulación fue el Geant4, incorporando el haz de radiación, un sistema de colimación y el maniquí. Al comparar los resultados de la simulación con los datos medidos, se obtuvo una diferencia de 3.50 % en la región definida entre el máximo de dosis y el 50 % en la curva de PDD, así como una desviación máxima de 0.96 mm en la zona de acumulación del PDD, además de una desviación de 1.42 mm en la zona de penumbra para el perfil dosis. Para la zona fuera del eje central en el perfil de dosis y en la zona de bajas dosis fueron del 2.41 % y del 2.56 %, respectivamente. Por los resultados obtenidos, se realizó la validación de las condiciones simuladas.

Resumen Se utilizan procedimientos no paramétricos para agregar flexibilidad a los modelos. Se propusieron tres modelos de respuesta al ítem no paramétricos, pero no se compararon directamente: Kernel smoothing (KS-IRT); la curva davidiana (DC-IRT); y el modelo bayesiano de Rasch semiparamétrico (SP-Rasch). El objetivo principal del presente estudio es comparar el desempeño de estos procedimientos en la recuperación de puntajes verdaderos simulados, utilizando puntajes de suma como puntos de referencia. El objetivo secundario es comparar su desempeño en términos de equivalencia práctica con datos reales. En general, los resultados muestran que, a excepción de DC-IRT, que es el modelo con peor desempeño, todos los otros modelos presentan resultados bastante similares a los obtenidos cuando se utilizan sumatorios. Estos resultados son seguidos por una discusión con implicaciones prácticas y recomendaciones para estudios futuros. directamente KSIRT KS IRT (KS-IRT) DCIRT DC (DC-IRT) SPRasch. SPRasch SP . (SP-Rasch) simulados referencia reales general DCIRT, IRT, DC-IRT sumatorios futuros (KS-IRT (DC-IRT (SP-Rasch

Abstract Nonparametric procedures are used to add flexibility to models. Three nonparametric item response models have been proposed, but not directly compared: the Kernel smoothing (KS-IRT); the Davidian-Curve (DC-IRT); and the Bayesian semiparametric Rasch model (SP-Rasch). The main aim of the present study is to compare the performance of these procedures in recovering simulated true scores, using sum scores as benchmarks. The secondary aim is to compare their performances in terms of practical equivalence with real data. Overall, the results show that, apart from the DC-IRT, which is the model that performs the worse, all the other models give results quite similar to those when sum scores are used. These results are followed by a discussion with practical implications and recommendations for future studies. proposed compared KSIRT KS IRT (KS-IRT) DavidianCurve Davidian Curve DCIRT DC (DC-IRT) SPRasch. SPRasch SP . (SP-Rasch) benchmarks data Overall DCIRT, IRT, DC-IRT worse studies (KS-IRT (DC-IRT (SP-Rasch

Resumo Procedimentos não paramétricos são usados para adicionar flexibilidade aos modelos. Três modelos não paramétricos de resposta ao item foram propostos, mas não comparados diretamente: o Kernel smoothing (KS-IRT); a Curva Davidiana (DC-IRT); e o modelo semiparamétrico Rasch Bayesiano (SP-Rasch). O objetivo principal do presente estudo é comparar o desempenho desses procedimentos na recuperação de escores verdadeiros simulados, utilizando escores de soma como benchmarks. O objetivo secundário é comparar seus desempenhos em termos de equivalência prática com dados reais. De forma geral, os resultados mostram que, além do DC-IRT, que é o modelo que apresenta o pior desempenho, todos os outros modelos apresentam resultados bastante semelhantes aos de quando se usam somatórios. Esses resultados são seguidos de uma discussão com implicações práticas e recomendações para estudos futuros. propostos diretamente KSIRT KS IRT (KS-IRT) DCIRT DC (DC-IRT) SPRasch. SPRasch SP . (SP-Rasch) simulados benchmarks reais geral DCIRT, IRT, DC-IRT somatórios futuros (KS-IRT (DC-IRT (SP-Rasch

Resumen Este artículo tiene por objetivos: (i) presentar, (ii) probar y (iii) validar un conjunto de técnicas matemáticas para predecir el número de incendios (No), la superficie incendiada (A), y la superficie promedio incendiada (AM), en bases anuales a escalas regionales. Una fuente de datos extensa proveniente de los bosques de coníferas de Durango se usó para ajustar (1970-2011) y validar (2012-2016) las técnicas de modelaje. La mayoría de los modelos probabilísticos y estocásticos explican menos del 70% de la varianza de los incendios. Sin embargo, el método de teleconecciones que emplea variables a escala oceánica y a escala local del hidroclima aumentó el nivel de precisión hasta cerca del 80% de la varianza total para ambas bases de datos. Los resultados muestran la complejidad de los factores que interactúan, incluyendo los procesos estocásticos y los físicos que hacen que las técnicas de predicción sean modificadas substancialmente. Por estas razones, se propone un modelo más físicamente basado y también un modelo conceptual más comprensivo. El modelo conceptual requiere del uso de todas las técnicas presentadas en este reporte, incluyendo los modelos probabilísticos, estocásticos y físicos en la parameterización de los diferentes submodelos. El uso del modelo conceptual además de sus técnicas de simulación Monte Carlo extraerían los escenarios futuros mas probables en la predicción del riesgo de incendios forestales severos en bosques templados del mundo.

Abstract This paper sets the following objectives: (i) presenting, (ii) testing, and (iii) evaluating a set of mathematical techniques to forecast the number of forest wildfires (No), the burned area (A), and the mean burned area (MA), on annual basis at regional scales. A comprehensive wildfire data set for coniferous forests of the State of Durango, Mexico was used to fit (1970-2011) and to validate (2012-2016) some modeling techniques. Most tested probabilistic and stochastic models hardly explain 70% of the wildfire variance. However, the teleconnection approach using a combination of large scale and local hydroclimate anomalies better predicted both data sets; explaining nearly 80% of the wildfire variance for fitting and for validating models. Results stress the complexity of interactive factors including the stochastic and underlying physical process that makes the prediction of wildfires losing precision and they should be further considered in future conceptual models. Therefore, proposing a more physical-based and conceptual models including Montecarlo models is an integral component of this paper; with the goal of increasing prediction capabilities and assisting decision-makers on the prevention activities inherent to better control wildfires. This proposed conceptual model stresses the need for using the probabilistic, stochastic and physical techniques to improve sub-model parameterization. Furthermore, the use of Monte Carlo simulation techniques would extract the most likely future scenarios for predicting the risk of high-severity wildfire regimes in temperate forests elsewhere.

Resumen Los mecanismos de transporte de radiación ionizante en presencia de campos magnéticos pueden ser descritos mediante formulaciones basadas en la aplicación de la ecuación principal de Boltzmann. Sin embargo, cuando el transporte no es en vacío sino en medios materiales dispersores la dificultad en el cálculo analítico de la trayectoria de los electrones se incrementa considerablemente debido al carácter aleatorio de los procesos de interacción de la radiación con el medio dispersor. En este contexto, el presente trabajo propone y describe una metodología para caracterizar mediante técnicas numéricas, como la simulación Monte Carlo, el efecto de campos magnéticos intensos, presentes típicamente en los dispositivos de MRI-LINAC, sobre la trayectoria de partículas cargadas, propagándose en un medio material comúnmente presente en situaciones de radioterapia. Los resultados obtenidos para el aire como medio dispersor, indican que las variaciones de las trayectorias de los electrones dependen conjuntamente de la intensidad y dirección del campo magnético y de las interacciones por colisión con el medio material. El aumento de la intensidad de campo magnético incrementa la curvatura de las trayectorias y disminuye el desplazamiento, o corrimiento, en la dirección radial.

Abstrac Ionizing radiation transport mechanisms in presence of magnetic fields can be described in terms of formalisms based on the Boltzmann radiation transport equation. However, when the transport occurs not in vacuum but in scattering media the difficulty of its analytical calculation considerably increases because of the random character of the radiation interaction processes with the scattering medium. In this context, the present work proposes and describe a methodology to characterize effects due to strong magnetic fields on electron trajectories by means of numerical techniques, such as Monte Carlo simulation, to be applied to scattering media, typically present in MRI-LINAC devices in radiotherapy situations. The obtained results for air as the scattering material show that variations of the electrons' trajectories depend jointly on the magnetic field and the interactions with the scattering medium. Increasing the magnetic field intensity enhances the curvature of the trajectories while decreasing the displacement in the radial direction.

Abstract For decades people near mine sites have complained of blasting-induced vibrations, because of that, tools of vibration prediction have a great significance. The idea of this work is to determine the probability of rock-blasting-induced vibrations exceed the threshold limits of international standards, related to structural damage and human discomfort. It is believed that a probabilistic analysis is a useful tool for avoiding harm to communities nearby. This paper shows how Monte Carlo simulation could be applied for blast-induced vibration from stochastic variables (maximum charge per delay and distance between blasting and the monitoring point) with the application of a representative attenuation law. For ground vibrations, it was observe 55% probability of structural damage in a community for a charge of 150 kg, 84% probability for human discomfort with a distance of 1000 m and 450 kg. However, for airblast overpressure it was observed no influence over the community studied.

Resumen Durante décadas, las personas cercanas a las minas se han quejado de vibraciones inducidas por voladuras, por eso, las herramientas de predicción de vibraciones tienen un gran significado. Este trabajo intenta a determinar la probabilidad de las vibraciones inducidas por voladuras de rocas superen los límites de los estándares internacionales de daños estructurales y malestar humano. Se cree que un análisis probabilístico es una herramienta útil para evitar daños a las comunidades cercanas. Este artículo muestra cómo se podría aplicar la simulación de Monte Carlo para las vibraciones a partir de variables estocásticas (carga explosiva y distancia) y una ley de atenuación representativa. Para vibraciones del suelo se observó 55% de probabilidad de daño estructural en una comunidad para una carga de 150 kg, 84% de probabilidad de malestar humano con una distancia de 1000 m y 450 kg. Considerando el ruido, no se observó influencia sobre la comunidad estudiada.

Abstract Many of the sectors in the economy were negatively affected, particularly insurance sector with the appearance of COVID-19. With the support of governments or reinsurers through the payment of a premium, insurance companies could receive a contingent resource in the face of excess infections caused by the pandemic. This paper calculates the premium to cover the excess affected population with a financial options model with a diffusion process without and with Poisson jumps and the Susceptible-Infected-Recovered (SIR) epidemiological model (this estimation is original). The obtained system is approximated with the Monte Carlo simulation method. The results show that there are important differences in the option premiums when Poisson jumps are included. Lastly, it is highlighted that the premium depends on the behavior trajectory of contagions and strike contagion value (K). This work has a limitation when applied to very particular cases, but a calibration of the parameters with more real information could be done in future research.

Resumen Con la aparición del COVID-19 muchos de los sectores de la economía resultaron afectados negativamente, particularmente el de seguros. Con el apoyo de los gobiernos o de las reaseguradoras mediante el pago de una prima las compañías aseguradoras podrían recibir un recurso contingente ante el exceso de contagios provocados por la pandemia. El presente trabajo calcula la prima para cubrir el exceso de población afectada con un modelo de opciones financieras con un proceso de difusión sin y con saltos de Poisson y el modelo epidemiológico Susceptibles-Infectados-Recuperados (SIR), esta propuesta no cuenta con muchos precedentes. El sistema obtenido se aproxima con el método de simulación Monte Carlo. Al final se presentan los resultados donde hay diferencias importantes en las primas de las opciones al incluir saltos de Poisson. Por último, dentro de las principales conclusiones se destaca que el valor de la prima depende de la trayectoria de comportamiento de los contagios y del valor de ejercicio del contagio (K). El trabajo tiene limitaciones al aplicarse a casos muy particulares, pero que podrían extenderse a casos más reales en investigaciones futuras.

Resumen Antecedentes: El diseño de medidas repetidas es uno de los más usados en ciencias sociales y de la salud. Aunque hay otras alternativas más avanzadas, el análisis de varianza de medidas repetidas (ANOVA-MR) sigue siendo el procedimiento más empleado para analizar las diferencias de medias. El impacto de la violación de la normalidad ha sido muy estudiado en el ANOVA intersujeto, pero los estudios son muy escasos en el ANOVA-MR. Por ello, el objetivo de este trabajo es realizar dos estudios de simulación Monte Carlo para analizar el error de Tipo I y la potencia cuando se incumple este supuesto bajo el cumplimiento de la esfericidad. Método: El estudio 1 incluye 20 distribuciones, tanto conocidas como desconocidas, manipulando el número de medidas repetidas (3, 4, 6 y 8) y el tamaño muestral (de 10 a 300). El estudio 2 incluye diferentes distribuciones en cada medida repetida. Las distribuciones analizadas representan desviación leve, moderada y severa de la normalidad. Resultados: En general, los resultados muestran que tanto el error Tipo I como la potencia del estadístico F no se alteran con la violación de la normalidad. Conclusiones: El ANOVA-MR es generalmente robusto a la no normalidad cuando la esfericidad se satisface.

Abstract Background: Repeated measures designs are commonly used in health and social sciences research. Although there are other, more advanced, statistical analyses, theF-statistic of repeated measures analysis of variance (RM-ANOVA) remains the most widely used procedure for analyzing differences in means. The impact of the violation of normality has been extensively studied for between-subjects ANOVA, but this is not the case for RM-ANOVA. Therefore, studies that extensively and systematically analyze the robustness of RM-ANOVA under the violation of normality are needed. This paper reports the results of two simulation studies aimed at analyzing the Type I error and power of RM-ANOVA when the normality assumption is violated but sphericity is fulfilled. Method: Study 1 considered 20 distributions, both known and unknown, and we manipulated the number of repeated measures (3, 4, 6, and 8) and sample size (from 10 to 300). Study 2 involved unequal distributions in each repeated measure. The distributions analyzed represent slight, moderate, and severe deviation from normality. Results: Overall, the results show that the Type I error and power of theF-statistic are not altered by the violation of normality. Conclusions: RM-ANOVA is generally robust to non-normality when the sphericity assumption is met.

Resumen Se realizó una evaluación de la rentabilidad económica del uso del método de regeneración de árboles padre para poblaciones arboladas en contraste con el uso de una plantación. La comparación proporciona información sobre la eficiencia económica de ambos métodos de regeneración que es útil para la toma de decisiones. La evaluación se basa en un algoritmo para estimar el número de árboles padre que se deben dejar en pie de acuerdo con la calidad de sitio, variable que integra factores ambientales e intrínsecos de cada especie. Se muestra un ejemplo de calibración del modelo y, a través de una simulación Monte Carlo, variando parámetros del modelo y costos, se evalúa el valor presente neto, tanto de la regeneración natural en un periodo de t años, como de la realización de una plantación. Los resultados muestran que, si no existen limitantes fisiológicas de la especie y de sitio para que se regenere a través de plantación, esta es más rentable que una repoblación natural. Se discuten algunas variables económicas y sociales que influyen en la decisión final de un productor de adoptar la regeneración asistida. Se concluye que la regeneración asistida es más rentable y potencialmente más eficiente que la regeneración natural.

Abstract A comparison of the profitability of the seed tree regeneration method versus planting is carried out. The comparison provides useful information on the economic efficiency of both methods of regeneration which is useful for decision making. The evaluation is based on an algorithm to estimate the number of seed trees to left standing based on-site quality, variable that accounts for environmental and species intrinsic factors. An example of model calibration is shown and through a Monte Carlo simulation varying several model and costs parameters, the net present value is evaluated for both, regeneration, and planting methods over a period of t years. Results show that if there are no physiological constraints for the species regeneration through planting, this method is more cost effective than natural regeneration. Some economic and social variables that influence the producer's final decision to adopt planting are discussed. It is concluded that assisted regeneration is more cost-effective and potentially more efficient than natural regeneration.

Resumen Debido a la incertidumbre existente en los procesos productivos, las organizaciones necesitan establecer en qué momento una variación en el costo con respecto a una estimación de referencia se considera aceptable y cuando se deben emprender acciones correctivas, preventivas o de mejora. Considerar valores tolerables genera la idea que el costo no es un parámetro determinístico, sino que puede oscilar en un rango definido, por lo que las estimaciones que pertenezcan a dicho rango se consideran normales y las que estén por fuera generan medidas de intervención por parte de los tomadores de decisiones. Este artículo presenta una metodología que permite establecer el intervalo de tolerancia. La metodología utiliza un modelo matricial de costeo basado en actividades que incorpora la simulación Monte Carlo para representar la variabilidad inherente en cada uno de los procesos. La metodología fue aplicada en un caso de estudio en el sector transporte obteniendo como resultado el rango de valores de control, el cual sirve como base para la toma de decisiones tácticas y operativas dentro de la organización.

Abstract Due to the uncertainty in production processes, organizations must establish when a variation in cost with respect to a reference estimate is considered acceptable and when corrective, preventive or improvement actions must be taken. Considering tolerable values generates the idea that the cost is not a deterministic parameter, but that it can oscillate in a defined range, so that the estimates that belong to said range are considered normal and those that are outside generate intervention measures by part of decision makers. This paper presents a methodology that allows establishing the tolerance interval. The methodology uses an activity-based costing matrix model that incorporates Monte Carlo simulation to represent the inherent variability in each of the processes. The methodology was applied in a case study in the transportation sector, obtaining as a result the range of control values, which serves as the basis for making tactical and operational decisions within the organization.

Resumen Las comparaciones por pares o métodos post-hoc se utilizan para determinar la fuente de la diferencia de medias de grupo en ANOVA unidireccional. Estos métodos dependen principalmente de la suposición de normalidad. Sin embargo, no normales distribuciones son más frecuentes que la distribución normal. Por lo tanto, los métodos robustos de estimación se convierten en herramientas muy importantes en el análisis estadístico. En este artículo, asumimos que la distribución de los términos de error es la de Azzalini sesgar t y obtener los estimadores robustos para realizar pruebas post-hoc en ANOVA de una vía. Utilizamos la metodología de máxima verosimilitud (ML) y comparamos esta metodología con algunos de los estimadores robustos como el estimador M, estimador de onda, media recortada y máxima verosimilitud modificada (MML) metodología con estudio de simulación Monte Cario. Los resultados de la simulación muestran que la metodología propuesta es más preferible. También comparamos potencia valores de las estadísticas de prueba y concluimos que las estadísticas de prueba basadas en los estimadores ML son más poderosos que las estadísticas de prueba basadas en otros métodos.

Abstract The pairwise comparisons or post-hoc methods are used for determining the source of the difference of group means in one-way ANOVA. These methods are mostly depend on normality assumption. However, nonnormal distributions are more prevalent than normal distribution. Therefore, robust estimation methods become very important tools in statistical analysis. In this paper, we assume that the distribution of the error terms is Azzalini's skew t and obtain the robust estimators in order to make post-hoc tests in one-way ANOVA. We use maximum likelihood (ML) methodology and compare this methodology with some of robust estimators like M estimator, Wave estimator, trimmed mean and modified maximum likelihood (MML) methodology with Monte Carlo simulation study. Simulation results show that the proposed methodology is more preferable. We also compare power values of the test statistics and conclude that the test statistics based on the ML estimators are more powerful than the test statistics based on other methods.

ABSTRACT This study aimed to determine the economic feasibility of the chlorophyll meters SPAD-502 and ClorofiLOG CFL1030 for nitrogen fertilization management in corn using deterministic and probabilistic models for small farms. To do so, representative cash flows for corn production were analyzed using three variable rate fertilization systems based on different nitrogen sufficiency indices. Thus, for each system, eight scenarios of different-sized cultivation areas and two sampling methodologies of the chlorophyll meter index (CMI) were simulated as follows: scenarios 1, 2, 3 and 4 used the CMI average obtained from 30 individual samplings in 1, 3, 5 and 7 ha areas; scenarios 5, 6, 7 and 8 used the CMI average in zones (30 individual samplings per zone), with 3, 5, 7 and 9 zones in 1, 3, 5 and 7 ha areas, respectively. The economic feasibility of the sensors was determined using different economic indicators, such as the internal rate of return (IRR) and the net present value (NPV), considering a minimum rate of attractiveness (MRA) of 6% per year with an 8-year project horizon. The risk was assessed by the Monte Carlo simulation. The results showed that regardless of the sampling approach, the SPAD-502 is economically feasible when acquired in Brazil or via importation through the USA in cultivation areas from 7 or 5 ha in size, respectively. Conversely, the ClorofiLOG CFL1030 could result in greater economic feasibility in areas from 3 ha.

RESUMEN El objetivo de este estudio fue precisar la viabilidad económica de los medidores de clorofila SPAD-502 y ClorofiLOG CFL1030para la gestión de la fertilización nitrogenada en maíz, utilizando modelos deterministas y probabilísticos para pequeñas explotaciones. Con este fin, se analizaron flujos de caja representativos de la producción de maíz mediante tres sistemas de fertilización de tasa variable basados en diferentes índices de suficiencia de nitrógeno. Para cada sistema se simularon ocho escenarios de áreas de cultivo de distinto tamaño y dos metodologías de muestreo del índice clorofílico (CMI). En los escenarios 1, 2, 3 y 4 se utilizó la media del CMI obtenida de 30 muestreos individuales en áreas de 1, 3, 5 y 7 ha. En los escenarios 5, 6, 7 y 8 se usó la media del CMI por zonas (30 muestreos individuales por zona), con 3, 5, 7 y 9 zonas en áreas de 1, 3, 5 y 7 ha, respectivamente. La viabilidad económica de los sensores se determinó a través de diferentes indicadores económicos, como la tasa interna de rendimiento (TIR) y el valor actual neto (VAN), considerando una tasa mínima de atracción (TMA) del 6% anual con un horizonte de proyecto de 8 años. El riesgo se evaluó mediante la simulación de Montecarlo. Los resultados evidenciaron que, independientemente del enfoque de muestreo, el SPAD-502 es económicamente viable cuando se adquiere en Brasil o mediante la importación a través de EE UU en áreas de cultivo a partir de 7 o 5 ha, respectivamente. Por el contrario, el ClorofiLOG CFL1030 podría resultar más factible económicamente en áreas a partir de 3 ha.