Search

The inverted (or inverse) distributions are sometimes very useful to explore additional properties of the phenomenons which non-inverted distributions cannot. We introduce a new inverted mo...

We study some mathematical properties of a new generator of continuous distributions with one extra parameter called the odd power Cauchy family including asymptotics, linear representation...

We introduce a new class of distributions called the Burr XII system of densities with two extra positive parameters. We provide a comprehensive treatment of some of its mathematical properties. We estimate the model parameters by maximum likelihood. We assess the performance of the maximum likelihood estimators in terms of biases and mean squared errors by means of a simulation study. We also introduce a new family of regression models based on this system of densities. The usefulness of the proposed models is illustrated by means of three real data sets.

We introduce a new class of continuous distributions named the Topp–Leone odd log-logistic family, which extends the one-parameter distribution pioneered by Topp and Leone [A family of J-shaped frequency functions. J Amer Statist Assoc. 1955;50:209–219]. We study some of its mathematical properties and describe two special cases. Further, we propose a regression model based on the new Topp–Leone odd log-logistic Weibull distribution. The usefulness and flexibility of the proposed family are illustrated by means of three real data sets.

We study a new family of continuous distributions with two extra shape parameters called the Burr generalized family of distributions. We investigate the shapes of the density and hazard rate function. We derive explicit expressions for some of its mathematical quantities. The estimation of the model parameters is performed by maximum likelihood. We prove the flexibility of the new family by means of applications to two real data sets. Furthermore, we propose a new extended regression model based on the logarithm of the Burr generalized distribution. This model can be very useful to the analysis of real data and provide more realistic fits than other special regression models.

Recently, Ristic and Nadarajah [A new lifetime distribution. J Stat Comput Simul. 2014;84:135–150] introduced the Poisson generated family of distributions and investigated the properties of a special case named the exponentiated-exponential Poisson distribution. In this paper, we study general mathematical properties of the Poisson-X family in the context of the T-X family of distributions pioneered by Alzaatreh et al. [A new method for generating families of continuous distributions. Metron. 2013;71:63–79], which include quantile, shapes of the density and hazard rate functions, asymptotics and Shannon entropy. We obtain a useful linear representation of the family density and explicit expressions for the ordinary and incomplete moments, mean deviations and generating function. One special lifetime model called the Poisson power-Cauchy is defined and some of its properties are investigated. This model can have flexible hazard rate shapes such as increasing, decreasing, bathtub and upside-down ba...

In applications of survival analysis, the failure rate function may frequently present a unimodal shape. In such cases, the log-normal and log-logistic distributions are used. In this paper, we shall be concerned only with parametric forms, so a location-scale regression model based on the odd log-logistic Weibull distribution is proposed for modelling data with a decreasing, increasing, unimodal and bathtub failure rate function as an alternative to the log-Weibull regression model. For censored data, we consider a classic method to estimate the parameters of the proposed model. We derive the appropriate matrices for assessing local influences on the parameter estimates under different perturbation schemes and present some ways to assess global influences. Further, for different parameter settings, sample sizes and censoring percentages, various simulations are performed. In addition, the empirical distribution of some modified residuals is determined and compared with the standard normal distribution. T...

In this paper, we formulate a simple latent cure rate model with repair mechanism for a cell exposed to radiation. This latent approach is a flexible alternative to the models proposed by Klebanov et al. [A stochastic model of radiation carcinogenesis: latent time distributions and their properties. Math Biosci. 1993;18:51–75], Kim et al. [A new threshold regression model for survival data with a cure fraction. Lifetime Data Anal. 2011;17:101–122], and is along the lines of the destructive cure rate model formulated recently by Rodrigues et al. [Destructive weighted Poisson cure rate model. Lifetime Data Anal. 2011b;17:333–346]. A new version of the modified Gompertz model and the promotion cure rate model that takes into account the first passage time of reaching the critical point are discussed, and the estimation of tumor size at detection is then addressed from the Bayesian viewpoint. In addition, a simulation study and an application to real data set illustrate the usefulness of the proposed cure rat...

We propose some statistical tools for diagnosing the class of generalized Weibull linear regression models [A.A. Prudente and G.M. Cordeiro, Generalized Weibull linear models, Comm. Statist. Theory Methods 39 (2010), pp. 3739–3755]. This class of models is an alternative means of analysing positive, continuous and skewed data and, due to its statistical properties, is very competitive with gamma regression models. First, we show that the Weibull model induces ma-ximum likelihood estimators asymptotically more efficient than the gamma model. Standardized residuals are defined, and their statistical properties are examined empirically. Some measures are derived based on the case-deletion model, including the generalized Cook's distance and measures for identifying influential observations on partial F-tests. The results of a simulation study conducted to assess behaviour of the global influence approach are also presented. Further, we perform a local influence analysis under the case-weights, response and e...

In this article, we derive explicit expansions for the moments of beta generalized distributions from power series expansions for the quantile functions of the baseline distributions. We apply our formula to the beta normal, beta Student t, beta gamma and beta beta generalized distributions. We propose a simple way to express the quantile function of any beta generalized distribution as a power series expansion with known coefficients.

For any continuous baseline G distribution, Zografos and Balakrishnan [On families of beta- and generalized gamma-generated distributions and associated inference. Statist Methodol. 2009;6:344–362] proposed a generalized gamma-generated distribution with an extra positive parameter. A new three-parameter continuous distribution called the gamma-Lomax distribution, which extends the Lomax distribution is proposed and studied. Various structural properties of the new distribution are derived including explicit expressions for the moments, generating and quantile functions, mean deviations and Renyi entropy. The estimation of the model parameters is performed by maximum likelihood. We also determine the observed information matrix. An application illustrates the usefulness of the proposed model.

For any continuous baseline G distribution, Zografos and Balakrishnan [On families of beta- and generalized gamma-generated distributions and associated inference. Statist Methodol. 2009;6:344–362] introduced the generalized gamma-generated distribution with an extra positive parameter. A new three-parameter continuous model called the gamma-linear failure rate (LFR) distribution, which extends the LFR model, is proposed and studied. Various structural properties of the new distribution are derived, including some explicit expressions for ordinary and incomplete moments, generating function, probability-weighted moments, mean deviations and Renyi and Shannon entropies. We estimate the model parameters by maximum likelihood and obtain the observed information matrix. The new model is modified to cope with possible long-term survivors in lifetime data. We illustrate the usefulness of the proposed model by means of two applications to real data.

In this paper, we study some mathematical properties of the beta Weibull (BW) distribution, which is a quite flexible model in analysing positive data. It contains the Weibull, exponentiated exponential, exponentiated Weibull and beta exponential distributions as special sub-models. We demonstrate that the BW density can be expressed as a mixture of Weibull densities. We provide their moments and two closed-form expressions for their moment-generating function. We examine the asymptotic distributions of the extreme values. Explicit expressions are derived for the mean deviations, Bonferroni and Lorenz curves, reliability and two entropies. The density of the BW-order statistics is a mixture of Weibull densities and two closed-form expressions are derived for their moments. The estimation of the parameters is approached by two methods: moments and maximum likelihood. We compare the performances of the estimates obtained from both the methods by simulation. The expected information matrix is derived. For th...

In this paper, we introduce the Birnbaum-Saunders power series class of distributions which is obtained by compounding Birnbaum-Saunders and power series distributions. The new class of distributions has as a particular case the two-parameter Birnbaum-Saunders distribution. The hazard rate function of the proposed class can be increasing and upside-down bathtub shaped. We provide important mathematical properties such as moments, order statistics, estimation of the parameters and inference for large sample. Three special cases of the new class are investigated with some details. We illustrate the usefulness of the new distributions by means of two applications to real data sets.

In this paper, we propose a new three-parameter model called the exponential–Weibull distribution, which includes as special models some widely known lifetime distributions. Some mathematical properties of the proposed distribution are investigated. We derive four explicit expressions for the generalized ordinary moments and a general formula for the incomplete moments based on infinite sums of Meijer's G functions. We also obtain explicit expressions for the generating function and mean deviations. We estimate the model parameters by maximum likelihood and determine the observed information matrix. Some simulations are run to assess the performance of the maximum likelihood estimators. The flexibility of the new distribution is illustrated by means of an application to real data.

We introduce an extended Burr III distribution as an important model for problems in survival analysis and reliability. The new distribution can be expressed as a linear combination of Burr III distributions and then it has tractable properties for the ordinary and incomplete moments, generating and quantile functions, mean deviations and reliability. The density of its order statistics can be given in terms of an infinite linear combination of Burr III densities. The estimation of the model parameters is approached by maximum likelihood and the observed information matrix is derived. The proposed model is applied to a real data set to illustrate its potentiality.

For the first time, we introduce the beta log-normal (LN) distribution for which the LN distribution is a special case. Various properties of the new distribution are discussed. Expansions for the cumulative distribution and density functions that do not involve complicated functions are derived. We obtain expressions for its moments and for the moments of order statistics. The estimation of parameters is approached by the method of maximum likelihood, and the expected information matrix is derived. The new model is quite flexible in analysing positive data as an important alternative to the gamma, Weibull, generalized exponential, beta exponential, and Birnbaum–Saunders distributions. The flexibility of the new distribution is illustrated in an application to a real data set.

The use of bivariate distributions plays a fundamental role in survival and reliability studies. In this paper, we introduce a location-scale model for bivariate survival times based on the copula to model the dependence of bivariate survival data with cure fraction. We create the correlation structure between the failure times using the Clayton family of copulas, which is assumed to have any distribution. It turns out that the model becomes very flexible with respect to the choice of the marginal distributions. For the proposed model, we consider inferential procedures based on constrained parameters under maximum likelihood. We derive the appropriate matrices for assessing local influence under different perturbation schemes and present some ways to perform global influence analysis. The relevance of the approach is illustrated using a real data set and a diagnostic analysis is performed to select an appropriate model.

The Weibull distribution is one of the most important distributions in reliability. For the first time, we introduce the beta exponentiated Weibull distribution which extends recent models by Lee et al. [Beta-Weibull distribution: some properties and applications to censored data, J. Mod. Appl. Statist. Meth. 6 (2007), pp. 173–186] and Barreto-Souza et al. [The beta generalized exponential distribution, J. Statist. Comput. Simul. 80 (2010), pp. 159–172]. The new distribution is an important competitive model to the Weibull, exponentiated exponential, exponentiated Weibull, beta exponential and beta Weibull distributions since it contains all these models as special cases. We demonstrate that the density of the new distribution can be expressed as a linear combination of Weibull densities. We provide the moments and two closed-form expressions for the moment-generating function. Explicit expressions are derived for the mean deviations, Bonferroni and Lorenz curves, reliability and entropies. The density of...

A five-parameter extension of the Weibull distribution capable of modelling a bathtub-shaped hazard rate function is introduced and studied. The beauty and importance of the new distribution lies in its ability to model both monotone and non-monotone failure rates that are quite common in lifetime problems and reliability. The proposed distribution has a number of well-known lifetime distributions as special sub-models, such as the Weibull, extreme value, exponentiated Weibull, generalized Rayleigh and modified Weibull (MW) distributions, among others. We obtain quantile and generating functions, mean deviations, Bonferroni and Lorenz curves and reliability. We provide explicit expressions for the density function of the order statistics and their moments. For the first time, we define the log-Kumaraswamy MW regression model to analyse censored data. The method of maximum likelihood is used for estimating the model parameters and the observed information matrix is determined. Two applications illustrate t...

We formulate and study a four-parameter lifetime model called the beta extended half-normal distribution. This model includes as sub-models the exponential, extended half-normal and half-normal distributions. We derive expansions for the new density function which do not depend on complicated functions. We obtain explicit expressions for the moments and incomplete moments, generating function, mean deviations, Bonferroni and Lorenz curves and Renyi entropy. In addition, the model parameters are estimated by maximum likelihood. We provide the observed information matrix. The new model is modified to cope with possible long-term survivors in the data. The usefulness of the new distribution is shown by means of two real data sets.

For the first time, we introduce a generalized form of the exponentiated generalized gamma distribution [Cordeiro et al. The exponentiated generalized gamma distribution with application to lifetime data, J. Statist. Comput. Simul. 81 (2011), pp. 827–842.] that is the baseline for the log-exponentiated generalized gamma regression model. The new distribution can accommodate increasing, decreasing, bathtub- and unimodal-shaped hazard functions. A second advantage is that it includes classical distributions reported in the lifetime literature as special cases. We obtain explicit expressions for the moments of the baseline distribution of the new regression model. The proposed model can be applied to censored data since it includes as sub-models several widely known regression models. It therefore can be used more effectively in the analysis of survival data. We obtain maximum likelihood estimates for the model parameters by considering censored data. We show that our extended regression model is very useful...

For any continuous baseline G distribution [G.M. Cordeiro and M. de Castro, A new family of generalized distributions, J. Statist. Comput. Simul. 81 (2011), pp. 883–898], proposed a new generalized distribution (denoted here with the prefix ‘Kw-G’ (Kumaraswamy-G)) with two extra positive parameters. They studied some of its mathematical properties and presented special sub-models. We derive a simple representation for the Kw-G density function as a linear combination of exponentiated-G distributions. Some new distributions are proposed as sub-models of this family, for example, the Kw-Chen [Z.A. Chen, A new two-parameter lifetime distribution with bathtub shape or increasing failure rate function, Statist. Probab. Lett. 49 (2000), pp. 155–161], Kw-XTG [M. Xie, Y. Tang, and T.N. Goh, A modified Weibull extension with bathtub failure rate function, Reliab. Eng. System Safety 76 (2002), pp. 279–285] and Kw-Flexible Weibull [M. Bebbington, C.D. Lai, and R. Zitikis, A flexible Weibull extension, Reliab. Eng. S...

Lemonte and Cordeiro [Birnbaum–Saunders nonlinear regression models, Comput. Stat. Data Anal. 53 (2009), pp. 4441–4452] introduced a class of Birnbaum–Saunders ( ) nonlinear regression models potentially useful in lifetime data analysis. We give a general matrix Bartlett correction formula to improve the likelihood ratio (LR) tests in these models. The formula is simple enough to be used analytically to obtain several closed-form expressions in special cases. Our results generalize those in Lemonte et al. [Improved likelihood inference in Birnbaum–Saunders regressions, Comput. Stat. Data Anal. 54 (2010), pp. 1307–1316], which hold only for the ( linear regression models. We consider Monte Carlo simulations to show that the corrected tests work better than the usual LR tests.

For the first time, a five-parameter distribution, called the Kumaraswamy Burr XII (KwBXII) distribution, is defined and studied. The new distribution contains as special models some well-known distributions discussed in lifetime literature, such as the logistic, Weibull and Burr XII distributions, among several others. We obtain the complete moments, incomplete moments, generating and quantile functions, mean deviations, Bonferroni and Lorenz curves and reliability of the KwBXII distribution. We provide two representations for the moments of the order statistics. The method of maximum likelihood and a Bayesian procedure are adopted for estimating the model parameters. For different parameter settings and sample sizes, various simulation studies are performed and compared to the performance of the KwBXII distribution. Three applications to real data sets demonstrate the usefulness of the proposed distribution and that it may attract wider applications in lifetime data analysis.

Interval-valued variables have become very common in data analysis. Up until now, symbolic regression mostly approaches this type of data from an optimization point of view, considering neither the probabilistic aspects of the models nor the nonlinear relationships between the interval response and the interval predictors. In this article, we formulate interval-valued variables as bivariate random vectors and introduce the bivariate symbolic regression model based on the generalized linear models theory which provides much-needed exibility in practice. Important inferential aspects are investigated. Applications to synthetic and real data illustrate the usefulness of the proposed approach.

Kumaraswamy [Generalized probability density-function for double-bounded random-processes, J. Hydrol. 462 (1980), pp. 79–88] introduced a distribution for double-bounded random processes with hydrological applications. For the first time, based on this distribution, we describe a new family of generalized distributions (denoted with the prefix ‘Kw’) to extend the normal, Weibull, gamma, Gumbel, inverse Gaussian distributions, among several well-known distributions. Some special distributions in the new family such as the Kw-normal, Kw-Weibull, Kw-gamma, Kw-Gumbel and Kw-inverse Gaussian distribution are discussed. We express the ordinary moments of any Kw generalized distribution as linear functions of probability weighted moments (PWMs) of the parent distribution. We also obtain the ordinary moments of order statistics as functions of PWMs of the baseline distribution. We use the method of maximum likelihood to fit the distributions in the new class and illustrate the potentiality of the new model with a...

A four-parameter extension of the generalized gamma distribution capable of modelling a bathtub-shaped hazard rate function is defined and studied. The beauty and importance of this distribution lies in its ability to model monotone and non-monotone failure rate functions, which are quite common in lifetime data analysis and reliability. The new distribution has a number of well-known lifetime special sub-models, such as the exponentiated Weibull, exponentiated generalized half-normal, exponentiated gamma and generalized Rayleigh, among others. We derive two infinite sum representations for its moments. We calculate the density of the order statistics and two expansions for their moments. The method of maximum likelihood is used for estimating the model parameters and the observed information matrix is obtained. Finally, a real data set from the medical area is analysed.

In this paper, we derive Bartlett and Bartlett-type corrections [G.M. Cordeiro and S.L.P. Ferrari 1991, A modified score test statistic having chi-squared distribution to order n −1 , Biometrika 78 (1991), pp. 573–582] to improve the likelihood ratio and Rao's score statistics for testing the mean parameter and the concentration parameter in the von Mises distribution. Simple formulae are suggested for the corrections valid for small and large values of the concentration parameter that do not depend on the modified Bessel functions and can be useful in practical applications.

In this article, we derive general matrix formulae for second-order biases of maximum likelihood estimators (MLEs) in a class of heteroscedastic symmetric nonlinear regression models, thus generalizing some results in the literature. This class of regression models includes all symmetric continuous distributions, and has a wide range of practical applications in various fields such as engineering, biology, medicine and economics, among others. The variety of distributions with different kurtosis coefficients than the normal may give more flexibility in the choice of an appropriate distribution, particularly to accommodate outlying and influential observations. We derive a joint iterative process for estimating the mean and dispersion parameters. We also present simulation studies for the biases of the MLEs.

In this work we present a source code for programs (scripts) that may be used with the symbolic computation software Maple and Mathematica, for generating analytic expressions for the second-order bias corrections of the maximum likelihood estimators in regular two parameter continuous distributions. The scripts are tested on more than 20 continuous distributions, and the results are compared with those published in earlier works, confirming all of the previously reported expressions.

In this paper, we give matrix formulae of order 𝒪(n −1), where n is the sample size, for the first two moments of Pearson residuals in exponential family nonlinear regression models [G.M. Cordeiro and G.A. Paula, Improved likelihood ratio statistic for exponential family nonlinear models, Biometrika 76 (1989), pp. 93–100.]. The formulae are applicable to many regression models in common use and generalize the results by Cordeiro [G.M. Cordeiro, On Pearson's residuals in generalized linear models, Statist. Prob. Lett. 66 (2004), pp. 213–219.] and Cook and Tsai [R.D. Cook and C.L. Tsai, Residuals in nonlinear regression, Biometrika 72(1985), pp. 23–29.]. We suggest adjusted Pearson residuals for these models having, to this order, the expected value zero and variance one. We show that the adjusted Pearson residuals can be easily computed by weighted linear regressions. Some numerical results from simulations indicate that the adjusted Pearson residuals are better approximated by the standard normal distribu...

In this work, we present source code for a program (script) that may be used with the algebraic manipulation software Maple, for generating analytic expressions for Bartlett and Bartlett-type corrections to improve likelihood ratio and score statistics on general two parameter distributions. The script is tested on normal, inverse Gaussian, gamma and Weibull distributions, and the results are compared with those published in earlier works, confirming all of the previously reported expressions. The presented script may be in principle applied for calculating the Bartlett and Bartlett-type corrections for any two parameter distribution under orthogonality of parameters.

The article derives Bartlett corrections for improving the chi-square approximation to the likelihood ratio statistics in a class of symmetric nonlinear regression models. This is a wide class of models which encompasses the t model and several other symmetric distributions with longer-than normal tails. In this paper we present, in matrix notation, Bartlett corrections to likelihood ratio statistics in nonlinear regression models with errors that follow a symmetric distribution. We generalize the results obtained by Ferrari, S. L. P. and Arellano-Valle, R. B. (1996). Modified likelihood ratio and score tests in linear regression models using the t distribution. Braz. J. Prob. Statist., 10, 15–33, who considered a t distribution for the errors, and by Ferrari, S. L. P. and Uribe-Opazo, M. A. (2001). Corrected likelihood ratio tests in a class of symmetric linear regression models. Braz. J. Prob. Statist., 15, 49–67, who considered a symmetric linear regression model. The formulae derived are simple enough...

We derive general formulae for second-order biases of maximum partial likelihood estimates for the Cox regression model that are easy to compute and yield bias-corrected maximum partial likelihood estimates to order n −1. Monte Carlo simulations indicate smaller biases without variance inflation. †On sabbatical leave from the Departamento de Estatistica, Universidade Federal de Minas Gerais, 31270-901 – Belo Horizonte – MG, Brazil.

In this paper we derive a general closed-form expression for the Bartlett correction for testing a scalar parameter of a two-parameter exponential family model. The correction has the advantage for algebraical and numerical purposes since it involves only trivial operations on suitably defined functions. The formula derived is general enough to cover many important and commonly used distributions. Simulations show that the corrected likelihood ratio tests have empirical sizes closer to the nominal size than the classical uncorrected tests.

We present correction formulae to improve likelihood ratio and score teats for testing simple and composite hypotheses on the parameters of the beta distribution. As a special case of our results we obtain improved tests for the hypothesis that a sample is drawn from a uniform distribution on (0, 1). We present some Monte Carlo investigations to show that both corrected tests have better performances than the classical likelihood ratio and score tests at least for small sample sizes.

In this paper, we derive general formulae for second-order biases of maximum likelihood estimates of the regression, dispersion and precision parameters in nonlinear regression models with t distributed errors. Our formulae are easy to compute, giving the biases by means of ordinary linear regressions. They generalize some previous results due to Cook, Tsai and Wei (1986); Cordeiro and McCullagh (1991) and Cordeiro and Vasconcellos (1997). We derive simple closed-form expressions for these biases in special models. We present some simulations that indicate that the bias-corrected estimates produce a reduction in bias without a given corresponding increase in variability. The bias correction achieves a second-order reduction in the mean square errors of the corrected estimates.

This paper gives closed-form expressions for bias-corrected maximum likelihood estimates of the parameters of the beta distribution that can be used to define bias-corrected estimates that are nearly unbiased. Some approximations based on asymptotic expansions for the bias corrections are given. We also present simulation results comparing the performances of the maximum likelihood estimates and corrected ones. The results suggest that bias-corrected estimates have better finite-sample performance than standard maximum likelihood estimates.

This paper addresses the issue of improving likelihood ratio (LR) statistics for exponential censored data. We give, in matrix notation, a simple formula to compute Bartlett corrections for LR statistics in exponential regression models. The formula derived is simple enough to be used analytically to obtain closed-form Bartlett corrections in special models. We also present Monte Carlo simulations to compare the performance of the usual LR test and its Bartlett modified version. The Bartlett corrections seem to be effective in pushing the true sizes of the tests towards the nominal levels.

In this paper we derive a general closed-form expression for the Bartlett correction for the test of Ho:λ = λ(0), where λ is a scalar parameter of a one-parameter exponential family model. Our results are general enough to cover many important and commonly used distributions. Several special cases and classes of variance functions of considerable importance are discussed, and some approximations based on asymptotic expansions are given. We also use a graphical analysis to examine how the correction varies with λ in some special cases. Stimulation results are also given.

Cordeiro (1983) has derived the expected value of the deviance for generalized linear models correct to terms of order n -1 being the sample size. Then a Bartlett-type factor is available for correcting the first moment of the deviance and for fitting its distribution. If the model is correct, the deviance is not, in general, distributed as chi-squared even asymptotically and very little is known about the adequacy of the X 2 approximation. This paper through simulation studies examines the behaviour of the deviance and a Bartlett adjusted deviance for testing the goodness-of-fit of a generalized linear model. The practical use of such adjustment is illustrated for some gamma and Poisson models. It is suggested that the null distribution of the adjusted deviance is better approximated by chi-square than the distribution of the deviance.

We give simple matrix formulae for computing Bartlett-type corrections for the score statistic for testing a general composite null hypothesis. These formulae, involving only simple operations on matrices and vectors, have computational advantages over the tensor notation of Harris' (1985) asymptotic expansion for the null distribution of the score statistic. Two situations are discussed. First, the test of a general composite null hypothesis is considered. Second, the test of a scalar parameter is studied under the assumption of global orthogonality of the parameter of interest and the vector of remaining parameters. The formulae can be implemented in algebraic computer systems in order to obtain closed-form expressions for the corrections in a wide range of problems. Two illustrative examples are presented.