Current Organizations / Institutes

The organizations/institutes from the scholar's latest article


Roger Tovar-Falón has a total of 27 co-authors. Here are the top co-authors by number of shared publications:


Latest Publications

  • Open Access

    A Bivariate Copula–Driven Multi-State Model for Statistical Analysis in Medical Research

    by Hugo Brango, Roger Tovar-Falón, Guillermo Martínez-Flórez

    We develop and evaluate a copula-based multistate model for illness–death processes with dependent transition times. The framework couples Cox proportional hazards models for the marginal transition intensities with Archimedean copulas to capture dependence, and it is estimated via the Inference Functions for Margins (IFM) approach under right censoring. A Monte Carlo study shows that assuming independence between transitions can severely underestimate joint survival, yielding coverage as low as 40% under strong dependence, compared with 92% to 97% when copulas are used. We apply the method to a large Colombian cohort of COVID-19 patients (2021 to 2022) that includes sociodemographic, clinical, and vaccination data. The Gumbel copula best captures the strong positive dependence between hospitalization and death, producing more accurate joint survival estimates than independence-based models. Model diagnostics, including proportional hazards tests, Kaplan-Meier comparisons, hazard rate functions, and TTT plots, support the adequacy of the Cox margins. We also discuss limitations and avenues for extension, such as parametric or cure-fraction margins, nested or vine copulas, and full-likelihood estimation. Overall, the results underscore the methodological and applied value of integrating copulas into multistate models, offering a robust framework for analyzing dependent event times in epidemiology and biomedicine.


  • Open Access

    Mathematical Formalization and Applications to Data with Excess of Zeros and Ones of the Unit-Proportional Hazard Inflated Models

    by Guillermo Martínez-Flórez, Roger Tovar-Falón, Héctor W. Gómez

    In this study, we model the rate or proportion of a specific phenomenon using a set of known covariates. To fit the regression model, which explains the phenomenon within the intervals (0,1), [0,1), (0,1], or [0,1], we employ a logit link function. This approach ensures that the model’s predictions remain within the appropriate range of zero to one. In cases of inflation at zero, one, or both, the logit link function is similarly applied to model the dichotomous Bernoulli-type variable with a multinomial response. The findings demonstrate that the model yields a non-singular information matrix, ensuring valid statistical inference. This ensures the invertibility of the information matrix, allowing for hypothesis testing based on likelihood statistics regarding the parameters in the model. This is not possible with other asymmetric models, such as those derived from the skew-normal distribution, which has a singular information matrix at the boundary of the skewness parameter. Finally, empirical results show the model’s effectiveness in analyzing proportion data with inflation at zero and one, proving its robustness and practicality for analyzing bounded data in various fields of research.


  • Open Access

    Skew-Normal Inflated Models: Mathematical Characterization and Applications to Medical Data with Excess of Zeros and Ones

    by Guillermo Martínez-Flórez, Roger Tovar-Falón, Víctor Leiva, Cecilia Castro

    The modeling of data involving proportions, confined to a unit interval, is crucial in diverse research fields. Such data, expressing part-to-whole relationships, span from the proportion of individuals affected by diseases to the allocation of resources in economic sectors and the survival rates of species in ecology. However, modeling these data and interpreting information obtained from them present challenges, particularly when there is high zero–one inflation at the extremes of the unit interval, which indicates the complete absence or full occurrence of a characteristic or event. This inflation limits traditional statistical models, which often fail to capture the underlying distribution, leading to biased or imprecise statistical inferences. To address these challenges, we propose and derive the skew-normal zero–one inflated (SNZOI) models, a novel class of asymmetric regression models specifically designed to accommodate zero–one inflation presented in the data. By integrating a continuous-discrete mixture distribution with covariates in both continuous and discrete parts, SNZOI models exhibit superior capability compared to traditional models when describing these complex data structures. The applicability and effectiveness of the proposed models are demonstrated through case studies, including the analysis of medical data. Precise modeling of inflated proportion data unveils insights representing advancements in the statistical analysis of such studies. The present investigation highlights the limitations of existing models and shows the potential of SNZOI models to provide more accurate and precise inferences in the presence of zero–one inflation.


  • Open Access

    Dichotomous Proportional Hazard Regression Model: A Case Study on Students’ Dropout

    by Guillermo Martínez-Flórez, Roger Tovar-Falón, Carlos Barrera-Causil

    In problems involving binary classification, researchers often encounter data suitable for modeling dichotomous responses. These scenarios include medical diagnostics, where outcomes are classified as “disease” or “no disease”, and credit scoring in finance, determining whether a loan applicant is “high risk” or “low risk”. Dichotomous response models are also useful in many other areas for estimating binary responses. The logistic regression model is one option for modeling dichotomous responses; however, other statistical models may be required to improve the quality of fits. In this paper, a new regression model is proposed for cases where the response variable is dichotomous. This novel, non-linear model is derived from the cumulative distribution function of the proportional hazard distribution, and is suitable for modeling binary responses. Statistical inference is performed using a classical approach with the maximum likelihood method for the proposed model. Additionally, it is demonstrated that the introduced model has a non-singular information matrix. The results of a simulation study, along with an application to student dropout data, show the great potential of the proposed model in practical and everyday situations.


  • Open Access

    Bivariate Unit-Weibull Distribution: Properties and Inference

    by Roger Tovar-Falón, Guillermo Martínez-Flórez, Luis Páez-Martínez

    In this article, we introduce a novel bivariate probability distribution that is absolutely continuous. Considering the Farlie–Gumbel–Morgenstern (FGM) copula and the unit-Weibull distribution, we can obtain a bivariate unit-Weibull distribution. We evaluate the main properties of the new proposal and use two estimation methods to estimate the parameter for the bivariate probability distribution. A brief Monte Carlo simulation study is conducted to assess the behavior of the employed estimation method and the characteristics of the estimators. Ultimately, as an illustration, a real-life application is presented, demonstrating the utility of the proposal.


  • Open Access

    The Log-Bimodal Asymmetric Generalized Gaussian Model with Application to Positive Data

    by Guillermo Martínez-Flórez, Roger Tovar-Falón, Heleno Bolfarine

    One of the most widely known probability distributions used to explain the probabilistic behavior of positive data is the log-normal (LN). Although the LN distribution is capable of adjusting data types, it is not always fully true that the model manages to adequately model the behavior of the response of interest since in some cases, the degree of skewness and/or kurtosis of the data are greater or less than those that the LN distribution can capture. Another peculiarity of the LN distribution is that it only fits unimodal positive data, which constitutes a limitation when dealing with data that present more than one mode (bimodality). On the other hand, the log-normal model only fits unimodal positive data and in reality there are multiple applications where the behavior of materials is bimodal. To fill this gap, this paper introduces a new probability distribution that is capable of fitting unimodal or bimodal positive data with a high or low degree of skewness and/or kurtosis. The new distribution is a generalization of the LN distribution. For the new proposal, its main properties are studied and the process of estimation of the parameters involved in the model is carried out from a classical perspective using the maximum likelihood method. An important feature of this distribution is the non-singularity of the Fisher information matrix, which guarantees the use of asymptotic theory to study the properties of the parameter estimators. A Monte Carlo type simulation study is carried out to evaluate the properties of the estimators and finally, an illustration is presented with a set of data related to the concentration of nickel in soil samples, allowing to show that the proposed distribution fits extremely well in certain situations.

  • Open Access

    Some Extensions of the Asymmetric Exponentiated Bimodal Normal Model for Modeling Data with Positive Support

    by Roger Tovar-Falón, Guillermo Martínez-Flórez, Isaías Ceña-Tapia

    It is common in many fields of knowledge to assume that the data under study have a normal distribution, which often generates mistakes in the results, since this assumption does not always coincide with the characteristics of the observations under analysis. In some cases, the data may have degrees of skewness and/or kurtosis greater than what the normal model can capture, and in others, they may present two or more modes. In this work, two new families of skewed distributions are presented that fit bimodal data with positive support. The new families were obtained from the extension of the bimodal normal distribution to the alpha-power family class. The proposed distributions were studied for their main properties, such as their probability density function, cumulative distribution function, survival function, and hazard function. The parameter estimation process was performed from a classical perspective using the maximum likelihood method. The non-singularity of Fisher’s information was demonstrated, which made it possible to find the stochastic convergence of the vector of the maximum likelihood estimators and, based on the latter, perform statistical inference via the likelihood ratio. The applicability of the proposed distributions was exemplified using real data sets.


  • Open Access

    The Multivariate Skewed Log-Birnbaum–Saunders Distribution and Its Associated Regression Model

    by Guillermo Martínez-Flórez, Sandra Vergara-Cardozo, Roger Tovar-Falón, Luisa Rodriguez-Quevedo

    In this article, a multivariate extension of the unit-sinh-normal (USHN) distribution is presented. The new distribution, which is obtained from the conditionally specified distributions methodology, is absolutely continuous, and its marginal distributions are univariate USHN. The properties of the multivariate USHN distribution are studied in detail, and statistical inference is carried out from a classical approach using the maximum likelihood method. The new multivariate USHN distribution is suitable for modeling bounded data, especially in the (0,1)p region. In addition, the proposed distribution is extended to the case of the regression model and, for the latter, the Fisher information matrix is derived. The numerical results of a small simulation study and two applications with real data sets allow us to conclude that the proposed distribution, as well as its extension to regression models, are potentially useful to analyze the data of proportions, rates, or indices when modeling them jointly considering different degrees of correlation that may exist in the study variables is of interest.


  • Open Access

    The Bivariate Unit-Sinh-Normal Distribution and Its Related Regression Model

    by Guillermo Martínez-Flórez, Artur J. Lemonte, Germán Moreno-Arenas, Roger Tovar-Falón

    In this paper, a new bivariate absolutely continuous probability distribution is introduced. The new distribution, which is called the bivariate unit-sinh-normal (BVUSHN) distribution, arises by applying a transformation to the bivariate Birnbaum–Saunders distribution (BVBS). The main properties of the new proposal are studied in detail. In addition, from the new distribution, the BVUSHN regression model is also introduced. For both the bivariate probability distribution and the respective associated regression model, parameter estimation is conducted from a classical approach by using the maximum likelihood method together with the two-step estimation method. A small Monte Carlo simulation study is carried out to evaluate the behavior of the used estimation method and the properties of the estimators. Finally, for illustrative purposes, two applications with real data are presented in which the usefulness of the proposals is evidenced.