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.     

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.     

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.     

Survival analysis for a new compounded bivariate failure time distribution in shock and competing risk models via an EM algorithm

Abstract

In this paper, two bivariate models based on the proposed methods of Marshall and Olkin are introduced. In the first model, the new bivariate distribution is presented based on the proposed method of Marshall and Olkin (1967 Marshall, A. W., and I. Olkin. 1967. A multivariate exponential distribution. Journal of the American Statistical Association 62 (317):30–44. doi: 10.2307/2282907.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) which has natural interpretations, and it can be applied in fatal shock models or in competing risks models. In the second model, the proposed method of Marshall and Olkin (1997 Marshall, A. W., and I. Olkin. 1997. A new method of adding a parameter to a family of distributions with application to the exponential and weibull families. Biometrika 84 (3):641–52. doi: 10.1093/biomet/84.3.641.[Crossref], [Web of Science ®] , [Google Scholar]) is generalized to bivariate case and a new bivariate distribution is introduced. We call these new distributions as the bivariate Gompertz (BGP) distribution and bivariate Gompertz-geometric (BGPG) distribution, respectively. Moreover, the BGP model can be obtained as a special case of the BGPG model. Then, we present various properties of the new bivariate models. In this regard, the joint and conditional density functions, the joint cumulative distribution function can be obtained in compact forms. Also, the aging properties and the bivariate hazard gradient are discussed. This model has five unknown parameters and the maximum likelihood estimators cannot be obtained in explicit form. We propose to use the EM algorithm to compute the maximum likelihood estimators of the unknown parameters, and it is computationally quite tractable. Also, Monte Carlo simulations are performed to investigate the effectiveness of the proposed algorithm. Finally, we analyze three real data sets for illustrative purposes.     

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.     

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.     

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.     

Adaptive estimation for varying coefficient models with nonstationary covariates

In this paper, the adaptive estimation for varying coefficient models proposed by Chen, Wang, and Yao (2015 Chen, Y., Q. Wang, and W. Yao. 2015. Adaptive estimation for varying coefficient models. Journal of Multivariate Analysis 137:17–31.[Crossref], [Web of Science ®] , [Google Scholar]) is extended to allowing for nonstationary covariates. The asymptotic properties of the estimator are obtained, showing different convergence rates for the integrated covariates and stationary covariates. The nonparametric estimator of the functional coefficient with integrated covariates has a faster convergence rate than the estimator with stationary covariates, and its asymptotic distribution is mixed normal. Moreover, the adaptive estimation is more efficient than the least square estimation for non normal errors. A simulation study is conducted to illustrate our theoretical results.     

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 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.     

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.     

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.     

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.     

A nonparametric specification test for the volatility functions of diffusion processes

This paper develops a new test for the parametric volatility function of a diffusion model based on nonparametric estimation techniques. The proposed test imposes no restriction on the functional form of the drift function and has an asymptotically standard normal distribution under the null hypothesis of correct specification. It is consistent against any fixed alternatives and has nontrivial asymptotic power against a class of local alternatives with proper rates. Monte Carlo simulations show that the test performs well in finite samples and generally has better power performance than the nonparametric test of Li (2007 Li, F. (2007). Testing the parametric specification of the diffusion function in a diffusion process. Econometric Theory 23(2):221–250.[Crossref], [Web of Science ®] , [Google Scholar]) and the stochastic process-based tests of Dette and Podolskij (2008 Dette, H., Podolskij, M. (2008). Testing the parametric form of the volatility in continuous time diffusion models–a stochastic process approach. Journal of Econometrics 143(1):56–73.[Crossref], [Web of Science ®] , [Google Scholar]). When applying the test to high frequency data of EUR/USD exchange rate, the empirical results show that the commonly used volatility functions fit more poorly when the data frequency becomes higher, and the general volatility functions fit relatively better than the constant volatility function.     

Estimation of nonseparable models with censored dependent variables and endogenous regressors

In this article we develop a nonparametric estimator for the local average response of a censored dependent variable to endogenous regressors in a nonseparable model where the unobservable error term is not restricted to be scalar and where the nonseparable function need not be monotone in the unobservables. We formalize the identification argument put forward in Altonji, Ichimura, and Otsu (2012 Altonji, J. G., Ichimura, H., Otsu, T. (2012). Estimating derivatives in nonseparable models with limited dependent variables. Econometrica 80:1701–1719.[Crossref], [Web of Science ®] , [Google Scholar]), construct a nonparametric estimator, characterize its asymptotic property, and conduct a Monte Carlo investigation to study its small sample properties. Identification is constructive and is achieved through a control function approach. We show that the estimator is consistent and asymptotically normally distributed. The Monte Carlo results are encouraging.     

Test for high dimensional regression coefficients of partially linear models

Abstract

Partially linear models attract much attention to investigate the association between predictors and the response variable when the dependency on some predictors may be nonlinear. However, the hypothesis test for significance of predictors is still challenging, especially when the number of predictors is larger than sample size. In this paper, we reconsider the test procedure of Zhong and Chen (2011 Zhong, P., and S. Chen. 2011. Tests for high-dimensional regression coefficients with factorial designs. Journal of the American Statistical Association 106 (493):260–74. doi:10.1198/jasa.2011.tm10284.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) when regression models have nonlinear components, and propose a generalized U-statistic for testing the linear components of the high dimensional partially linear models. The asymptotic properties of test statistic are obtained under null and alternative hypotheses, where the effect of nonlinear components should be considered and thus is different from that in linear models. Through simulation studies, we demonstrate good finite-sample performance of the proposed test in comparison with the existing methods. The practical utility of our proposed method is illustrated by a real data example.     

Bias corrected empirical Bayes confidence intervals for the selected mean

Empirical Bayes (EB) methods are very useful for post selection inference. Following Datta et al. (2002 Datta, G. S., M. Ghosh, D. D. Smith, and P. Lahiri. 2002. On an asymptotic theory of conditional and unconditional coverage probabilities of empirical Bayes confidence intervals. Scandinavian Journal of Statistics 29:139–52.[Crossref], [Web of Science ®] , [Google Scholar]), we construct EB confidence intervals for the selected population mean. The EB intervals are adjusted to achieve the target coverage probabilities asymptotically up to the second order. Both unconditional coverage probabilities of EB intervals and corresponding probabilities conditional on ancillary statistics are found.     

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.     

Nonparametric test for stratum effects in the cox model with interval-censored data

The testing of the stratum effects in the Cox model is an important and commonly asked question in medical research as well as in many other fields. In this paper, we will discuss the problem where one observes interval-censored failure time data and generalize the procedure given in Sun and Yang (2000 Sun, J., and I. Yang. 2000. Nonparametric test for stratum effects in the cox model. Lifetime Data Analysis 6:321–30.[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]) for right-censored data. The asymptotic distribution of the new test statistic is established and the simulation study conducted for the evaluation of the finite sample properties of the method suggests that the generalized procedure seems to work well for practical situations. An application is provided.     

Letter to the Editor

The aim of this letter to acknowledge of priority on calibration estimation. There are numerous studies on calibration estimation in literature. The studies on calibration estimation are reviewed and it is found out that an existing calibration estimator is reprocessed in the recent paper published by Nidhi et al. (2007 Nidhi, B. V. S. Sisodia, Subedar Singh, and Sanjay K. Singh. 2017. Calibration approach estimation of the mean in stratified sampling and stratified double sampling. Commun.Statist.Theor.Meth. 46 (10):4932–4942.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]).