Next Article in Journal
Fourier Series Related to p-Trigonometric Functions
Next Article in Special Issue
On the Conflation of Negative Binomial and Logarithmic Distributions
Previous Article in Journal
The Invariant Subspace Problem for Separable Hilbert Spaces
Previous Article in Special Issue
New Flexible Asymmetric Log-Birnbaum–Saunders Nonlinear Regression Model with Diagnostic Analysis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Unit-Power Half-Normal Distribution Including Quantile Regression with Applications to Medical Data

by
Karol I. Santoro
1,
Yolanda M. Gómez
2,*,
Darlin Soto
2 and
Inmaculada Barranco-Chamorro
3
1
Departamento de Estadística y Ciencia de Datos, Facultad de Ciencias Básicas, Universidad de Antofagasta, Antofagasta 1240000, Chile
2
Departamento de Estadística, Facultad de Ciencias, Universidad del Bío-Bío, Concepción 4081112, Chile
3
Departamento de Estadística e Investigación Operativa, Facultad de Matemáticas, Universidad de Sevilla, 41012 Sevilla, Spain
*
Author to whom correspondence should be addressed.
Axioms 2024, 13(9), 599; https://doi.org/10.3390/axioms13090599
Submission received: 28 June 2024 / Revised: 20 August 2024 / Accepted: 28 August 2024 / Published: 2 September 2024
(This article belongs to the Special Issue Probability, Statistics and Estimations, 2nd Edition)

Abstract

:
In this paper, we present the unit-power half-normal distribution, derived from the power half-normal distribution, for data analysis in the open unit interval. The statistical properties of the unit-power half-normal model are described in detail. Simulation studies are carried out to evaluate the performance of the parameter estimators. Additionally, we implement the quantile regression for this model, which is applied to two real healthcare data sets. Our findings suggest that the unit power half-normal distribution provides a robust and flexible alternative for existing models for proportion data.

1. Introduction

To study and interpret real events, new statistical models must continuously be built. In recent times, models limited to the interval ( 0 , 1 ) have generated a lot of interest. This type of models is mainly used to study proportion data, such as survival rate data and product failure, among others. Different areas, such as health, actuarial and financial sciences, require this type of distributions. So, recently, several distributions with positive support have been transformed to distributions with unit support, such as the cases of Jones [1] in the Kumaraswamy distribution, Gómez-Déniz et al. [2] in the log-Lindley distribution, Mazucheli et al. [3], in the Lindley distribution, Abd El-Monsef et al. [4], who proposed a new two-parameter omega unitary distribution, Altun et al. [5], who studied a distribution called enhanced second-order Lindley distribution modeling data on the interval (0,1), and, recently, Ahmad et al. [6], who studied the exponential pareto distribution.
One of the distributions that is mainly used for this type of data is the Beta distribution. Interest in this type of model has been increasing, and many researchers have transformed known distributions with positive support into distributions with unitary support, an example of which can be found in articles such as Grassia [7], based on the Gamma distribution, Ghitany et al. [8], based on the Inverse Gamma distribution, Mazucheli et al. [9], based on the Birnbaum-Saunders distribution, Modi et al. [10], based on Burr distribution III, Korkmaz and Chesneau [11], based on Burr distribution XII, Haq et al. [12], based on the modified Burr III distribution, and, more recently, Bakouch et al. [13], based on half-normal distribution.
The aim of this paper is to introduce a new distribution with a bounded domain on (0,1). This distribution is originated by modifying the representation of the power half-normal (PHN) distribution proposed by Gómez et al. [14]. One of the motivations of this work is to generate alternatives to well-known distributions used in the statistical analysis of certain type of data. This work is based on the PHN distribution whose probability density function f (pdf), cumulative distribution function F (cdf), and quantile function Q are:
f Y ( y ; σ , α ) = 2 α σ ϕ y σ 2 Φ y σ 1 α 1 , y > 0 ,
F Y ( y ; σ , α ) = 2 Φ y σ 1 α , y > 0 ,
Q ( p ; σ , α ) = σ Φ 1 1 + p 1 / α 2 , 0 < p < 1 ,
where the parameters are σ > 0 and α > 0 . ϕ ( · ) and Φ ( · ) denote the pdf and cdf of the standard normal distribution, respectively.
At present, this model continues to be studied. For example, Barrios et al. [15] performed an extension of the distribution using the slash process, resulting in a model with higher kurtosis, i.e., with heavier tails. In addition, Pallini [16] introduced a discretized model based on the PHN distribution. The interest in further investigating models with support on the interval (0,1) is of great importance, since in many areas, such as finance, actuarial sciences, engineering, and health, databases that belong to this type of interval are handled. To consider a new model with unit support based on the PHN model is a promising alternative for more accurate inferential studies.
The rest of the paper is organized as follows. In Section 2, we introduce our proposal, the unit-power half-normal (UPHN) distribution. Several important properties of this new model are presented. In Section 3, inference is performed, and maximum likelihood (ML) estimators are obtained. In Section 4, the reparametrized model in terms of a quantile is presented. In Section 5, a simulation study is carried out to analyse the performance of ML estimators in finite samples for the proposed model without and with covariates. In Section 6, two real data applications are presented, and we are again dealing with cases without and with covariates. Finally, in Section 7, some concluding comments are presented.

2. Unit Power Half-Normal Distribution

In this section, we will discuss the stochastic representation of the UPHN model, including its pdf, cdf, and some properties of the model.

2.1. Stochastic Representation

Definition 1.
A random variable X follows a U P H N distribution with parameters σ > 0 and α > 0 if its stochastic representation is given by:
X = 1 1 + Y
where Y P H N ( σ , α ) . The notation X U P H N ( σ , α ) will be used.

2.2. Pdf, Cdf, Survival, and Hazard Functions

Proposition 1.
Let X U P H N ( σ , α ) . Then, the pdf of X is given by:
f X ( x ; σ , α ) = 2 α σ x 2 ϕ 1 x σ x 2 Φ 1 x σ x 1 α 1 , 0 < x < 1 , σ , α > 0 .
Proof. 
By considering the stochastic representation given in (4), x = g ( y ) = 1 / ( 1 + y ) , we have that
f X ( x ; σ , α ) = f Y ( g 1 ( x ) ; σ , α ) d g 1 ( x ) d x = 2 α σ ϕ g 1 ( x ) σ 2 Φ g 1 ( x ) σ 1 α 1 1 x 2 .
Taking into account that g 1 ( x ) = ( 1 x ) / x , (5) follows.    □
Proposition 2.
Let X U P H N ( σ , α ) . Then, its cdf is given by:
F ( x ; σ , α ) = 0 , x 0 , 1 2 Φ 1 x σ x 1 α , 0 < x < 1 , 1 , x 1 .
Proof. 
  • If x 0 , then f ( x ) = 0 . Therefore, F ( x , σ , α ) = 0 .
  • It follows from the fact that for 0 < x < 1 ,
    F ( x ; σ , α ) = 0 x 2 α σ t 2 ϕ 1 t σ t 2 Φ 1 t σ t 1 α 1 d t .
    Making u = 2 Φ 1 t σ t 1 , we obtain the result.
  • If x 1 , then:
    F ( x ; σ , α ) = x f ( t ) d t = 0 0 d t + 0 1 2 α σ t 2 ϕ 1 t σ t 2 Φ 1 t σ t 1 α 1 d t + 1 x 0 d t = 1 .
   □
Figure 1 shows the pdf and cdf of the UPHN ( σ , α ) model of σ = 1 with different values for α . Figure 2 shows the pdf and cdf of the UPHN ( σ , α ) model, for different values of ( σ , α ) .
Proposition 3.
Let X U P H N ( σ , α ) . Then, the survival and hazard function of X are given by:
S ( x ; σ , α ) = 2 Φ 1 x σ x 1 α a n d h ( x ; σ , α ) = 2 α σ x 2 ϕ 1 x σ x 2 Φ 1 x σ x 1 1 ,
for 0 < x < 1 .
Proof. 
Both functions are immediately obtained from their definitions, since S ( x ) = 1 F ( x ) and h ( x ) = f ( x ) 1 F ( x ) .    □
Figure 3 shows the survival and the hazard function for the UPHN ( σ , α ) model, considering several values of σ and α . We highlight that the hazard function is increasing.

2.3. Moments

Proposition 4.
Let X U P H N ( σ , α ) . Then, the n-th moment of X, where n is a positive integer, is given by:
μ X = α b n ( σ , α ) ,
where b n ( σ , α ) = 0 1 σ Φ 1 u + 1 2 + 1 n u α 1 d u , are computed numerically.
Proof. 
The moments are immediately obtained from their definition.
E [ X n ] = 2 α σ 0 1 x n 2 ϕ 1 x σ x 2 Φ 1 x σ x 1 α 1 d x .
Making the change of variable u = 2 Φ 1 x σ x 1 , we have that:
E [ X n ] = α 0 1 σ Φ 1 u + 1 2 + 1 n u α + 1 d u ,
and defining b n ( σ , α ) as the integral, the result is obtained.    □
Corollary 1.
Let X U P H N ( σ , α ) . Then, the skewness β 1 and kurtosis ( β 2 ) coefficients are
β 1 = b 3 3 b 1 b 2 + 2 b 1 3 ( b 2 b 1 2 ) 3 / 2 a n d β 2 = b 4 4 b 1 b 3 + 6 b 1 2 b 2 3 b 1 4 ( b 2 b 1 2 ) 2 ,
where b n = b n ( σ , α ) .
Proof. 
It is immediately obtained from the definition of these coefficients.    □
Table 1 and Table 2 provide us β 1 and ( β 2 ) for several values of parameters σ and α . We can observe that for low values of σ and α , the skewness coefficient is negative and for high values, positive skewness is obtained. As for the kurtosis coefficient, the observed effect is more dispersed.

2.4. Incomplete Moments

The k-th incomplete moments of the UPHN distribution are given by:
T k ( t ) = E X k | x < t = 2 α σ 0 t x k 2 ϕ 1 x σ x 2 Φ 1 x σ x 1 α 1 d x , k = 0 , 1 , 2 ,
This integral cannot be solved analytically. In particular, for k = 0 , 1 , 2 , we have:
T 0 ( t ) = F ( t ; σ , α ) , T 1 ( t ) = α γ 1 ( a ( t ) ) , T 2 ( t ) = α γ 2 ( a ( t ) ) ,
where γ j ( a ( t ) ) = a ( t ) 1 α Φ 1 u + 1 2 + 1 j u α 1 d u and a ( y ) = 2 Φ 1 y σ y 1 . γ j ( a ( t ) ) , are computed numerically.

2.5. The Lorenz Curve and the Gini Index

The standard definition of the Lorenz curve [17] is provided in terms of the first incomplete moment and the expected value of X. Specifically, for the UPHN model, the following closed form expression is obtained
L ( p ) = 1 μ X 0 t x f ( x ) d x = 1 μ X T 1 ( t ) ; = α μ X γ 1 ( a ( t ) ) .
The Gini index, also known as the Gini coefficient (see [18,19]), is a statistical dispersion metric associated with the Lorenz curve, intended to represent income inequality, wealth inequality, or consumption inequality within a nation or social group. The Gini index is defined as:
G ( σ , α ) = 1 1 μ X 0 1 T 0 ( t ) 2 d t .
Proposition 5.
Let X U P H N ( σ , α ) . Then, the Gini index is given by:
G ( σ , α ) = 1 4 α μ X 0 Φ 1 t σ t 2 α d t .
Proof. 
By definition, the proof is direct.    □
Proposition 6.
Let X U P H N ( σ , α ) ; the Renyi entropy of order δ for X is given by:
R δ ( x ) = 1 1 δ δ log ( α ) log ( ( α 1 ) δ + 1 ) + log 1 1 2 Φ ( 1 / σ ) ( α 1 ) δ + 1 ,
where δ > 0 and δ 1 .
Proof. 
The Renyi entropy is defined as R δ ( x ) = 1 1 δ log 0 f ( x ) δ d x . Therefore, we have:
R δ ( y ) = 1 1 δ log 0 2 α σ x 2 ϕ 1 x σ x 2 Φ 1 x σ x 1 α 1 δ d x ,
Making the change of variable u = 2 Φ 1 x σ x 1 , the result is obtained.    □

2.6. Shannon Entropy

The entropy of a random variable X is a measure of its uncertainty. The Shannon entropy measure is defined by:
J S = E log ( f ( x ) ) .
If follows, after extensive algebraic manipulations, that the Shannon entropy for the UPHN model is:
J S = log σ 2 π 2 α + 2 m 1 + 1 2 σ 2 α b 2 ( σ , α ) 2 α b 1 ( σ , α ) + 1 + ( α + 1 ) m 2 .
where m 1 = E ( log ( x ) ) , m 2 = E ( log 2 Φ 1 x σ x 1 ) , and with b 1 ( σ , α ) and b 2 ( σ , α ) as given above.

2.7. Quantiles

Proposition 7.
Let X U P H N ( σ , α ) . Then, the quantile function is given by:
Q ( p ; σ , α ) = σ Φ 1 ( 1 p ) 1 / α 2 + 1 2 + 1 1 , 0 < p < 1 .
Proof. 
It follows from a direct computation, by applying the definition of quantile function.    □
Corollary 2.
The quartiles for the U P H N distribution are:
1. 
(First quartile) Q ( 0.25 ; σ , α ) = σ Φ 1 0 . 75 1 / α 2 + 1 2 + 1 1 .
2. 
(Median) Q ( 0.5 ; σ , α ) = σ Φ 1 0 . 5 1 / α 2 + 1 2 + 1 1 .
3. 
(Third quartile) Q ( 0.75 ; σ , α ) = σ Φ 1 0 . 25 1 / α 2 + 1 2 + 1 1 .
The quantile function can also be used to study the kurtosis and skewness in the UPHN model. A classical measurement for skewness was introduced by MacGillivray [20], which is given by:
ρ X i ( p ) = Q X i ( 1 p ; σ , α ) + Q X i ( p ; σ , α ) 2 Q X i ( 0.5 ; σ , α ) Q X i ( 1 p ; σ , α ) Q X i ( p ; σ , α ) , p ( 0 , 1 ) .
In particular, the MacGuillevray skewness measurement can efficiently describe the effect of the parameters ( σ , α ) on asymmetry. In Figure 4, the behavior of the asymmetry coefficient is plotted with respect to the values of the parameters.
The kurtosis of the UPHN ( σ , α ) distribution can also be studied using the Moors kurtosis coefficient [21], usually given by:
K ( σ , α ) = Q ( 7 / 8 ; σ , α ) Q ( 5 / 8 ; σ , α ) + Q ( 3 / 8 ; σ , α ) Q ( 1 / 8 ; σ , α ) Q ( 3 / 4 ; σ , α ) Q ( 1 / 4 ; σ , α ) .
It can be seen in [21] that for large values of (11), the distribution has heavy tails and for small values, the model has lighter tails. Figure 5 shows the behavior of the Moors kurtosis coefficient for the UPHN distribution.

2.8. Order Statistics

Order statistics have a wide range of applications in physical and life sciences (see, for instance, Balakrishnan and Cohen [22]). From a statistical perspective, they allow the computation of useful functions such as the sample range and the sample median. The following result states the pdf of k-th order statistic from a UPHN random sample of size n, which is arranged in a non-decreasing order.
Proposition 8.
Let X 1 , X 2 , , X n be independent and identically distributed U P H N ( σ , α ) random variables. Then, the pdf of the k-th order statistic, X k : n , is given by:
f X k : n ( x ) = 2 α n ! ( k 1 ) ! ( n k ) ! σ x 2 ϕ 1 x σ x 2 Φ 1 x σ x 1 α ( n k + 1 ) 1 1 2 Φ 1 x σ x 1 α k 1 .
Corollary 3.
Let X 1 , X 2 , , X n be a random sample from a U P H N ( σ , α ) distribution. Then,
1. 
The pdf of the minimum, X 1 : n , is
f X 1 : n ( x ) = 2 α n σ x 2 ϕ 1 x σ x 2 Φ 1 x σ x 1 n α 1 ,
2. 
The pdf of the maximum, X n : n , is
f X n : n ( x ) = 2 α n σ x 2 ϕ 1 x σ x 2 Φ 1 x σ x 1 α 1 1 2 Φ 1 x σ x 1 α n 1 .

2.9. Bonferroni Curve

In various disciplines, such as socioeconomics and public health sciences, there is a need to compare and analyze inequality in distributions. Bonferroni curves are tools useful to reach this aim. It is worth mentioning that these curves have applications not only in economics to study income and poverty, but also in medicine or reliability. A detailed discussion can be seen in Bonferroni [23] or Arcagni and Porro [24]. The expressions of these curves for the UPHN model are presented below.
Proposition 9.
Let X U P H N ( σ , α ) . Then, for 0 < p < 1 , the Bonferroni curve, say B ( p ) , is given by:
B ( p ) = α μ X p j = 0 ( 1 ) j σ j γ a ( q ) j ,
where μ X = E ( X ) .
Proof. 
The expression above is obtained by using the definition of the Bonferroni curve, that is,
B ( p ) = 1 μ X p 0 p F 1 ( t ) d t = 1 μ X p 0 q x f ( x ) d x , 0 < p < 1 ,
where μ X is the expected value of the corresponding non-negative random variable and q = F 1 ( p ) .    □

3. Inference

In this section, the inference about the parameters in the UPHN distribution is carried out from a classical point of view. Let us consider X 1 , X 2 , , X n a random sample of X U P H N ( σ , α ) . The maximum likelihood (ML) estimation method is discussed below. Given a random sample x 1 , x 2 , , x n of size n from U P H N ( σ , α ) , the log-likelihood function is given by:
( θ ) = n log 2 α σ n 2 log ( 2 π ) 2 i = 1 n log ( x i ) 1 2 i = 1 n v i 2 + ( α 1 ) i = 1 n log 2 Φ v i 1 ,
where v i = 1 x i σ x i . Therefore, the score assumes the form S ( θ ) = ( S σ ( θ ) , S α ( θ ) ) where:
S σ ( θ ) = n σ 2 + 1 σ i = 1 n v i 2 2 ( α 1 ) σ i = 1 n ϕ v i v i 2 Φ v i 1 ,
S α ( θ ) = n α + i = 1 n log 2 Φ v i 1 .
The ML estimators (MLEs) can be obtained by solving the likelihood Equations (13) and (14). Numerical methods, such as the Newton–Raphson procedure, must be used. Alternative maximization techniques could also be applied, for instance, those proposed in MacDonald [25].
The asymptotic variance of the MLEs, say θ ^ = ( σ ^ , α ^ ) , can be estimated by the Fisher information matrix defined as I ( θ ) = E 2 ( θ ) / θ θ , where ( θ ) is the log-likelihood function of the UPHN model given in (12). Recall that, under the regularity conditions,
I ( θ ) 1 / 2 θ ^ θ D N 2 ( 0 2 , I 2 ) , a s n + ,
where D stands for convergence in distribution and N 2 ( 0 2 , I 2 ) denotes the standard bivariate normal distribution. The elements of the matrix 2 ( θ ) / θ θ are given by I σ σ = 2 ( θ ) / σ 2 , I σ α = 2 ( θ ) / σ α and I α α = 2 ( θ ) / α 2 . Explicitly, we have:
I σ σ = n σ 2 + 3 σ 2 i = 1 n v i 2 4 ( α 1 ) σ 2 i = 1 n ϕ ( v i ) v i 2 Φ ( v i ) 1 + 2 ( α 1 ) σ 2 i = 1 n v i 2 ϕ ( v i ) v i ( 2 Φ ( v i ) 1 ) + 2 ϕ ( v i ) ( 2 ϕ ( v i ) 1 ) 2 , I σ α = 2 σ 2 i = 1 n v i ϕ v i 2 Φ v i 1 , I α α = n α 2 .
In practice, it is not possible to obtain in a closed form the expected value of previous expressions. However, the covariance matrix of the MLEs, I ( θ ) 1 , can be estimated consistently by I ( θ ^ ) 1 , where I ( θ ^ ) denotes the observed information matrix, which is obtained as
I ( θ ^ ) = 2 ( θ ) θ θ | θ = θ ^ .
The asymptotic variances of σ ^ and α ^ are estimated by the diagonal elements of I ( θ ^ ) 1 , and their standard errors by the square root of asymptotic variances.

4. Quantile Regression

We now introduce a regression study, in which a quantile regression model is proposed to describe the conditional quantile of the response variable. Given the quantile function of the U P H N distribution, the pdf of the UPHN distribution can be reparameterized in terms of its p t h quantile, denoted as ρ = Q ( p ; σ , α ) .
Let σ = 1 ρ ρ Φ 1 1 2 ( 1 p ) 1 / α + 1 . Then, the reparameterized pdf for the UPHN model is:
f X ( x ; ρ , α , p ) = 2 α ( 1 α ) x 2 h ( α , p ) ϕ ( 1 x ) ρ ( 1 ρ ) x h ( α , p ) 2 Φ ( 1 x ) ρ ( 1 ρ ) x h ( α , p ) 1 α 1 ,
where h ( α , p ) = Φ 1 1 2 ( 1 p ) 1 / α + 1 , 0 < x < 1 , 0 < p < 1 , 0 < ρ < 1 , and α > 0 .
The cdf of the reparameterized UPHN distribution is:
F ( x ; ρ , α , p ) = 1 2 Φ ( 1 x ) ρ Φ 1 1 2 [ ( 1 p ) 1 / α + 1 ] ( 1 ρ ) x 1 α .

4.1. Model Formulation

Let us suppose that 0 < p < 1 is fixed and we have a random sample x i ( p ) U P H N ( ρ i , α ) , i = 1 , 2 , , n . Here, the parameters ρ i and α are unknown. As a result, the suggested quantile regression model can be defined as:
g ( ρ i ( p ) ) = z i β ( p ) ,
where β ( p ) = ( β 0 ( p ) , β 1 ( p ) , , β k ( p ) ) and z i = ( 1 , z i 1 , , z i k ) act as the regression coefficient vector and the ith vector of covariate values, respectively. It is important to note that g is the link function that associates the covariates with the response variable. We can obtain a median quantile regression for a unit response variable when p = 0.5 . Taking into account that the UPHN distribution is a probability model with support on ( 0 , 1 ) , we use and define the logit-link function as:
g ( ρ i ( p ) ) = log ρ i ( p ) 1 ρ i ( p ) , i = 1 , , n .

4.2. ML Method for the Regression Coefficients

Now, our objective is to estimate the unknown parameters of the UPHN quantile regression model using the ML technique. For this purpose, the logit link function is considered:
g ( ρ i ( p ) ) = log ρ i ( p ) 1 ρ i ( p ) = z i β .
From the above formula, the inverse link function is:
ρ i ( p ) = exp ( z i β ( p ) ) 1 + exp ( z i β ( p ) ) .
Let us now have a random sample X 1 , , X n with X i UPNH ( ρ i ( p ) , α ) , and observed values x 1 , , x n . Then, using the pdf introduced in (15), the associated log-likelihood function is:
l ( Θ ( p ) ) = n log ( 2 ) + log ( α ) + n log [ h ( α , p ) ] + i = 1 n log ϕ [ 1 x i ] ρ i [ 1 ρ i ] x i h ( α , p ) + i = 1 n log ρ i [ 1 ρ i ] x i 2 + ( α 1 ) i = 1 n log 2 Φ [ 1 x i ] ρ i [ 1 ρ i ] x i h ( α , p ) 1 ,
where Θ = ( α , β ) denotes the unknown vector of parameters. The ML estimators of Θ , denoted as Θ ^ = ( α ^ , β ^ ) is achieved by maximizing l ( Θ ) with respect to Θ . The optim function of the R software, 4.3.1 version, [26] can be used to maximize l ( Θ ) .

4.3. Model Checking

Once we have fitted the model, it is crucial to assess whether the regression model is appropriate for the data. In this sense, the analysis of residuals is key to check or validate the fitted model. We focus on the Cox–Snell residuals (Cox and Snell 1968) [27], which are calculated for the i t h observation as:
η ^ i = log 1 F ( x i ; ρ i ( p ) , α ) , i = 1 , , n .
Here, F ( x i , ρ i ( p ) , α ) stands for the estimated cdf of the UPHN distribution reparameterized in terms of quantiles. If the fitted model provides a good fit for the data, the residual η ^ i will be an observation of a random variable from an exponential distribution with a scale parameter of one.

5. Simulation

In this section, a simulation study is carried out to evaluate the performance of the ML estimators. We used R software 4.3.1 version, for our calculations, developing customized code that integrates the optim function from R package stats [28].

5.1. Simulation 1

Next, an algorithm to generate samples from the UPHN ( σ , α ) is provided. Algorithm 1 is based on the inverse transformation method, where the inverse of the cdf is used.
Algorithm 1 Simulating values from the UPHN ( σ , α ) distribution.
  • 1: Simulate U U ( 0 , 1 ) .
  • 2: Calculate X = σ Φ 1 ( 1 U ) 1 / α 2 + 1 2 + 1 1
As parameter values in our simulation, we consider σ 1 , 2 and α 0.8 , 1 , 2 , 3 , 5 . For the sample size, we consider n 150 , 300 , 600 . For each sample size and every combination of σ , α , we perform 1000 repetitions and the corresponding ML estimates are calculated. The results are given in Table 3. As summaries, we provide the estimated bias (bias) for the ML estimators of σ and α , standard errors (SE), the root of the estimated mean square error (RMSE), and the empirical coverage probability (CP) for the asymptotic intervals based on MLEs to 95%. For the ML estimators of σ and α , it should be noted that, as the sample size increases, the bias, SE, and RMSE decrease. It should also be noted that as the sample size increases, the SE and RMSE are closer, which suggests that the standard errors of the ML estimators are well estimated. As for the CPs, we highlight that they are close to the nominal level 0.95.

5.2. Simulation 2

In this section, data are generated from the quantile regression model with sample sizes n { 100 , 200 , 300 , 400 } and quantile levels p { 0.10 , 0.25 , 0.50 , 0.75 , 0.90 } under two different scenarios. In the first scenario, three covariates are considered, that is,
l o g i t ( ρ i ( p ) ) = β 0 ( p ) + β 1 ( p ) z i 1 + β 2 ( p ) z i 2 + β 3 ( p ) z i 3 , i = 1 , 2 , , n ,
where β 0 ( p ) = 0.275 , β 1 ( p ) = 0.098 , β 2 ( p ) = 0.015 , β 3 ( p ) = 0.397 , z i 1 , z i 2 , z i 3 iid N ( 0 , 1 ) , and α = 1.285 , whereas in the second scenario, five covariates are proposed
l o g i t ( ρ i ( p ) ) = β 0 ( p ) + β 1 ( p ) z i 1 + β 2 ( p ) z i 2 + β 3 ( p ) z i 3 + β 4 ( p ) z i 4 + β 5 ( p ) z i 5 ,
where β 0 ( p ) = 2.975 , β 1 ( p ) = 0.005 , β 2 ( p ) = 0.092 , β 3 ( p ) = 0.958 , β 4 ( p ) = 0.122 , β 5 ( p ) = 0.289 , z i 1 , z i 2 , z i 3 , z i 4 , z i 5 iid N ( 0 , 1 ) , and α = 7.425 .
For fixed values of n and p, the response variable is generated as
x i = σ i Φ 1 ( 1 u i ) 1 / α 2 + 1 2 + 1 1 ,
where u i is an observation generated from a continuous uniform distribution on ( 0 , 1 ) , that is, u i U ( 0 , 1 ) and σ i is calculated as
σ i = 1 ρ i ρ i Φ 1 1 2 [ ( 1 p ) 1 / α + 1 ] .
For both scenarios, and fixed ( n , p , α ) , M = 500 Monte Carlo replicates were performed. Then, the empirical bias, RMSE, and CP are calculated by
Bias ( γ ^ ) = 1 M i = 1 M ( γ ^ i γ ) , RMSE ( γ ^ ) = 1 M i = 1 M ( γ ^ i γ ) 2 1 / 2 , CP 95 % = 1 M i = 1 M 1 [ γ ^ i ± 1.96 × SE ( γ ^ i ) ] ( γ ) ,
respectively, where γ = α , β j ( p ) , j = 1 , , 5 , 1 [ γ ^ i ± 1.96 × SE ( γ ^ i ) ] ( γ ) is the indicator function defined by
1 [ γ ^ i ± 1.96 × SE ( γ ^ i ) ] ( γ ) = 1 , γ ( γ ^ i ± 1.96 × SE ( γ ^ i ) ) 0 , γ ( γ ^ i ± 1.96 × SE ( γ ^ i ) ) ,
and SE ( γ ^ i ) is the estimated standard error of γ ^ i .
These summaries are provided in Table 4 and Table 5. In both settings, the good performance of our estimators is observed.

6. Applications

In this section, the unit-power half-normal quantile regression model is applied to two real-world health data sets. The first one involves patient health outcomes and the second one examines survival rates under a specific treatment. We compare the results with those obtained from the Kumaraswamy (kum), the unit-Birnbaum–Saunders (ubs) [9], the unit-generalized half-normal-X (ughnx) [13], the unit-Gompertz (ugompertz) [29], and the unit-Weibull (uweibull) [30] quantile regression models.
We performed our calculations by using the R software. To fit the unit-power half-normal quantile regression model, we implemented our own code, combining existing functions available in R packages. For the Kumaraswamy, the unit-Birnbaum–Saunders, the unit-generalized half-normal-X, the unit-Gompertz, and the unit-Weibull quantile regression models, we used the R package unitquantreg [31]. The codes are available in a public GitHub repository (https://github.com/DarlinSoto/UPHN (accessed on 12 August 2024)).

6.1. Application 1: Body Fat Data Set

This data set was reported and studied in [32] and consists of 298 observations about the body fat percentage of individuals treated at a public hospital in Curitiba, Paraná, Brazil. The data set includes the variables: fat percentage of arms, legs, body, android, gynecoid, the body mass index ( k g / m 2 ), age (years), gender (female or male), and level of physical activity (sedentary, insufficiently active, or active). The goal is to explore the functional relationship between the covariates gender, age, body mass index, and the level of physical activity, with the response variable body fat in arms. The assumed regression model for ρ i is
l o g i t ( ρ i ( p ) ) = β 0 ( p ) + β 1 ( p ) Age i + β 2 ( p ) BodyMassIndex i + β 3 ( p ) Gender_M i + β 4 ( p ) IPAQ_suffactive + β 5 ( p ) IPAQ_active , i = 1 , , 298 ,
where
Gender_M i = 1 , if the person i is a male 0 , if the person i is a female
IPAQ_suffactive i = 1 , if the physical activity of the person i is insufficiently active 0 , if the physical activity of the person i is sedentary or active
IPAQ_active i = 1 , if the physical activity of the person i is active 0 , if the physical activity of the person i is sedentary or insufficiently active .
To explore the relationship between gender, age, body mass index, and the level of physical activity with the body fat in arms, the unit-power half-normal, the Kumaraswamy, the unit-generalized half-normal-X, the unit-Gompertz, and the unit-Weibull quantile regression models were fitted to the data and the fits were compared by using goodness-of-fit measurements. Table 6 presents the maximum likelihood estimates and their standard errors, Table 7 shows their p-values, and Table 8 shows the negative value of the log-likelihood, Akaike information criterion, Bayesian information criterion, and Hannan and Quinn information criterion given by
Neg 2 LogLike = 2 log ( L )
AIC = 2 log ( L ) + 2 k
BIC = 2 log ( L ) + k log ( n )
HQIC = 2 log ( L ) + 2 k log [ l o g ( n ) ]
where n is the sample size, k is the number of parameters, and L is the likelihood function evaluated in the estimated parameters.
According to the results presented in Table 6 and Table 7, and for all quantile levels p, the unit-power half-normal, the unit-generalized half-normal-X, the unit-Gompertz, and the unit-Weibull quantile regression models show that all covariates are statistically significant at the 0.05 level. On the other hand, the Kumaraswamy quantile regression model shows that all variables are significant except for the IPAQ_suffactive variable. Upon analyzing the signs of β ^ j , we observe that an increase in covariates such as Age and BodyMassIndex corresponds to an increase in arm fat percentage. Moreover, males exhibit lower arm fat percentages than females, and individuals with a sedentary lifestyle demonstrate an increase in arm fat percentages compared to those with insufficiently active or active lifestyles.
Based on the results given in Table 8, the unit-power half-normal quantile regression model consistently shows higher log-likelihood values across all quantile levels p compared to the Kumaraswamy, the unit-generalized half-normal-X, the unit-Gompertz, and the unit-Weibull quantile regression models. Consequently, all goodness-of-fit measurements indicate that the unit-power half-normal quantile regression model provides the best fit to these data, while the unit-Gompertz model provides the worst one.
To validate these inferences, we conducted a diagnostic analysis for the fitted models. The Cox–Snell quantile residuals are considered, QQ plots are carried out by using their simulated envelopes. The Cramér–von Mises test is applied, where the null hypothesis is that the residuals come from an exponential distribution with a scale parameter of one. These plots are given in Figure 6. This figure suggests that the unit-power half-normal and Kumaraswamy quantile regression models provide a more substantial agreement with the data set compared to the unit-generalized half-normal-X, unit-Gompertz, and unit-Weibull quantile regression models, mainly because most of the residuals envelopes fail to cover the grey points, and the null hypothesis that the residuals come from an exponential distribution is rejected at the 5% significance level.
Therefore, the unit-power half-normal quantile regression model provides the best fit for the data, as it consistently shows the minimum goodness-of-fit measures, with all coefficients being statistically significant. The model’s assumptions are validated through diagnostic analyses, which confirm its reliability to describe the relationships between the covariates and arm fat percentage.

6.2. Application 2: Autologous Stem Cell Transplants Data Set

We consider another application linked to autologous peripheral blood stem cell (PBSC) transplants, which have been extensively utilized to accelerate hematologic recovery after myeloablative therapy for diverse malignant hematological disorders. The data set comprises a study involving 239 patients who agreed to undergo autologous PBSC transplant following myeloablative chemotherapy doses between 2003 and 2008 at the Edmonton Hematopoietic Stem Cell Lab within the Cross Cancer Institute of Alberta Health Services. This data set has been studied in [33] and contains information about the patients’ age, gender, and clinical characteristics. We aim to explain the response variable recovery rates for viable CD34+ cells (rcd) in terms of age, gender, and chemotherapy. The regression model assumed for ρ i is given by
l o g i t ( ρ i ( p ) ) = β 0 ( p ) + β 1 ( p ) Sexmale_M i + β 2 ( p ) Age i + β 3 ( p ) Chemo i , i = 1 , , 239 ,
where
Sexmale_M i = 1 , if the gender of the patient i is male 0 , if the gender of the patient i is female
Chemo i = 1 , if a patient i receives a chemotherapy on a 3 - day protocol 0 , if a patient i receives a chemotherapy on a 1 - day protocol .
We examined the relationship between gender, age, chemotherapy, and the recovery rate for viable CD34+ cells by fitting the unit-power half-normal, the Kumaraswamy, the unit-generalized half-normal-X, the unit-Gompertz, and the unit-Birnbaum–Saunders quantile regression models, and evaluated the fits using goodness-of-fit measures.
Table 9 and Table 10 present the maximum likelihood estimates, their standard errors, and their p-values. For the unit-power half-normal and the unit-Gompertz quantile regression models, all quantile levels show statistically significant regression coefficients for the variables age and chemo at the 5% level. On the other hand, for the unit-generalized half-normal-X and the unit-Birnbaum–Saunders quantile regression models, the variable chemo is the only one statistically significant at the 5% level, while for the Kumaraswamy quantile regression, none of the three variables is statistically significant at the 5% level. Thus, based on the unit-power half-normal and the unit-Gompertz quantile regressions, we may conclude that the covariates Age and Chemo significantly affect the response variable recovery rate for viable CD34+ cells. Specifically, as Age increases, the recovery rate for viable CD34+ cells also increases. Additionally, patients receiving chemotherapy on a three-day protocol exhibit a higher recovery rate for viable CD34+ cells than those on a one-day protocol.
Table 11 displays the model selection criteria for the unit-power half-normal, the Kumaraswamy, the unit-generalized half-normal-X, the unit-Gompertz, and the unit-Birnbaum–Saunders quantile regressions. From this table, it can be said that the log-likelihood values of the unit-power half-normal quantile regression surpass those of the Kumaraswamy, the unit-generalized half-normal-X, the unit-Gompertz, and the unit-Birnbaum–Saunders quantile regressions. Therefore, based on these criteria, we may conclude that the unit-power half-normal model has a superior performance.
Finally, Figure 7 displays the QQ plots of the Cox–Snell quantile residuals and the p-value of the Cramér–von Mises test for the null hypothesis that the residuals come from an exponential distribution with a scale parameter of one. This figure indicates that most of the Cox–Snell quantile residuals obtained from the unit-power half-normal and the Kumaraswamy quantile regressions are inside the simulated envelopes. However, this is not the case for the Cox–Snell quantile residuals obtained from the unit-generalized half-normal-X, the unit-Gompertz, and the unit-Birnbaum–Saunders quantile regressions. These facts suggest that the assumption about Cox–Snell quantile residuals may not be right in those cases. This is validated by the p-value of the Cramér–von Mises test, where for the unit-generalized half-normal, the unit-Gompertz, and the unit-Birnbaum–Saunders quantile regressions, the p-value is less than 0.05, rejecting the null hypothesis that residuals come from an exponential distribution.
Therefore, the unit-power half-normal quantile regression model is the best fit for the autologous stem cell transplants data. It consistently shows higher log-likelihood values and better model selection criteria compared to the other models. The significant coefficients for age and chemotherapy at the 5% level highlight their strong association with recovery rates for viable CD34+ cells. Furthermore, diagnostic analyses, including QQ plots and the Cramér–von Mises test, validate the model’s assumptions, confirming its superior performance.

7. Conclusions

In this paper, we have put forth a new quantile regression model. The predictor covariates were linked to the quantile of the dependent variable through the logit link function. We examined various aspects such as estimation, inference relying on the maximum likelihood method, and diagnoses via the quantile residual. We evaluated two real data sets drawn from medical data using our proposed quantile regression model and the previously introduced Kumaraswamy quantile regression models. The outcomes showed that our proposed quantile regression fitted well to the two real-world data sets. In the first one, the covariates age, body mass index, gender, and level of physical activity are statistically significant to evaluate the level of physical activity with the body fat in arms. For the second data set, the covariates age and chemo are statistically significant to explain the response variable recovery rates for viable CD34+ cells (rcd), using our proposed quantile regression model as an alternative to other rival models outlined in the related literature. We carried out a comparison based on different information criteria and the model’s goodness of fit, showing its suitability.

Author Contributions

Conceptualization, K.I.S. and Y.M.G.; software, D.S.; formal analysis, K.I.S., Y.M.G. and I.B.-C.; investigation, K.I.S., Y.M.G. and I.B.-C.; writing—original draft preparation, K.I.S., Y.M.G., D.S. and I.B.-C.; writing—review and editing, K.I.S., Y.M.G., D.S., I.B.-C. and D.S.; supervision, Y.M.G. and I.B.-C.; funding acquisition, Y.M.G. and I.B.-C. All authors have read and agreed to the published version of the manuscript.

Funding

The research of I. Barranco-Chamorro was supported by the IOAP of the University of Seville, Spain.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data sets are available in the references given in Section 6.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Jones, M.C. Kumaraswamy’s Distribution: A Beta-Type Distribution with Some Tractability Advantages. Stat. Methodol. 2009, 6, 70–81. [Google Scholar] [CrossRef]
  2. Gómez-Déniz, E.; Sordo, M.A.; Calderín-Ojeda, E. The Log-Lindley distribution as an alternative to the beta regression model with applications in insurance. Insur. Math. Econ. 2014, 54, 49–57. [Google Scholar] [CrossRef]
  3. Mazucheli, J.; Bapat, S.R.; Menezes, A.F.B. A new one-parameter unit Lindley distribution. Chil. J. Stat. 2019, 1, 53–67. [Google Scholar]
  4. Abd El-Monsef, M.M.E.; El-Awady, M.M.; Seyam, M.M. A new quantile regression model for modelling child mortality. Int. J. Biomath. 2022, 10, 142–149. [Google Scholar]
  5. Altun, E.; Cordeiro, G.M. The unit-improved second-degree Lindley distribution: Inference and regression modelling. Comput. Stat. 2020, 35, 259–279. [Google Scholar] [CrossRef]
  6. Ahmad, H.H.; Almetwally, E.M.; Elgarhy, M.; Ramadan, D.A. On Unit Exponential Pareto Distribution for Modeling the Recovery Rate of COVID-19. Processes 2023, 11, 232. [Google Scholar] [CrossRef]
  7. Grassia, A. On a Family of Distributions with Argument between 0 and 1 Obtained by Transformation of the Gamma Distribution and Derived Compound Distributions. J. Stat. 1977, 19, 108–114. [Google Scholar] [CrossRef]
  8. Ghitany, M.; Mazucheli, J.; Menezes, A.; Alqallaf, F. The unit-inverse Gaussian distribution: A new alternative to two-parameter distributions on the unit interval. Commun. Stat. Theory Methods 2018, 48, 3423–3438. [Google Scholar] [CrossRef]
  9. Mazucheli, J.; Menezes, A.F.B.; Dey, S. The unit-Birnbaum-Saunders distribution with applications. Chil. J. Stat. 1977, 1, 47–57. [Google Scholar]
  10. Modi, K.; Gill, V. Unit Burr III distribution with application. J. Stat. Manag. Syst. 2019, 23, 579–592. [Google Scholar] [CrossRef]
  11. Korkmaz, M.; Chesneau, C. On the unit Burr-XII distribution with the quantile regression modeling and applications. Comput. Appl. Math. 2021, 40, 29. [Google Scholar] [CrossRef]
  12. Haq, M.; Hashmi, S.; Aidi, K.; Ramos, P.F.L. Unit Modified Burr-III Distribution: Estimation, Characterizations and Validation Test. Ann. Data Sci. 2023, 10, 415–449. [Google Scholar] [CrossRef]
  13. Bakouch, H.S.; Nikb, A.S.; Asgharzadehb, A.; Salinas, H.S. A flexible probability model for proportion data: Unit-half-normal distribution. Commun. Stat. Case Stud. Data Anal. 2021, 7, 271–288. [Google Scholar] [CrossRef]
  14. Gómez, Y.M.; Bolfarine, H. Likelihood-based inference for the power half-normal distribution. J. Stat. Theory Appl. 2015, 14, 383–398. [Google Scholar] [CrossRef]
  15. Barrios, L.; Gómez, Y.M.; Venegas, O.; Barranco-Chamorro, I.; Gómez, H.W. The Slashed Power Half-Normal Distribution with Applications. Mathematics 2022, 10, 1528. [Google Scholar] [CrossRef]
  16. Pallini, A. The discrete Powe Half-Normal Distribution. Statistica 2022, 82, 229–242. [Google Scholar]
  17. Lorenz, M.O. Methods of measuring the concentration of wealth. J. Am. Stat. Assoc. 1905, 9, 209–219. [Google Scholar] [CrossRef]
  18. Gini, C. On the measurement of concentration and variability of characters. Metron 2005, 1, 1–38. [Google Scholar]
  19. Gini, C. Measurement of inequality of incomes. Econ. J. 1921, 31, 124–126. [Google Scholar] [CrossRef]
  20. MacGillivray, H.L. Skewness and Asymmetry: Measures and Orderings. Ann. Stat. 1986, 14, 994–1011. [Google Scholar] [CrossRef]
  21. Moors, J.J. A quantile alternative for kurtosis. J. R. Stat. Soc. Ser. D 1988, 37, 25–32. [Google Scholar] [CrossRef]
  22. Balakrishnan, N.; Cohen, C.A. Order Statistics and Inference: Estimation Methods; Statistical Modeling and Decision Science; Elsevier Science: Amsterdam, The Netherlands, 1991. [Google Scholar]
  23. Bonferroni, C.E. Elementi di Statistica Generale; Libreria Seber: Firenze, Italy, 1930. [Google Scholar]
  24. Arcagni, A.; Porro, F. The Graphical Representation of Inequality. Rev. Colomb. Estad. 2014, 37, 419–436. [Google Scholar] [CrossRef]
  25. MacDonald, I.L. Does Newton-Raphson really fail? Stat. Methods Med Res. 2014, 23, 308–311. [Google Scholar] [CrossRef] [PubMed]
  26. R Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2024; Available online: https://www.R-project.org/ (accessed on 31 July 2024).
  27. Cox, D.R.; Snell, E.J. A general definition of residuals. J. Royal Stat. Soc. Ser. B 1968, 30, 248–265. [Google Scholar] [CrossRef]
  28. Bolar, K. STAT: Interactive Document for Working with Basic Statistical Analysis. 2019. Available online: https://cran.r-project.org/web/packages/STAT/index.html (accessed on 1 April 2024).
  29. Mazucheli, J.; Menezes, A.F.; Dey, S. Unit-Gompertz Distribution with Applications. Statistica 2019, 79, 25–43. [Google Scholar]
  30. Mazucheli, J.; Menezes, A.F.B.; Alqallaf, F.; Ghitany, M.E. Bias-Corrected Maximum Likelihood Estimators of the Parameters of the Unit-Weibull Distribution. Austrian J. Stat. 2021, 50, 41–53. [Google Scholar]
  31. Menezes, A.F.B. unitquantreg: Parametric Quantile Regression Models for Bounded Data. 2023. Available online: https://cran.r-project.org/web/packages/unitquantreg/index.html (accessed on 15 April 2024).
  32. Petterle, R.; Bonat, W.; Scarpin, C.; Jonasson, T.; Borba, V. Multivariate quasi-beta regression models for continuous bounded data. Int. J. Biostat. 2021, 17, 39–53. [Google Scholar] [CrossRef]
  33. Yang, H.; Acker, J.P.; Cabuhat, M.; Letcher, B.; Larratt, L.; McGann, L.E. Association of post-thaw viable CD34+ cells and CFU-GM with time to hematopoietic engraftment. Bone Marrow Transplant. 2005, 35, 881–887. [Google Scholar] [CrossRef]
Figure 1. Pdf and cdf for UPHN ( σ = 1 , α ) model with different values of α .
Figure 1. Pdf and cdf for UPHN ( σ = 1 , α ) model with different values of α .
Axioms 13 00599 g001
Figure 2. Pdf and cdf for UPHN ( σ , α ) model for some values of parameters σ and α .
Figure 2. Pdf and cdf for UPHN ( σ , α ) model for some values of parameters σ and α .
Axioms 13 00599 g002
Figure 3. Survival and hazard functions for U P H N ( σ , α ) model for some values of parameters σ and α .
Figure 3. Survival and hazard functions for U P H N ( σ , α ) model for some values of parameters σ and α .
Axioms 13 00599 g003
Figure 4. Plots of the MacGillivray skewness coefficient in UPHN model.
Figure 4. Plots of the MacGillivray skewness coefficient in UPHN model.
Axioms 13 00599 g004
Figure 5. Plots of the Moors kurtosis coefficient for the U P H N model.
Figure 5. Plots of the Moors kurtosis coefficient for the U P H N model.
Axioms 13 00599 g005
Figure 6. QQ plots with envelopes of Cox–Snell residuals for the indicated 100 p-th quantile level and Cramér-von Mises test p-value for residuals obtained from different distributions (uphn, kum, ughnx, ugompertz, and uweibull) with body fat data.
Figure 6. QQ plots with envelopes of Cox–Snell residuals for the indicated 100 p-th quantile level and Cramér-von Mises test p-value for residuals obtained from different distributions (uphn, kum, ughnx, ugompertz, and uweibull) with body fat data.
Axioms 13 00599 g006
Figure 7. QQ plots with envelopes of Cox–Snell residuals for the proposed 100 p-th quantile level and Cramér–von Mises test p-value for residuals obtained from different distributions (uphn, kum, ughnx, ugompertz, and uweibull) with PBSC data.
Figure 7. QQ plots with envelopes of Cox–Snell residuals for the proposed 100 p-th quantile level and Cramér–von Mises test p-value for residuals obtained from different distributions (uphn, kum, ughnx, ugompertz, and uweibull) with PBSC data.
Axioms 13 00599 g007
Table 1. Skewness coefficient of U P H N ( σ , α ) model for different values of σ and α .
Table 1. Skewness coefficient of U P H N ( σ , α ) model for different values of σ and α .
α
σ 0.30.71.535
0.1−1.7390−0.9557−0.5196−0.3301−0.2892
0.6−1.1043−0.31910.17640.38170.3805
1.4−0.70970.12180.71150.93000.8517
2.5−0.43810.45971.16691.39151.2040
5−0.12280.90351.84452.05451.6279
100.18051.40332.74192.86752.0197
Table 2. Kurtosis coefficient of U P H N ( σ , α ) model for different values of σ and α .
Table 2. Kurtosis coefficient of U P H N ( σ , α ) model for different values of σ and α .
α
σ 0.30.71.535
0.15.83563.42392.85392.87062.9599
0.63.18962.07882.32582.84543.0395
1.42.18771.85422.91894.00654.1003
2.51.74922.03034.07745.83785.5157
51.47772.69566.929810.05058.0999
101.44354.017912.864818.333411.6586
Table 3. Estimated bias, SE, and RMSE for ML estimators in finite samples from the U P H N model.
Table 3. Estimated bias, SE, and RMSE for ML estimators in finite samples from the U P H N model.
True Value n = 150 n = 300 n = 600
σ α EstimatorBiasseRMSECPBiasseRMSECPBiasseRMSECP
10.8 σ −0.00750.07340.07470.930−0.00360.05220.05200.948−0.00130.03700.03640.947
α 0.01600.07740.07930.9530.00840.05410.05490.9550.00370.03800.03750.954
1 σ −0.00840.06800.06900.932−0.00480.04820.04790.944−0.00120.03430.03350.956
α 0.01810.09930.10280.9460.01270.06970.07240.9460.00170.04870.04900.949
2 σ −0.00270.05580.05540.9440.00120.03970.04030.9480.00010.02800.02930.933
α 0.04410.22290.23120.9530.00830.15420.15550.9510.00800.10900.11350.944
3 σ −0.00380.05070.05060.951−0.00220.03590.03600.940−0.00120.02540.02630.931
α 0.07180.36300.36010.9710.02300.25140.24510.9530.02110.17760.18110.942
5 σ −0.00500.04590.04630.945−0.00200.03260.03340.938−0.00040.02310.02400.937
α 0.14460.67990.72060.9580.06780.47150.49760.9480.02850.33000.33990.948
20.8 σ −0.00890.14740.14610.954−0.00800.10440.10440.948−0.00290.07410.07480.941
α 0.01240.07700.07730.9520.00600.05390.05530.9520.00200.03790.03670.956
1 σ −0.01940.13580.13200.949−0.00610.09670.09890.946−0.00110.06870.06890.944
α 0.01690.09910.10210.9530.01090.06960.07330.9470.00140.04870.04910.948
2 σ −0.00500.11170.10970.946−0.00490.07890.07710.953−0.00220.05590.05700.942
α 0.03350.22150.22440.9570.02180.15540.15480.9510.00920.10910.11210.949
3 σ −0.01250.10090.09970.946−0.00150.07190.07600.9340.00070.05090.05110.942
α 0.07890.36380.37050.9550.02380.25140.25340.9550.00980.17670.17410.949
5 σ −0.00590.09210.09290.949−0.00220.06520.06670.942−0.00020.04620.04670.945
α 0.11660.67550.70790.9470.05800.47020.47610.9580.02040.32930.33920.943
Table 4. Empirical bias, RMSE, and CP 95 % for the true values β 0 ( p ) = 0.275 , β 1 ( p ) = 0.098 , β 2 ( p ) = 0.015 , β 3 ( p ) = 0.397 , and α = 1.285 with simulated data.
Table 4. Empirical bias, RMSE, and CP 95 % for the true values β 0 ( p ) = 0.275 , β 1 ( p ) = 0.098 , β 2 ( p ) = 0.015 , β 3 ( p ) = 0.397 , and α = 1.285 with simulated data.
pnBiasRMSECP95%
β 0 ( p ) β 1 ( p ) β 2 ( p ) β 3 ( p ) α β 0 ( p ) β 1 ( p ) β 2 ( p ) β 3 ( p ) α β 0 ( p ) β 1 ( p ) β 2 ( p ) β 3 ( p ) α
0.11000.0180.0000.000−0.0050.0700.0690.0720.0640.0670.1840.9400.9500.9560.9480.960
2000.011−0.0030.0010.0020.0220.0480.0460.0430.0430.1160.9240.9320.9420.9540.952
3000.005−0.0020.003−0.0010.0150.0370.0350.0370.0390.0960.9440.9640.9480.9400.958
4000.0050.0000.000−0.0010.0140.0340.0340.0320.0340.0810.9400.9500.9500.9520.958
0.251000.012−0.003−0.0040.0000.0730.0720.0600.0590.0710.1900.9180.9420.9360.9440.952
2000.006−0.0020.000−0.0010.0260.0460.0460.0440.0440.1180.9380.9460.9560.9480.956
3000.0050.001−0.0010.0020.0160.0370.0400.0360.0390.0900.9380.9560.9540.9440.962
4000.0010.001−0.0010.0020.0210.0320.0330.0340.0310.0840.9440.9420.9400.9600.956
0.51000.010−0.006−0.002−0.0060.0570.0790.0650.0640.0750.1880.9400.9440.9440.9340.950
2000.004−0.0010.001−0.0030.0380.0510.0480.0480.0460.1180.9560.9460.9300.9400.962
3000.0020.0050.001−0.0020.0190.0460.0370.0350.0360.0970.9260.9560.9480.9560.952
400−0.001−0.0010.000−0.0010.0220.0350.030.0340.0310.0810.9560.9600.9400.9540.954
0.75100−0.0070.001−0.0020.0020.0590.1170.0680.0680.0640.1830.9400.9460.9420.9400.954
2000.0010.002−0.001−0.0020.0270.0820.0410.0440.0470.1220.9500.9440.9460.9380.958
300−0.004−0.0010.000−0.0020.0240.0660.0360.0350.0400.0940.9540.9580.9660.9260.956
4000.0020.0010.0010.0010.0140.0530.0330.0320.0310.0790.9660.9360.9580.9500.964
0.9100−0.0310.0000.0010.0040.0670.1960.0700.0660.0650.1980.9200.9440.9520.9360.940
200−0.0200.0010.0000.0040.0350.1400.0500.0430.0440.1230.9280.9440.9420.9640.954
300−0.008−0.001−0.0010.0000.0190.1130.0380.0390.0390.0970.9420.9480.9580.9360.948
400−0.004−0.0010.0000.0020.0130.0950.0340.0310.0330.0830.9520.9420.9440.9480.942
Table 5. Empirical bias, RMSE, and CP 95 % for the true values β 0 ( p ) = 2.975 , β 1 ( p ) = 0.005 , β 2 ( p ) = 0.092 , β 3 ( p ) = 0.958 , β 4 ( p ) = 0.122 , β 5 ( p ) = 0.289 , and α = 7.425 with simulated data.
Table 5. Empirical bias, RMSE, and CP 95 % for the true values β 0 ( p ) = 2.975 , β 1 ( p ) = 0.005 , β 2 ( p ) = 0.092 , β 3 ( p ) = 0.958 , β 4 ( p ) = 0.122 , β 5 ( p ) = 0.289 , and α = 7.425 with simulated data.
pn Bias RMSE CP95%
β 0 ( p ) β 1 ( p ) β 2 ( p ) β 3 ( p ) β 4 ( p ) β 5 ( p ) α β 0 ( p ) β 1 ( p ) β 2 ( p ) β 3 ( p ) β 4 ( p ) β 5 ( p ) α β 0 ( p ) β 1 ( p ) β 2 ( p ) β 3 ( p ) β 4 ( p ) β 5 ( p ) α
0.11000.0090.0010.000−0.003−0.0010.0010.8440.0380.0340.0300.0310.0290.0301.7450.9280.9360.9600.9460.9440.9560.964
2000.0060.0010.0000.0010.0000.0010.3440.0270.0200.0210.0210.0200.0221.0690.9260.9440.9480.9500.9680.9380.952
3000.0030.0000.0000.0000.0000.0000.2100.0200.0180.0190.0180.0170.0170.8290.9420.9360.9240.9640.9580.9540.958
4000.0030.0010.0010.000−0.001−0.0010.1870.0180.0150.0150.0150.0140.0170.6990.9520.9520.9620.9300.9600.9260.952
0.251000.0060.0020.0020.0020.0010.0000.8510.0350.0350.0320.030.0310.0321.8060.9260.9300.9520.9580.9380.9480.956
2000.0020.0000.0000.0010.0010.0010.3410.0230.0210.0240.0240.0230.0231.0530.9460.9440.9260.9420.9380.9420.968
3000.0030.000−0.001−0.0010.0010.0000.2390.0190.0180.0170.0180.0180.0170.7930.9400.9480.9420.9460.9420.9480.956
4000.0010.0000.0010.000−0.0010.0010.2230.0160.0150.0160.0160.0160.0150.7010.9500.950.9380.9320.9500.9560.96
0.51000.0030.0010.0000.0000.000−0.0020.7830.0320.0310.0350.0280.0330.0271.6980.9360.9480.9160.9540.9380.9420.974
2000.0010.001−0.0010.0020.001−0.0010.4410.0220.0210.0200.0200.0210.0201.0650.9500.9440.9400.9540.9280.9500.964
3000.0010.0000.0000.000−0.0010.0020.2490.0180.0170.0170.0180.0180.0180.840.9360.9460.9520.9400.9460.9400.964
400−0.0010.0000.0010.0000.0000.0000.2380.0150.0150.0160.0150.0150.0150.6880.9620.9500.9620.9420.9520.9480.95
0.750100−0.003−0.0010.0010.0020.0040.0000.7300.0370.0300.0310.0290.0330.0301.7410.9300.9420.9460.9560.9260.9300.956
200−0.0020.000−0.001−0.0010.0000.0000.4150.0260.0220.0220.0220.0220.0211.0880.9480.9540.9440.9640.9420.9420.956
300−0.0020.000−0.0020.0000.0010.0010.3020.0210.0170.0170.0170.0170.0200.8440.9460.9540.9580.9600.9620.9320.948
4000.0000.002−0.0010.0010.000−0.0010.2120.0170.0150.0160.0160.0150.0150.6950.9600.9480.9300.9560.9420.9520.964
0.9100−0.0170.0010.0000.0000.0000.0020.9190.0530.0310.0290.0320.0270.0311.9490.8960.9480.9420.9540.9300.9400.954
200−0.0060.0010.0000.0000.000−0.0010.3480.0360.0230.0220.0210.0240.0221.0120.9440.9600.9260.9500.9480.9520.968
300−0.005−0.0010.000−0.0010.0000.0000.2500.0280.0190.0160.0180.0190.0180.8300.9560.9360.9440.9420.9400.9340.942
400−0.0030.0000.0010.000−0.0010.0000.1790.0250.0150.0150.0160.0170.0160.7060.9380.9520.9500.9440.9300.9600.948
Table 6. Maximum likelihood estimates and their standard errors for different distributions (uphn, kum, ughnx, ugompertz, and uweibull) in body fat data.
Table 6. Maximum likelihood estimates and their standard errors for different distributions (uphn, kum, ughnx, ugompertz, and uweibull) in body fat data.
pDistribution β ^ 0 ( p ) SE ( β ^ 0 ( p ) ) β ^ 1 ( p ) SE ( β ^ 1 ( p ) ) β ^ 2 ( p ) SE ( β ^ 2 ( p ) ) β ^ 3 ( p ) SE ( β ^ 3 ( p ) ) β ^ 4 ( p ) SE ( β ^ 4 ( p ) ) β ^ 5 ( p ) SE ( β ^ 5 ( p ) ) α ^ SE ( α ^ ) θ ^ SE ( θ ^ )
0.1uphn−3.3150.1430.0050.0010.0910.006−0.9550.035−0.1180.052−0.2830.0507.9210.867
kum−3.0700.1370.0040.0010.0730.006−0.7300.033−0.0730.048−0.1960.047 4.7050.210
ughnx−3.2940.1550.0040.0010.0800.006−0.7990.035−0.0990.049−0.2280.050 2.6190.116
ugompertz−3.4360.1650.0050.0010.0940.007−0.9670.041−0.1280.062−0.3520.06 3.6570.155
uweibull−3.2140.1680.0060.0010.0860.007−0.9660.042−0.1310.055−0.3550.055 5.9760.243
0.25uphn−3.1740.1470.0050.0010.0920.007−0.9580.036−0.1220.053−0.2890.0517.4250.789
kum−2.8940.1420.0040.0010.0770.006−0.7710.035−0.0740.050−0.2060.050 4.7150.210
ughnx−2.9390.1490.0040.0010.0800.006−0.7990.035−0.0990.049−0.2280.050 2.6190.116
ugompertz−3.3270.1690.0060.0010.0980.007−1.0000.042−0.1360.064−0.3640.062 3.6740.155
uweibull−2.9600.1590.0050.0010.0820.007−0.9190.040−0.1270.053−0.3390.053 5.9800.243
0.5uphn−2.9750.1460.0050.0010.0920.007−0.9580.036−0.1220.053−0.2890.0517.4250.789
kum−2.7540.1490.0040.0010.0820.006−0.8200.037−0.0760.054−0.2160.053 4.7240.210
ughnx−2.6530.1460.0040.0010.0800.006−0.7990.035−0.0990.049−0.2280.050 2.6190.116
ugompertz−3.1900.1780.0060.0010.1040.008−1.0570.044−0.1480.068−0.3850.066 3.7000.155
uweibull−2.6400.1480.0050.0010.0770.006−0.8630.037−0.1210.050−0.3190.050 5.9850.243
0.75uphn−2.7720.1470.0050.0010.0930.007−0.9580.036−0.1220.053−0.2870.0517.4400.791
kum−2.6580.1570.0040.0010.0870.006−0.8700.039−0.0780.057−0.2270.057 4.7310.210
ughnx−2.4490.1440.0040.0010.0800.006−0.7990.034−0.0990.049−0.2280.050 2.6190.116
ugompertz−3.0400.1910.0060.0020.1130.008−1.1510.048−0.1670.074−0.4170.071 3.7370.154
uweibull−2.2720.1380.0050.0010.0720.006−0.8030.035−0.1150.047−0.2980.046 5.9920.243
0.9uphn−2.5380.1480.0050.0010.0910.007−0.9570.036−0.1210.053−0.2850.0517.5010.799
kum−2.5990.1650.0050.0010.0920.007−0.9150.042−0.0780.06−0.2360.060 4.7360.210
ughnx−2.3120.1430.0040.0010.080.006−0.8000.034−0.0990.049−0.2280.050 2.6190.116
ugompertz−2.9020.2120.0070.0020.1270.009−1.2960.054−0.1950.083−0.4660.079 3.7890.153
uweibull−1.8960.1300.0040.0010.0670.006−0.7460.032−0.1090.044−0.2780.043 6.0000.243
Table 7. p-values of Wald test for β j ( p ) , j = 0 , , 5 , coefficients for different distributions (uphn, kum, ughnx, ugompertz, and uweibull) in body fat data.
Table 7. p-values of Wald test for β j ( p ) , j = 0 , , 5 , coefficients for different distributions (uphn, kum, ughnx, ugompertz, and uweibull) in body fat data.
pDistribution p-Value
β 0 ( p ) β 1 ( p ) β 2 ( p ) β 3 ( p ) β 4 ( p ) β 5 ( p )
0.1uphn< 10 3 < 10 3 < 10 3 < 10 3 0.023< 10 3
kum< 10 3 < 10 3 < 10 3 < 10 3 0.127< 10 3
ughnx< 10 3 < 10 3 < 10 3 < 10 3 0.046< 10 3
ugompertz< 10 3 < 10 3 < 10 3 < 10 3 0.038< 10 3
uweibull< 10 3 < 10 3 < 10 3 < 10 3 0.017< 10 3
0.25uphn< 10 3 < 10 3 < 10 3 < 10 3 0.021< 10 3
kum< 10 3 < 10 3 < 10 3 < 10 3 0.139< 10 3
ughnx< 10 3 < 10 3 < 10 3 < 10 3 0.045< 10 3
ugompertz< 10 3 < 10 3 < 10 3 < 10 3 0.034< 10 3
uweibull< 10 3 < 10 3 < 10 3 < 10 3 0.016< 10 3
0.5uphn< 10 3 < 10 3 < 10 3 < 10 3 0.021< 10 3
kum< 10 3 < 10 3 < 10 3 < 10 3 0.154< 10 3
ughnx< 10 3 < 10 3 < 10 3 < 10 3 0.045< 10 3
ugompertz< 10 3 < 10 3 < 10 3 < 10 3 0.029< 10 3
uweibull< 10 3 < 10 3 < 10 3 < 10 3 0.015< 10 3
0.75uphn< 10 3 < 10 3 < 10 3 < 10 3 0.022< 10 3
kum< 10 3 < 10 3 < 10 3 < 10 3 0.173< 10 3
ughnx< 10 3 < 10 3 < 10 3 < 10 3 0.045< 10 3
ugompertz< 10 3 < 10 3 < 10 3 < 10 3 0.023< 10 3
uweibull< 10 3 < 10 3 < 10 3 < 10 3 0.014< 10 3
0.9uphn< 10 3 < 10 3 < 10 3 < 10 3 0.022< 10 3
kum< 10 3 < 10 3 < 10 3 < 10 3 0.193< 10 3
ughnx< 10 3 < 10 3 < 10 3 < 10 3 0.044< 10 3
ugompertz< 10 3 < 10 3 < 10 3 < 10 3 0.019< 10 3
uweibull< 10 3 < 10 3 < 10 3 < 10 3 0.012< 10 3
Table 8. Model selection criteria for different distributions (uphn, kum, ughnx, ugompertz, and uweibull) for body fat data.
Table 8. Model selection criteria for different distributions (uphn, kum, ughnx, ugompertz, and uweibull) for body fat data.
pDistributionNeg2LogLikeAICBICHQIC
0.1uphn−876.795−862.795−836.916−852.436
kum−853.660−839.660−813.781−829.301
ughnx−858.131−844.131−818.251−833.771
ugompertz−763.903−749.903−724.023−739.544
uweibull−852.151−838.151−812.271−827.791
0.25uphn−877.216−863.216−837.336−852.856
kum−856.224−842.224−816.344−831.865
ughnx−858.131−844.131−818.251−833.771
ugompertz−766.845−752.845−726.965−742.486
uweibull−853.364−839.364−813.484−829.004
0.5uphn−877.215−863.215−837.336−852.856
kum−859.03−845.03−819.151−834.671
ughnx−858.131−844.131−818.251−833.771
ugompertz−771.300−757.300−731.420−746.941
uweibull−854.932−840.932−815.052−830.573
0.75uphn−877.222−863.222−837.343−852.863
kum−861.657−847.657−821.778−837.298
ughnx−858.131−844.131−818.251−833.771
ugompertz−777.425−763.425−737.546−753.066
uweibull−856.757−842.757−816.877−832.397
0.9uphn−877.185−863.185−837.306−852.826
kum−863.813−849.813−823.933−839.454
ughnx−858.131−844.131−818.251−833.771
ugompertz−784.259−770.259−744.379−759.899
uweibull−858.602−844.602−818.723−834.243
Table 9. Maximum likelihood estimates with their standard errors across different distributions (uphn, kum, ughnx, ugompertz, and ubs) with PBSC data.
Table 9. Maximum likelihood estimates with their standard errors across different distributions (uphn, kum, ughnx, ugompertz, and ubs) with PBSC data.
pDistribution β ^ 0 SE ( β ^ 0 ( p ) ) β ^ 1 SE ( β ^ 1 ( p ) ) β ^ 2 SE ( β ^ 2 ( p ) ) β ^ 3 SE ( β ^ 3 ( p ) ) α ^ SE ( α ^ ) θ ^ SE ( θ ^ )
0.1uphn−0.4980.2110.0970.0900.0150.0030.3970.0941.2840.104
kum0.2810.1560.0230.0550.0030.0020.0850.065 6.7070.454
ughnx−1.1220.457−0.1520.1670.0010.0070.3800.192 0.6450.028
ugompertz−0.4560.2310.0940.0980.0140.0040.3370.106 2.7890.342
ubs−0.1500.350−0.0050.1430.0080.0060.3440.155 0.8830.040
0.25uphn−0.1800.2100.0970.0900.0150.0030.3970.0941.2850.104
kum0.6530.1910.0290.0710.0040.0030.1130.085 6.7080.454
ughnx0.3200.422−0.1520.1600.0010.0070.3800.179 0.6450.029
ugompertz−0.1690.2520.1020.1060.0150.0040.3650.115 2.7890.342
ubs0.4910.314−0.0060.1300.0070.0050.3160.140 0.8830.040
0.5uphn0.2750.2110.0980.0900.0150.0030.3970.0941.2850.104
kum1.0530.2530.0380.0960.0060.0040.1590.115 6.7080.454
ughnx1.4820.401−0.1520.1580.0010.0060.380.174 0.6450.029
ugompertz0.2770.2850.1120.1180.0170.0050.4010.126 2.7880.342
ubs1.1790.291−0.0060.1220.0070.0050.2980.131 0.8830.040
0.75uphn0.8880.2180.0970.0900.0150.0030.3970.0941.2850.104
kum1.4430.3430.0500.1320.0080.0060.2260.158 6.7080.454
ughnx2.3090.393−0.1520.1570.0010.0060.3800.172 0.6450.029
ugompertz0.9450.3230.1220.1290.0190.0050.4360.138 2.7850.342
ubs1.8200.282−0.0070.1170.0070.0050.2890.126 0.8830.040
0.9uphn1.6260.2380.0970.090.0150.0030.3970.0941.2850.104
kum1.7860.4520.0630.1750.0120.0070.3130.208 6.7080.454
ughnx2.8640.389−0.1520.1570.0010.0060.3800.171 0.6450.029
ugompertz1.8250.3500.1280.1370.0200.0050.4570.146 2.7810.341
ubs2.3380.279−0.0070.1150.0070.0050.2840.123 0.8820.040
Table 10. p-values of Wald test for coefficients β j ( p ) , j = 0 , 1 , 2 , 3 , across different distributions (uphn, kum, ughnx, ugompertz, and ubs) for PBSC data.
Table 10. p-values of Wald test for coefficients β j ( p ) , j = 0 , 1 , 2 , 3 , across different distributions (uphn, kum, ughnx, ugompertz, and ubs) for PBSC data.
pDistribution p-Value
β 0 ( p ) β 1 ( p ) β 2 ( p ) β 3 ( p )
0.1uphn0.0180.277< 10 3 < 10 3
kum0.0710.6770.1640.189
ughnx0.0140.3650.9320.048
ugompertz0.0480.338< 10 3 0.001
ubs0.6690.9740.1650.027
0.25uphn0.3920.277< 10 3 < 10 3
kum< 10 3 0.6850.1570.181
ughnx0.4480.3430.9280.034
ugompertz0.5020.338< 10 3 0.001
ubs0.1180.9650.1560.024
0.5uphn0.1930.276< 10 3 < 10 3
kum< 10 3 0.6950.1470.168
ughnx< 10 3 0.3370.9260.028
ugompertz0.3310.341< 10 3 0.001
ubs< 10 3 0.9580.150.023
0.75uphn< 10 3 0.278< 10 3 < 10 3
kum< 10 3 0.7060.1340.151
ughnx< 10 3 0.3350.9260.027
ugompertz0.0030.346< 10 3 0.002
ubs< 10 3 0.9540.1460.022
0.9uphn< 10 3 0.277< 10 3 < 10 3
kum< 10 3 0.7180.1190.132
ughnx< 10 3 0.3350.9250.026
ugompertz< 10 3 0.350< 10 3 0.002
ubs< 10 3 0.9510.1440.022
Table 11. Model selection criteria for different distributions (uphn, kum, ughnx, ugompertz, and ubs) for PBSC data.
Table 11. Model selection criteria for different distributions (uphn, kum, ughnx, ugompertz, and ubs) for PBSC data.
pDistributionNeg2LogLikeAICBICHQIC
0.1uphn−384.643−374.643−357.26−367.638
kum−384.207−374.207−356.824−367.202
ughnx−190.055−180.055−162.673−173.051
ugompertz−369.921−359.921−342.539−352.916
ubs−313.824−303.824−286.442−296.82
0.25uphn−384.643−374.643−357.26−367.638
kum−384.234−374.234−356.851−367.229
ughnx−190.055−180.055−162.673−173.051
ugompertz−369.908−359.908−342.526−352.904
ubs−313.861−303.861−286.479−296.856
0.5uphn−384.643−374.643−357.26−367.638
kum−384.271−374.271−356.889−367.267
ughnx−190.055−180.055−162.673−173.051
ugompertz−369.865−359.865−342.483−352.861
ubs−313.889−303.889−286.507−296.885
0.75uphn−384.643−374.643−357.26−367.638
kum−384.318−374.318−356.936−367.313
ughnx−190.055−180.055−162.673−173.051
ugompertz−369.793−359.793−342.41−352.788
ubs−313.906−303.906−286.524−296.902
0.9uphn−384.643−374.643−357.26−367.638
kum−384.367−374.367−356.985−367.362
ughnx−190.055−180.055−162.673−173.051
ugompertz−369.728−359.728−342.345−352.723
ubs−313.915−303.915−286.532−296.910
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

Santoro, K.I.; Gómez, Y.M.; Soto, D.; Barranco-Chamorro, I. Unit-Power Half-Normal Distribution Including Quantile Regression with Applications to Medical Data. Axioms 2024, 13, 599. https://doi.org/10.3390/axioms13090599

AMA Style

Santoro KI, Gómez YM, Soto D, Barranco-Chamorro I. Unit-Power Half-Normal Distribution Including Quantile Regression with Applications to Medical Data. Axioms. 2024; 13(9):599. https://doi.org/10.3390/axioms13090599

Chicago/Turabian Style

Santoro, Karol I., Yolanda M. Gómez, Darlin Soto, and Inmaculada Barranco-Chamorro. 2024. "Unit-Power Half-Normal Distribution Including Quantile Regression with Applications to Medical Data" Axioms 13, no. 9: 599. https://doi.org/10.3390/axioms13090599

APA Style

Santoro, K. I., Gómez, Y. M., Soto, D., & Barranco-Chamorro, I. (2024). Unit-Power Half-Normal Distribution Including Quantile Regression with Applications to Medical Data. Axioms, 13(9), 599. https://doi.org/10.3390/axioms13090599

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