Roger Tovar-Falón
Current Organizations / Institutes
The organizations/institutes from the scholar's latest article
Research Subjects
Roger Tovar-Falón has a total of 27 co-authors. Here are the top co-authors by number of shared publications:
- Guillermo Martínez-Flórez
Departamento de Matemática y Estadística, Facultad de Ciencias Básicas, Universidad de Córdoba, Montería 230001, Colombia
25 shared publications - Hector W. Gomez
Departamento de Ciencias Matemáticas y Físicas, Facultad de Ingeniería, Universidad Católica de Temuco, Temuco 4813302, Chile
4 shared publications - Heleno Bolfarine
Universidade de São Paulo
3 shared publications - Carlos Barrera-Causil
Instituto Tecnológico Metropolitano
2 shared publications - Sandra Vergara-Cardozo
Departamento de Estadística, Facultad de Ciencias, Universidad Nacional de Colombia, Bogotá 111321, Colombia
2 shared publications - Artur J. Lemonte
Departamento de Estatística, Universidade Federal do Rio Grande do Norte, Natal, Brazil
1 shared publication
Latest Publications
- Open Access
- 8 August 2026
- Published by Elsevier in IJID Regions
- Vol. 20
Clinical, parasitological, and hematologic profile of pregnant women with malaria in Colombia: a retrospective cross-sectional study
- 8 August 2026
- Published by Elsevier in IJID Regions
- Vol. 20
Open AccessClinical, parasitological, and hematologic profile of pregnant women with malaria in Colombia: a retrospective cross-sectional study
- Open Access
- 24 September 2025
- Published by MDPI AG in Mathematics
- Vol. 13 (19) , 3072
A Bivariate Copula–Driven Multi-State Model for Statistical Analysis in Medical Research
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 under strong dependence, compared with to 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.- 24 September 2025
- Published by MDPI AG in Mathematics
- Vol. 13 (19) , 3072
Open AccessA Bivariate Copula–Driven Multi-State Model for Statistical Analysis in Medical Research
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 under strong dependence, compared with to 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
- 15 November 2024
- Published by MDPI AG in Mathematics
- Vol. 12 (22) , 3566
Mathematical Formalization and Applications to Data with Excess of Zeros and Ones of the Unit-Proportional Hazard Inflated Models
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 , , , or , 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.- 15 November 2024
- Published by MDPI AG in Mathematics
- Vol. 12 (22) , 3566
Open AccessMathematical Formalization and Applications to Data with Excess of Zeros and Ones of the Unit-Proportional Hazard Inflated Models
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 , , , or , 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
- 12 August 2024
- Published by MDPI AG in Mathematics
- Vol. 12 (16) , 2486
Skew-Normal Inflated Models: Mathematical Characterization and Applications to Medical Data with Excess of Zeros and Ones
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.
- 12 August 2024
- Published by MDPI AG in Mathematics
- Vol. 12 (16) , 2486
Open AccessSkew-Normal Inflated Models: Mathematical Characterization and Applications to Medical Data with Excess of Zeros and Ones
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
- 11 July 2024
- Published by MDPI AG in Mathematics
- Vol. 12 (14) , 2170
Dichotomous Proportional Hazard Regression Model: A Case Study on Students’ Dropout
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.
- 11 July 2024
- Published by MDPI AG in Mathematics
- Vol. 12 (14) , 2170
Open AccessDichotomous Proportional Hazard Regression Model: A Case Study on Students’ Dropout
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
- 1 September 2023
- Published by MDPI AG in Mathematics
- Vol. 11 (17) , 3760
Bivariate Unit-Weibull Distribution: Properties and Inference
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.
- 1 September 2023
- Published by MDPI AG in Mathematics
- Vol. 11 (17) , 3760
Open AccessBivariate Unit-Weibull Distribution: Properties and Inference
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
- 19 August 2023
- Published by MDPI AG in Mathematics
- Vol. 11 (16) , 3587
The Log-Bimodal Asymmetric Generalized Gaussian Model with Application to Positive Data
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.
- 19 August 2023
- Published by MDPI AG in Mathematics
- Vol. 11 (16) , 3587
Open AccessThe Log-Bimodal Asymmetric Generalized Gaussian Model with Application to Positive Data
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
- 23 March 2023
- Published by MDPI AG in Mathematics
- Vol. 11 (7) , 1563
Some Extensions of the Asymmetric Exponentiated Bimodal Normal Model for Modeling Data with Positive Support
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.
- 23 March 2023
- Published by MDPI AG in Mathematics
- Vol. 11 (7) , 1563
Open AccessSome Extensions of the Asymmetric Exponentiated Bimodal Normal Model for Modeling Data with Positive Support
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
- 22 February 2023
- Published by MDPI AG in Mathematics
- Vol. 11 (5) , 1095
The Multivariate Skewed Log-Birnbaum–Saunders Distribution and Its Associated Regression Model
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 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.- 22 February 2023
- Published by MDPI AG in Mathematics
- Vol. 11 (5) , 1095
Open AccessThe Multivariate Skewed Log-Birnbaum–Saunders Distribution and Its Associated Regression Model
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 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
- 31 August 2022
- Published by MDPI AG in Mathematics
- Vol. 10 (17) , 3125
The Bivariate Unit-Sinh-Normal Distribution and Its Related Regression Model
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.
- 31 August 2022
- Published by MDPI AG in Mathematics
- Vol. 10 (17) , 3125
Open AccessThe Bivariate Unit-Sinh-Normal Distribution and Its Related Regression Model
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.