On the type-B integral equation and the distribution of wilks statist for testing independence of several groups of variables

In this paper, WILKS'type-B integral equation is solved in the general form of a series of beta functions and a series of weighted gamma functions as proposed by WALD and BROOKNER 1941. The coefficients in both representations can be obtained by explicit recurrence relartions, therefore the results solve many distributional problems and have the fewest computational difficulties of any representation that has surfaced to date. The radius of convergence of the second series representation is given, whereas the convergence property of the first series representation is given, whereas the convergence property of the first series representation was studied by WALD and Brookner. The exact null distributions of WILKS' statistic A for testing the independence of several groups of variables and of V = -log A are given. The coefficients in all the series representration can be computed recursilvely and hence can be obtained easily with the help of modern computatinal facilities     

A Study of Interval Censoring in Parametric Regression Models

Parametric models for interval censored data can now easily be fitted with minimal programming in certain standard statistical software packages. Regression equations can be introduced, both for the location and for the dispersion parameters. Finite mixture models can also be fitted, with a point mass on right (or left) censored observations, to allow for individuals who cannot have the event (or already have it). This mixing probability can also be allowed to follow a regression equation.Here, models based on nine different distributions are compared for three examples of heavily censored data as well as a set of simulated data. We find that, for parametric models, interval censoring can often be ignored and that the density, at centres of intervals, can be used instead in the likelihood function, although the approximation is not always reliable. In the context of heavily interval censored data, the conclusions from parametric models are remarkably robust with changing distributional assumptions and generally more informative than the corresponding non-parametric models.     

The S-system computation of non-central gamma distribution

The non-central gamma distribution can be regarded as a general form of non-central χ² distributions whose computations were thoroughly investigated (Ruben, H., 1974, Non-central chi-square and gamma revisited. Communications in Statistics, 3(7), 607–633; Knüsel, L., 1986, Computation of the chi-square and Poisson distribution. SIAM Journal on Scientific and Statistical Computing, 7, 1022–1036; Voit, E.O. and Rust, P.F., 1987, Noncentral chi-square distributions computed by S-system differential equations. Proceedings of the Statistical Computing Section, ASA, pp. 118–121; Rust, P.F. and Voit, E.O., 1990, Statistical densities, cumulatives, quantiles, and power obtained by S-systems differential equations. Journal of the American Statistical Association, 85, 572–578; Chattamvelli, R., 1994, Another derivation of two algorithms for the noncentral χ² and F distributions. Journal of Statistical Computation and Simulation, 49, 207–214; Johnson, N.J., Kotz, S. and Balakrishnan, N., 1995, Continuous Univariate Distributions, Vol. 2 (2nd edn) (New York: Wiley). Both distributional function forms are usually in terms of weighted infinite series of the central one. The ad hoc approximations to cumulative probabilities of non-central gamma were extended or discussed by Chattamvelli, Knüsel and Bablok (Knüsel, L. and Bablok, B., 1996, Computation of the noncentral gamma distribution. SIAM Journal on Scientific Computing, 17, 1224–1231), and Ruben (Ruben, H., 1974, Non-central chi-square and gamma revisited. Communications in Statistics, 3(7), 607–633). However, they did not implement and demonstrate proposed numerical procedures. Approximations to non-central densities and quantiles are not available. In addition, its S-system formulation has not been derived. Here, approximations to cumulative probabilities, density, and quantiles based on the method of Knüsel and Bablok are derived and implemented in R codes. Furthermore, two alternate S-system forms are recast on the basis of techniques of Savageau and Voit (Savageau, M.A. and Voit, E.O., 1987, Recasting nonlinear differential equations as S-systems: A canonical nonlinear form. Mathematical Biosciences, 87, 83–115) as well as Chen (Chen, Z.-Y., 2003, Computing the distribution of the squared sample multiple correlation coefficient with S-Systems. Communications in Statistics—Simulation and Computation, 32(3), 873–898.) and Chen and Chou (Chen, Z.-Y. and Chou, Y.-C., 2000, Computing the noncentral beta distribution with S-system. Computational Statistics and Data Analysis, 33, 343–360.). Statistical densities, cumulative probabilities, quantiles can be evaluated by only one numerical solver power low analysis and simulation (PLAS). With the newly derived S-systems of non-central gamma, the specialized non-central χ² distributions are demonstrated under five cases in the same three situations studied by Rust and Voit. Both numerical values in pairs are almost equal. Based on these, nine cases in three similar situations are designed for demonstration and evaluation. In addition, exact values in finite significant digits are provided for comparison. Demonstrations are conducted by R package and PLAS solver in the same PC system. By doing these, very accurate and consistent numerical results are obtained by three methods in two groups. On the other hand, these three methods are performed competitively with respect to speed of computation. Numerical advantages of S-systems over the ad hoc approximation and related properties are also discussed.     

The log-exponentiated generalized gamma regression model for censored data

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 by means of two applications to real data.     

Consistent estimation of parameters and quantiles of the three-parameter gamma distribution based on Type-II right-censored data

In this paper, we propose a method of estimation of parameters and quantiles of the three-parameter gamma distribution based on Type-II right-censored data. In the proposed method, under mild conditions, the estimates always exist uniquely, and the estimators have consistency over the entire parameter space. Through Monte Carlo simulations, we further show that the proposed method performs well compared with another prominent method of estimation in terms of bias and root mean-squared error in small-sample situations. Finally, two real data sets are used for illustrating the proposed method.     

Linear estimation for the extended exponential power distribution

Tiao and Lund [The use of OLUMV estimators in inference robustness studies of the location parameter of a class of symmetric distributions. J Amer Statist Assoc. 1970;65(329):370–386] tabulated the coefficients of the best linear unbiased estimators (BLUEs) of location and scale for a particular family of symmetric distributions. This family was a reparameterization of the extended exponential power distribution (EEPD) with the shape parameter restricted to be greater than or equal to one. In this work, we consider the BLU estimation of the location and scale parameters of the EEPD when the shape parameter is one-third and one-half. We obtain closed-form expressions for the single and product moments of the order statistics when the shape parameter is in general in the form of a reciprocal of an integer. These expressions are then used to determine the BLUEs and the corresponding variances for complete samples of size 20 and less. We consider some other linear estimators of the location and scale parameters and then compare them with the BLUEs. Finally, we present a numerical example to illustrate the developed results.     

The cross-product ratio in bivariate lognormal and gamma distributions,with an application to non-randomized trials

Non-randomized trials can give a biased impression of the effectiveness of any intervention. We consider trials in which incidence rates are compared in two areas over two periods. Typically, one area receives an intervention, whereas the other does not. We outline and illustrate a method to estimate the bias in such trials under two different bivariate models. The illustrations use data in which no particular intervention is operating. The purpose is to illustrate the size of the bias that could be observed purely due to regression towards the mean (RTM). The illustrations show that the bias can be appreciably different from zero, and even when centred on zero, the variance of the bias can be large. We conclude that the results of non-randomized trials should be treated with caution, as interventions which show small effects could be explained as artefacts of RTM.     

Three strings of inequalities among six Bayes estimators

We discover three interesting strings of inequalities among six Bayes estimators, where for the parameter space (0, 1), (0, ∞), and ( ? ∞, ∞), each case has a string of inequalities. The three strings of inequalities only depend on the loss functions, and the inequalities are independent of the chosen models and the used priors provided the Bayes estimators exist. Therefore, they exist in a general setting which makes them quite interesting. Finally, the numerical simulations exemplify the two strings of inequalities defined on (0, 1) and (0, ∞), and that there does not exist a string of inequalities among the six smallest posterior expected losses.     

Generalized method of moments for an extended gamma process

In reliability theory, a widely used process to model the phenomena of the cumulative deterioration of a system over time is the standard gamma process (SGP). Based on several restrictions, such as a constant variance-to-mean ratio, this process is not always a suitable choice to describe the deterioration. A way to overcome these restrictions is to use an extended version of the gamma process introduced by Cinlar (1980), which is characterized by shape and scale functions. In this article, the aim is to propose statistical methods to estimate the unknown parameters of parametric forms of the shape and scale functions. We here develop two generalized methods of moments (Hansen 1982 Hansen, L. P. 1982. Large sample properties of generalized method of moments estimators. Econometrica 50 (4):1029–54.[Crossref], [Web of Science ®] , [Google Scholar]), based either on the moments or on the Laplace transform of an extended gamma process. Asymptotic properties are provided and a Wald-type test is derived, which allows to test SGPs against extended ones with a specific parametric shape function. Also, the performance of the proposed estimation methods is illustrated on simulated and real data.