An efficient new “Partial” randomized response model

ABSTRACT

In this article, a new “Partial” randomized response model has been proposed. Its properties are studied both theoretically and empirically. The proposed model is proved to be more efficient than the randomized response models studied by Eichhorn and Hayre (1983 Eichhorn, B.H., Hayre, L.S. (1983). Scrambled randomized response methods for obtaining sensitive quantitative dada. J. Stat. Plann. Inference. 7:307–316.[Crossref], [Web of Science ®] , [Google Scholar]) and the “Partial” randomized response model.     

POOLED SLICING METHODS VERSUS SLICING METHODS

We consider the semiparametric regression model introduced by [1] Duan, N. and Li, K. C. 1991. Slicing regression: a link-free regression method. The Annals of Statistics, 19: 505–530. [Crossref], [Web of Science ®] , [Google Scholar]. The dependent variable y is linked to the index x′ β through an unknown link function. [1] Duan, N. and Li, K. C. 1991. Slicing regression: a link-free regression method. The Annals of Statistics, 19: 505–530. [Crossref], [Web of Science ®] , [Google Scholar] and [2] Li, K. C. 1991. Sliced inverse regression for dimension reduction, with discussions. Journal of the American Statistical Association, 86: 316–342. [Taylor & Francis Online], [Web of Science ®] , [Google Scholar] present Slicing methods (the Sliced Inverse Regression methods SIR-I, SIR-II and SIR_α) in order to estimate the direction of the unknown slope parameter β. These methods are computationally simple and fast but depend on the choice of an arbitrary slicing fixed by the user. When the sample size is small, the number and the position of slices have an influence on the estimated direction. In this paper, we suggest to use the corresponding Pooled Slicing methods: PSIR-I (proposed by [3] Aragon, Y. and Saracco, J. 1997. Sliced Inverse Regression (SIR): an appraisal of small sample alternatives to slicing. Computational Statistics, 12: 109–130. [Web of Science ®] , [Google Scholar]), PSIR-II and PSIR_α. These methods combine the results from a number of slicings. We compare the sample behaviour of Slicing and Pooled Slicing methods on simulations. We also propose a practical choice of α in SIR_α and PSIR_α methods.     

Complete convergence and complete moment convergence for arrays of rowwise negatively superadditive dependent random variables

In this paper, some complete convergence and complete moment convergence results for arrays of rowwise negatively superadditive dependent (NSD, in short) random variables are studied. The obtained theorems not only extend the result of Gan and Chen (2007 Gan, S. X., and P. Y. Chen. 2007. On the limiting behavior of the maximum partial sums for arrays of rowwise NA random variables. Acta Mathematica Scientia. Series B 27 (2):283–90.[Crossref], [Web of Science ®] , [Google Scholar]) to the case of NSD random variables, but also improve them.     

Equivalent Conditions of Complete Moment Convergence of Weighted Sums for ϕ-Mixing Sequence of Random Variables

In this article, the complete moment convergence of weighted sums for ?-mixing sequence of random variables is investigated. By applying moment inequality and truncation methods, the equivalent conditions of complete moment convergence of weighted sums for ?-mixing sequence of random variables are established. These results promote and improve the corresponding results obtained by Li et al. (1995 Li, D.L., Rao, M.B., Jiang, T.F., Wang, X.C. (1995). Complete convergence and almost sure convergence of weighted sums of random variables. J. Theoret. Probab. 8:49–76.[Crossref], [Web of Science ®] , [Google Scholar]) and Gut (1993 Gut, A. (1993). Complete convergence and Cesàro summation for i.i.d. random variables. Probab. Theory Related Fields 97:169–178.[Crossref], [Web of Science ®] , [Google Scholar]) from i.i.d. to ?-mixing setting. Moreover, we obtain the complete moment convergence of moving average processes based on ?-mixing random variables, which extends the result of Kim et al. (2008 Kim, T.S., Ko, M.H. (2008). Complete moment convergence of moving average processes under dependence assumptions. Statist. Probab. Lett. 78:839–846.[Crossref], [Web of Science ®] , [Google Scholar]) in the sense that it does not require a specific mixing rate.     

Quantile-based Chernoff distance for truncated random variables

Several probability distributions such as power-Pareto distribution (see Gilchrist 2000 Gilchrist, W. 2000. Statistical modelling with quantile functions. Boca Raton, FL: Chapman and Hall/CRC.[Crossref] , [Google Scholar] and Hankin and Lee 2006 Hankin, R. K. S., and A. Lee. 2006. A new family of non-negative distributions. Australian and New Zealand Journal of Statistics 48:67–78.[Crossref], [Web of Science ®] , [Google Scholar]), various forms of lambda distributions (see Ramberg and Schmeiser 1974 Ramberg, J. S., and B. W. Schmeiser. 1974. An appropriate method for generating asymmetric random variables. Communications of the ACM 17:78–82.[Crossref], [Web of Science ®] , [Google Scholar] and Freimer et al. 1988 Freimer, M., S. Mudholkar, G. Kollia, and C. T. Lin. 1988. A study of the generalized lambda family. Communications in Statistics - Theory and Methods 17:3547–67.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]), Govindarajulu distribution (see Nair, Sankaran, and Vineshkumar 2012 Nair, U. N., P. G. Sankaran, and B. Vineshkumar. 2012. The Govindarajulu distribution: some properties and applications. Communications in Statistics—Theory and Methods 41:4391–406.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]), etc., do not have manageable distribution functions, though they have tractable quantile functions. Hence, analytical study of the properties of Chernoff distance of two random variables associated with these distributions via traditional distribution function-based tool becomes difficult. To make this simple, in this paper, we introduce quantile-based Chernoff distance for (left or right) truncated random variables and study its various properties. Some useful bounds as well as characterization results are obtained.     

An adaptive multivariate CUSUM control chart for signaling a range of location shifts

ABSTRACT

In this work, we proposed an adaptive multivariate cumulative sum (CUSUM) statistical process control chart for signaling a range of location shifts. This method was based on the multivariate CUSUM control chart proposed by Pignatiello and Runger (1990 Pignatiello, J.J., Runger, G.C. (1990). Comparisons of multivariate CUSUM charts. J. Qual. Technol. 22(3):173–186.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]), but we adopted the adaptive approach similar to that discussed by Dai et al. (2011 Dai, Y., Luo, Y., Li, Z., Wang, Z. (2011). A new adaptive CUSUM control chart for detecting the multivariate process mean. Qual. Reliab. Eng. Int. 27(7):877–884.[Crossref], [Web of Science ®] , [Google Scholar]), which was based on a different CUSUM method introduced by Crosier (1988 Crosier, R.B. (1988). Multivariate generalizations of cumulative sum quality-control schemes. Technometrics 30(3):291–303.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]). The reference value in this proposed procedure was changed adaptively in each run, with the current mean shift estimated by exponentially weighted moving average (EWMA) statistic. By specifying the minimal magnitude of the mean shift, our proposed control chart achieved a good overall performance for detecting a range of shifts rather than a single value. We compared our adaptive multivariate CUSUM method with that of Dai et al. (2001 Dai, Y., Luo, Y., Li, Z., Wang, Z. (2011). A new adaptive CUSUM control chart for detecting the multivariate process mean. Qual. Reliab. Eng. Int. 27(7):877–884.[Crossref], [Web of Science ®] , [Google Scholar]) and the non adaptive versions of these two methods, by evaluating both the steady state and zero state average run length (ARL) values. The detection efficiency of our method showed improvements over the comparative methods when the location shift is unknown but falls within an expected range.     

A quantile-based study on ageing intensity function

Jiang, Ji, and Xiao (2003 Jiang, R., P. Ji, and X. Xiao. 2003. Aging property of unimodal failure rate models. Reliability Engineering and System Safety 79(1):113–6.[Crossref], [Web of Science ®] , [Google Scholar]) has introduced a quantitative measure known as the ageing intensity function for evaluating the ageing properties of a component/system. In recent years, there has been a great interest on the study of quantile function, an equivalent alternative to the distribution function approach. Unlike the distribution function approach, the quantile method possess some unique properties (see Gilchrist 2000 Gilchrist, W. 2000. Statistical modelling with quantile functions. Boca Raton, Florida: Chapman and Hall/CRC.[Crossref] , [Google Scholar], Nair, Sankaran, and Balakrishnan 2013 Nair N. U., P. G. Sankaran, N. Balakrishnan. 2013. Quantile-based reliability concepts. In: Quantile-Based Reliability Analysis. Statistics for Industry and Technology. New York, NY: Birkhäuser.[Crossref] , [Google Scholar]). Motivated with this, in the present paper we introduce a quantile-based ageing intensity function and study its various ageing properties. We also study some stochastic comparison of random variables based on the proposed measure.     

A study of uniformity pattern for extended designs

The objective of this paper is to investigate the issue of projection discrepancy for extended U-type or nearly U-type extended designs along the line of Fang and Qin (2005 Fang, K. T., and H. Qin. 2005. Uniformity pattern and related criteria for two-level factorials. Science China Series A 48:1–11.[Crossref] , [Google Scholar]) based on the centered L₂-discrepancy proposed in Hickernell (1998 Hickernell, F. J. 1998. A generalized discrepancy and quadrature error bound. Mathematics of Computation 67:299–322.[Crossref], [Web of Science ®] , [Google Scholar]). Extended designs are obtained through augmenting optimally few runs (or points) to an optimal U-type design. Lower bounds to projection discrepancy with reference to the centered L₂-discrepancy of extended designs have been obtained. Some illustrative examples are also provided.     

Bayes minimax ridge regression estimators

The problem of estimating of the vector β of the linear regression model y = Aβ + ? with ? ～ N_p(0, σ²I_p) under quadratic loss function is considered when common variance σ² is unknown. We first find a class of minimax estimators for this problem which extends a class given by Maruyama and Strawderman (2005 Maruyama, Y., and W. E. Strawderman. 2005. A new class of generalized Bayes minimax ridge regression estimators. Annals of Statistics 33:1753–70.[Crossref], [Web of Science ®] , [Google Scholar]) and using these estimators, we obtain a large class of (proper and generalized) Bayes minimax estimators and show that the result of Maruyama and Strawderman (2005 Maruyama, Y., and W. E. Strawderman. 2005. A new class of generalized Bayes minimax ridge regression estimators. Annals of Statistics 33:1753–70.[Crossref], [Web of Science ®] , [Google Scholar]) is a special case of our result. We also show that under certain conditions, these generalized Bayes minimax estimators have greater numerical stability (i.e., smaller condition number) than the least-squares estimator.     

A stratified Tracy and Osahan's two-stage randomized response model

ABSTRACT

This paper proposes an alternative two-stage stratified randomized response model based on Tracy and Osahan (1999 Tracy, D.S., Osahan, S.S. (1999). An improved randomized response technique. Pak. J. Stat. 15(1):1–6. [Google Scholar]) model that has an optimal allocation and large gain in precision. It is also shown that the proposed model is more efficient than Kim and Warde (2004 Kim, J., Warde, W. (2004). A stratified Warner randomized response model. J. Stat. Plan. Infer. 120:155–165.[Crossref], [Web of Science ®] , [Google Scholar]) and Kim and Elam (2005 Kim, J.M., Elam, M.E. (2005). A two-stage stratified Warner's randomized response model using optimal allocation. Metrika 61:1–7.[Crossref], [Web of Science ®] , [Google Scholar]) stratified randomized response models under the conditions presented in both the cases of completely truthful reporting and that of not completely truthful reporting by the respondents. Numerical illustrations and graphs are also given in support of the present study.     

Generalized linear regression models incorporating original outcome distributions

ABSTRACT

In this article, we propose an approach for incorporating continuous and discrete original outcome distributions into the usual exponential family regression models. The new approach is an extension of the works of Suissa (1991 Suissa, S. (1991). Binary methods for continuous outcomes: A parametric alternative. J. Clin. Epidemiol. 44:241–248.[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]) and Suissa and Blais (1995 Suissa, S., Blais, L. (1995). Binary regression with continuous outcomes. Stat. Med. 14:247–255.[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]), which present methods to estimate the risk of an event defined in a sample subspace of an original continuous outcome variable. Simulation studies are presented in order to illustrate the performance of the developed methodology. Real data sets are analyzed by using the proposed models.     

On simulation of normal records

ABSTRACT

In the present paper, we discuss algorithms of record generation when records are taken from a normal population. We propose three new generation algorithms, compare their efficiency and find the most efficient algorithm (Algorithm 2.1). We then compare these algorithms with known generation algorithms presented in the work of Balakrishnan, So, and Zhu (2016 Balakrishnan, N., H. Y. So, and X. J. Zhu. 2016. On Box-Muller transformation and simulation of normal record data. Communication in Statistics – Simulation and Computations 45(10):3670–82.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]).     

EVALUATION OF POWER OF TESTS INVOLVING COVARIATES UNDER NONCONSTANT VARIANCE

In many experiments where pre-treatment and post-treatment measurements are taken, investigators wish to determine if there is a difference between two treatment groups. For this type of data, the post-treatment variable is used as the primary comparison variable and the pre-treatment variable is used as a covariate. Although most of the discussion in this paper is written with the pre-treatment variable as the covariate the results are applicable to other choices of a covariate. Tests based on residuals have been proposed as alternatives to the usual covariance methods. Our objective is to investigate how the powers of these tests are affected when the conditional variance of the post-treatment variable depends on the magnitude of the pre-treatment variable. In particular, we investigate two cases. [1] Crager, Michael R. 1987. Analysis of Covariance in Parallel-Group Clinical Trials With Pretreatment Baselines. Biometrics, 43: 895–901. [Crossref], [PubMed], [Web of Science ®] , [Google Scholar] The conditional variance of the post-treatment variable gradually increases as the magnitude of the pre-treatment variable increases. (In many biological models this is the case.) [2] Knoke, James D. 1991. Nonparametric Analysis of Covariance for Comparing Change in Randomized Studies with Baseline Values Subject to Error. Biometrics, 47: 523–533. [Crossref], [PubMed], [Web of Science ®] , [Google Scholar] The conditional variance of the post-treatment variable is dependent upon natural or imposed subgroups contained within the pre-treatment variable. Power comparisons are made using Monte Carlo techniques.     

ESTIMATION OF PARAMETERS IN ARFIMA PROCESSES: A SIMULATION STUDY

It is known that, in the presence of short memory components, the estimation of the fractional parameter d in an Autoregressive Fractionally Integrated Moving Average, ARFIMA(p, d, q), process has some difficulties (see [1] Smith, J., Taylor, N. and Yadav, S. 1997. Comparing the bias and misspecification in ARFIMA models. Journal of Time Series Analysis, 18(5): 507–527. [Crossref] , [Google Scholar]). In this paper, we continue the efforts made by Smith et al. [1] Smith, J., Taylor, N. and Yadav, S. 1997. Comparing the bias and misspecification in ARFIMA models. Journal of Time Series Analysis, 18(5): 507–527. [Crossref] , [Google Scholar] and Beveridge and Oickle [2] Beveridge, S. and Oickle, C. 1993. Estimating fractionally integrated time series models. Economics Letters, 43: 137–142.  [Google Scholar] by conducting a simulation study to evaluate the convergence properties of the iterative estimation procedure suggested by Hosking [3] Hosking, J. 1981. Fractional differencing. Biometrika, 68(1): 165–176. [Crossref], [Web of Science ®] , [Google Scholar]. In this context we consider some semiparametric approaches and a parametric method proposed by Fox-Taqqu[4] Fox, R. and Taqqu, M. S. 1986. Large-sample properties of parameter estimates for strongly dependent stationary gaussian time series. The Annals of Statistics, 14(2): 517–532. [Crossref], [Web of Science ®] , [Google Scholar]. We also investigate the method proposed by Robinson [5] Robinson, P. M. 1995a. Log-periodogram regression of time series with long range dependence. The Annals of Statistics, 23(3): 1048–1072. [Crossref], [Web of Science ®] , [Google Scholar] and a modification using the smoothed periodogram function.     

Generalized linear failure rate power series distribution

Abstract

In this article, we introduce a new class of lifetime distributions. This new class includes several previously known distributions such as those of Chahkandi and Ganjali (2009 Chahkandi, M., Ganjali, M. (2009). On some lifetime distributions with decreasing failure rate. Computat. Statist. Data Anal. 53:4433–4440.[Crossref], [Web of Science ®] , [Google Scholar]), Mahmoudi and Jafari (2012 Mahmoudi, E., Jafari, A.A. (2012). Generalized exponential power series distributions. Comput. Statist. Data Anal. 56(12):4047–4066.[Crossref], [Web of Science ®] , [Google Scholar]), and Nadarajah et al. (2012 Nadarajah, S., Shahsanaei, F., Rezaei, S. (2012). A new four-parameter lifetime distribution. J. Statist. Computat. Simul.. ifirst, 1–16. [Google Scholar]). This new class of four-parameter distributions allows for flexible failure rate behavior. Indeed, the failure rate function here can be increasing, decreasing, bathtub-shaped or upside-down bathtub-shaped. Several distributional properties of the new class including moments, quantiles and order statistics are studied. An EM algorithm for computing the estimates of the parameters involved is proposed and some maximum entropy characterizations are discussed. Finally, to show the flexibility and potential of the new class of distributions, applications to two real data sets are provided.     

First- and Second-order Asymptotics for the Tail Distortion Risk Measure of Extreme Risks

The tail distortion risk measure at level p was first introduced in Zhu and Li (2012 Zhu, L., Li, H. (2012). Tail distortion risk and its asymptotic analysis. Insur. Math. Econ. 51(1):115–121.[Crossref], [Web of Science ®] , [Google Scholar]), where the parameter p ∈ (0, 1) indicates the confidence level. They established first-order asymptotics for this risk measure, as p↑1, for the Fréchet case. In this article, we extend their work by establishing both first-order and second-order asymptotics for the Fréchet, Weibull, and Gumbel cases. Numerical studies are also carried out to examine the accuracy of both asymptotics.     

Bootstrap test of hypothesis for the multi-state models in survival analysis

Abstract

In the present paper we develop bootstrap tests of hypothesis, based on simulation, for the transition probability matrix arising in the context of a multi-state model. The bootstrap test statistic is based on the paper of Tattar and Vaman (2008 Tattar, P. N., Vaman, H. J. (2008). Testing transition probability matrix of a multi-state model with censored data. Lifetime Data Anal. 14(2):216–230.[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]), which develops a statistic for the testing problems concerning the transition probability matrix of the non homogeneous Markov process.     

A MONTE CARLO COMPARISON OF VARIOUS ASYMPTOTIC APPROXIMATIONS TO THE DISTRIBUTION OF INSTRUMENTAL VARIABLES ESTIMATORS

ABSTRACT

We examine empirical relevance of three alternative asymptotic approximations to the distribution of instrumental variables estimators by Monte Carlo experiments. We find that conventional asymptotics provides a reasonable approximation to the actual distribution of instrumental variables estimators when the sample size is reasonably large. For most sample sizes, we find Bekker[11] Bekker, P. A. 1994. Alternative Approximations to the Distributions of Instrumental Variable Estimators. Econometrica, 62: 657–681. [Crossref], [Web of Science ®] , [Google Scholar] asymptotics provides reasonably good approximation even when the first stage R ² is very small. We conclude that reporting Bekker[11] Bekker, P. A. 1994. Alternative Approximations to the Distributions of Instrumental Variable Estimators. Econometrica, 62: 657–681. [Crossref], [Web of Science ®] , [Google Scholar] confidence interval would suffice for most microeconometric (cross-sectional) applications, and the comparative advantage of Staiger and Stock[5] Staiger, D. and Stock, J. H. 1997. Instrumental Variables Regression with Weak Instruments. Econometrica, 65: 556–586. [Crossref], [Web of Science ®] , [Google Scholar] asymptotic approximation is in applications with sample sizes typical in macroeconometric (time series) applications.     

Group Variable Selection with Oracle Property by Weight-Fused Adaptive Elastic Net Model for Strongly Correlated Data

This paper is the generalization of weight-fused elastic net (Fu and Xu, 2012 Fu, G., Xu, Q. (2012). Grouping variable selection by weight fused elastic net for multi-collinear data. Communications in Statistics-Simulation and Computation 41(2):205–221.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]), which performs group variable selection by combining weight-fused LASSO(wfLasso) and elastic net (Zou and Hastie, 2005 Zou, H., Hastie, T. (2005). Regularization and variable selection via the elastic net. Journal of the Royal Statistical Society: Series B (Statistical Methodology) 67(2):301–320.[Crossref], [Web of Science ®] , [Google Scholar]) penalties. In this study, the elastic net penalty is replaced by adaptive elastic net penalty (AdaEnet) (Zou and Zhang, 2009 Zou, H., Zhang, H. (2009). On the adaptive elastic-net with a diverging number of parameters. Annals of Statistics 37(4):1733–1751.[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]), and a new group variable selection algorithm with oracle property (Fan and Li, 2001 Fan, J., Li, R. (2001). Variable selection via nonconcave penalized likelihood and its oracle properties. Journal of the American Statistical Association 96(456):1348–1360.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]; Zou, 2006 Zou, H. (2006). The adaptive lasso and its oracle properties. Journal of the American Statistical Association 101(476):1418–1429.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) is obtained.