Finite-sample distributions of the Wald,likelihood ratio,and Lagrange multiplier test statistics in the classical linear model

In this note, it is shown that the finite-sample distributions of the Wald, likelihood ratio, and Lagrange multiplier statistics in the classical linear regression model are members of the generalized beta model introduced by McDonald and Xu (1995a McDonald, J.B., Xu, Y.J. (1995a). A generalization of the beta distribution with applications. J. Econom. 66:133–152.[Crossref], [Web of Science ®] , [Google Scholar]). This is useful for examining the properties of these test statistics. For example, this characterization makes it easy to find distribution, quantile, and density functions for each test statistic, makes it clear why Wald tests may overreject the null hypothesis using asymptotic critical values, and formalizes the fact that the Lagrange multiplier statistic follows a distribution with bounded support.     

The Marshall–Olkin extended generalized Rayleigh distribution: Properties and applications

In this paper, a new extension for the generalized Rayleigh distribution is introduced. The proposed model, called Marshall–Olkin extended generalized Rayleigh distribution, arises based on the scheme introduced by Marshall and Olkin (1997) Marshall, A.W., Olkin, I. (1997). A new method for adding a parameter to a family of distributions with application to the exponential and Weibull families. Biometrika 84:641–652.[Crossref], [Web of Science ®] , [Google Scholar]. A comprehensive account of the mathematical properties of the new distribution is provided. We discuss about the estimation of the model parameters based on two estimation methods. Empirical applications of the new model to real data are presented for illustrative purposes.     

Generalized Point Estimation with Application to Small Response Estimation

Motivated by a number of drawbacks of classical methods of point estimation, we generalize the definitions of point estimation, and address such notions as unbiasedness and estimation under constraints. The utility of the extension is shown by deriving more reliable estimates for small coefficients of regression models, and for variance components and random effects of mixed models. The extension is in the spirit of generalized confidence intervals introduced by Weerahandi (1993 Weerahandi , S. ( 1993 ). Generalized confidence intervals . J. Amer. Statist. Assoc. 88 : 899 – 905 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) and should encourage much needed further research in point estimation in unbalanced models, multi-variate models, non normal models, and nonlinear models.     

A Simple Method for Testing Variance Components in Unbalanced Nested Model

In this article, we consider the unbalanced case of the three fold nested random effects model under partial balance. The distributions of unweighted sums of squares are obtained first. Using the method of generalized p value introduced in Tsui and Weerahandi (1989 Tsui , K. , Weerahandi , S. ( 1989 ). Generalized p-values in significance testing of hypotheses in the presence of nuisance parameters . Journal of the American Statistical Association 84 : 602 – 607 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]), a new method is proposed for hypothesis tests involving functions of variance components. To evaluate the sizes of the generalized p value, a simulation study is conducted. The results indicate that the proposed method performs well under all examined conditions.     

Gamma-Generalized Inverse Gaussian Class of Distributions with Applications

In this article, a new family of probability distributions with domain in ?⁺ is introduced. This class can be considered as a natural extension of the exponential-inverse Gaussian distribution in Bhattacharya and Kumar (1986 Bhattacharya , S. K. , Kumar , S. ( 1986 ). E-IG model in life testing . Calcutta Statist. Assoc. Bull. 35 : 85 – 90 . [Google Scholar]) and Frangos and Karlis (2004 Frangos , N. , Karlis , D. ( 2004 ). Modelling losses using an exponential-inverse Gaussian distribution . Insur. Math. Econo. 35 : 53 – 67 .[Crossref], [Web of Science ®] , [Google Scholar]). This new family is obtained through the mixture of gamma distribution with generalized inverse Gaussian distribution. We also show some important features such as expressions of probability density function, moments, etc. Special attention is paid to the mixture with the inverse Gaussian distribution, as a particular case of the generalized inverse Gaussian distribution. From the exponential-inverse Gaussian distribution two one-parameter family of distributions are obtained to derive risk measures and credibility expressions. The versatility of this family has been proven in numerical examples.     

On Some New α-Modified Binomial and Poisson Distributions and Their Applications

A new class of α-modified binomial distribution has been proposed, and its distributional properties like probability generating function (pgf), moments, and their interrelations have been studied. Two new α-modified Poisson distributions and Poisson distribution have been obtained as limiting distributions. Modified binomial and Poisson distributions introduced by Berg and Jaworski (1988 Berg , S. , Jaworski , J. ( 1988 ). Modified binomial and Poisson distributions with application in random mapping theory . J. Statist. Plann. Infer. 18 : 313 – 322 . [Google Scholar]) have been seen as particular cases. Mixture distributions of α-modified binomial distributions have been derived. A new distributions called α-modified binomial distributions of type j, their moment properties, limiting distributions as α-modified Poisson distribution of type j, their different convolution properties, pgf, parameter estimators have been studied. Two more new distributions namely Doubly α-modified binomial distributions of type (i, j) and α-modified weighted generalized Poisson distributions of type (j ? 1) have also been studied. Various α-modified binomial and Poisson distributions of Berg and Mutafchiev (1990 Berg , S. , Mutafchiev , L. ( 1990 ). Random mapping with an attracting center: Lagrangian distributions and a regression function . J. Appl. Probab. 27 : 622 – 636 . [Google Scholar]) and Berg and Nowicki (1991 Berg , S. , Nowicki , K. ( 1991 ). Statistical inference for a class of modified power series distributions with applications to random mapping theory . J. Statist. Plann. Infer. 28 : 247 – 261 . [Google Scholar]) have been seen as special cases. Application of some of these proposed distributions have been identified.     

New Generalizations of Cauchy Distribution

The generalized skew-normal distribution introduced by Balakrishnan (2002 Balakrishnan , N. ( 2002 ). Discussion on ‘Skew multivariate models related to hidden truncation and/or selective reporting’ by B. C. Arnold and R. J. Beaver . Test 11 : 37 – 39 .[Web of Science ®] , [Google Scholar]) is used to obtain new generalizations of univariate Cauchy distribution with two parameters, denoted by GC _{m, n} (a, b) with m and n non-negative integer numbers and a, b ∈ R. For cases (m, n) = (1, 2), (m, n) = (2, 1), (m, n) = (0, 3) and (m, n) = (3, 0) explicit forms of the density functions are derived and compared to previous generalizations of Cauchy and skew-Cauchy distributions.     

Fisher Information Matrix for the Generalized Feller-Pareto Distribution

In this article, the exact form of Fisher information matrix for the generalized Feller-Pareto (GFP) distribution is determined. The GFP family is a general distribution which includes a variety of distributions as special cases. For example:

??generalized Singh-Maddala distribution which in turn includes Burr, Fisk, and Lomax distribution (see Kleiber and Kotz, 2003 Kleiber, C., Kotz, S. (2003). Statistical Size Distributions in Economics and Actuarial Sciences. Hoboken, NJ: John Wiley.[Crossref] , [Google Scholar]);

??a Pareto IV distribution which includes a hierarchy of Pareto models, omitted an additional location parameter (see Arnold, 1983 Arnold, B.C. (1983). Pareto Distributions. Fairland, MD: International Cooperative Publishing House. [Google Scholar], 2008 Arnold, B.C. (2008). Pareto and generalized pareto distributions. In: Modeling Income Distributions and Lorenz Curves, Economic Studies in Equality, Social Exclusion and Well-Being, Chotikapanich, D. (Ed.), New York: Springer. pp. 119–145.[Crossref] , [Google Scholar]); and

??beta Lomax distribution which includes, for example, beta II and Lomax distributions.

Application of these distributions covers a wide spectrum of areas ranging from actuarial science, economics, finance to bioscience, telecommunications, and medicine.     

Skew generalized extreme value distribution: Probability-weighted moments estimation and application to block maxima procedure

ABSTRACT

Following the work of Azzalini (1985 Azzalini, A. (1985). A class of distributions which includes the normal ones. Scand. J. Stat. 12:171–178.[Web of Science ®] , [Google Scholar] and 1986 Azzalini, A. (1986). Further results on a class of distributions which includes the normal ones. Statistica 46:199–208. [Google Scholar]) on the skew-normal distribution, we propose an extension of the generalized extreme value (GEV) distribution, the SGEV. This new distribution allows for a better fit of maxima and can be interpreted as both the distribution of maxima when maxima are taken on dependent data and when maxima are taken over a random block size. We propose to estimate the parameters of the SGEV distribution via the probability-weighted moment method. A simulation study is presented to provide an application of the SGEV on block maxima procedure and return level estimation. The proposed method is also implemented on a real-life data.     

An approximation to the convolution of gamma distributions

In general, the exact distribution of a convolution of independent gamma random variables is quite complicated and does not admit a closed form. Of all the distributions proposed, the gamma-series representation of Moschopoulos (1985 Moschopoulos, P. G. (1985). The distribution of the sum of independent gamma random variables. Annals of the Institute of Statistical Mathematics 37Part A:541–544. [Google Scholar]) is relatively simple to implement but for particular combinations of scale and/or shape parameters the computation of the weights of the series can result in complications with too much time consuming to allow a large-scale application. Recently, a compact random parameter representation of the convolution has been proposed by Vellaisamy and Upadhye (2009 Vellaisamy, P., Upadhye, N. S. (2009). On the sums of compound negative binomial and gamma random variables. Journal of Applied Probability 46:272–283.[Crossref], [Web of Science ®] , [Google Scholar]) and it allows to give an exact interpretation to the weights of the series. They describe an infinite discrete probability distribution. This result suggested to approximate Moschopoulos’s expression looking for an approximating theoretical discrete distribution for the weights of the series. More precisely, we propose a general negative binomial distribution. The result is an “excellent” approximation, fast and simple to implement for any parameter combination.     

Joint Moment Generating Functions of Nonadjacent Dual Generalized Order Statistics from Reflected Generalized Pareto Distributions

In this article, a class of reflected generalized Pareto distributions (cf. Burkschat et al., 2003 Burkschat , M. , Cramer , E. , Kamps , U. ( 2003 ). Dual generalized order statistics . Metron LXI ( 1 ): 13 – 26 . [Google Scholar]) is considered. Recurrence relations for joint moment generating functions of higher non adjacent dual generalized order statistics based on a random sample drawn from the considered class are derived. Higher joint moments of non adjacent dual generalized order statistics (reversed ordered order statistics and lower k-records as special cases) are obtained. Recurrence relations for single and product moment generating functions and moments of higher non adjacent dual generalized order statistics are derived. Some results of higher moments of non adjacent generalized order statistics from generalized Pareto distributions (cf. Johnson et al., 1995 Johnson , N. L. , Kotz , S. , Balakrishnan , N. ( 1995 ). Continuous Univariate Distributions. , 2nd ed. Vol. 2. New York : Wiley & Sons . [Google Scholar]), are obtained by using a relation connecting higher moments of generalized order statistics and its dual.     

On the generalized cumulative residual entropy weighted distributions

Recently, Feizjavadian and Hashemi (2015 Feizjavadian, S.H., Hashemi, R. (2015). Mean residual weighted versus the length-biased Rayleigh distribution. J. Stat. Comput. Simul. 85:2823–2838.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]) introduced and studied the mean residual weighted (MRW) distribution as an alternative to the length-biased distribution, by using the concepts of the mean residual lifetime and the cumulative residual entropy (CRE). In this article, a new sequence of weighted distributions is introduced based on the generalized CRE. This sequence includes the MRW distribution. Properties of this sequence are obtained generalizing and extending previous results on the MRW distribution. Moreover, expressions for some known distributions are given, and finite mixtures between the new sequence of weighted distributions and the length-biased distribution are studied. Numerical examples are given to illustrate the new results.     

A new lifetime model with different types of failure rate

A new class of lifetime distributions, which can exhibit with upside-down bathtub-shaped, bathtub-shaped, decreasing, and increasing failure rates, is introduced. The new distribution is constructed by compounding generalized Weibull and logarithmic distributions, leading to improvement on the lifetime distribution considered in Dimitrakopoulou et al. (2007 Dimitrakopoulou, T., K. Adamidis, and S. Loukas. 2007. A lifetime distribution with an upside-down bathtub-shaped hazard function. IEEE Transactions on Reliability 56:308–11.[Crossref], [Web of Science ®] , [Google Scholar]) by having no restriction on the shape parameter and extending the result studied by Tahmasbi and Rezaei (2008 Tahmasbi, R., and S. Rezaei. 2008. A two-parameter lifetime distribution with decreasing failure rate. Computational Statistics and Data Analysis 52:3889–901.[Crossref], [Web of Science ®] , [Google Scholar]) in the general form. The proposed model includes the exponential–logarithmic and Weibull–logarithmic distributions as special cases. Various statistical properties of the proposed class are discussed. Furthermore, estimation via the maximum likelihood method and the Fisher information matrix are discussed. Applications to real data demonstrate that the new class of distributions is more flexible than other recently proposed classes.     

Moments for Some Kumaraswamy Generalized Distributions

Explicit expansions for the moments of some Kumaraswamy generalized (Kw-G) distributions (Cordeiro and de Castro, 2011 Cordeiro, G.M., de Castro, M. (2011). A new family of generalized distributions. J. Statist. Computat. Simul. 81:883–898.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) are derived using special functions. We explore the Kw-normal, Kw-gamma, Kw-beta, Kw-t, and Kw-F distributions. These expressions are given as infinite weighted linear combinations of well-known special functions for which numerical routines are readily available.     

Recurrence relations for single and product moments of generalized order statistics from Marshall–Olkin extended Pareto distributions

Marshall and Olkin (1997 Marshall, A.W., Olkin, I. (1997). A new method for adding a parameter to a family of distributions with application to the exponential and Weibull families. Biometrika 84(3):641–652.[Crossref], [Web of Science ®] , [Google Scholar]) introduced a new method of adding parameter to expand a family of distributions. Using this concept, in this article, the Marshall–Olkin extended Pareto distribution is introduced and some recurrence relations for single and product moments of generalized order statistics are studied. Also the results are deduced for record values and order statistics.     

Shared inverse Gaussian frailty models based on additive hazards

Frailty models are used in the survival analysis to account for the unobserved heterogeneity in individual risks to disease and death. To analyze the bivariate data on related survival times (e.g., matched pairs experiments, twin, or family data), the shared frailty models were suggested. These models are based on the assumption that frailty acts multiplicatively to hazard rate. In this article, we assume that frailty acts additively to hazard rate. We introduce the shared inverse Gaussian frailty models with three different baseline distributions, namely the generalized log-logistic, the generalized Weibull, and exponential power distribution. We introduce the Bayesian estimation procedure using Markov chain Monte Carlo technique to estimate the parameters involved in these models. We apply these models to a real-life bivariate survival dataset of McGilchrist and Aisbett (1991 McGilchrist, C.A., Aisbett, C.W. (1991). Regression with frailty in survival analysis. Biometrics 47:461–466.[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]) related to the kidney infection data, and a better model is suggested for the data.     

The Beta Exponential-Geometric Distribution

A new four-parameter distribution with decreasing, increasing, and upside-down bathtub failure rate called the beta exponential-geometric distribution is proposed. The new distribution, generated from the logit of a beta random variable, extends the exponential-geometric distribution of Adamidis and Loukas (1998 Adamidis , K. , Loukas , S. ( 1998 ). A lifetime distribution with decreasing failure rate . Statistics and Probability Letters 39 : 35 – 42 .[Crossref], [Web of Science ®] , [Google Scholar]) and some other distributions. A comprehensive mathematical treatment of this distribution is provided. Some expressions for the moment generating function, moments, order statistics, and Rényi entropy of the new distribution are derived. Estimation of the stress-strength parameter is also obtained. The model parameters are estimated by the maximum likelihood method and Fisher information matrix is discussed. Finally, an application to a real data set is illustrated.     

Confidence Intervals for Linear Function of Variance Components in Two-Fold Nested Random Effects Model

In the article, we consider the unbalanced case of the two-way nested random effects model under partial balance. Using the method of generalized confidence intervals (GCIs) introduced in Weeranhandi (1993 Weeranhandi , S. ( 1993 ). Generalized confidence intervals . J. Amer. Statist. Assoc. 88 : 899 – 905 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar] 1995 Weeranhandi , S. ( 1995 ). Exact Statistical Methods for Data Analysis . New York : Springer-Verlag . [Google Scholar]), a new method is proposed for constructing confidence intervals on linear function of variance components. To compare the resulted intervals with the Modified Large Sample (MLS) intervals by Hernandez and Burdick (1993 Hernandez , R. P. , Burdick , R. K. ( 1993 ). Confidence intervals on the total variance in unbalanced two-fold nested designs . Biom. J. 35 : 515 – 522 .[Crossref] , [Google Scholar]), a simulation study is conducted. The results indicate that the proposed method performs better than the MLS method, especially for very unbalanced designs.     

Algorithmic Inference of Two-Parameter Gamma Distribution

We provide an estimation procedure of the two-parameter Gamma distribution based on the Algorithmic Inference approach. As a key feature of this approach, we compute the joint probability distribution of these parameters without assuming any prior. To this end, we propose a numerical algorithm which is often beneficial of a highly efficient speed up based on an approximate analytical expression of the probability distribution. We contrast our interval and point estimates with those recently obtained in Son and Oh (2006 Son , Y. S. , Oh , M. ( 2006 ). Bayesian estimation of the two-parameter gamma distribution . Communications in Statistics – Simulation and Computation 35 : 285 – 293 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) for the same problem. From this benchmark we realize that our estimates are both unbiased and more accurate, albeit more dispersed, in some cases, than the competitors' methods, where the dispersion drawback is notably mitigated w.r.t. Bayesian methods by a greater estimate decorrelation. We also briefly discuss the theoretical novelty of the adopted inference paradigm which actually represents a brush up on a Fisher perspective dating to almost a century, made feasible today by the available computational tools.