Linear Subspace Spanned by Principal Points of a Mixture of Spherically Symmetric Distributions

For each positive integer k, a set of k-principal points of a distribution is the set of k points that optimally represent the distribution in terms of mean squared distance. However, explicit form of k-principal points is often difficult to obtain. Hence a theorem established by Tarpey et al. (1995 Tarpey , T. , Li , L. , Flury , B. D. ( 1995 ). Principal points and self-consistent points of elliptical distributions . Ann. Statist. 23 : 102 – 112 .[Crossref], [Web of Science ®] , [Google Scholar]) has been influential in the literature, which states that when the distribution is elliptically symmetric, any set of k-principal points is in the linear subspace spanned by some principal eigenvectors of the covariance matrix. This theorem is called a “principal subspace theorem”. Recently, Yamamoto and Shinozaki (2000b Yamamoto , W. , Shinozaki , N. ( 2000b ). Two principal points for multivariate location mixtures of spherically symmetric distributions . J. Japan Statist. Soc. 30 : 53 – 63 .[Crossref] , [Google Scholar]) derived a principal subspace theorem for 2-principal points of a location mixture of spherically symmetric distributions. In their article, the ratio of mixture was set to be equal. This article derives a further result by considering a location mixture with unequal mixture ratio.     

Relations for Moments of Progressively Type-II Censored Order Statistics from Log-Logistic Distribution with Applications to Inference

In this article, we establish several recurrence relations for the single and product moments of progressively Type-II right censored order statistics from a log-logistic distribution. The use of these relations in a systematic recursive manner would enable the computation of all the means, variances and covariances of progressively Type-II right censored order statistics from the log-logistic distribution for all sample sizes n, effective sample sizes m, and all progressive censoring schemes (R ₁,…, R _m ). The results established here generalize the corresponding results for the usual order statistics due to Balakrishnan and Malik (1987 Balakrishnan , N. , Malik , H. J. ( 1987 ). Moments of order statistics from truncated log-logistic distribution . J. Statist. Plann. Infer. 17 : 251 – 267 .[Crossref], [Web of Science ®] , [Google Scholar]) and Balakrishnan et al. (1987 Balakrishnan , N. , Malik , H. J. , Puthenpura , S. ( 1987 ). Best linear unbiased estimation of location and scale parameters of the log-logistic distribution . Commun. Statist. Theor. Meth. 16 : 3477 – 3495 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]). The moments so determined are then utilized to derive best linear unbiased estimators for the scale- and location-scale log-logistic distributions. A comparison of these estimates with the maximum likelihood estimates is made through Monte Carlo simulation. The best linear unbiased predictors of progressively censored failure times is then discussed briefly. Finally, a numerical example is presented to illustrate all the methods of inference developed here.     

The Poisson Generalized Linear Failure Rate Model

We formulate a new cure rate survival model by assuming that the number of competing causes of the event of interest has the Poisson distribution, and the time to this event has the generalized linear failure rate distribution. A new distribution to analyze lifetime data is defined from the proposed cure rate model, and its quantile function as well as a general expansion for the moments is derived. We estimate the parameters of the model with cure rate in the presence of covariates for censored observations using maximum likelihood and derive the observed information matrix. We obtain the appropriate matrices for assessing local influence on the parameter estimates under different perturbation schemes and present some ways to perform global influence analysis. The usefulness of the proposed cure rate survival model is illustrated in an application to real data.     

Interval Estimation of the Stress-Strength Reliability with Independent Normal Random Variables

This article develops a procedure to obtain highly accurate confidence interval estimates for the stress-strength reliability R = P(X > Y) where X and Y are data from independent normal distributions of unknown means and variances. Our method is based on third-order likelihood analysis and is compared to the conventional first-order likelihood ratio procedure as well as the approximate methods of Reiser and Guttman (1986 Reiser, B., Guttman, I. (1986). Statistical inference for Pr(Y < X): the normal case. Technometrics 28: 253257.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) and Guo and Krishnamoorthy (2004 Guo, H., Krishnamoorthy, K. (2004). New approximate inferential methods for the reliability parameter in a stress-strength model: the normal case. Commun. Statist. Theor. Meth. 33: 17151731.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]). The use of our proposed method is illustrated by an empirical example and its superior accuracy in terms of coverage probability and error rate are examined through Monte Carlo simulation studies.     

Singular Ridge Regression With Stochastic Constraints

This article considers the estimation of the restricted ridge regression parameter in singular models. The problem is commenced with considering elliptically contoured equality constrained and then followed by proposing the preliminary test estimator. Along with proposing some important properties of this estimator, a real example satisfying the elliptical assumption is also given to bring the problem into a noticeable issue.     

Asymptotic Expansions for i.i.d. Sums Via Lower-order Convolutions

In this article, we introduce new asymptotic expansions for probability functions of sums of independent and identically distributed random variables. Results are obtained by efficiently employing information provided by lower-order convolutions. In comparison with Edgeworth-type theorems, advantages include improved asymptotic results in the case of symmetric random variables and ease of computation of main error terms and asymptotic crossing points. The first-order estimate can perform quite well against the corresponding renormalized saddlepoint approximation and, pointwise, requires evaluation of only a single convolution integral. While the new expansions are fairly straightforward, the implications are fortuitous and may spur further related work.     

Upper Bound for p-Value of the Test of Multivariate Normal Ordered Mean Vectors Against all Alternatives

Consider the problem of testing the isotonic of several p-variate normal mean vectors against all alternatives. It is difficult to compute the exact p-value for this problem of testing with the classical method when the covariance matrices are completely unknown. In the present paper, a test statistic is proposed for this problem of testing. A reformulation of the test statistic is given based on the orthogonal projections on the closed convex cones and then the upper bound for p-value of the test statistic is computed.     

Modeling bivariate survival data using shared inverse Gaussian frailty model

ABSTRACT

The shared frailty models are often used to model heterogeneity in survival analysis. The most common shared frailty model is a model in which hazard function is a product of a random factor (frailty) and the baseline hazard function which is common to all individuals. There are certain assumptions about the baseline distribution and the distribution of frailty. In this paper, we consider inverse Gaussian distribution as frailty distribution and three different baseline distributions, namely the generalized Rayleigh, the weighted exponential, and the extended Weibull distributions. With these three baseline distributions, we propose three different inverse Gaussian shared frailty models. We also compare these models with the models where the above-mentioned distributions are considered without frailty. We develop the Bayesian estimation procedure using Markov Chain Monte Carlo (MCMC) technique to estimate the parameters involved in these models. We present a simulation study to compare the true values of the parameters with the estimated values. A search of the literature suggests that currently no work has been done for these three baseline distributions with a shared inverse Gaussian frailty so far. We also apply these three models by using a real-life bivariate survival data set 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 using the Bayesian model selection criteria.     

Tail behavior of the generalized Maxwell distribution

Abstract

In this article, we investigate the tail propertiesof the generalized Maxwell distribution and gain an asymptoticbehavior of its Mills-type ratio. Meanwhile, we show two applications. The first application thinks about the asymptotic property of the ratio of density functions and the ratio of the tails of the generalized Maxwell and classical Maxwell distributions. Another application obtains the asymptotic distribution of the partial maximum of an independent and identically distributed sequence from the distribution.