Next Article in Journal
Investigating Oscillatory Behavior in Third-Order Neutral Differential Equations with Canonical Operators
Next Article in Special Issue
OUCH: Oversampling and Undersampling Cannot Help Improve Accuracy in Our Bayesian Classifiers That Predict Preeclampsia
Previous Article in Journal
Child Health Dataset Publishing and Mining Based on Differential Privacy Preservation
Previous Article in Special Issue
Modeling Residential Energy Consumption Patterns with Machine Learning Methods Based on a Case Study in Brazil
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

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

by
Guillermo Martínez-Flórez
1,
Roger Tovar-Falón
1,
Víctor Leiva
2,* and
Cecilia Castro
3,*
1
Department of Mathematics and Statistics, Universidad de Córdoba, Montería 230002, Colombia
2
School of Industrial Engineering, Pontificia Universidad Católica de Valparaíso, Valparaíso 2362807, Chile
3
Centre of Mathematics, Universidade do Minho, 4710-057 Braga, Portugal
*
Authors to whom correspondence should be addressed.
Mathematics 2024, 12(16), 2486; https://doi.org/10.3390/math12162486
Submission received: 20 June 2024 / Revised: 20 July 2024 / Accepted: 8 August 2024 / Published: 12 August 2024
(This article belongs to the Special Issue Applied Statistics in Real-World Problems)

Abstract

:
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.

1. Introduction

Analysis of phenomena represented by proportion data, restricted to values between zero and one, is crucial in various scientific fields. These data describe part-to-whole relationships and appear in diverse applications, such as disease prevalence, resource distribution in economics, species survival rates, and habitat usage in ecology [1]. Accurately modeling these data poses challenges, especially when there is high zero–one inflation, which can reflect the absence or complete manifestation of a trait or event. Such inflation creates obstacles for traditional statistical models, like the tobit model [2,3,4], which, despite its prevalent application in economics, medicine, and social sciences, often struggles to accurately represent the underlying distribution of proportion data with inflated extremes.
To enhance the robustness of the tobit models, recent efforts have incorporated the Student-t [5,6], power-normal [7,8,9], and Birnbaum–Saunders [10,11] distributions, addressing asymmetry and kurtosis in the data distribution. Similarly, the adoption of the skew-normal (SN) distribution [12,13,14] has provided a framework for fitting data with inherent asymmetry, demonstrating the versatility required for contemporary datasets [10,15,16,17].
Recent advancements include the use of SN distributions in various applications, such as quantile estimation in child weight data and classification methods for serological status [18,19]. Furthermore, issues of inference and analyses utilizing SN distributions have been extensively studied, addressing challenges in modeling near-normality data [20,21,22,23].
The beta distribution, used for unit interval data when asymmetry is apparent, necessitates beta regression when independent variables (covariates) are considered [24]. However, standard approaches often fall short when zero–one inflation is present, limiting the existing models [25,26,27]. Recent studies have proposed extensions to the beta distribution to handle these limitations, such as the zero–one inflated beta models, which have shown promise in dealing with highly skew distributed data [28,29,30,31,32,33].
Despite strides in modeling data with inflation and asymmetry, a gap remains in adequately addressing zero–one inflation within data of proportions. Existing models do not fully capture the unique distributional characteristics and complexities introduced by these inflations, which can lead to biased estimators and imprecise inferences [34,35].
The main objective of the present study is to propose and derive skew-normal zero–one inflated (SNZOI) models, a novel class of asymmetric regressions specifically designed to address the challenges of zero–one inflation in proportion data confined to the unit interval. The SNZOI models integrate a continuous-discrete mixture distribution with covariates, enabling them to effectively capture the complex dynamics of such data. This integration accounts for asymmetry and mixed discrete-continuous characteristics of the data distribution, providing more precise and unbiased statistical inferences.
Traditional formulations, such as the tobit and beta regressions, often struggle with zero–one inflation, leading to imprecise model fits and biased inferences. Through detailed case studies, including the analysis of disease prevalence, we demonstrate the limitations of these traditional models. In contrast, the SNZOI models exhibit superior performance, offering deeper insights and more accurate fits to the data. To achieve precise parameter estimators, we employ the maximum likelihood (ML) method supported by iterative procedures like the Newton–Raphson algorithm. The accuracy of these ML estimators is rigorously assessed using the corresponding Fisher information matrix. The present study addresses a critical gap in the field of statistical modeling by providing a robust alternative for researchers dealing with inflated data in the unit interval. Our proposal details the development, application, and comparative analysis of SNZOI models, showing their potential to enhance statistical analysis in areas where zero–one inflation is common.
The remainder of this article is structured as follows, as illustrated in Figure 1. Section 2 provides background on skew models, focusing on the SN distribution and extensions to handle inflated data. In Section 3, we present regressions for data in the unit interval with inflation, detailing the complementary log–log SN, probit SN, and logit doubly censored models, along with their information matrices. In Section 4, the new SNZOI regression model is introduced, including its formulation, parameter estimation, and its particular case based on the probit structure. In Section 5, we demonstrate the application of these models through empirical case studies on doubly censored data and data inflated at one. Section 6 states a discussion of the key findings, comparisons with previous studies, limitations, and future research directions. The article concludes with Section 7.

2. Modeling with Skew Distributions

In this section, we provide background on general concepts on statistical modeling, link functions, and skew models, with emphasis on the SN distribution. Extensions to the single and doubly censored cases are also considered. Also, in this section we indicate how the censoring thresholds are determined.

2.1. Fundamental Concepts of Regression and Link Functions

Various statistical models have been proposed to model data of rates and proportions within the closed unit interval denoted by [ 0 , 1 ] . Let x i = ( 1 , x 1 i , , x p i ) be the vector of p covariates for observation i, where the first element is a one to account for the intercept term, and  β = ( β 0 , β 1 , , β p ) be the vector of regression coefficients. The response (dependent) variable Y i is modeled as
Y i = x i β + ε i = μ i + ε i , i { 1 , , n } ,
where ε i is the error term.
The estimated regression coefficients are denoted by β ^ = ( β ^ 0 , β ^ 1 , , β ^ p ) . A common issue with some of these models is that the predicted value of the mean response, namely, E ^ ( Y i ) = μ ^ i = x i β ^ , may not always fall within the interval [ 0 , 1 ] . To address this issue, continuous, strictly monotonic, and twice differentiable link functions are used to transform the mean response μ i , ensuring that the predicted values fall within the unit interval.
Several link functions are commonly employed in practice, and they can be categorized as either symmetrical or asymmetrical [36]. Symmetrical link functions include the following:
  • Logit link function g ( μ ) = log μ / ( 1 μ ) .
  • Probit link function g ( μ ) = Φ 1 ( μ ) , where Φ 1 is the inverse of the standard normal cumulative distribution function (CDF).
These two links functions are characterized by their symmetry, that is, the response curve, which plots g ( μ ) against μ , is symmetric around μ , meaning they approach to the extremes at the same rate. Both the logit and probit link functions often yield similar results, except in extreme cases where predicted probabilities are close to 0 or 1. In such cases, the logit link function tends to compress extreme values more than the probit link function, leading to slightly different predictions at the boundaries of the probability scale.
In addition to symmetrical link functions, asymmetrical link functions are also utilized when the data are not symmetrically distributed within [ 0 , 1 ] . Specifically, we have the following:
  • Log–log link function g ( μ ) = log ( log ( 1 μ ) ) .
  • Complementary log–log link function g ( μ ) = log ( log ( μ ) ) .
The log–log link function is suitable for modeling data where early occurrences are rare but become more frequent over time, effectively capturing the increasing hazard rate. Conversely, the complementary log–log link function is ideal for situations where the data belonging to [0, 1] increase slowly from small to moderate values and then rises sharply near one. This function approaches to zero infinitely slower than others, making it particularly helpful for modeling skew distributed data with a heavy upper tail.
In summary, the rationale for selecting specific link functions is based on the characteristics of the data and the underlying assumptions of the model. For symmetrically distributed data, logit and probit link functions are often appropriate due to their symmetric nature around the mean μ . However, in cases where the response variable belonging to [0, 1] has a distribution that demonstrates high skewness, particularly with a heavy tail, the complementary log–log link function provides a better fit due to its ability to handle such skewness effectively. For example, in survival analysis, where the event of interest might occur at the end of the observation period (resulting in a high proportion of censored data), the complementary log–log link function is commonly used.
Another scenario is in dose-response models, where the probability of a response might sharply increase after a certain dose level, making the complementary log–log link function more suitable. Additionally, the log–log link function may be advantageous when modeling time-to-event data, where early events are rare but become more frequent as time progresses. This function captures the increasing hazard over time, which is often seen in failure time data.

2.2. Skew-Normal Distribution and Its Modeling

The probability density function (PDF) of the SN distribution is given by
ϕ SN ( y ; θ ) = 2 ω ϕ y ξ ω Φ λ y ξ ω , y R ,
where ξ R is a location parameter, ω R + is a scale parameter, λ R is an asymmetry parameter, ϕ is the standard normal PDF, and  θ = ( ξ , ω , λ ) is the parameter vector.
The notation Y SN ( ξ , ω , λ ) indicates that Y follows an SN distribution with parameters ξ , ω , and  λ . For  λ = 0 , the PDF stated in (2) reduces to the N ( ξ , ω 2 ) PDF.
The SN distribution can be extended to regression by setting ξ i = x i β , where x i is the observed covariate vector and β is the regression coefficient vector. Then, the response Y i is modeled using the formulation stated in (1), where ε i SN ( 0 , ω , λ ) , for  i { 1 , , n } .
Under the SN model, the expectation of the error term ε i is presented as
E ( ε i ) = ω 2 π λ 1 + λ 2 0 , i { 1 , , n } ,
which implies that E ( Y i ) x i β . To correct this, an adjustment to the intercept is defined as
β 0 * = β 0 + ω 2 π λ 1 + λ 2 ,
leading to E ( Y i ) = x i β * . Thus, we redefine the coefficients as β * = ( β 0 * , β 1 , , β p ) .
The log-SN (LSN) model extends the SN model to positive data and has PDF defined as
ϕ log SN ( y ; ξ , ω , λ ) = 2 ω y ϕ log ( y ) ξ ω Φ λ log ( y ) ξ ω , y R + ,
where ξ and ω are stated in (2). For  λ = 0 , the LSN model is the log-normal model.

2.3. Single and Doubly Censored Data

In practical situations, the response may be censored, meaning that, for some observations, we only know that their values fall above or below a certain threshold. To handle such situations, researchers commonly use the tobit model, in which the true value of the response variable is only partially observed due to censoring.
To relax the constraint on tail probabilities inherent in the tobit specification, a two-part model was proposed in [37]. The PDF for the random variable Y i in the two-part model, as proposed in [38], is given by
f ( y i ) = p i , if y i 0 ; ( 1 p i ) f ( y i ) , if y i > 0 ;
where f is the PDF of the model continuous part with positive support and p i = P ( Y i 0 ) = P ( Y i = 0 ) is the probability of observation i being exactly zero due to censoring.
The model presented in (3) was generalized in [39], allowing censored (limited) responses to result from interval censoring and described by
f ( y i ) = p i + ( 1 p i ) F ( c ) , if y i c ; ( 1 p i ) f ( y i ) , if y i > c ;
where F is the CDF corresponding to the PDF f.
The standard tobit model is a special case of the Moulton–Halsey model stated in [39] when p i = 0 , indicating no inflation at zero. The Moulton–Halsey model was later extended to a two-part form in [37], where log ( Y i ) SN ( x i β , ω , λ ) , with a probit link function defined as p i = Φ ( x i β ) . This extension is particularly advantageous for modeling skew distributed data with a heavy upper tail, as the probit function provides a smooth CDF-based transformation that effectively handles skewness. Furthermore, the Moulton–Halsey model, in its two-part form, accommodates limited responses resulting from interval censoring. This accommodation makes it suitable for scenarios where the response variable can take on a broader range of values, thereby improving the fit for data with complex distributions. Other probit and Heckman models can be found in [40,41,42,43].
In contrast to interval censoring, which limits responses to within specific ranges, in some cases data may be doubly censored, meaning that observations are only known to fall within a specific interval, while those outside this interval are recorded at the censoring points. Assume that Y * follows an SN distribution with parameters ξ , ω , and  λ with PDF defined in (2). Given a sample of size n, { Y 1 * , , Y n * } say, with observations { y 1 * , , y n * } , only values within the interval [ c 0 , c 1 ] are recorded. For  Y i * c 0 and Y i * c 1 , we record c 0 and c 1 , respectively. Therefore, for  i { 1 , , n } , the obtained values are defined as
y i = c 0 , if y i * c 0 ; y i * , if c 0 < y i * < c 1 ; c 1 , if y i * c 1 .
The sample { Y 1 * , , Y n * } is considered to be collected from a doubly censored skew-normal (DCSN) distribution, denoted by Y i * DCSN ( ξ , ω , λ ) , with PDF defined in (2) for c 0 < Y i * < c 1 . For  Y i * c 0 and Y i * c 1 , we have
P ( Y i = c 0 ) = P ( Y i * c 0 ) = Φ SN ( z 0 ) , P ( Y i = c 1 ) = P ( Y i * c 1 ) = 1 Φ SN ( z 1 ) ,
respectively, where z 0 = ( c 0 ξ ) / ω , z 1 = ( c 1 ξ ) / ω , and  Φ SN is the SN CDF. When λ = 0 , the usual doubly censored normal tobit model follows. Considering U = log ( Y ) , with Y DCSN ( ξ , ω , λ ) , the doubly censored log-skew-normal (DCLSN) model is obtained, denoted by DCLSN ( ξ , ω , λ ) . Similarly, when λ = 0 , we obtain the doubly censored log-normal tobit model.

2.4. Determination of Censoring Thresholds

The censoring thresholds c 0 and c 1 are determined based on prior knowledge about the data distribution or empirical observations. Typically, c 0 and c 1 are chosen to reflect cut-off points in the data. For example, c 0 might represent the minimum detectable limit of a measurement instrument, while c 1 could represent a saturation point beyond which measurements are no longer reliable.
The impact of censoring thresholds on parameter estimation may be important. If  c 0 and c 1 are not appropriately chosen, they can lead to biased parameter estimators and misleading inference. To mitigate this, sensitivity analyses are often conducted to assess the robustness of the parameter estimators to different choices of c 0 and c 1 . Additionally, graphical methods, such as inspecting the distribution of the censored and uncensored observations, can provide insights into the appropriateness of the chosen thresholds.
For instance, one may plot the empirical CDF of the observed data and identify natural breakpoints or plateaus that suggest appropriate values for c 0 and c 1 . Furthermore, domain-specific knowledge can guide the selection of these thresholds, ensuring they align with practical considerations and the specific context of the study. By carefully determining and justifying the censoring thresholds, researchers may improve the reliability and validity of their parameter estimators and subsequent inferences.

3. Regression Models for Unit Interval Data with Inflation

In this section, we present three regression models for unit interval (proportion) data that account for inflation at specific values, denoted by c 0 and c 1 . The values c 0 and c 1 represent the points of inflation where the response variable can take the values zero, one, or any value between zero and one. These models are helpful in scenarios where there is a high probability of the response variable being exactly at these specific values.

3.1. Complementary Log–Log Unit Skew-Normal Regression Model

We construct the unit regression model inflated at the values c 0 and c 1 , where the link for the location parameter μ i is the complementary log–log function log ( log ( μ i ) ) . We begin this construction by introducing the complementary log–log unit SN regression model truncated (SNT) at c 0 and c 1 , whose PDF is stated as
ϕ SNT ( y i ; θ ) = 1 ω ϕ SN y i μ i ω Φ SN c 1 μ i ω Φ SN c 0 μ i ω , c 0 < y i < c 1 , i { 1 , , n } ,
where ϕ SN and Φ SN are defined in (2). The standardization terms, which help to normalize the data within the given bounds, are established as
z 0 i = c 0 μ i ω , z i = y i μ i ω , z 1 i = c 1 μ i ω , i { 1 , , n } ,
where μ i is the location parameter, ω is the scale parameter, and  x i = ( 1 , x i 1 , , x i p ) is a 1 × ( p + 1 ) observed covariate vector associated with the parameter vector β = ( β 0 , β 1 , , β p ) . The full parameter vector is θ = ( β , ω , λ ) .
The log-likelihood function for a sample of size n, Y 1 , , Y n say, with observed values y 1 , , y n collected from the PDF ϕ SNT presented in (4), is obtained by taking the natural logarithm of the joint likelihood function defined as L ( θ ; y ) = i = 1 n ϕ SNT ( y i ; θ ) , where now θ = ( μ , ω , λ ) . Taking the natural logarithm in the above expression, we obtain the log-likelihood function established as
( θ ) = i = 1 n log 1 ω ϕ SN y i μ i ω Φ SN c 1 μ i ω Φ SN c 0 μ i ω .
To simplify the logarithmic terms presented in (6), we break down the expression inside the logarithm as
log 1 ω ϕ SN y i μ i ω Φ SN c 1 μ i ω Φ SN c 0 μ i ω = log 1 ω + log ϕ SN y i μ i ω log Φ SN c 1 μ i ω Φ SN c 0 μ i ω .
Note that the first term formulated in (7) simplifies to i = 1 n log ( 1 / ω ) = n log ( ω ) , the second term remains as i = 1 n log ( ϕ SN ( ( y i μ i ) / ω ) ) , and the third term simplifies to i = 1 n log ( Φ SN ( ( c 1 μ i ) / ω ) Φ SN ( ( c 0 μ i ) / ω ) ) . Thus, the log-likelihood function (with no constant term and using the symbol ∝ indicating proportional to) is now given by
( θ ) n log ( ω ) + i = 1 n log ϕ SN y i μ i ω i = 1 n log Φ SN c 1 μ i ω Φ SN c 0 μ i ω .
To further simplify and generate the score equations, we introduce the following notations: δ i = ( 1 μ i ) log ( 1 μ i ) ,   w i λ = ϕ SN ( λ z i ) / Φ SN ( λ z i ) ,
Δ j k i = z 1 i j ϕ SN k ( z 1 i ; λ ) z 0 i j ϕ SN k ( z 0 i ; λ ) Φ SN ( z 1 i ; λ ) Φ SN ( z 0 i ; λ ) , Γ j k i = z 1 i j ϕ SN k ( 1 + λ 2 z 1 i ) z 0 i j ϕ SN k ( 1 + λ 2 z 0 i ) Φ SN ( z 1 i ; λ ) Φ SN ( z 0 i ; λ ) ,
where δ i is a term involving the mean response μ i related to the complementary log–log link function and w i λ is the ratio of the SN PDF and CDF, whereas Δ j k i and Γ j k i are terms involving the standardized variables z i , z 1 i , and  z 0 i .
By taking the first derivatives of the function presented in (8) with respect to the parameters, we obtain the score elements. For the coefficient β j , the element is stated as
˙ ( β j ) = 1 ω i = 1 n x i j z i δ i + λ ω i = 1 n x i j w i λ δ i 1 ω i = 1 n x i j Δ 01 i δ i , j { 0 , 1 , , p } .
For the scale parameter ω , the score function is formulated as
˙ ( ω ) = n ω + 1 ω i = 1 n z i 2 λ ω i = 1 n z i w i λ + 1 ω i = 1 n Δ 11 i .
For the shape parameter λ , the score is defined as
˙ ( λ ) = i = 1 n w i λ z i 2 π 1 1 + λ 2 i = 1 n Γ 01 i .
The solution of the system of equations obtained by setting (9), (10), and (11) to zero yields the ML estimates of the parameters. This system is typically resolved using iterative numerical algorithms such as Newton–Raphson or quasi-Newton, which iteratively update parameter estimates to maximize the likelihood function.
The Newton–Raphson algorithm is known for its quadratic convergence properties, provided that the initial estimates are close to the true parameter values and the likelihood function is well-behaved. Specifically, the Newton–Raphson update rule is given by θ ( k + 1 ) = θ ( k ) ¨ ( θ ( k ) ) 1 ˙ ( θ ( k ) ) , where ¨ ( θ ) = 2 ( θ ) / θ θ is the Hessian matrix, from which the observed Fisher information matrix can be obtained, and  ˙ ( θ ) = ( θ ) / θ is its score function. The observed information matrix provides a local approximation to the curvature of the log-likelihood surface, which is used to iteratively update the parameter estimates. In cases where the likelihood surface is complex or the initial estimates are not close to the true values, quasi-Newton algorithms, which approximate the Hessian matrix, can provide more robust convergence. These algorithms do not require the explicit computation of the Hessian matrix at each iteration. Instead, they utilize rank-one updates to approximate the inverse Hessian matrix.
The popular quasi-Newton algorithm includes the Broyden–Fletcher–Goldfarb–Shanno (BFGS) and limited-memory BFGS (L-BFGS) approaches. The BFGS algorithm iteratively updates the parameter estimates as θ ( k + 1 ) = θ ( k ) + ϱ k p ( k ) , where ϱ k is a step length determined by a line search and p ( k ) = H k ˙ ( θ ( k ) ) is the search direction. Here, H k is the approximation to the inverse Hessian matrix at iteration k updated as H k + 1 = I ρ k s k y k H k I ρ k y k s k + ρ k s k s k , where I is the identity matrix, s k = θ ( k + 1 ) θ ( k ) , y k = ˙ ( θ ( k + 1 ) ) ˙ ( θ ( k ) ) , and  ρ k = 1 / y k s k . The L-BFGS algorithm, a variant of the BFGS algorithm, is particularly helpful for large-scale optimization problems as it maintains only a limited number of vectors to approximate the inverse Hessian matrix, so reducing computer memory requirements.
The update rule for L-BFGS is similar to BFGS but uses a limited memory representation, as stated by H k + 1 i = k m + 1 k I ρ i s i y i H i I ρ i y i s i + ρ i s i s i , where m is the number of past updates stored. This rule ensures efficiency in the optimization, making it suitable for maximizing the likelihood function in complex models [44].

3.2. Probit Unit Skew-Normal Regression Model

We now examine the unit SN model with the probit link defined as μ i = Φ ( x i β ) . This is similar to the case of the complementary log–log link in the unit SN model. The PDF ϕ SNT used to construct the likelihood function of the model is represented in (4), where z i , z 1 i , and  z 0 i are stated in (5), with  μ i being the location parameter, ω the scale parameter, and  x i = ( 1 , x i 1 , , x i p ) , a 1 × ( p + 1 ) observed covariate vector associated with β = ( β 0 , β 1 , , β p ) . For parameter estimation, the score function for β j is given by
˙ ( β j ) = 1 ω i = 1 n x i j z i ϕ ( x i β ) λ ω i = 1 n x i j w i λ ϕ ( x i β ) + 1 ω i = 1 n x i j Δ 01 i ϕ ( x i β ) , j { 0 , 1 , , p } .
The elements of the score equation for the parameters ω and λ are defined in (10) and (11), respectively, with the standardization terms z i , z 1 i , and  z 0 i stated in (5), as mentioned.

3.3. Logit Unit Skew-Normal Regression Model

Now, we explore the unit SN model with the logit link function expressed as
μ i = exp ( x i β ) 1 + exp ( x i β ) , i { 1 , , n } .
For this model, the parameters are interpreted from odds ratios between the prediction or response odds when one covariate is increased by m units (with the rest of the covariates remaining constant) and the odds without this increase. Note that this ratio is given by exp ( m β k ) , where β k is the parameter associated with the covariate k increased by m units. As in the two link functions previously analyzed, the distribution of the response variable is Y i SN ( μ i , ω , λ ) , for  i { 1 , , n } .
From the logit unit SN regression, the logit normal case is obtained with λ = 0 . The parameter estimates of the SN regression model for the response in the unit interval with logit link function can be obtained using the ML method. The corresponding log-likelihood function for a sample of n observations is obtained similarly from the PDF ϕ SNT defined in (4), where z i , z 1 i , and  z 0 i are given in (5) with μ i being stated in (12). The elements of the score function for ω and λ are similar to those of the complementary log–log unit SN model, but with μ i being formulated in (12), while for the parameter β j , we have that
˙ ( β j ) = 1 ω i = 1 n x i j ( z i λ w i λ + Δ 01 i ) μ i ( 1 μ i ) , j { 0 , 1 , , p } .
As mentioned, the score equation system is reached by setting the elements of the score function to zero, obtained from the first derivative of the function with respect to the parameters β 0 , β 1 , , β p , ω , and  λ . Solving this system of equations provides the ML estimates. To maximize this function, iterative algorithmic methods are necessary.

3.4. Information Matrices in Skew-Normal Models

To analyze the properties and behavior of SN regression models, it is essential to examine the information matrices. These matrices provide insights into the variance and covariance of the parameter estimators, as well as their overall performance. Specifically, the observed information matrix can be approximated by the negative of the Hessian matrix, which is obtained from the second derivatives of the log-likelihood function.
To obtain the Hessian matrix elements in the complementary log–log case, we define
τ j k i = z 1 i j ϕ SN k ( z 1 i ) z 0 i j ϕ SN k ( z 0 i ) Φ SN ( z 1 i ; λ ) Φ SN ( z 0 i ; λ ) .
The second derivatives of the log-likelihood function for β j β j , β j ω , β j λ , ω ω , ω λ , and  λ λ are, respectively, expressed as
¨ ( β j β j ) = 1 ω i = 1 n x i j x i j δ i ( log ( 1 μ i ) + 1 ) ( z i + λ w i λ Δ 01 i ) 1 ω 2 i = 1 n x i j x i j δ i 2 τ 01 i ( ϕ ( z 1 i ) + ϕ ( z 0 i ) ) 1 ω 2 i = 1 n x i j x i j δ i 2 ( 1 + λ 2 w i λ ( λ z i + w i λ ) ( Δ 11 i + Δ 01 i 2 ) ) , ¨ ( β j ω ) = 2 ω 2 i = 1 n x i j δ i z i λ ω i = 1 n x i j δ i w i λ + λ 2 ω 2 i = 1 n x i j δ i w i λ z i ( λ z i + w i λ ) 1 ω 2 i = 1 n x i j δ i ( Δ 01 i ( Δ 11 i 1 ) + ( Δ 21 i 2 τ 12 i ) ) , ¨ ( β j λ ) = 1 ω i = 1 n x i j δ i w i λ λ ω i = 1 n x i j δ i w i λ z i ( λ z i + w i λ ) 2 π i = 1 n x i j δ i Γ 11 i + 1 1 + λ 2 Δ 01 i Γ 01 i , ¨ ( ω ω ) = 1 ω 2 i = 1 n ( 1 + 3 z i 2 2 λ z i w i λ + λ 2 z i 2 w i λ ( λ z i + w i λ ) ) 1 ω 2 i = 1 n Δ 11 i 1 ω i = 1 n ( Δ 11 i Δ 21 i Δ 11 i 2 2 τ 12 i ) , ¨ ( ω λ ) = 1 ω i = 1 n w i λ z i ( 1 λ z i ( λ z i + w i λ ) ) + 2 π 1 ω i = 1 n Γ 21 i + 1 1 + λ 2 Δ 11 i Γ 01 i , ¨ ( λ λ ) = i = 1 n w i λ z i 2 ( λ z i + w i λ ) 2 π 1 ( 1 + λ 2 ) 2 i = 1 n 2 π Γ 01 i 2 λ Γ 01 i + 2 2 π λ 1 + λ 2 i = 1 n Γ 21 i .
The elements of the Hessian matrix in the probit case are given by
¨ ( β j β j ) = 1 ω i = 1 n x i j x i j ϕ ( x i β ) ϕ ( x i β ) + ( x i β ) z i + λ w i λ x i β + λ ϕ ( x i β ) ( λ z i + w i λ ) + 1 ω 2 i = 1 n x i j x i j ϕ ( x i β ) Δ 01 i x i β + ϕ ( x i β ) Δ 01 i 1 ω 2 i = 1 n x i j x i j ϕ 2 ( x i β ) ( Δ 11 i + τ 01 i ( ϕ ( z 1 i ) + ϕ ( z 0 i ) ) ) , ¨ ( β j ω ) = 1 ω 2 i = 1 n x i j ϕ ( x i β ) ( 2 λ w i λ + λ 2 w i λ z i ( λ z i + w i λ ) ) 1 ω 2 i = 1 n x i j ϕ ( x i β ) Δ 01 i ( 1 Δ 11 i ) ( Δ 21 i 2 τ 12 i ) , ¨ ( β j λ ) = 1 ω i = 1 n x i j ϕ ( x i β ) w i λ ( 1 λ z i ( λ z i + w i λ ) ) + 2 π i = 1 n x i j ϕ ( x i β ) Γ 11 i + 1 1 + λ 2 Δ 01 i Γ 01 i .
The Hessian matrix elements for the parameters ω and λ in the logit case follow the same structure as those for the complementary log–log unit SN model, with the standardization terms z i , z 1 i , and  z 0 i defined in (5) for μ i = Φ ( x i β ) .
For the probit case, the elements of the Hessian matrix are similarly obtained as for the complementary log–log and logit cases and are represented as
¨ ( β j β l ) = 1 ω 2 i = 1 n x i j x i l μ i 2 ( 1 μ i ) 2 1 + λ 2 w i λ ( λ z i + w i λ ) 1 ω i = 1 n x i j x i l μ i ( 1 μ i ) ( 1 2 μ i ) z i + λ w i λ + 1 ω 2 i = 1 n x i j x i l μ i ( 1 μ i ) ω ( 1 2 μ i ) + μ i ( 1 μ i ) ( Δ 11 i Δ 01 i 2 τ 02 i ) , ¨ ( β j ω ) = 1 ω 2 i = 1 n x i j μ i ( 1 μ i ) 2 z i λ w i λ + λ 2 z i w i λ ( λ z i + w i λ ) 1 ω 2 i = 1 n x i j μ i ( 1 μ i ) ( Δ 01 i ( 1 Δ 11 i ) Δ 21 i 2 τ 12 i ) , ¨ ( β j λ ) = 1 ω i = 1 n x i j μ i ( 1 μ i ) w i λ 1 λ z i ( λ z i + w i λ ) Γ 11 i + 1 1 + λ 2 Δ 01 i Γ 01 i .
The Hessian matrix elements for the parameters ω and λ in the probit case are the same as the complementary log–log and logit unit SN models, with its respective link function.
The elements of the observed information matrix can be written as ¨ θ j θ l = i = 1 n j θ j θ l , where j θ j θ l = 2 ( θ ) / θ j θ l = ¨ ( θ j θ l ) and the sum is over the set of observations i { 1 , , n } .
To study the conditions of the information matrix of the models under consideration, we analyze this matrix for the general case, that is, the non-inflated case at c 0 and c 1 . Thus, when λ = 0 , the information matrix is stated as
I SN β , ω , λ = 1 ω 2 ( M x ) ( M x ) 0 p + 1 2 π ω x M 1 n 0 p + 1 2 n ω 2 0 p + 1 2 π ω 1 n M x 0 p + 1 2 n π 2 ,
where 0 p + 1 is a ( p + 1 ) × 1 vector of zeros, 1 n is an n × 1 vector of ones, x = ( 1 n , x ) , where x = ( x 1 , , x p ) is an n × p matrix, with  x j = ( x j 1 , , x j n ) , for  j { 1 , , p } , and  M being a diagonal matrix stated as M = diag { μ 1 ( 1 μ 1 ) , , μ n ( 1 μ n ) } for the logit link, whereas
M = diag { δ 1 , , δ n } , M = diag { ϕ ( x 1 β ) , , ϕ ( x n β ) } ,
for the complementary log–log and probit links, respectively. The determinant of I SN is presented as
I SN β , ω , λ = 2 n 2 / π 2 ω 2 ( p + 1 ) ( M x ) ( M x ) n 1 n M x ( x M M x ) 1 x M 1 n .
Letting k = 1 n M M 1 n , l = 1 n M 1 n and writing x in partitioned form as x = ( 1 , x ) , it follows that the matrix x M M x can be expressed as a partitioned matrix.
We may use the existing expressions of matrix algebra to find the inverse of a partitioned matrix, since M is a diagonal matrix. Therefore, the rank of the matrix M x is the same rank of x , while M x M x is also full rank and then invertible. After some calculations, we obtain
1 n M x x M M x 1 x M 1 n = l 2 k + 1 n R V 1 R 1 n ,
where
R = x M ( k l M I ) , V = x M ( I k 1 M J M ) ,
with I n being the n × n identity matrix and J = 1 n 1 n the matrix of ones. Noting that 1 n R V 1 R 1 n > 0 is a quadratic form, it is clear that
n 1 n M x ( x M M x ) 1 x M 1 n 0 .
Note that, for our models, this determinant is a function of δ i , μ i ( 1 μ i ) , and  ϕ ( x i β ) for i { 1 , , n } . In particular, for the logit model, assuming that 0 < x i β < 1 , it follows that
16 n e 2 ( 1 + e ) 2 < l 2 k = i = 1 n μ i ( 1 μ i ) 2 i = 1 n μ i 2 ( 1 μ i ) 2 < e 2 ( 1 + e ) 4 n 16 ,
where e is the Euler or exponential number. For the complementary log–log and probit models, l 2 / k is given by i = 1 n δ i 2 / i = 1 n δ i 2 and ( i = 1 n ϕ ( x i β ) ) 2 / i = 1 n ϕ 2 ( x i β ) , respectively. Then, the information matrix is non-singular for each studied case of the unit SN model.

4. Unit Skew-Normal Zero–One Inflated Regression Models

In this section, we describe zero–one inflated data and explain how the ML method is used to estimate the parameters of the SNZOI model inflated at c 0 = 0 and c 1 = 1 .

4.1. Formulation of the Skew-Normal Zero–One Inflated Model

With the inflations set at c 0 = 0 and c 1 = 1 , the PDF for the SNZOI model is stated as
f ( y i ) = p 0 i , if   y i = 0 ; 1 p 0 i p 1 i ω ϕ SN ( y i μ i ) / ω ( Φ SN ( 1 μ i ) / ω Φ SN μ i / ω ) , if   0 < y i < 1 ; p 1 i , if   y i = 1 ;
with μ i = g 1 ( x i β ) for i { 1 , , n } , where p 0 i = P ( Y i = 0 ) and p 1 i = P ( Y i = 1 ) are defined for the SN model. We assume that the responses inflated at c 0 = 0 and c 1 = 1 are explained by the observed covariates x 0 i = ( 1 , x 0 i 1 , , x 0 i l ) and x 1 i = ( 1 , x 1 i 1 , , x 1 i r ) , respectively, and data in 0 < y i < 1 correspond to the part of the SN model stated by a 1 × ( p + 1 ) observed covariate vector x i = ( 1 , x i 1 , , x i p ) . Hence, the regression parameters for the SNZOI model are β 0 = ( β 00 , β 01 , , β 0 l ) , β = ( β 0 , β 1 , , β p ) , and  β 1 = ( β 10 , β 11 , , β 1 r ) .
Consider Y i SN ( μ i , ω , λ ) , for  0 < y i < 1 and i { 1 , , n } . Assuming a polychotomous random variable with a logit link function, we have
p 0 i = P ( Y i = 0 ) = exp ( x 0 i β 0 ) 1 + exp ( x 0 i β 0 ) + exp ( x 1 i β 1 ) , p 1 i = P ( Y i = 1 ) = exp ( x 1 i β 1 ) 1 + exp ( x 0 i β 0 ) + exp ( x 1 i β 1 ) , p i = P ( 0 < Y i < 1 ) = 1 p 0 i p 1 i = 1 1 + exp ( x 0 i β 0 ) + exp ( x 1 i β 1 ) ,
where x 0 i = ( 1 , x 0 i 1 , , x 0 i l ) and x 1 i = ( 1 , x 1 i 1 , , x 1 i r ) are 1 × ( l + 1 ) and 1 × ( r + 1 ) observed covariate vectors associated with the parameter vectors β 0 and β 1 , respectively. A particular case of this model addresses zero–one response probabilities using β 0 and β 1 , and responses within the unit interval using β with the SN distribution. This particular case serves as an alternative to the beta regression model studied in [25].

4.2. Parameter Estimation in the SNZOI Model

Estimation of the parameters for the discrete part of the SNZOI model can be carried out using the ML method. Additionally, with the MNP library of the R software (version 4.3.2) [45], the estimates of the probit model can be obtained using the procedure described in [43]. For the continuous part, ϕ SN ( y i ; λ ) namely, the score function and the elements of the Fisher information matrices were derived in the previous section for the complementary log–log, logit, and probit unit SN models. Similar to the previous case, the estimates of the parameters for the discrete and continuous parts can be obtained separately due to orthogonality. The estimation of the parameters ( β 0 , β , β 1 , ω , λ ) for the SNZOI model encompasses the parameter vectors β 0 , β , β 1 , along with the parameters ω and λ .
Estimation of the SNZOI parameters can be performed using the ML method by maximizing the log-likelihood function, decomposed into parts for the discrete and continuous components of the model given by ( θ ) = ( β 0 , β 1 ) + ( β , ω , λ ) , where
( β 0 , β 1 ) i L x 0 i β 0 + i R x 1 i β 1 i = 1 n log 1 + exp ( x 0 i β 0 ) + exp ( x 1 i β 1 ) , ( β , ω , λ ) i S 1 2 y i μ i ω 2 + log Φ λ ( y i μ i ) ω i S log Φ SN 1 μ i ω Φ SN μ i ω + n 01 ( log ( 2 ) log ( ω ) ) ,
with n 01 being the number of observations between zero and one. The set L refers to the zero-inflated (left) data with y i = 0 , the set R refers to the one-inflated (right) data with y i = 1 , and the set S includes the data 0 < y i < 1 , for  i { 1 , , n } .
The score function for ( β 0 , β 1 ) , leading to the ML estimates, is given by
˙ ( β 0 k ) = i L x 0 i k i = 1 n x 0 i k exp ( x 0 i β 0 ) 1 + exp ( x 0 i β 0 ) + exp ( x 1 i β 1 ) , ˙ ( β 1 m ) = i R x 1 i m i = 1 n x 1 i m exp ( x 1 i β 1 ) 1 + exp ( x 0 i β 0 ) + exp ( x 1 i β 1 ) ,
for k L { 1 , , l } and m R { 1 , , r } . The score function for ( β , ω , λ ) depends on the complementary log–log, logit, or probit link function used. These scores were already obtained previously for these three models. The system of equations reached by equating the scores to zero has no closed-form solution and needs to be solved numerically. The Fisher information matrix is block diagonal, that is, it can be written as
I ( θ ) = diag { I ( β 0 , β 1 ) , I ( β , ω , λ ) } ,
where I ( β 0 , β 1 ) is the component of the information matrix for the discrete part of the model and I ( β , ω , λ ) is its continuous part. The elements of I ( β 0 , β 1 ) were described in Section 4, while the elements of I ( β , ω , λ ) for the complementary log–log, logit, and probit links were detailed in Section 3. Hence, for large n, we have that θ ^ · N p + l + r + 5 ( θ , Σ θ θ ) , meaning that θ ^ is consistent and asymptotically normal distributed with Σ θ θ = I ( θ ) 1 = diag { I ( β 0 , β 1 ) 1 , I ( β , ω , λ ) 1 } .
Due to orthogonality, the parameters of the two models can be estimated separately. Estimation methods for when the likelihood function is orthogonal relative to the partition of interest have been presented in [46]. The approximated distribution N p + l + r + 5 ( θ , Σ θ θ ) can be used for constructing confidence intervals for θ j , presented as θ ^ j ± z 1 α / 2 ( σ ^ θ ^ j ) 1 / 2 , where σ ^ θ ^ j is the j-th element in the diagonal of the matrix Σ θ ^ θ ^ corresponding to θ j , and  z 1 α / 2 is the ( 1 α / 2 ) 100 % -th quantile of the standard normal distribution.
The zero and one inflated cases can be obtained from the SNZOI model. Taking p 1 i = 0 , we obtain a model with zero inflation while treating values above zero as continuous between 0 and 1. This case is referred to as left-censored because all observations at zero are treated as a single inflated point. For  p 0 i = 0 , we obtain a model with one inflation while treating values below one as continuous between 0 and 1. This case is referred to as right-censored because all observations at one are treated as a single inflated point [25].

5. Empirical Applications

In this section, we illustrate the application of the SNZOI model and compare it with other models using real data. We demonstrate that the SNZOI model can be a valid alternative to the zero–one inflated beta regression and tobit models. Also, in this section, we indicate how the models can be selected using some statistical criteria.

5.1. Model Selection Criteria

To evaluate and compare the models, we use the Akaike information criterion (AIC) and the corrected AIC (AICc), which balance model fit and complexity [16,47]. These criteria are particularly helpful in our study as we compare multiple models with varying levels of complexity to identify the best-fitting model. The AIC is defined as
AIC = 2 p 2 ,
where p is the number of parameters in the model and is the maximum value of the log-likelihood function for the model under analysis.
The AICc is a corrected version of AIC that adjusts for small sample sizes, given by
AICc = AIC + 2 p ( p + 1 ) n p 1 ,
where n is the sample size, providing a more accurate indicator of model quality [48]. Both AIC and AICc are obtained from the likelihood function of the observed data under the specified model. By using these criteria, we ensure a fair comparison across models, considering both their goodness of fit and their complexity.

5.2. Case Study 1: Doubly Censored Data

The data for this first application come from a clinical study previously analyzed in [49], which assessed the status and progression of periodontal disease among Gullah-speaking African-Americans with type-2 diabetes. The study used a questionnaire focused on demographical, dental, medical, and social aspects. Ethical approval was obtained, and all procedures followed the guidelines of the Helsinki Declaration, with informed consent provided by all participants. Data were collected between January 2018 and December 2019.
In this study, the clinical attachment level, a marker of periodontal disease, was measured at each of the six sites of a subject tooth. We model the dependence of the proportion of diseased sites corresponding to specific tooth types (canines, first molars, incisors, and premolars) with respect to the following covariates: gender ( X 1 ), age ( X 2 ), glycosylated hemoglobin ( X 3 ), and smoking status ( X 4 ). The dataset includes records from 290 individuals, with the response variable being the proportion of diseased tooth sites specifically for the premolars (Y). Notably, the data exhibit a high frequency of zeros ( Y = 0 ) and, for some cases, ones ( Y = 1 ), indicating high inflation at these points.
In Case Studies 1 and 2, we fitted several models to these data, including the new SNZOI model derived in the present article and the following models: Bernoulli one inflated (BOI), Bernoulli zero inflated (BZI), complementary log–log doubly censored log-normal (CLL-DCLN), complementary log–log doubly censored skew-normal (CLL-DCLSN), complementary log–log log-normal zero–one inflated (CLL-LNZOI), complementary log–log skew-normal zero–one inflated (CLL-SNZOI), DCLSN, logit doubly censored skew-normal (L-DCSN), logit doubly censored log-skew-normal (L-DCLSN), logit log-skew-normal zero–one inflated (L-LSNZOI), probit doubly censored log-normal (P-DCLN), probit log-normal zero–one inflated (P-LNZOI), probit log-skew-normal zero–one inflated (P-LSNZOI), and tobit.
From our analysis, we found that only covariates X 1 and X 2 are statistically significant at 5%. For the L-DCLSN model, the parameter estimates for the discrete part are reported in Table 1.
Table 2 presents the ML estimates of the model parameters for the tooth data using different models, including the detailed performance metrics AIC and AICc stated in (13) and (14) for a comprehensive comparison between the SNZOI model and other commonly employed models, such as the zero–one inflated beta and tobit models. The results demonstrate that the SNZOI model achieves competitive performance, often outperforming the other models in terms of lower AIC and AICc values.
In all cases, it can be noted that the two best models are CLL-LSNZOI and P-LSNZOI, followed by the BZI and P-LNZOI models. The CLL-LNZOI and L-LSNZOI models, as well as the doubly censored models (LN and LSN) follow, respectively. Additionally, we report that there are no significant differences statistically at 5% between the CLL-LSNZOI models with logit and probit link functions in the discrete part.
We estimate the model using the MNP library, which employs a Bayesian method for this estimation. This method performs better than the ML method of the traditional case. The estimated models for the discrete part are y ^ 1 i = 2.8527 + 0.0440 × age and y ^ 0 1 i = 0.4405 + 0.0221 × age . The CLL-LSNZOI and P-LSNZOI models can be compared against the CLL-LNZOI and P-LNZOI models using the hypotheses H 0 : λ = 0 versus H 1 : λ 0 , for each link function g, based on the test statistic Λ g = L LN g ( θ ^ ) / L LSN g ( θ ^ ) . After substituting the estimated values, we obtain 2 log ( Λ CLL ) = 8.13 and 2 log ( Λ P ) = 6.96 , which are both greater than the critical value of 5% of the chi-square distribution, χ 1 , 0.95 2 = 3.8414 . This leads to the conclusion that the CLL-LSNZOI and P-LSNZOI models fit the data better than the CLL-LNZOI and P-LNZOI models.
We estimate the percentage of the dataset with low response by fitting the considered model. An estimate of the percentage of observations at or below the detection limit is 100 × 1 / ( 1 + exp ( 1.407 ) ) = 19.67 % , which, compared to the observed 21.98 % , indicates good agreement with the considered model.
To uncover atypical outcomes and/or model misspecifications, we examined the transformed martingale residuals as proposed in [3] and defined as
r i = sign ( r i ) 2 ( r i + ν i log ( ν i r i ) ) , i { 1 , , n } ,
where r i is the martingale residual proposed in [50], ν i = 0 or 1 indicates whether observation i is censored or not, respectively, and  sign ( a ) denotes the sign of a. The plots of r i , for  i { 1 , , n } , with simulated envelopes, are presented in Figure 2 and Figure 3. From these figures, it is evident that the P-LSNZOI model fits the data better than the other models. To illustrate the process followed for Case Study 1, we present Algorithm 1 and the flowchart in Figure 4.
Algorithm 1 Case Study 1: Doubly censored data
  1:
Begin
  2:
Load the clinical study data.
  3:
Define covariates X 1 (gender), X 2 (age), X 3 (glycosylated hemoglobin), X 4 (smoking status), and response variable Y (proportion of diseased tooth sites for premolars).
  4:
Check for inflation at Y = 0 and Y = 1 .
  5:
Fit the BZI, L-DCLSN, L-LSNZOI, and P-LSNZOI models.
  6:
Compare models using AIC and AICc to select the best model.
  7:
Determine statistical significance of the covariates at 5% in the best model.
  8:
Calculate and plot martingale residuals to evaluate model fit.
  9:
Identify the model to be used for predicting based on results and residual analysis.
10:
End

5.3. Case Study 2: One-Inflated Data

This case study focuses on analyzing clinical data to model the proportion of damaged incisor teeth among an individual cohort. The data are characterized by a high inflation at the value Y = 1 , indicating full damage, which necessitates the use of specialized models to handle this inflation. We want to identify significant covariates that influence the proportion of damage on incisor teeth. To illustrate the process for Case Study 2, we present Algorithm 2 and the flowchart in Figure 4. This variable exhibits significant inflation at Y = 1 , with 65 cases of complete damage and only one case of no damage ( Y = 0 ). For this analysis, we exclude the observation where Y = 0 , resulting in a dataset of 289 one-inflated observations with Y = 1 .
Algorithm 2 Case Study 2: One-inflated data
  1:
Begin
  2:
Load the clinical study data.
  3:
Define covariates X 1 (gender), X 2 (age), X 3 (glycosylated hemoglobin), X 4 (smoking status), and response variable Y (proportion of damage on incisor teeth).
  4:
Exclude observations where Y = 0 .
  5:
Check for inflation at Y = 1 .
  6:
Fit the BOI, L-DCLSN, L-LSNZOI, P-LNSZOI, and tobit models.
  7:
Compare models using AIC and AICc to select the best model.
  8:
Determine statistical significance of the covariates at 5% in the best model.
  9:
Calculate and plot martingale residuals to evaluate model fit.
10:
Identify the model to be used for predicting based on results and residual analysis.
11:
End
To identify the most suitable model for these one-inflated data, we employed the L-LSNZOI and P-LSNZOI models, selecting the model with the lowest AIC and/or AICc from the combinations studied. We applied the same selection criteria to other models, including the normal tobit model (on a logarithmic scale), and the L-DCLSN and P-DCLSN models inflated at Y = 1 . The results are reported in Table 3.
As observed in Table 3, the L-LSNZOI and L-DCLSN models showed the lowest AIC and AICc values, indicating their superior fit to the data compared to the other models. Specifically, the L-LSNZOI model had the lowest AIC (167.74) and AICc (170.08), suggesting it as the best model. The tobit model had the highest AIC (174.29) and AICc (176.50), making it the least suitable among the models evaluated. The BOI and P-LSNZOI models also performed well but were outperformed by the L-LSNZOI model. These results show the importance of considering model selection criteria like AIC and AICc in determining the best model for one-inflated data. The significant parameters across these models also provide insights into the factors influencing the proportion of damaged incisor teeth, aiding in targeted clinical interventions.

5.4. Computational Costs of Algorithms

The most computationally intensive steps in both Algorithms 1 and 2 are the model fitting and residual analysis due to the iterative nature of parameter estimation and the subsequent evaluation of model adequacy. For example, the fitting process for the LSNZOI models involves estimating several parameters iteratively until the likelihood function converges, which can be time-consuming for large datasets. The residual analysis, especially when involving bootstrapping for simulated envelopes or generating multiple diagnostic plots, further adds to the computational load.
In our implementation, the following details highlight the computational efforts:
  • Model fitting—The fitting process for models, such as L-LSNZOI and P-LSNZOI structures, involves maximizing the log-likelihood function through iterative numerical optimization methods. This process typically requires multiple iterations (often ranging from 50 to 200 iterations) to achieve convergence, especially in the presence of complex likelihood surfaces and multiple parameters. Each iteration involves recalculating the likelihood function and its gradient, adding to the computational time. The computational cost can be roughly estimated as O ( n × p × r ) , where n is the sample size, p is the number of parameters, and r is the number of iterations.
  • Residual analysis—Calculating martingale residuals and their diagnostic plots involve extensive computations. Generating simulated envelopes through bootstrapping, for instance, necessitates multiple re-samplings of the data and re-fitting the model to each bootstrap sample, which is computationally expensive. If B bootstrap samples are used, the computational cost for this step is roughly O ( B × n × p × r ) .
  • Model selection criteria—Calculating AIC and AICc requires obtaining the log-likelihood value at the estimated parameters, which involves evaluating the likelihood function at these estimates. While this step is less computationally intensive than the fitting process, it still adds to the overall computational burden. The computational cost for this step is approximately O ( n × p ) .
Overall, the total computational cost is a combination of these intensive steps, making the entire process demanding, especially for large datasets or complex models.

6. Discussion

In this article, we proposed a general class of skew regression models for response variables taking values in the unit interval, possibly with an excess of zeros or ones. The proposed models were developed from a continuous-discrete mixture distribution with covariates in both discrete and continuous parts. As demonstrated by applications with real data, the proposed models offer a valid alternative for describing data of rates and proportions that are inflated at zero or one.

6.1. Key Findings and Insights

Our results indicate that the SNZOI model, particularly with the logit and probit link functions (L-SNZOI and P-SNZOI), consistently outperformed other models in terms of AIC and AICc values. These models provided a better fit to the data from the clinical study on periodontal disease, where the response variable was the proportion of diseased tooth sites.
In the first application involving doubly censored data, the models effectively captured the significant effects of gender and age on the proportion of diseased premolar sites. The L-SNZOI and P-SNZOI models outperformed the other models, suggesting that the SN distribution is a more suitable assumption for this type of data. In the second application, which involved inflation at the one value, the SNZOI models again outperformed traditional beta regression and tobit models. The L-DCLSN, L-LSNZOI, and P-LSNZOI models, adjusted for right-censored data and inflation at one, demonstrated superior performance.
To further validate the model adequacy, we analyzed martingale residuals. The transformation of the martingale residuals and the generation of simulated envelopes, as proposed in [3], provided a robust method to identify model goodness of fit and outliers.
Specifically, we considered transformed martingale residuals, which help to detect discrepancies between the observed and expected values under the model. The simulated envelopes are generated to visually assess the distribution of these residuals and identify any points that fall outside the expected range, indicating potential model misfits or outliers. This comprehensive analysis of martingale residuals, along with the simulated envelopes, reinforces the robustness of our chosen models and provides a thorough evaluation of their performance, ensuring that the models accurately capture the underlying data structure.

6.2. Comparison with Previous Studies

Our findings contribute to the literature on regression models for bounded and inflated data. The superior performance of the SN models aligns with previous studies, such as [25,46], which demonstrated the advantages of using SN distributions for modeling asymmetrically distributed data with bounds. The practical implications of our study are important for clinical research where proportion data with inflation at specific values is common. The ability of the SNZOI models to handle such data accurately allows researchers to obtain more precise estimators and make better-informed decisions.

6.3. Model Limitations

Despite the promising results, there are several limitations to our study. Firstly, the complexity of the models and the need for iterative numerical methods for parameter estimation can be computationally intensive. Secondly, while the models performed well on the datasets used in this study, further validation on other types of data is necessary to confirm their generalizability.

6.4. Directions for Future Research

Future research could focus on several areas. One potential direction is the development of more efficient algorithms to reduce the computational burden of fitting these models. Additionally, exploring the application of these models to other fields, such as economics or environmental studies, could provide further validation and potentially uncover new areas of application.
While the present article focuses on the case study of periodontal disease prevalence to demonstrate the practical implications of the SNZOI model, future research could include additional case studies from diverse fields. These additional studies would provide more comprehensive evidence of the model versatility and practical utility across different domains.
As model performance is a critical aspect of our analysis, while the methods used, such as AIC, AICc, and martingale residuals, are robust for evaluating model adequacy, there is potential for further enhancement. Future work could explore additional goodness-of-fit tests that are specifically tailored to the characteristics of bounded and inflated data. Incorporating these tests could provide a more comprehensive assessment of model performance and robustness.
Another interesting direction for future research is the extension of these models to handle longitudinal data or hierarchical data structures. This would involve developing methods to account for correlations within subjects or groups, which is a common feature in many practical datasets. Additionally, investigating the robustness of these models under various misspecification scenarios can help in developing more resilient modeling strategies. Understanding how these models perform when the underlying assumptions are violated is crucial for practical applications.

7. Conclusions

The analysis of proportion data, especially when values are inflated at zero and one, poses challenges in various scientific fields. Traditional models, such as the beta and tobit regression models, often fall short in accurately capturing the complexities introduced by such data. This highlights the necessity for more advanced modeling techniques that can handle the unique distributional characteristics of zero–one inflation.
In this study, we addressed these challenges by proposing the skew-normal zero–one inflated models, particularly focusing on logit and probit link functions. These models integrate a continuous-discrete mixture distribution with covariates in both components, providing a sophisticated framework for analyzing proportion data with specific inflation points. By doing so, the skew-normal zero–one inflated models offer a robust and flexible approach for capturing asymmetrically distributed data and mixed discrete-continuous characteristics, which are common in fields such as ecology, economics, and medicine.
Our applications, including a detailed case study on periodontal disease prevalence, demonstrated that the skew-normal zero–one inflated models, particularly with the logit and probit link functions, consistently outperformed traditional models. These models showed superior performance in terms of the Akaike information criterion and its correction, offering more precise and unbiased statistical inferences. The transformation of martingale residuals and the generation of simulated envelopes further validated the robustness of our models, highlighting their ability to identify model misfits and outliers effectively. The proposed models fill a critical void in the statistical modeling and provide valuable insights and precise estimators for dealing with bounded and inflated data. The flexibility and robustness of the skew-normal zero–one inflated models make them a valid alternative for describing data of proportions that are inflated at zero or one.
Despite the promising results, several limitations remain. The complexity of the models and the need for numerical methods for parameter estimation can be computationally intensive. Additionally, while the models performed well on the datasets used in this study, further validation on other types of data is necessary to confirm their generalizability.
Future research should focus on developing more efficient algorithms to reduce the computational burden of fitting these models. Exploring the application of these models to other fields, such as economics or environmental studies, could provide further validation and potentially uncover new areas of application. Additionally, incorporating more comprehensive goodness-of-fit tests tailored to bounded and inflated data could enhance the evaluation of model performance and robustness.
In summary, the skew-normal zero–one inflated models represent an advancement in statistical modeling for proportion data with zero–one inflation. These models offer a robust and flexible framework for analyzing such data, providing deeper insights and more precise estimators. Continued research and development in this area holds promise for further advancements in statistical methodologies for bounded and inflated data.

Author Contributions

Conceptualization: G.M.-F., R.T.-F., V.L. and C.C.; data curation: G.M.-F., R.T.-F. and C.C.; formal analysis: G.M.-F., R.T.-F., V.L. and C.C.; investigation: G.M.-F., R.T.-F., V.L. and C.C.; methodology: G.M.-F., R.T.-F., V.L. and C.C.; writing—original draft: G.M.-F. and R.T.-F.; writing—review and editing: V.L. and C.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was partially supported by Vice-rectorate for Research of the Universidad de Córdoba, Colombia, project grant FCB-06-22 (G.M.-F. and R.T.-F); by Vice-rectorate for Research, Creation, and Innovation (VINCI) of the Pontificia Universidad Católica de Valparaíso (PUCV), Chile, grants VINCI 039.470/2024 (regular research), VINCI 039.493/2024 (interdisciplinary associative research), VINCI 039.309/2024 (PUCV centenary), and FONDECYT 1200525 (V.L.), from the National Agency for Research and Development (ANID) of the Chilean government under the Ministry of Science, Technology, Knowledge, and Innovation; and by Portuguese funds through the CMAT—Research Centre of Mathematics of University of Minho, Portugal, within projects UIDB/00013/2020 (https://doi.org/10.54499/UIDB/00013/2020) and UIDP/00013/2020 (https://doi.org/10.54499/UIDP/00013/2020) (C.C.).

Data Availability Statement

The data and codes used in this study are available under request to the authors.

Acknowledgments

The authors would like to thank the editors and anonymous reviewers for their valuable comments and suggestions, which helped us to improve the quality of this article.

Conflicts of Interest

The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this article.

References

  1. Schaminée, J.H.; Hennekens, S.M.; Chytry, M.; Rodwell, J.S. Vegetation-plot data and databases in Europe: An overview. Preslia 2009, 81, 173–185. [Google Scholar]
  2. Tobin, J. Estimation of relationships for limited dependent variables. Econometrica 1958, 26, 24–36. [Google Scholar] [CrossRef]
  3. Barros, M.; Galea, M.; González, M.; Leiva, V. Influence diagnostics in the tobit censored response model. Stat. Methods Appl. 2010, 19, 379–397. [Google Scholar] [CrossRef]
  4. Ferreira, P.H.; Shimizu, T.K.; Suzuki, A.K.; Louzada, F. On an asymmetric extension of the tobit model based on the tilted-normal distribution. Chil. J. Stat. 2019, 10, 99–122. [Google Scholar]
  5. Barros, M.; Galea, M.; Leiva, V.; Santos-Neto, M. Generalized tobit models: Diagnostics and application in econometrics. J. Appl. Stat. 2018, 45, 145–167. [Google Scholar] [CrossRef]
  6. Arellano-Valle, R.B.; Gómez, H.W.; Quintana, F. Statistical inference for a general class of asymmetric distributions. J. Stat. Plan. Inference 2005, 128, 427–443. [Google Scholar] [CrossRef]
  7. Gallardo, D.I.; Bourguignon, M.; Galarza, C.E.; Gómez, H.W. A parametric quantile regression model for asymmetric response variables on the real line. Symmetry 2022, 14, 1938. [Google Scholar] [CrossRef]
  8. Gupta, R.D.; Gupta, R.C. Analyzing skewed data by power normal model. Test 2008, 17, 197–210. [Google Scholar] [CrossRef]
  9. Pewsey, A.; Gómez, H.W.; Bolfarine, H. Developments in skew-symmetric distributions and their applications. Symmetry 2022, 14, 567. [Google Scholar]
  10. Desousa, M.; Saulo, H.; Leiva, V.; Scalco, P. On a tobit-Birnbaum–Saunders model with an application to medical data. J. Appl. Stat. 2018, 45, 932–955. [Google Scholar] [CrossRef]
  11. Sanchez, L.; Leiva, V.; Galea, M.; Saulo, H. Birnbaum–Saunders quantile regression and its diagnostics with application to economic data. Appl. Stoch. Model. Bus. Ind. 2021, 37, 53–73. [Google Scholar] [CrossRef]
  12. Azzalini, A. A class of distributions which includes the normal ones. Scand. J. Stat. 1985, 12, 171–178. [Google Scholar]
  13. Azzalini, A. Further results on a class of distributions which includes the normal ones. Statistica 1986, 46, 199–208. [Google Scholar]
  14. Henze, N. A probabilistic representation of the skew-normal distribution. Scand. J. Stat. 1986, 13, 271–275. [Google Scholar]
  15. Castillo, N.O.; Gómez, H.W.; Leiva, V.; Sanhueza, A. On the Fernández-Steel distribution: Inference and application. Comput. Stat. Data Anal. 2011, 55, 2951–2961. [Google Scholar] [CrossRef]
  16. Ventura, M.; Saulo, H.; Leiva, V.; Monsueto, S. Log-symmetric regression models: Information criteria, application to movie business and industry data with economic implications. Appl. Stoch. Model. Bus. Ind. 2019, 3, 963–977. [Google Scholar] [CrossRef]
  17. Massuia, M.B.; Garay, A.M.; Cabral, C.R.; Lachos, V.H. Bayesian analysis of censored linear regression models with scale mixtures of skew-normal distributions. Stat. Its Interface 2017, 10, 425–439. [Google Scholar] [CrossRef]
  18. Morán-Vásquez, R.A.; Giraldo-Melo, A.D.; Mazo-Lopera, M.A. Quantile estimation using the log-skew-normal linear regression model with application to children’s weight data. Mathematics 2023, 11, 3736. [Google Scholar] [CrossRef]
  19. Dias-Domingues, T.; Mouriño, H.; Sepúlveda, N. Classification methods for the serological status based on mixtures of skew-normal and skew-t distributions. Mathematics 2024, 12, 217. [Google Scholar] [CrossRef]
  20. Mudholkar, G.S.; Hutson, A.D. The epsilon-skew-normal distribution for analyzing near-normal data. J. Stat. Plan. Inference 2000, 83, 291–309. [Google Scholar] [CrossRef]
  21. Gómez, H.W.; Venegas, O.; Bolfarine, H. Skew-symmetric distributions generated by the distribution function of the normal distribution. Environmetrics 2007, 18, 395–407. [Google Scholar] [CrossRef]
  22. Arrué, J.; Arellano-Valle, R.; Gómez, H.W.; Leiva, V. On a new type of Birnbaum–Saunders models and its inference and application to fatigue data. J. Appl. Stat. 2020, 47, 2690–2710. [Google Scholar] [CrossRef] [PubMed]
  23. Pewsey, A. Problems of inference for Azzalini’s skew-normal distribution. J. Appl. Stat. 2000, 27, 859–870. [Google Scholar] [CrossRef]
  24. Ferrari, S.; Cribari-Neto, F. Beta regression for modelling rates and proportions. J. Appl. Stat. 2004, 31, 799–815. [Google Scholar] [CrossRef]
  25. Ospina, R.; Ferrari, S. Inflated beta distributions. Stat. Pap. 2010, 51, 111–126. [Google Scholar] [CrossRef]
  26. Ospina, R.; Ferrari, S. A general class of zero-or-one inflated beta regression models. Comput. Stat. Data Anal. 2012, 56, 1609–1623. [Google Scholar] [CrossRef]
  27. Couri, L.; Ospina, R.; da Silva, G.; Leiva, V.; Figueroa-Zúñiga, J. A study on computational algorithms in the estimation of parameters for a class of beta regression models. Mathematics 2022, 10, 299. [Google Scholar] [CrossRef]
  28. Mohammadi, Z.; Sajjadnia, Z.; Bakouch, H.S.; Sharafi, M. Zero-and-one inflated Poisson-Lindley INAR (1) process for modelling count time series with extra zeros and ones. J. Stat. Comput. Simul. 2022, 92, 2018–2040. [Google Scholar] [CrossRef]
  29. Lee, B.S.; Haran, M. A class of models for large zero-inflated spatial data. J. Agric. Biol. Environ. Stat. 2024. [Google Scholar] [CrossRef]
  30. Figueroa-Zúñiga, J.; Niklitschek, S.; Leiva, V.; Liu, S. Modeling heavy-tailed bounded data by the trapezoidal beta distribution with applications. REVSTAT Stat. J. 2022, 20, 387–404. [Google Scholar]
  31. Jornsatian, C.; Bodhisuwan, W. Zero-one inflated negative binomial-beta exponential distribution for count data with many zeros and ones. Commun. Stat. Theory Methods 2022, 51, 8517–8531. [Google Scholar] [CrossRef]
  32. Keim, J.L.; DeWitt, P.D.; Fitzpatrick, J.J.; Jenni, N.S. Estimating plant abundance using inflated beta distributions: Applied learnings from a Lichen-Caribou ecosystem. Ecol. Evol. 2017, 7, 486–493. [Google Scholar] [CrossRef] [PubMed]
  33. Benites, L.; Maehara, R.; Lachos, V.H.; Bolfarine, H. Linear regression models using finite mixtures of skew heavy-tailed distributions. Chil. J. Stat. 2019, 10, 21–41. [Google Scholar]
  34. Desousa, M.; Saulo, H.; Santos-Neto, M.; Leiva, V. On a new mixture-based regression model: Simulation and application to data with high censoring. J. Stat. Comput. Simul. 2020, 90, 2861–2877. [Google Scholar] [CrossRef]
  35. Arellano-Valle, R.B.; Gómez, H.W.; Quintana, F. A new class of skew-normal distributions. Commun. Stat. Theory Methods 2004, 33, 1465–1480. [Google Scholar] [CrossRef]
  36. Saulo, H.; Dasilva, A.; Leiva, V.; Sanchez, L.; de la Fuente, H. Log-symmetric quantile regression models. Stat. Neerl. 2022, 76, 124–163. [Google Scholar] [CrossRef]
  37. Chai, H.S.; Bailey, K.R. Use of log-skew-normal distribution in analysis of continuous data with a discrete component at zero. Stat. Med. 2008, 27, 3643–3655. [Google Scholar] [CrossRef] [PubMed]
  38. Cragg, J.G. Some statistical models for limited dependent variables with application to the demand for durable goods. Econometrica 1971, 39, 829–844. [Google Scholar] [CrossRef]
  39. Moulton, L.H.; Halsey, N.A. A mixture model with detection limits for regression analyses of antibody response to vaccine. Biometrics 1995, 51, 1570–1578. [Google Scholar] [CrossRef]
  40. McCulloch, R.; Rossi, P.E. An exact likelihood analysis of the multinomial probit model. J. Econom. 1994, 64, 207–240. [Google Scholar] [CrossRef]
  41. Keane, M.P. A note on identification in the multinomial probit model. J. Bus. Econ. Stat. 1992, 10, 193–200. [Google Scholar] [CrossRef]
  42. Heckman, J.; Sedlacek, G. Heterogeneity, aggregation, and market wage functions: An empirical model of self-selection in the labor market. J. Political Econ. 1985, 93, 1077–1125. [Google Scholar] [CrossRef]
  43. Imai, K.; Van-Dyk, D. A Bayesian analysis of the multinomial probit model using marginal data augmentation. J. Econom. 2005, 124, 311–334. [Google Scholar] [CrossRef]
  44. Nocedal, J.; Wright, S. Numerical Optimization; Springer: New York, NY, USA, 2006. [Google Scholar]
  45. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2022. [Google Scholar]
  46. Farias, R.; Moreno-Arenas, G.; Patriota, A. Reduction of models in the presence of nuisance parameters. Colomb. J. Stat. 2009, 32, 99–121. [Google Scholar]
  47. Burnham, K.P.; Anderson, D.R. Model Selection and Multimodel Inference: A Practical Information-Theoretic Approach; Springer: New York, NY, USA, 2004. [Google Scholar]
  48. Hurvich, C.M.; Tsai, C.L. Regression and time series model selection in small samples. Biometrika 1989, 76, 297–307. [Google Scholar] [CrossRef]
  49. Galvis, D.M.; Bandyopadhyay, D.; Lachos, V.H. Augmented mixed beta regression models for periodontal proportion data. Stat. Med. 2014, 33, 3759–3771. [Google Scholar] [CrossRef]
  50. Ortega, E.M.; Bolfarine, H.; Paula, G.A. Influence diagnostics in generalized log-gamma regression models. Comput. Stat. Data Anal. 2003, 42, 165–186. [Google Scholar] [CrossRef]
Figure 1. Flowchart of the article structure.
Figure 1. Flowchart of the article structure.
Mathematics 12 02486 g001
Figure 2. Plots of envelopes for r i using (a) BZI and (b) L-DCLSN models and tooth data.
Figure 2. Plots of envelopes for r i using (a) BZI and (b) L-DCLSN models and tooth data.
Mathematics 12 02486 g002
Figure 3. Plots of envelopes for r i using (a) L-LSNZOI and (b) P-LSNZOI models and tooth data.
Figure 3. Plots of envelopes for r i using (a) L-LSNZOI and (b) P-LSNZOI models and tooth data.
Mathematics 12 02486 g003
Figure 4. Overlayed flowchart for Case Studies 1 and 2.
Figure 4. Overlayed flowchart for Case Studies 1 and 2.
Mathematics 12 02486 g004
Table 1. ML estimates of the indicated parameter and model for the tooth data and their AIC and AICc.
Table 1. ML estimates of the indicated parameter and model for the tooth data and their AIC and AICc.
ParameterCLL-SNZOICLL-LSNZOICLL-DCLSNP-LSNZOIL-LSNZOIL-DCLSN
β 00 −0.44050.6337−3.91220.63370.63370.7301
(0.4174)(0.7408)(0.8863)(0.7408)(0.7408)(0.8914)
β 02 0.0221−0.0376−0.0764−0.0376−0.0376−0.0957
(0.0075)(0.0135)(0.0160)(0.0135)(0.0135)(0.0161)
β 10 −2.2137−2.2137−2.1137−2.1123−1.2683−1.7604
(1.2252)(1.2252)(0.5710)(1.0552)(1.2515)(0.7648)
β 11 −1.5892−1.5892−0.4207−1.4316−1.1768−0.4990
(0.7523)(0.7523)(0.1805)(0.5520)(0.4923)(0.2381)
β 12 0.05090.05090.02740.05330.05170.0291
(0.0260)(0.0260)(0.0101)(0.0230)(0.0218)(0.0131)
β 20 −2.8527−8.0316−8.0886−8.0316−8.0316−8.3703
(1.9217)(2.3145)(1.4958)(2.3145)(2.3145)(1.4941)
β 22 0.04400.0788−0.06880.07880.07880.0324
(0.0274)(0.0358)(0.0236)(0.0358)(0.0358)(0.0236)
ω 0.76710.64530.30230.61610.31770.3198
(0.2731)(0.2279)(0.0545)(0.2362)(0.0477)(0.0667)
λ −1.9835−1.9835−0.7232−1.7782−1.6450−0.9430
(0.9456)(0.9456)(0.0523)(0.7634)(0.4274)(0.1660)
AIC310.13309.94325.34308.84313.53326.40
AICc312.91312.73328.12311.63316.32329.47
Where CLL-SNZOI is complementary log–log skew-normal zero–one inflated, CLL-LSNZOI is complementary log–log log-skew-normal zero–one inflated, CLL-DCLSN is complementary log–log doubly censored skew-normal, P-LSNZOI is probit log-skew-normal zero–one inflated, L-LSNZOI is logit log-skew-normal zero–one inflated, L-DCLSN is logit doubly censored log-skew-normal, AIC is Akaike information criterion, and  AICc is corrected Akaike information criterion.
Table 2. ML estimates of the indicated parameter and model for the tooth data and their AIC and AICc.
Table 2. ML estimates of the indicated parameter and model for the tooth data and their AIC and AICc.
ParameterBZIP-LNZOICLL-DCLNP-DCLNCLL-LNZOI
β 00 0.63370.63370.77681.17310.6337
(0.7408)(0.7408)(0.7821)(2.1798)(0.7408)
β 02 −0.0376−0.0376−0.1791−0.2042−0.0376
(0.0135)(0.0135)(0.0139)(0.0450)(0.0135)
β 10 −1.3885−2.7994−3.2301−2.0330−2.8949
(0.3957)(1.0573)(0.5903)(0.4353)(1.1453)
β 11 −0.5366−1.1257−0.5926−0.4060−1.3134
(0.1613)(0.4304)(0.2087)(0.4905)(0.4387)
β 12 0.02170.04520.03800.02630.0393
(0.0068)(0.0186)(0.0106)(0.0091)(0.0194)
β 20 −8.0316−8.0316−8.4042−11.0378−8.0316
(2.3153)(2.3145)(2.8866)(1.3108)(2.3145)
β 22 0.07880.0788−0.10490.06110.0788
(0.0358)(0.0358)(0.0449)(0.0210)(0.0358)
ω 0.09030.32500.26490.26460.3096
(0.0652)(0.0418)(0.1012)(0.0126)(0.0796)
AIC311.70313.80324.33323.51316.07
AICc314.58316.45326.97326.15318.71
Where BZI is Bernoulli zero inflated; P-LNZOI is probit log-normal zero–one inflated; CLL-DCLN is complementary log–log doubly censored log-normal; P-DCLN is probit doubly censored log-normal; CLL-LNZOI is complementary log–log log-normal zero–one inflated, AIC is Akaike information criterion, and AICc is corrected Akaike information criterion.
Table 3. ML estimates of the indicated parameter and model with tooth data and their AIC and AICc.
Table 3. ML estimates of the indicated parameter and model with tooth data and their AIC and AICc.
ParameterBOIL-LSNZOIP-LSNZOIL-DCLSNTobit
β 10 −0.0086−3.3573−2.8368−3.35730.2968
(0.3339)(1.0938)(1.5888)(1.0938)(0.0897)
β 11 −0.28280.12030.10360.1203−0.1003
(0.1513)(0.0351)(0.0489)(0.0351)(0.0393)
β 12 0.0182 0.0104
(0.0061) (0.0015)
β 20 −6.1599−6.1599−5.73503.6472
(1.0003)(0.9992)(1.0661)(0.5480)
β 22 0.08530.08530.0756−0.0503
(0.0166)(0.0165)(0.0174)(0.0092)
ω 1.31120.38150.36660.3815
(0.0868)(0.0215)(0.0180)(0.0215)−1.3046
λ −7.7660−8.1582−7.76600.0493
(3.1311)(3.4725)(3.1311)
AIC170.67167.74172.65166.85174.29
AICc173.07170.08175.04169.24176.50
Where BOI is Bernoulli one inflated; L-LSNZOI is logit log-skew-normal zero–one inflated, P-LSNZOI is probit log-skew-normal zero–one inflated, L-DCLSN is logit doubly censored log-skew-normal, AIC is Akaike information criterion; and AICc is corrected Akaike information criterion.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Martínez-Flórez, G.; Tovar-Falón, R.; Leiva, V.; Castro, C. Skew-Normal Inflated Models: Mathematical Characterization and Applications to Medical Data with Excess of Zeros and Ones. Mathematics 2024, 12, 2486. https://doi.org/10.3390/math12162486

AMA Style

Martínez-Flórez G, Tovar-Falón R, Leiva V, Castro C. Skew-Normal Inflated Models: Mathematical Characterization and Applications to Medical Data with Excess of Zeros and Ones. Mathematics. 2024; 12(16):2486. https://doi.org/10.3390/math12162486

Chicago/Turabian Style

Martínez-Flórez, Guillermo, Roger Tovar-Falón, Víctor Leiva, and Cecilia Castro. 2024. "Skew-Normal Inflated Models: Mathematical Characterization and Applications to Medical Data with Excess of Zeros and Ones" Mathematics 12, no. 16: 2486. https://doi.org/10.3390/math12162486

APA Style

Martínez-Flórez, G., Tovar-Falón, R., Leiva, V., & Castro, C. (2024). Skew-Normal Inflated Models: Mathematical Characterization and Applications to Medical Data with Excess of Zeros and Ones. Mathematics, 12(16), 2486. https://doi.org/10.3390/math12162486

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop