Bootstrap approach to test the homogeneity of order restricted mean vectors when the covariance matrices are unknown

Testing homogeneity of multivariate normal mean vectors under an order restriction when the covariance matrices are unknown, arbitrary positive definite and unequal are considered. This problem of testing has been studied to some extent, for example, by Kulatunga and Sasabuchi (1984 Kulatunga, D. D. S., Sasabuchi, S. (1984). A test of homogeneity of mean vectors against multivariate isotonic alternatives. Mem Fac Sci, Kyushu Univ Ser A Mathemat 38:151–161. [Google Scholar]) when the covariance matrices are known and also Sasabuchi et al. (2003 Sasabuchi, S., Tanaka, K., Tsukamodo, T. (2003). Testing homogeneity of multivariate normal mean vectors under an order restriction when the covariance matrices are common but unknown. Annals of Statistics. 31(5):1517–1536.[Web of Science ®] , [Google Scholar]) and Sasabuchi (2007 Sasabuchi, S. (2007). More powerful tests for homogeneity of multivariate normal mean vectors under an order restriction. Sankhya 69(4):700–716. [Google Scholar]) when the covariance matrices are unknown but common. In this paper, a test statistic is proposed and because of the main advantage of the bootstrap test is that it avoids the derivation of the complex null distribution analytically, a bootstrap test statistic is derived and since the proposed test statistic is location invariance the bootstrap p-value defined logical and some steps are presented to estimate it. Our numerical studies via Monte Carlo simulation show that the proposed bootstrap test can correctly control the type I error rates. The power of the test for some of the p-dimensional normal distributions is computed by Monte Carlo simulation. Also, the null distribution of test statistic is estimated using kernel density. Finally, the bootstrap test is illustrated using a real data.     

Estimation of finite population mean using two auxiliary variables in stratified two-phase sampling

We propose a new ratio type estimator for estimating the finite population mean using two auxiliary variables in stratified two-phase sampling. Expressions for bias and mean squared error of the proposed estimator are derived up to the first order of approximation. The proposed estimator is more efficient than the usual stratified sample mean estimator, traditional stratified ratio estimator and some other stratified estimators including Bahl and Tuteja (1991 Bahl, S., Tuteja, R. K. (1991). Ratio and product type exponential estimators. Information and Optimization Sciences 12:159–163. [Google Scholar]), Chami et al. (2012 Chami, P. S., Singh, B., Thomas, D. (2012). A two-prameter ratio-product-ratio estimator using auxiliary information. ISRN Probability and Statistics 2012:1–15, doi: 10.5402/2012/103860.[Crossref] , [Google Scholar]), Chand (1975 Chand, L. (1975) Some Ratio Type Estimator Based on two or more Auxiliary Variables, Ph.D. dissertation, Iowa State University, Ames, Iowa (unpublished). [Google Scholar]), Choudhury and Singh (2012 Choudhury, S., Singh, B. K. (2012). A class of chain ratio-product type estimators with two auxiliary variables under double sampling scheme. Journal of the Korean Statistical Society 41:247–256. [Google Scholar]), Hamad et al. (2013 Hamad, N., Hanif, M., Haider, N. (2013). A regression type estimator with two auxiliary variables for two-phase sampling. Open Journal of Statistics, 3:74–78. [Google Scholar]), Vishwakarma and Gangele (2014 Vishwakarma, G. K., Gangele, R. K. (2014). A class of chain ratio-type exponential estimators in double sampling using two auxiliary variates. Applied Mathematics and Computation 227:171–175. [Google Scholar]), Sanaullah et al. (2014 Sanaullah, A., Ali, H. M., Noor ul Amin, M., Hanif, M. (2014). Generalized exponential chain ratio estimators under stratified two-phase random sampling. Applied Mathematics and Computation 226:541–547. [Google Scholar]), and Chanu and Singh (2014 Chanu, W. K., Singh, B. K. (2014). Improved class of ratio-cum-product estimators of finite population mean in two phase sampling. Global Journal of Science Frontier Research: F Mathematics and Decision Sciences 14(2):1–14. [Google Scholar]).     

Wavelets and estimation of long memory in nonstationary models: Does anything beat the exact local whittle estimator?

In this article, we analyze the performance of five estimation methods for the long memory parameter d. The goal of our article is to construct a wavelet estimate for the fractional differencing parameter in nonstationary long memory processes that dominate the well-known estimate of Shimotsu and Phillips (2005) Shimotsu, K., Phillips, P. (2005). Exact local whittle estimation of fractional integration. Annals of statistics 20:87–127. [Google Scholar]. The simulation results show that the wavelet estimation method of Lee (2005) Lee, J. (2005). Estimating memory parameter in the US inflation rate. Economics Letters 87:207–210. [Google Scholar] with several tapering techniques performs better under most cases in nonstationary long memory. The comparison is based on the empirical root mean squared error of each estimate.     

A new estimator of finite population mean based on the dual use of the auxiliary information

When a sufficient correlation between the study variable and the auxiliary variable exists, the ranks of the auxiliary variable are also correlated with the study variable, and thus, these ranks can be used as an effective tool in increasing the precision of an estimator. In this paper, we propose a new improved estimator of the finite population mean that incorporates the supplementary information in forms of: (i) the auxiliary variable and (ii) ranks of the auxiliary variable. Mathematical expressions for the bias and the mean-squared error of the proposed estimator are derived under the first order of approximation. The theoretical and empirical studies reveal that the proposed estimator always performs better than the usual mean, ratio, product, exponential-ratio and -product, classical regression estimators, and Rao (1991 Rao, T.J. (1991). On certail methods of improving ration and regression estimators. Commun. Stat. Theory Methods 20(10):3325–3340.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]), Singh et al. (2009 Singh, R., Chauhan, P., Sawan, N., Smarandache, F. (2009). Improvement in estimating the population mean using exponential estimator in simple random sampling. Int. J. Stat. Econ. 3(A09):13–18. [Google Scholar]), Shabbir and Gupta (2010 Shabbir, J., Gupta, S. (2010). On estimating finite population mean in simple and stratified random sampling. Commun. Stat. Theory Methods 40(2):199–212.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]), Grover and Kaur (2011 Grover, L.K., Kaur, P. (2011). An improved estimator of the finite population mean in simple random sampling. Model Assisted Stat. Appl. 6(1):47–55. [Google Scholar], 2014) estimators.     

On evaluating the performance of different forecasters

In this article, we have evaluated the performance of different forecasters and tested association between their performances for different pairs of variables. We have used three data sets of track records of professional U.S. economic forecasters participating in the Blue Chip consensus forecasting service (the data sets contain the root mean square errors (RMSE) of different forecasters for different years). To evaluate the performance of forecasters we have covered three well-known tests, namely the usual F test (cf. Fisher (1923 Fisher, R. A., Mackenzie, M. A. (1923). Studied in crop variation II. The manurial response of different potato. Journal of Agricultural Science 13:311–320. [Google Scholar])), Kruskal Wallis test (cf. Kruskal and Wallis (1952 Kruskall, W. H., Wallis, W. A. (1952). Use of ranks in one-criterion variance analysis. Journal of American Statistical Association 47:583–621. [Google Scholar])), and Extension of Median test (cf. Daniel (1990 Daniel, W. W. (1990). Applied Nonparametric Statistics. Duxbury Classic Series. (2nd Ed.), Boston. [Google Scholar])). To test the association between the forecaster's performances for different pairs of variables, we have considered Gini mean correlation coefficient r_g1 (cf. Yitzhaki, S., and Olkin, I. (1991 Yitzhaki, S., Olkin, I. (1991). Concentration indices and concentration curves, in K. Mosler and M. Scarsini (eds.), Stochastic Orders and Decisions under Risk, Institute of Mathematical Statistics: Lecture-Notes Monograph Series, 19, 1991, 380–392. [Google Scholar]) and Yitzhaki (2003 Yitzaki, S. (2003). Gini mean difference: A superior measure of variability for non normal distribution. Metron-International Journal of Statistics, LXI:285–316. [Google Scholar])), Modified rank correlation coefficient (cf. Zimmerman (1994 Zimmerman, D. W. (1994). A Note on modified rank correlation. Journal of educational and Behavioral Statistics 19:357–362. [Google Scholar])) and three modifications of Spearman rank correlation coefficient. We have observed that different forecasters do not necessarily offer same average performance. Moreover, an evidence of association between two criteria does not always lead us reaching at the same decision. The outcomes of the study may help the practitioners in selecting the best forecaster(s) for policymaking purposes.     

A two-stage algorithm for computation of the exact distribution of sums of spacings

Based on the recursions in Huffer (1988 Huffer, F. (1988). Divided differences and the joint distribution of linear combinations of spacings. Journal of Applied Probability 25:346–354. [Google Scholar]) and Huffer and Lin (2001 Huffer, F. W., Lin, C. T. (2001). Computing the joint distribution of general linear combinations of spacings or exponential variates. Statistica Sinica 11:1141–1157. [Google Scholar]), we present a two-stage algorithm and two specialized methods for evaluating the probabilities involving linear combination of spacings of special forms. The two-stage algorithm combines the advantages of marking algorithm in Huffer and Lin (1997 Huffer, F. W., Lin, C. T. (1997). Computing the exact distribution of the extremes of sums of consecutive spacings. Computational Statistics and Data Analysis 26:117–132. [Google Scholar]) and general algorithm in Huffer and Lin (2001 Huffer, F. W., Lin, C. T. (2001). Computing the joint distribution of general linear combinations of spacings or exponential variates. Statistica Sinica 11:1141–1157. [Google Scholar]). The proposed methods can analytically derive the exact expressions for some specific problems, and efficiently handle problems such as the distribution of the circular scan statistic and multiple coverage probabilities.     

An approximation to the convolution of gamma distributions

In general, the exact distribution of a convolution of independent gamma random variables is quite complicated and does not admit a closed form. Of all the distributions proposed, the gamma-series representation of Moschopoulos (1985 Moschopoulos, P. G. (1985). The distribution of the sum of independent gamma random variables. Annals of the Institute of Statistical Mathematics 37Part A:541–544. [Google Scholar]) is relatively simple to implement but for particular combinations of scale and/or shape parameters the computation of the weights of the series can result in complications with too much time consuming to allow a large-scale application. Recently, a compact random parameter representation of the convolution has been proposed by Vellaisamy and Upadhye (2009 Vellaisamy, P., Upadhye, N. S. (2009). On the sums of compound negative binomial and gamma random variables. Journal of Applied Probability 46:272–283.[Crossref], [Web of Science ®] , [Google Scholar]) and it allows to give an exact interpretation to the weights of the series. They describe an infinite discrete probability distribution. This result suggested to approximate Moschopoulos’s expression looking for an approximating theoretical discrete distribution for the weights of the series. More precisely, we propose a general negative binomial distribution. The result is an “excellent” approximation, fast and simple to implement for any parameter combination.     

Some new results on allocation of coverage limits and deductibles to mutually independent risks

This article further investigates the allocation of coverage limits and deductibles to multiple independent risks from the viewpoint of policyholders with increasing utility functions. In a more general setup, we develop the usual stochastic orders on the retained loss, which either generalize or supplement the corresponding results due to Lu and Meng (2011 Lu, Z., Meng, L. (2011). Stochastic comparisons for allocations of policy limits and deductibles with applications. Insur. Math. Econ. 48:338–343. [Google Scholar]) and Hu and Wang (2014 Hu, S., Wang, R. (2014). Stochastic comparisons and optimal allocation for policy limits and deductibles. Commun. Stat. Theory Methods 43:151–164. [Google Scholar]). Also, the most unfavorable and favorable allocations of coverage limits and deductibles are developed for multiple risks with dominated reversed hazard rates and hazard rates, respectively.     

MULTI-LEVEL FACTORIAL DESIGNS WITH MINIMUM NUMBERS OF LEVEL CHANGES

The order of experimental runs in a fractional factorial experiment is essential when the cost of level changes in factors is considered. The generalized foldover scheme given by [1] Coster, D. C. and Cheng, C. S. 1988. Minimum cost trend free run orders of fractional factorial designs. The Annals of Statistics, 16: 1188–1205. [Crossref], [Web of Science ®] , [Google Scholar]gives an optimal order to experimental runs in an experiment with specified defining contrasts. An experiment can be specified by a design requirement such as resolution or estimation of some interactions. To meet such a requirement, we can find several sets of defining contrasts. Applying the generalized foldover scheme to these sets of defining contrasts, we obtain designs with different numbers of level changes and then the design with minimum number of level changes. The difficulty is to find all the sets of defining contrasts. An alternative approach is investigated by [2] Cheng, C. S., Martin, R. J. and Tang, B. 1998. Two-level factorial designs with extreme numbers of level changes. The Annals of Statistics, 26: 1522–1539. [Crossref], [Web of Science ®] , [Google Scholar]for two-level fractional factorial experiments. In this paper, we investigate experiments with all factors in slevels.     

A modified chi-square test for Bertholon model with censored data

In this work, we propose the construction of a chi-squared goodness-of-fit test in censored data case, for Bertholon model which can analyse various competing risks of failure or death. This test is based on a modification of the Nikulin-Rao-Robson (NRR) statistic proposed by Bagdonavicius and Nikulin (2011a Bagdonavicius, V., Nikulin, M. (2011a). Chi-squared tests for general composite hypotheses from censored samples. Comptes Rendus Mathématiques: Series I 349(3–4):219–223. [Google Scholar], 2011b Bagdonavicius, V., Nikulin, M. (2011b). Chi-squared goodness-of-fit test for right censored data. International Journal of Applied Mathematics and Statistics 24:30–50. [Google Scholar]) for censored data. We applied this test to numerical examples from simulated samples and real data.     

Almost sure convergence of vertex degree densities in the vertex splitting model

We study the limiting degree distribution of the vertex splitting model introduced in Ref.^[3 David, F.; Dukes, M.; Jonsson, T.; Stefansson, S.Ö. Random tree growth by vertex splitting. J. Statist. Mech. Theory Exp. 2009, 04. doi:10.1088/1742-5468/2009/04/P04009. [Google Scholar]^]. This is a model of randomly growing ordered trees, where in each time step the tree is separated into two components by splitting a vertex into two, and then inserting an edge between the two new vertices. Under some assumptions on the parameters, related to the growth of the maximal degree of the tree, we prove that the vertex degree densities converge almost surely to constants which satisfy a system of equations. Using this, we are also able to strengthen and prove some previously non-rigorous results mentioned in the literature.     

On the performance of the poisson and the negative binomial ridge predictors

Mansson and Shukur (2011 Mansson, K., Shukur, G. (2011). A Poisson ridge regression estimator. Economic Modelling 28:1475–1481. [Google Scholar]) investigated the performance of the Poisson ridge regression (PRR) estimator in terms of the mean square error (MSE) criterion. Similarly, Mansson (2012 Mansson, K. (2012). On ridge estimators for the negative binomial regression model. Economic Modelling 29:178–184. [Google Scholar]) investigated the performance of the Negative binomial ridge regression (NBRR) according to the MSE criterion. But there is no any analysis of the predictive performance of the PRR and NBRR estimators. Therefore, we define the PRR and the NBRR predictors to evaluate their predictive performances according to the prediction mean squared error under the target function. The Monte Carlo simulations and the real life numerical example are conducted to investigate the defined predictors' performance.     

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.     

Bootstrap methods for multivariate hypothesis testing

The nonparametric and parametric bootstrap methods for multivariate hypothesis testing are developed. They are used to approximate the null distribution of the test statistics proposed by Duchesne and Francq (2015 Duchesne, P., Francq, C. (2015). Multivariate hypothesis testing using generalized and {2}-inverses—with applications. Statistics 49:475–496.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]), resulting in bootstrap testing procedures. In the problem of testing for the mean vector of a multivariate distribution, the asymptotic validity of the bootstrap methods is proved. The finite sample performance of the new solutions is demonstrated by means of Monte Carlo simulation studies. They indicate that for small-sample size, the bootstrap tests provide a better finite sample properties than the asymptotic tests considered by Duchesne and Francq (2015 Duchesne, P., Francq, C. (2015). Multivariate hypothesis testing using generalized and {2}-inverses—with applications. Statistics 49:475–496.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]).     

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.     

Control charts for generalized exponential distribution percentiles

The generalized exponential (GE) distribution, which was introduced by Mudholkar and Srivastava in 1993 Mudholkar, G. S., Srivastava, D. K. (1993). Exponentiated Weibull family for analyzing bathtub failure data. IEEE Transactions on Reliability 42:299–302. [Google Scholar], has been studied for various applications of lifetime modelings. In this article, five control charts, that comprise the Shewhart-type chart and four parametric bootstrap charts based on maximum likelihood estimation method, the moment estimation method, probability plot method, and least-square error method for the GE percentiles, are investigated. An extensive Monte Carlo simulation study is conducted to compare the performance among all five control charts in terms of average run length. Finally, an example is given for illustration.     

ROBUST TESTS FOR THE SIGNIFICANCE OF ORTHANT RESTRICTED MEAN VECTOR

Mudholkar and Srivastava [1] Mudholkar, G. S. and Srivastava, D. K. A class of robust stepwise alternatives to Hotelling's T²tests. Submitted to the Journal of Applied Statistics 1999 [Google Scholar]adapted Mudholkar and Subbaiah's [2] Mudholkar, G. S. and Subbaiah, P. 1980. Testing significance of a mean vector–a possible alternative to Hotelling's T². Ann. Inst. Statist. Math., 32(A): 43–52.  [Google Scholar]modified stepwise procedure, using the trimmed means in place of the means and appropriate studentization, to construct robust tests for the significance of a mean vector. They concluded that the robust alternatives provide excellent type I error control, and a substantial gain in power over Hotelling's T ²test in case of heavy tailed populations without significant loss of power when the population is normal. In this paper we adapt the modified stepwise approach to construct simple tests for the significance of the orthant constrained mean vector of a p-variate normal population with unknown covariance matrix, and also for constructing robust tests without assuming normality. The simple normal theory tests have exact type I error, whereas the robust tests provide a reasonably type I error control and substantial power advantage over Perlman's [3] Perlman, M. D. 1969. One-sided testing problems in multivariate analysis. Annals of Mathematical Statistics, 40: 549–567. [Crossref] , [Google Scholar]likelihood ratio test.     

Optimality of ratio type estimation methods for population mean in the presence of missing data

This article proposes various Searls-type ratio imputation methods (STRIM) on the lines of Ahmed et al. (2006 Ahmed, M. S., O. Al-Titi, Z. Al-Rawi, and W. Abu-Dayyeh. 2006. Estimation of a population mean using different imputation methods. Stat. Trans. 7 (6):1247–1264. [Google Scholar]). It is a well-known fact that the optimal ratio type estimator attains the MSE of regression estimator (or optimal difference estimator) but while using Searls-type transformation (STT) (Searls (1964 Searls, D. T. 1964. The utilization of a known coefficient of variation in the estimation procedure. J. Am. Stat. Assoc. 59:1225–1226.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar])) this may not always happen. These STRIM are shown to perform better than the imputation procedures of Ahmed et al. (2006 Ahmed, M. S., O. Al-Titi, Z. Al-Rawi, and W. Abu-Dayyeh. 2006. Estimation of a population mean using different imputation methods. Stat. Trans. 7 (6):1247–1264. [Google Scholar]). The STRIM may even outperform the Searls type difference imputation methods (STDIM) proposed by us in our earlier work, Bhushan and Pandey (2016 Bhushan, S., and A. P. Pandey. 2016. Optimal imputation of the missing data for estimation of population mean. Journal of Statistics and Management System 19 (6):755–69.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]). This study is concluded with the numerical study along with the theoretical comparison.     

Some extensions of circular balanced and circular strongly balanced repeated measurements designs

Repeated measurement designs are widely used in medicine, pharmacology, animal sciences, and psychology. In this paper the works of Iqbal and Tahir (2009 Iqbal, I., and M. H. Tahir. 2009. Circular strongly balanced repeated measurements designs. Communications in Statistics—Theory and Methods 38:3686–96.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]) and Iqbal, Tahir, and Ghazali (2010 Iqbal, I., M. H. Tahir, and S. S. A. Ghazali. 2010. Circular first- and second-order balanced repeated measurements designs. Communications in Statistics—Theory and Methods 39:228–40.[Taylor &; Francis Online], [Web of Science ®] , [Google Scholar]) are generalized for the construction of circular-balanced and circular strongly balanced repeated measurements designs through the method of cyclic shifts for three periods.