Home The Unit-Gompertz Quantile Regression Model for the Bounded Responses
Article
Licensed
Unlicensed Requires Authentication

The Unit-Gompertz Quantile Regression Model for the Bounded Responses

  • Josmar Mazucheli , Bruna Alves and Mustafa Ç. Korkmaz EMAIL logo
Published/Copyright: August 4, 2023
Become an author with De Gruyter Brill

ABSTRACT

This paper proposes a regression model for the continuous responses bounded to the unit interval which is based on the unit-Gompertz distribution as an alternative to quantile regression models based on the unit-Birnbaum-Saunders, unit-Weibull, L-Logistic, Kumaraswamy and Johnson SB distributions. Re-parameterizing the unit-Gompertz distribution as a function of its quantile allows us to model the effect of covariates across the entire response distribution, rather than only at the mean. Our proposal sometimes outperforms the other distributions available in the literature. These discoveries are provided by Monte Carlo simulations and one application using a real data set. An R package, including parameter estimation, model checking as well as density, cumulative distribution, quantile and random number generating functions of the unit-Gompertz distribution are developed and can be readily used in applications.

2020 Mathematics Subject Classification: Primary 60E05; 62E10; 62N05

(Communicated by Gejza Wimmer)


REFERENCES

[1] ALTUN, E.: The log-weighted exponential regression model: alternative to the beta regression model, Comm. Statist. Theory Methods 50 (2021), 2306–2321.10.1080/03610926.2019.1664586Search in Google Scholar

[2] ALTUN, E.—CORDEIRO, G. M.: The unit-improved second-degree Lindley distribution: inference and re- gression modeling, Comput. Statist. 35 (2020), 259–279.Search in Google Scholar

[3] ALTUN, E.—EL-MORSHEDY, M.—ELIWA, M. S.: A new regression model for bounded response variable: An alternative to the beta and unit-Lindley regression models, PLoS ONE 16 (2021), e0245627.10.1371/journal.pone.0245627Search in Google Scholar PubMed PubMed Central

[4] ATKINSON, A. C.: Two graphical displays for outlying and influential observations in regression, Biometrika 68 (1981), 13–20.10.1093/biomet/68.1.13Search in Google Scholar

[5] ALTUN, E.—CORDEIRO, G. M.: The unit-improved second-degree Lindley distribution: inference and re- gression modeling, Comput. Statist. 35 (2020), 259–279.10.1007/s00180-019-00921-ySearch in Google Scholar

[6] BALAKRISHNAN, N.: Handbook of the Logistic Distribution, Marcel Dekker, New York, 1992.10.1201/9781482277098Search in Google Scholar

[7] BIRNBAUM, Z. W.—SAUNDERS, S. C.: A new family of life distributions, J. Appl. Probab. 6 (1969), 319–327.10.1017/S0021900200032848Search in Google Scholar

[8] CANCHO, V. G.—BAZÁN, J. L.—DEY, D. K.: A new class of regression model for a bounded response with application in the study of the incidence rate of colorectal cancer, Stat. Methods Med. Res. 29 (2020), 2015–2033.10.1177/0962280219881470Search in Google Scholar PubMed

[9] COX, D. R.—SNELL, E. J.: A general definition of residuals, J. R. Stat. Soc. Ser. B. Stat. Methodol. 30 (1968), 248–265.10.1111/j.2517-6161.1968.tb00724.xSearch in Google Scholar

[10] DA PAZ, R.—BALAKRISHNAN, N.—BAZÁN, J. L.: L-Logistic regression models: Prior sensitivity analysis, robustness to outliers and applications, Braz. J. Probab. Stat. 33 (2019), 455–479.10.1214/18-BJPS397Search in Google Scholar

[11] DUNN, P. K.—SMYTH, G. K.: Randomized quantile residuals, J. Comput. Graph. Stat. 5 (1996), 236–244.10.1080/10618600.1996.10474708Search in Google Scholar

[12] FERRARI, S.—CRIBARI–NETO, F.: Beta regression for modelling rates and proportions, J. Appl. Stat. 31 (2004), 799–815.10.1080/0266476042000214501Search in Google Scholar

[13] GOMPERTZ, B.: On the nature of the function expressive of the law of human mortality, and on a new mode of determining the value of life contingencies, Philos. Trans. Roy. Soc.Search in Google Scholar

[14] GUPTA, R. D.—KUNDU, D.: Theory & methods: Generalized exponential distributions, Aust. N. Z. J. Stat. 41 (1999), 173–188.10.1111/1467-842X.00072Search in Google Scholar

[15] JODRÁ, P.—JIME´NEZ–GAMERO, M. D.: A quantile regression model for bounded responses based on the exponential-geometric distribution, REVSTAT 18 (2020), 415–436.Search in Google Scholar

[16] JOHNSON, N. L.: Systems of frequency curves generated by methods of translation, Biometrika 36 (1949), 149–176.10.1093/biomet/36.1-2.149Search in Google Scholar

[17] KOENKER, R.—BASSETT, J. R. G.: Regression quantiles, Econometrica 46 (1978), 33–50.10.2307/1913643Search in Google Scholar

[18] KORKMAZ, M. Ç.—ALTUN, E.—CHESNEAU, C.—YOUSOF, H. M.: On the unit-Chen distribution with associated quantile regression and applications, Math. Slovaca 72(3) (2022), 765–786.10.1515/ms-2022-0052Search in Google Scholar

[19] KORKMAZ, M. Ç.—CHESNEAU, C.: On the unit Burr-XII distribution with the quantile regression modeling and applications, Comput. Appl. Math. 40 (2021), 1–26.10.1007/s40314-021-01418-5Search in Google Scholar

[20] KORKMAZ, M. Ç.—CHESNEAU, C.—KORKMAZ, Z. S.: A new alternative quantile regression model for the bounded response with educational measurements applications of OECD countries, J. Appl. Stat. 50 (2022).10.1080/02664763.2021.1981834Search in Google Scholar PubMed PubMed Central

[21] KORKMAZ, M. Ç.—CHESNEAU, C.—KORKMAZ, Z. S.: On the arcsecant hyperbolic normal distribution. Properties, quantile regression modeling and applications, Symmetry 13 (2021), 1–24.10.3390/sym13010117Search in Google Scholar

[22] KORKMAZ, M. Ç.—CHESNEAU, C.—KORKMAZ, Z. S.: Transmuted unit Rayleigh quantile regression model: Alternative to beta and Kumaraswamy quantile regression models, U. Politeh. Buch. Ser. A 83 (2021), 149–158.Search in Google Scholar

[23] KORKMAZ, M. Ç.—KORKMAZ, Z. S.: The unit log-log distribution: A new unit distribution with alternative quantile regression modeling and educational measurements applications, J. Appl. Stat. 50 (2023), 889–908.10.1080/02664763.2021.2001442Search in Google Scholar PubMed PubMed Central

[24] KUMARASWAMY, P.: A generalized probability density function for double-bounded random processes, J. Hydrol. 46 (1980), 79–88.10.1016/0022-1694(80)90036-0Search in Google Scholar

[25] LEMONTE, A.—MORENA-ARENAS, G.: On a heavy-tailed parametric quantile regression model for limited range response variables, Comput. Statist. 35 (2020), 379–398.10.1007/s00180-019-00898-8Search in Google Scholar

[26] MAZUCHELI, J.—ALVES, B.: Quantile Regression Modeling for Unit-Gompertz Responses, https://cran. r-project.org/web/packages/ugomquantreg/index.html, R package version 1.0.0, 2021.10.32614/CRAN.package.ugomquantregSearch in Google Scholar

[27] MAZUCHELI, J.—LEIVA, V.—ALVES, B.—MENEZES, A. F. B.: A new quantile regression for modeling bounded data under a unit BirnbaumSaunders distribution with applications in medicine and politics, Symmetry 13(4) (2021), Art. No. 682.10.3390/sym13040682Search in Google Scholar

[28] MAZUCHELI, J.—MENEZES, A. F. B.—FERNANDES, L. B.—DE OLIVEIRA, R. P.—GHITANY, M. E.: The unit-Weibull distribution as an alternative to the Kumaraswamy distribution for the modeling of quantiles conditional on covariates, J. Appl. Stat. 47 (2020), 954–974.10.1080/02664763.2019.1657813Search in Google Scholar PubMed PubMed Central

[29] MAZUCHELI, J.—MENEZES, A. F. B.—CHAKRABORTY, S.: On the one parameter unit-Lindley distribu- tion and its associated regression model for proportion data, J. Appl. Stat. 46 (2019), 700–714.10.1080/02664763.2018.1511774Search in Google Scholar

[30] MAZUCHELI, J.—MENEZES, A. F. B.—DEY, S.: Unit-Gompertz distribution with applications, Statistica 79 (2019), 25–43.Search in Google Scholar

[31] MITNIK, P. A.—BAEK, S. The Kumaraswamy distribution: median-dispersion re-parameterizations for re- gression modeling and simulation-based estimation, Stat. Pap. 54 (2013), 177–192.10.1007/s00362-011-0417-ySearch in Google Scholar

[32] PARZEN, E.: Quantile probability and statistical data modeling, Stat. Sci. 4 (2004), 652–662.10.1214/088342304000000387Search in Google Scholar

[33] RIGBY, R. A.—STASINOPOULOS, D. M.: Generalized additive models for location, scale and shape, J. R. Stat. Soc. Ser. C. Appl. Stat. 54 (2005), 507–554.10.1111/j.1467-9876.2005.00510.xSearch in Google Scholar

[34] SAS Institute Inc.: SAS/STAT 15.1 User’s Guide. The NLMIXED Procedure, SAS Institute Inc., Cary, NC: SAS Institute Inc., 2018.Search in Google Scholar

[35] TADIKAMALLA, P. R.—JOHNSON, N. L.: Systems of frequency curves generated by transformations of logistic variables, Biometrika 69 (1982), 461–465.10.1093/biomet/69.2.461Search in Google Scholar

[36] WEIBULL, W.: A statistical distribution function of wide applicability, J. Appl. Mech. 18 (1951), 293–297.10.1115/1.4010337Search in Google Scholar

Appendix

  1. The UBSA distribution [27] is obtained from the transformation Y = exp(–X), where X ~ Birnbaum-Saunders(α, θ) denotes a Birnbaum-Saunders [7] distributed random variable. The corresponding p.d.f., c.d.f. and q.f. are written as

    f(yα,θ)=12yαθ2π[(αlog(y))12+(αlog(y))32]exp[12θ2(2+log(y)α+αlog(y))],F(yα,θ)=1Φ{1θ[(log(y)α)12(αlog(y))12]}

    and

    Q(τα,θ)=exp{2α2+[θΦ1(1τ)]2θΦ1(1τ)4+[θΦ1(1τ)]2},

    respectively, where 0 < y < 1, θ > 0, α > 0 and Φ(·) is the c.d.f. of the standard normal distribution.

  2. The ULOG distribution [35] is obtained from the transformation Y=exp(Xαθ)1+exp(Xαθ)) , where X ~ Log (0,1) denotes a standard Logistic distributed random variable [6]. The corresponding p.d.f., c.d.f. and q.f. are written as

    f(yα,θ)=θexp(α)(y1y)θ1[1+exp(α)(y1y)θ]2,F(yα,θ)=exp(α)(y1y)θ1+exp(α)(y1y)θ,

    and

    Q(τα,θ)=exp(αθ)(τ1τ)1θ1+exp(αθ)(τ1τ)1θ,

    respectively, where 0 < y, τ < 1, α > 0 and θ > 0. From above q.f., the parameter α can be re-parameterized as α=g1(μ)=logτ1τθlogμ1μ .

  3. The JOSB distribution [16] can be obtained from the transformation Y=exp(Xαθ)1+exp (Xαθ) ) , where X ~ N (0, 1) denotes a standard Normal random variable. The corresponding p.d.f., c.d.f. and q.f. are written as

    f(yα,θ)=θ2π1y(1y)exp{12[α+θlog(y1y)]2},F(yα,θ)=Φ[α+θlog(y1y)],Q(τα,θ)=exp[Φ1(τ)αθ]1+exp[Φ1(τ)αθ],

    respectively, where 0 < y, τ < 1, α ∈ ℝ and θ > 0 and Φ–1(·) is the qf of the standard Normal distribution. From above q.f., the parameter α can be re-parameterized as α=g1(μ)=Φ1(τ)θlog(μ1μ) .

  4. The KUMA distribution [24] can be obtained from the transformation Y = exp(–X), where X ~ EE(α, θ) denotes a Exponentiated-Exponential distributed random variable [14]. The corresponding p.d.f., c.d.f. and q.f. are written as

    f(y|α,θ)=αθyθ1(1yθ)α1,F(y|α,θ)=1(1yθ)α,Q(τ|α,θ)=[1(1τ)1α]1θ,

    respectively, where 0 < y, τ < 1, α > 0 and θ > 0 are shapes parameters. Prom above q.f., the parameter α can be re-parameterized as α=g1(μ)=log(1τ)log(1μθ) .

  5. The UWEI distribution [28] is obtained from the transformation Y = exp(–X), where X ~ WEI(α, θ) denotes a Weibull distributed random variable [36]. The corresponding p.d.f., c.d.f. and q.f. are written as

    f(yα,θ)=αθylog(y)θ1expαlog(y)θ,F(yα,θ)=expαlog(y)θ,Qτα,θ=explog(τ)α1θ,

    respectively, where 0 < y, τ < 1, α > 0 and θ > 0. Prom above q.f., the parameter α can be re-parameterized as α=g1(μ)=log(τ)[log(μ)]θ .

Received: 2021-09-14
Accepted: 2022-01-26
Published Online: 2023-08-04

© 2023 Mathematical Institute Slovak Academy of Sciences

Downloaded on 14.9.2025 from https://www.degruyterbrill.com/document/doi/10.1515/ms-2023-0077/html
Scroll to top button