Test for comparing complete expectations of life of two groups

Abstract

Complete expectation of life of an individual gives an intuitive and interesting perspective on the ageing process and is an important concept in the insurance sector for determination of premium. We propose a new test for testing equality of complete expectations of life of two groups/populations. Power of the new test is calculated through simulations and compared with the power of the tests given by Berger, Boos, and Guess (1988 Berger, R. L., D. D. Boos, and F. M. Guess. 1988. Tests and confidence sets for comparing two mean residual life functions. Biometrics 44 (1):103–15.[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]) and Aly (1997 Aly, E. E. A. A. 1997. Nonparametric tests for comparing two mean residual life functions. Lifetime Data Analysis 3 (4):353–66.[Crossref], [PubMed] , [Google Scholar]). It is observed that the proposed test statistic is more powerful than the competing tests for the cases considered in this paper. A real life illustration is included.     

Testing economic “genetic pleiotropy” for Box-Cox linear model

Abstract

Genetic pleiotropy occurs when a single gene influences two or more seemingly unrelated phenotypic traits. It is significant to detect pleiotropy and understand its causes. However, most current statistical methods to discover pleiotropy mainly test the null hypothesis that none of the traits is associated with a variant, which departures from the null to test just one associated trait or k associated traits. Schaid et al. (2016 Schaid, D. J., X. Tong, B. Larrabee, R. B. Kennedy, G. A. Poland, and J. P. Sinnwell. 2016. Statistical methods for testing genetic pleiotropy. Genetics 204 (2):483–97. doi:10.1534/genetics.116.189308.[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]) first proposed a sequential testing framework to analyze pleiotropy based on a linear model and a multivariate normal distribution. In this paper, we analyze the Economic pleiotropy which occurs when an economic action or policy influences two or more economic phenomena. In this paper, we extend the linear model to Box-Cox transformation model and proposed a new decision method. It improves the efficiency of hypothesis test and controls the Type I error. We then apply the method using economic data to multivariate sectoral employments in response to governmental expenditures and provide a quantitative assessment and some insights of different impacts from economic policy.     

Use of scrambled response for estimating mean of the sensitivity variable

Abstract

In this article, we propose new efficient and more generalized difference-cum-exponential type estimator and generalized-difference-cum-generalized exponential type estimators for estimating the mean of sensitivity variable using the auxiliary information. We also discuss theoretically that proposed generalized estimators are more efficient than Sousa et al. (2010 Sousa, R., J. Shabbir, P. C. Real, and S. Gupta. 2010. Ratio estimation of the mean of a sensitive variable in the presence of auxiliary information. Journal of Statistical Theory and Practice 4 (3):495–507.[Taylor & Francis Online] , [Google Scholar]), Gupta et al. (2012 Gupta, S., J. Shabbir, R. Sousa, and P. C. Real. 2012. Estimation of the mean of a sensitive variable in the presence of auxiliary information. Communications in Statistics-Theory and Methods 41:1–12.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) and Koyuncu, Gupta, and Sousa (2014 Koyuncu, N., S. Gupta, and R. Sousa. 2014. Exponential-type estimators of the mean of a sensitive variable in the presence of non sensitive auxiliary information. Communications in Statistics-Simulation and Computation 43 (7):1583–94. doi: 10.1080/03610918.2012.737492.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) estimators. Results from a real life application and simulation study are presented to demonstrate the performance of the proposed mean estimators in relation to some of the existing mean estimators.     

Variance estimation in adaptive cluster sampling

Abstract

In this paper, a class of variance estimator is proposed of a finite population variance under an adaptive cluster sampling design in the presence of information on an auxiliary variable. We obtain expressions for the mean square error and bias for the developed estimators and their performance is evaluated on a Poisson clustered process and a real data set. The simulation study evaluates the efficiency of the suggested estimators for an adaptive cluster sampling (ACS) design and the Isaki (1983 Isaki, C. T. 1983. Variance estimation using auxiliary information. Journal of the American Statistical Association 78 (381):117–23. doi: 10.1080/01621459.1983.10477939.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) estimator of the variance for SRSWOR over the sample variance for SRSWOR.     

Optimal nonparametric estimator of the area under ROC curve based on clustered data

Abstract

In diagnostic trials, clustered data are obtained when several subunits of the same patient are observed. Intracluster correlations need to be taken into account when analyzing such clustered data. A nonparametric method has been proposed by Obuchowski (1997 Obuchowski, N. A. 1997. Nonparametric analysis of clustered ROC curve data. Biometrics 53 (2):567–78.[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]) to estimate the Receiver Operating Characteristic curve area (AUC) for such clustered data. However, Obuchowski’s estimator is not efficient as it gives equal weight to all pairwise rankings within and between cluster. In this paper, we propose a more efficient nonparametric AUC estimator with two sets of optimal weights. Simulation results show that the loss of efficiency of Obuchowski’s estimator for a single AUC or the AUC difference can be substantial when there is a moderate intracluster test correlation and the cluster size is large. The efficiency gain of our weighted AUC estimator for a single AUC or the AUC difference is further illustrated using the data from a study of screening tests for neonatal hearing.     

The exponential distribution analog of the Grubbs–Weaver method

Abstract

Grubbs and Weaver (1947 Grubbs, F. E., and C. L. Weaver. 1947. The best unbiased estimate of population standard deviation based on group ranges. Journal of the American Statistical Association 42 (238):224–41. doi: 10.2307/2280652.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) suggest a minimum-variance unbiased estimator for the population standard deviation of a normal random variable, where a random sample is drawn and a weighted sum of the ranges of subsamples is calculated. The optimal choice involves using as many subsamples of size eight as possible. They verified their results numerically for samples of size up to 100, and conjectured that their “rule of eights” is valid for all sample sizes. Here we examine the analogous problem where the underlying distribution is exponential and find that a “rule of fours” yields optimality and prove the result rigorously.     

On high-dimensional tests for mutual independence based on Pearson’s correlation coefficient

Abstract

This paper investigates a class of statistics based on Pearson’s correlation coefficient for testing the mutual independence of a random vector in high dimensions. Two existing statistics, proposed by Schott (2005 Schott, J. R. 2005. Testing for complete independence in high dimensions. Biometrika 92 (4):951–6. doi:10.1093/biomet/92.4.951.[Crossref], [Web of Science ®] , [Google Scholar]) and Mao (2014 Mao, G. 2014. A new test of independence for high-dimensional data. Statistics & Probability Letters 93:14–8. doi:10.1016/j.spl.2014.05.024.[Crossref], [Web of Science ®] , [Google Scholar]) respectively, are special cases of the class. A generic testing theory for the class of statistics is developed, which clarifies under what conditions the class of statistics can be employed for the testing purpose. By virtue of the theory, three new tests are introduced, and related statistical properties are discussed. To examine our theoretical findings and check the performance of the new tests, simulation studies are applied. The simulation results justify the theoretical findings and show that the newly introduced tests perform well, as long as both the dimension and the sample size of the data are moderately large.     

Complete integral convergence for arrays of row-wise extended independent random variables under Sub-linear expectations

Abstract

In this paper, we establish the complete convergence and complete integral convergence for arrays of row-wise extended independent random variables under sub-linear expectation space with some conditions. At the same time we extend some complete convergence and complete integral convergence theorems from the classical probability space to the sub-linear expectation space. The results generalize corresponding results obtained by Wu et al. (2017 Zhang, L. X. 2016b. Exponential inequalities under the sub-linear expectations with applications to laws of the iterated logarithm. Science China Mathematics 59 (12):2503–26. doi: 10.1007/s11425-016-0079-1.[Crossref], [Web of Science ®] , [Google Scholar]).     

On a class exhibiting trend change in average failure rate

Abstract

Weak convergence and moment convergence issues are investigated for the New Better than Average Failure Rate (NBAFR) family (introduced by Loh (1984 Loh, W. Y. 1984. A new generalization of the class of NBU distributions. IEEE Transactions on Reliability R-33 :97–113[Crossref], [Web of Science ®] , [Google Scholar])). We explore the validity of these results in the context of a more general ageing class that we introduce. We prove some new properties of this class and derive its interrelationships with other non-monotonic ageing families. Reliability and moment bounds are obtained and an interesting characterization of exponentiality is proved. Special cases of our results lead to new theorems for the NBAFR class. Finally weak convergence and related issues are established for this class.     

Bayesian Estimation of Parameters and Comparison of Shared Gamma Frailty Models

In this article, we introduce shared gamma frailty models with three different baseline distributions namely, Weibull, generalized exponential and exponential power distributions. We develop 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. Also we apply these three models to a real life bivariate survival dataset of McGilchrist and Aisbett (1991 McGilchrist, C. A. and Aisbett, C. W. 1991. Regression with frailty in survival analysis. Biometrics, 47: 461–466. [Crossref], [PubMed], [Web of Science ®] , [Google Scholar]) related to kidney infection data and a better model is suggested for the data.     

The Johnson System of Frequency Curves—Historical,Graphical, and Limiting Perspectives

Abstract

The idea of transforming one random variate to another with a more convenient density has been developed in the first half of the 20th century. In his thesis, Norman L. Johnson (1917–2004) developed a pioneering system of transformations of the standard normal distribution which gained substantial popularity in the second half of the 20th century and beyond. In Johnson’s 1949 Johnson, N. L. (1949), “Systems of Frequency Curves Generated by Methods of Translation,” Biometrika, 36, 149–176.[Crossref], [PubMed], [Web of Science ®] , [Google Scholar] Biometrika paper entitled Systems of frequency curves generated by methods of translation, summarizing that thesis, one of his primary interests was the behavior of the shape of the probability density functions as their parameter values change. Herein, we attempt to further elucidate this behavior through a series of geometric expositions of that transformation process. In these expositions insight is obtained into the behavior of Johnson’s density functions, and their skewness and kurtosis, as they converge to their limiting distributions, a topic which received little attention.     

Is Each NPMLE of a Continuous Bivariate Distribution Function with Singly Right-Censored Data Really Inconsistent?

We consider non-parametric estimation of a continuous cdf of a random vector (X ₁, X ₂). With bivariate RC data, it is stated in van der Laan (1996 Van der Laan , M. J. ( 1996 ) Efficient estimation in the bivariate censoring model and repairing NPMLE . Ann. Statist. 24 : 596 – 627 .[Crossref], [Web of Science ®] , [Google Scholar], p. 598¹⁰, Ann. Statist.), Quale et al. (2006 Quale , C. M. , van der Laan , M. J. , Robins , J. R. ( 2006 ). Locally efficient estimation with bivariate right-censored data . JASA. 101 : 1076 – 1084 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar], JASA) etc. that “it is well known that the NPMLE for continuous data is inconsistent (Tsai et al. (1986 Tsai , W. Y. , Leurgans , S. , Crowley , J. ( 1986 ). Nonparametric estimation of a bivariate survival function in the presence of censoring . Ann. Statist. 14 : 1351 – 1365 .[Crossref], [Web of Science ®] , [Google Scholar])).” The claim is based on a result in Tsai et al. (1986 Tsai , W. Y. , Leurgans , S. , Crowley , J. ( 1986 ). Nonparametric estimation of a bivariate survival function in the presence of censoring . Ann. Statist. 14 : 1351 – 1365 .[Crossref], [Web of Science ®] , [Google Scholar], p.1352, Ann. Statist.) that if X ₁ is right censored but not X ₂, then common ways for defining one NPMLE lead to inconsistency. If X ₁ is right censored and X ₂ is type I right-censored (which includes the case in Tsai et al.), we present a consistent NPMLE. The result corrects a common misinterpretation of Tsai's example (Tsai et al., 1986 Tsai , W. Y. , Leurgans , S. , Crowley , J. ( 1986 ). Nonparametric estimation of a bivariate survival function in the presence of censoring . Ann. Statist. 14 : 1351 – 1365 .[Crossref], [Web of Science ®] , [Google Scholar], Ann. Statist.).     

Log-epsilon-skew normal: A generalization of the log-normal distribution

Abstract

The log-normal distribution is widely used to model non-negative data in many areas of applied research. In this paper, we introduce and study a family of distributions with non-negative reals as support and termed the log-epsilon-skew normal (LESN) which includes the log-normal distributions as a special case. It is related to the epsilon-skew normal developed in Mudholkar and Hutson (2000 Mudholkar, G. S., and A. D. Hutson. 2000. The epsilon-skew-normal distribution for analyzing near-normal data. Journal of Statistical Planning and Inference 83 (2):291–309. doi:10.1016/S0378-3758(99)00096-8.[Crossref], [Web of Science ®] , [Google Scholar]) the way the log-normal is related to the normal distribution. We study its main properties, hazard function, moments, skewness and kurtosis coefficients, and discuss maximum likelihood estimation of model parameters. We summarize the results of a simulation study to examine the behavior of the maximum likelihood estimates, and we illustrate the maximum likelihood estimation of the LESN distribution parameters to two real world data sets.     

Additive Mixed Effect Model for Clustered Doubly Censored Data

Cai and Zeng (2011 Cai, J. and Zeng, D. 2011. Additive mixed effect model for clustered failure time data. Biometrics, 67(4): 1340–1351. [Crossref], [PubMed] , [Google Scholar]) proposed an additive mixed effect model to analyze clustered right-censored data. In this article, we demonstrate that the approach of Cai and Zeng (2011 Cai, J. and Zeng, D. 2011. Additive mixed effect model for clustered failure time data. Biometrics, 67(4): 1340–1351. [Crossref], [PubMed] , [Google Scholar]) can be extended to clustered doubly censored data. Furthermore, when both left- and right-censoring variables are always observed, we propose alternative estimators using the approach of Cai and Cheng (2004 Cai, T. and Cheng, S. C. 2004. Semiparametric regression analysis for doubly censored data. Biometrika, 91: 277–290. [Crossref], [Web of Science ®] , [Google Scholar]). A simulation study is conducted to investigate the performance of the proposed estimators.     

Second-Order Covariance Matrix Formula for Heteroskedastic Generalized Linear Models

For the first time, we provide a matrix formula for second-order covariances of maximum likelihood estimates in heteroskedastic generalized linear models, thus generalizing the results of Cordeiro (2004 Cordeiro , G. M. ( 2004 ). Second-order covariance matrix of maximum likelihood estimates in generalized linear models . Statist. Probab. Lett. 66 : 153 – 160 .[Crossref], [Web of Science ®] , [Google Scholar]) and Cordeiro et al. (2006 Cordeiro , G. M. , Barroso , L. P. , Botter , D. A. (2006). Covariance matrix formula for generalized linear models with unknown dispersion. Commun. Statist. Theor. Meth. 35:113–120.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) related to the generalized linear models with known and unknown dispersion parameter, respectively. The covariance matrix formula does not involve cumulants of log-likelihood derivatives and can be easily obtained using simple matrix operations. We apply our main result to a simple model. Some simulations show that the second-order covariances can be quite pronounced in small to moderate samples. The usual covariances of the maximum likelihood estimates can be corrected by these second-order covariances.     

Elementary renewal theorems for widely dependent random variables with applications to precise large deviations

Abstract

On the basis of Wang and Cheng (J. Math. Anal. Appl. 384 (2011) 597–606), this paper further investigates elementary renewal theorems for counting processes generated by random walks with widely orthant dependent increments. The obtained results improve the corresponding ones of the above-mentioned paper mainly in the sense of weakening the moment conditions on the positive parts of the increments. Meanwhile, a revised version of strong law of large numbers for random walks with widely orthant dependent increments is established, which improves Theorem 1.4 of Wang and Cheng (2011 Wang, Y., and D. Cheng. 2011. Basic renewal theorems for a random walk with widely dependent increments and their applications. Journal of Mathematical Analysis and Applications 384 (2):597–606. doi:10.1016/j.jmaa.2011.06.010.[Crossref], [Web of Science ®] , [Google Scholar]) by enlarging the regions of dominating coefficients. Finally, by using the above results, some precise large deviation results for a nonstandard renewal risk model are established, in which the innovations are widely orthant dependent random variables with common heavy tails, and the inter-arrival times are also widely orthant dependent.     

On Estimation of the Bivariate Poisson INAR Process

In a recent article, Pedeli and Karlis (2010 Pedeli, X. and Karlis, D. 2010. A Bivariate INAR(1) Process with Application. Statistical Modelling: An international Journal, 11: 325–349. [Crossref], [Web of Science ®] , [Google Scholar]) examined the extension of the classical Integer–valued Autoregressive (INAR) model to the bivariate case. In the present article, we examine estimation methods for the case of bivariate Poisson innovations. This is a simple extension of the classical INAR model allowing for two discrete valued time series to be correlated. Properties of different estimators are given. We also compare their properties via a small simulation experiment. Extensions to incorporate covariate information is discussed. A real data application is also provided.     

Quantile estimation for a progressively censored exponential distribution

Abstract

In this paper, we consider the problem of estimating the quantile of a two-parameter exponential distribution with respect to an arbitrary strictly convex loss function under progressive type II censoring. Inadmissibility of the best affine equivariant (BAE) estimator is established through a conditional risk analysis. In particular we provide dominance results for quadratic, linex and absolute value loss functions. Further, a class of dominating estimators is derived using the IERD (integral expression of risk difference) approach of Kubokawa (1994 Kubokawa, T. 1994. A unified approach to improving equivariant estimators. The Annals of Statistics 22 (1):290–9. doi:10.1214/aos/1176325369.[Crossref], [Web of Science ®] , [Google Scholar]). In sequel the generalized Bayes estimator is shown to improve the BAE estimator.     

A side sensitive modified group runs control chart to detect shifts in the process mean

Gadre and Rattihalli [5 Gadre, M. P. and Rattihalli, R. N. 2006. Modified group runs control charts to detect increases in fraction non-conforming and shifts in the process mean. Commun. Stat. Simul. Comput., 35: 225–240. [Taylor & Francis Online], [Web of Science ®] [Google Scholar]] have introduced the Modified Group Runs (MGR) control chart to identify the increases in fraction non-conforming and to detect shifts in the process mean. The MGR chart reduces the out-of-control average time-to-signal (ATS), as compared with most of the well-known control charts. In this article, we develop the Side Sensitive Modified Group Runs (SSMGR) chart to detect shifts in the process mean. With the help of numerical examples, it is illustrated that the SSMGR chart performs better than the Shewhart's X¯ chart, the synthetic chart [12 Wu, Z. and Spedding, T. A. 2000. A synthetic control chart for detecting small shifts in the process mean. J. Qual. Technol., 32: 32–38. [Taylor & Francis Online], [Web of Science ®] [Google Scholar]], the Group Runs chart [4 Gadre, M. P. and Rattihalli, R. N. 2004. A group runs control chart for detecting shifts in the process mean. Econ. Qual. Control, 19: 29–43. [Crossref] [Google Scholar]], the Side Sensitive Group Runs chart [6 Gadre, M. P. and Rattihalli, R. N. 2007. A side sensitive group runs control chart for detecting shifts in the process mean. Stat. Methods Appl., 16: 27–37. [Crossref], [Web of Science ®] [Google Scholar]], as well as the MGR chart [5 Gadre, M. P. and Rattihalli, R. N. 2006. Modified group runs control charts to detect increases in fraction non-conforming and shifts in the process mean. Commun. Stat. Simul. Comput., 35: 225–240. [Taylor & Francis Online], [Web of Science ®] [Google Scholar]]. In some situations it is also superior to the Cumulative Sum chart p9 Page, E. S. 1954. Continuous inspection schemes. Biometrika, 41: 100–114. [Crossref], [Web of Science ®] [Google Scholar]] and the exponentially weighed moving average chart [10 Roberts, S. W. 1959. Control chart tests based on geometric moving averages. Technometrics, 1: 239–250. [Taylor & Francis Online] [Google Scholar]]. In the steady state also, its performance is better than the above charts.