Analyzing Supersaturated Designs by Means of an Information Based Criterion

The cost and time consumption of many industrial experimentations can be reduced using the class of supersaturated designs since this can be used for screening out the important factors from a large set of potentially active variables. A supersaturated design is a design for which there are fewer runs than effects to be estimated. Although there exists a wide study of construction methods for supersaturated designs, their analysis methods are yet in an early research stage. In this article, we propose a method for analyzing data using a correlation-based measure, named as symmetrical uncertainty. This method combines measures from the information theory field and is used as the main idea of variable selection algorithms developed in data mining. In this work, the symmetrical uncertainty is used from another viewpoint in order to determine more directly the important factors. The specific method enables us to use supersaturated designs for analyzing data of generalized linear models for a Bernoulli response. We evaluate our method by using some of the existing supersaturated designs, obtained according to methods proposed by Tang and Wu (1997 Tang , B. , Wu , C. F. J. (1997). A method for constructing supersaturated designs and its E(s ²)-optimality. Canadian Journal of Statistics 25:191–201.[Crossref], [Web of Science ®] , [Google Scholar]) as well as by Koukouvinos et al. (2008 Koukouvinos , C. , Mylona , K. , Simos , D. E. ( 2008 ). E(s ²)-optimal and minimax-optimal cyclic supersaturated designs via multi-objective simulated annealing . Journal of Statistical Planning and Inference 138 : 1639 – 1646 .[Crossref], [Web of Science ®] , [Google Scholar]). The comparison is performed by some simulating experiments and the Type I and Type II error rates are calculated. Additionally, Receiver Operating Characteristics (ROC) curves methodology is applied as an additional statistical tool for performance evaluation.     

Tests of Multivariate Independence for Ordinal Data

Population and sample versions of Kendall and Spearman measures of association suitable for multivariate ordinal data are defined. The latter generalize the indices of dependence of Ruymgaart and van Zuijlen (1978 Ruymgaart , F. H. , van Zuijlen , M. C. A. ( 1978 ). Asymptotic normality of multivariate linear rank statistics in the non-i.i.d. case . Ann. Statist. 6 : 588 – 602 .[Crossref], [Web of Science ®] , [Google Scholar]), Joe (1990 Joe , H. ( 1990 ). Multivariate concordance . J. Multivariate Anal. 35 : 12 – 30 .[Crossref], [Web of Science ®] , [Google Scholar]), and Schmid and Schmidt (2007 Schmid , F. , Schmidt , R. ( 2007 ). Multivariate extensions of Spearman's rho and related statistics . Statist. Probab. Lett. 77 : 407 – 416 .[Crossref], [Web of Science ®] , [Google Scholar]) by allowing atoms in the underlying distribution. The representation of the proposed empirical measures as U-statistics enables to establish their asymptotic normality under general distributions. A special attention is given to tests of independence for multivariate ordinal data, where the power of the new methodologies are investigated under fixed and contiguous alternatives.     

The Identification of Potential Longitudinal Biomarkers and Measurements of Effectiveness for Biomarkers as Surrogates in Multivariate Survival Data

A proposed method based on frailty models is used to identify longitudinal biomarkers or surrogates for a multivariate survival. This method is an extention of earlier models by Wulfsohn and Tsiatis (1997 Wulfsohn , M. S. , Tsiatis , A. A. ( 1997 ). A joint model for survival and longitudinal data measured with error . Biometrics 53 ( 1 ): 330 – 339 .[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]) and Song et al. (2002 Song , X. , Davidian , M. , Tsiatis , A. A. ( 2002 ). A Semiparametric likelihood approach to joint modeling of longitudinal and time-to-event data . Biometrics 58 ( 4 ): 742 – 753 .[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]). In this article, similar to Henderson et al. (2002 Henderson , R. , Diggle , P. J. , Dobson , A. ( 2002 ). Identification and efficacy of longitudinal markers for survival . Biostatistics 3 ( 1 ): 33 – 50 .[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]), a joint likelihood function combines the likelihood functions of the longitudinal biomarkers and the multivariate survival times. We use simulations to explore how the number of individuals, the number of time points per individual and the functional form of the random effects from the longitudianl biomarkers influence the power to detect the association of a longitudinal biomarker and the multivariate survival time. The proposed method is illustrate by using the gastric cancer data.     

Designs Robust Against the Presence of an Aberration in a Diallel Cross Design Setup

In this article, we find designs insensitive to the presence of an outlier in a diallel cross design setup for estimating a complete set of orthonormal contrasts among the effects of the general combining abilities of a set of parental lines. The criterion of robustness, suggested by Mandal (1989 Mandal , N. K. ( 1989 ). On robust designs . Cal. Stat. Assoc. Bull. 38 : 115 – 119 . [Google Scholar]) in block design setup and used by Biswas (2012 Biswas , A. ( 2012 ). Block designs robust against the presence of an aberration in a treatment-control setup . Commun. Statist. Theor. Meth. 41 : 920 – 933 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) in treatment-control setup, is adapted here. Complete diallel cross designs, suggested by Gupta and Kageyama (1994 Gupta , S. , Kageyama , S. ( 1994 ). Optimal complete diallel crosses . Biometrika 81 : 420 – 424 .[Crossref], [Web of Science ®] , [Google Scholar]), and partial diallel cross designs, suggested by Gupta et al. (1995 Gupta , S. , Das , A. , Kageyama , S. ( 1995 ). Single replicate orthogonal block designs for circulant partial diallel crosses . Commun. Statist. A 24 : 2601 – 2607 .[Taylor & Francis Online] , [Google Scholar]) and Mukerjee (1997 Mukerjee , R. ( 1997 ). Optimal partial diallel crosses . Biometrika 84 : 939 – 948 .[Crossref], [Web of Science ®] , [Google Scholar]), are found to be robust under certain conditions.     

Spurious Rejections with Endogenous Break Unit Root Tests in the Presence of Outliers and Breaks

Unit root tests with structural break developed by Zivot and Andrews (1992 Zivot , E. , Andrews , D. W. K. ( 1992 ). Further evidence on the great crash, the oil price shock and the unit root hypothesis . Journal of Business and Economic Statistics 10 : 251 – 270 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) and Perron and Rodriguez (2003 Perron , P. , Rodriguez , G. ( 2003 ). GLS detrending, efficient unit root tests and structural change . Journal of Econometrics 115 : 1 – 27 .[Crossref], [Web of Science ®] , [Google Scholar]) in the presence of additive outliers and breaks are studied by simulation experiments. The results show that the Zivot–Andrews test appears to have size distortions due to the additive outliers whereas the Perron–Rodriguez test exhibits good properties of size and power. However, the two tests are biased when a second break is present but not taken into account. Furthermore, these endogenous break unit root tests tend to determine the break point incorrectly at one period behind the true break point, leading to spurious rejections of the unit root null hypothesis.     

On Estimating the Finite Population Mean Using Two Auxiliary Variables in Two Phase Sampling in the Presence of Non Response

This article considers the problem of estimating the population mean using information on two auxiliary variables in the presence of non response under two-phase sampling. Some improved ratio-in-regression type estimators have been proposed in four different situations of non response along with their properties under large sample approximation. Efficiency comparisons of the proposed estimators with the usual unbiased estimator by Hansen and Hurwitz (1946 Hansen , M. H. , Hurwitz , W. N. ( 1946 ). The problem of non response in sample surveys . J. Amer. Statist. Assoc. 41 : 517 – 529 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]), conventional ratio and regression estimators using single auxiliary variable and Singh and Kumar (2010b Singh , H. P. , Kumar , S. ( 2010b ). Improved estimation of population mean under double sampling with sub-sampling the non-respondents . J. Statist. Plann. Infer. 140 ( 9 ): 2536 – 2550 .[Crossref], [Web of Science ®] , [Google Scholar]) estimators using two auxiliary variables have been made. Finally, these theoretical findings are illustrated by a numerical example.     

Adaptive Estimation of Causal Periodic Autoregressive Model

This article deals with the adaptive estimation of a periodic autoregressive model, with unspecified innovation density satisfying only some general technical assumptions. We first establish, while verifying the adapted sufficient conditions of Swensen (1985 Swensen , A. R. ( 1985 ). The asymptotic distribution of the likelihood ratio for autoregressive time series with a regression trend . Journal of Multivariate Analysis 16 : 54 – 70 .[Crossref], [Web of Science ®] , [Google Scholar]) to our model, the Local Asymptotic Normality (LAN), the Local Asymptotic Quadratic (LAQ), and the Local Asymptotic properties satisfied by its central sequence. Secondly, the Locally Asymptotically Minimax (LAM) estimators are constructed. Using these results, we construct the adaptive estimators of the unknown autoregressive parameters. The performances of the established estimators are shown, via simulation studies.     

The Modified Pulak and Al-Sultan's Model for Determining the Optimum Process Parameters

Economic selection of process parameters has been an important topic in modern statistical process control. The optimum process parameters setting have a major effect on the expected profit/cost per item. There are some concerns on the problem of setting process parameters. Boucher and Jafari (1991 Boucher , T. O. , Jafari , M. A. ( 1991 ). The optimum target value for single filling operations with quality sampling plans . J. Qual. Technol. 23 : 44 – 47 . [CSA] [Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) first considered the attribute single sampling plan applied in the selection of process target. Pulak and Al-Sultan (1996 Pulak , M. F. S. , Al-Sultan , K. S. ( 1996 ). The optimum targeting for a single filling operation with rectifying inspection . Omega 24 : 727 – 733 . [CSA] [CROSSREF] [Crossref], [Web of Science ®] , [Google Scholar]) extended Boucher and Jafari's model and presented the rectifying inspection plan for determining the optimum process mean. In this article, we further propose a modified Pulak and Al-Sultan model for determining the optimum process mean and standard deviation under the rectifying inspection plan with the average outgoing quality limit (AOQL) protection. Taguchi's (1986 Taguchi , G. ( 1986 ). Introduction to Quality Engineering . Asian Productivity Organization . [Google Scholar]) symmetric quadratic quality loss function is adopted for evaluating the product quality. By solving the modified model, we can obtain the optimum process parameters with the maximum expected profit per item and the specified quality level can be reached.     

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.     

Generalized E(s 2) Criterion for Multilevel Supersaturated Designs

Several extensions of the popular E(s ²) criterion of Booth and Cox (1962 Booth , K. H. V. , Cox , D. R. ( 1962 ). Some systematic supersaturated designs . Technometrics 4 : 489 – 495 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) to multilevel supersaturated designs have been advanced in literature. These extensions are not unique due to different ways they measure overall nonorthogonality between all pairs of the columns of the model matrix. We exploit the connection of the E(s ²) criterion with A- and D-optimality that naturally lends itself to a generalized criterion for the multilevel situation in a unified way. The extensions provided in literature follow as special cases of the generalized criterion. A lower bound to the generalized criterion is derived for a wide class of designs, and a method of construction for the symmetrical case is discussed.     

Some Properties of Semi-Folding Designs

Suppose that some information is available for the particular factor. The experimenter may apply the technique of foldover to isolate the factor and the two-factor interactions involving this factor. In fact, under some situations this can be done by the method of semi-folding. We will discuss this property in detail in this article. Furthermore, we use the computer to search the corresponding optimal semi-folding design for the given 2 ^k?p designs that are tabulated in Chen et al. (1993 Chen , J. , Sun , D. X. , Wu , C. F. J. ( 1993 ). A catalogue of two-level and three-level fractional factorial designs with small runs . International Statistical Review 61 : 131 – 145 .[Crossref], [Web of Science ®] , [Google Scholar]) research.     

A Comparison of Procedures for Controlling the False Discovery Rate in the Presence of Small Variance Genes: A Simulation Study

The Significance Analysis of Microarrays (SAM; Tusher et al., 2001 Tusher , V. G. , Tibshirani , R. , Chu , G. ( 2001 ). Significance analysis of microarrys applied to the ionizing radiation response . Proceedings of the National Academy of Sciences 98 : 5116 – 5121 .[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]) method is widely used in analyzing gene expression data while controlling the FDR by using resampling-based procedure in the microarray setting. One of the main components of the SAM procedure is the adjustment of the test statistic. The introduction of the fudge factor to the test statistic aims at deflating the large value of test statistics due to the small standard error of gene-expression. Lin et al. (2008 Lin , D. , Shkedy , Z. , Burzykowski , T. , Göhlmann , H. W. H. , De Bondt , A. , Perera , T. , Geerts , T. , Bijnens , L. ( 2008 ). Significance analysis of microarray (SAM) for comparisons of several treatments with one control . Biometric Journal, MCP 50 ( 5 ): 801 – 823 .[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]) pointed out that the fudge factor does not effectively improve the power and the control of the FDR as compared to the SAM procedure without the fudge factor in the presence of small variance genes. Motivated by the simulation results presented in Lin et al. (2008 Lin , D. , Shkedy , Z. , Burzykowski , T. , Göhlmann , H. W. H. , De Bondt , A. , Perera , T. , Geerts , T. , Bijnens , L. ( 2008 ). Significance analysis of microarray (SAM) for comparisons of several treatments with one control . Biometric Journal, MCP 50 ( 5 ): 801 – 823 .[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]), in this article, we extend our study to compare several methods for choosing the fudge factor in the modified t-type test statistics and use simulation studies to investigate the power and the control of the FDR of the considered methods.     

Inverse Sampling with Unequal Selection Probabilities

In some real situations the population of interest is divided into two groups, of which one contains only a few units. In other cases, the population may be considered as subdivided into two group', for example, if only a few units display a value of the variable of interest which is highly different from zero, while all the other units show a value equal to or near zero. In both cases, inverse sampling is more efficient than classical fixed sample-size designs to obtain the parameter estimators for the whole population as well as for its groups (e.g., Salehi and Seber, 2004 Salehi , M. M. , Seber , G. A. F. ( 2004 ). A general inverse sampling scheme and its application to adaptive cluster sampling . Austral. NZ J. Statist. 46 : 483 – 494 .[Crossref], [Web of Science ®] , [Google Scholar]). In fact, in this design the procedure selection continues until a prefixed number of units with the characteristic of interest is sampled. Since it is not known a priori to which group the population units belong, the sample size is a random variable. Christman and Lan (2001 Christman , M. C. , Lan , F. ( 2001 ). Inverse adaptive cluster sampling . Biometrics 57 : 1096 – 1105 .[Crossref], [PubMed], [Web of Science ®] , [Google Scholar]) and Salehi and Seber (2001 Salehi , M. M. , Seber , G. A. F. ( 2001 ). A new proof of Murthy's estimator which applies to sequential sampling . Austral. NZ J. Statist. 43 : 281 – 286 . [Google Scholar] 2004 Salehi , M. M. , Seber , G. A. F. ( 2004 ). A general inverse sampling scheme and its application to adaptive cluster sampling . Austral. NZ J. Statist. 46 : 483 – 494 .[Crossref], [Web of Science ®] , [Google Scholar]) considered inverse sampling designs when all the population units have equal selection probabilities. In this article, we consider the general case in which the units may have unequal probabilities of being included in the sample. In fact, in many real situations different units may have different selection probabilities because of some inherent features of the sampling procedure, or in order to obtain better estimates. We derive unbiased estimators of the totals of the two groups, their variance and the corresponding unbiased variance estimators in inverse sampling with replacement. Finally, we derive similar results for more complex designs, where the selection procedure stops before observing the prefixed number of units from the rare group.     

Circular Strongly Balanced Repeated Measurements Designs

Magda (1980 Magda , C. G. ( 1980 ). Circular balanced repeated measurements designs . Commun. Statist. Theor. Meth. 9 : 1901 – 1918 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) and Hedayat (1981 Hedayat , A. S. ( 1981 ). Repeated measurements designs-IV: recent advances (with discussion) . Bull. Int. Statist. Inst. 1 : 591 – 610 . [Google Scholar]) first considered the construction of circular strongly balanced repeated measurements designs. Sen and Mukerjee (1987 Sen , M. , Mukerjee , R. ( 1987 ). Optimal repeated measurements designs under interaction . J. Statist. Plann. Infer. 17 : 81 – 91 .[Crossref], [Web of Science ®] , [Google Scholar]) and Roy (1988 Roy , B. K. ( 1988 ). Construction of strongly balanced uniform repeated measurements designs . J. Statist. Plann. Infer. 19 : 341 – 348 .[Crossref], [Web of Science ®] , [Google Scholar]) considered the optimality and existence of circular strongly balanced repeated measurements designs based on the method of differences and Hamiltonian decomposition of lexicographic product of two graphs. In this article, we consider the construction of circular strongly balanced repeated measurements designs using the newly proposed method called cyclic shifts, and propose some new designs for p < v.     

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

Orthogonal Arrays for the Estimation of Global Sensitivity Indices Based on ANOVA High-Dimensional Model Representation

We present an efficient new approach to estimate global sensitivity indices based on ANOVA high-dimensional model representation. The method makes use of the properties of orthogonal arrays and extend the results of Wang et al. (2010 Wang , X. , Tang , Y. , Chen , X. , Zhang , Y. ( 2010 ). Design of experiment in global sensitivity analysis based on ANOVA high-dimensional model representation . Communications in Statistics—Simulation and Computation 39 : 1183 – 1195 .[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) to a more general situation. A theorem is given to illustrate that both the non influential and significant sensitivity indices are well estimated. This new method offers higher sampling efficiency than those of Sobol and Saltelli. We test its performance on two different models, which are widely used in engineering and statistical areas.     

D-Optimal Chemical Balance Weighing Designs with N ≡ 0 (mod 4) and 3 Objects

In this article, the problem of estimation of the individual weights of three objects using a chemical balance weighing design is considered. We use the criterion of D-optimality. We assume that the covariance matrix of errors is the matrix of first-order autoregressive process. Such problems were discussed in Li and Yang (2005 Li , C. H. , Yang , S. Y. ( 2005 ). On a conjecture in D-optimal designs with n ≡ 0 (mod 4) . Lin. Alg. Applic. 400 : 279 – 290 .[Crossref], [Web of Science ®] , [Google Scholar]) and also in Yeh and Lo Huang (2005 Yeh , H. G. , Lo Huang , M. N. ( 2005 ). On exact D-optimal designs with 2 two-level factors and n autocorrelated observations . Metrika 61 : 261 – 275 .[Crossref], [Web of Science ®] , [Google Scholar]). We present some results of D-optimal designs in certain class of designs with the design matrix X  ∈ M _n×3(±1) such that each column of matrix X has at least one 1 and one ?1.     

Generalized Scrambling in Quantitative Optional Randomized Response Models

Huang (2010 Huang , K. C. ( 2010 ). Unbiased estimators of mean, variance and sensitivity level for quantitative characteristics in finite population sampling . Metrika 71 : 341 – 352 .[Crossref], [Web of Science ®] , [Google Scholar]) proposed an optional randomized response model using a linear combination scrambling which is a generalization of the multiplicative scrambling of Eichhorn and Hayre (1983 Eichhorn , B. H. , Hayre , L. S. ( 1983 ). Scrambled randomized response methods for obtaining sensitive quantitative data . J. Statist. Plann. Infer. 7 : 307 – 316 .[Crossref], [Web of Science ®] , [Google Scholar]) and the additive scrambling of Gupta et al. (2006, 2010). In this article, we discuss two main issues. (1) Can the Huang (2010 Huang , K. C. ( 2010 ). Unbiased estimators of mean, variance and sensitivity level for quantitative characteristics in finite population sampling . Metrika 71 : 341 – 352 .[Crossref], [Web of Science ®] , [Google Scholar]) model be improved further by using a two-stage approach?; (2) Does the linear combination scrambling provide any benefit over the additive scrambling of Gupta et al. (2010 Gupta , S. N. , Shabbir , J. , Sehra , S. ( 2010 ). Mean and sensitivity estimation in optional randomized response models . J. Statist. Plann. Infer. 140 : 2870 – 2874 .[Crossref], [Web of Science ®] , [Google Scholar])? We will note that the answer to the first question is “yes” but the answer to the second question is “no.”     

A Note on Importance Resampling for Multi-Dimensional Statistics

Johns (1988 Johns , M. V. (1988). Importance sampling for bootstrap confidence intervals. Journal of the American Statistical Association 83:709–714.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]), Davison (1988 Davison , A. C. ( 1988 ). Discussion of paper by D. V. Hinkley . Journal of the Royal Statistical Society Series B 50 : 356 – 357 . [Google Scholar]), and Do and Hall (1991 Do , K. A. , Hall , P. ( 1991 ). On importance sampling for the bootstrap . Biometrika 78 : 161 – 167 .[Crossref], [Web of Science ®] , [Google Scholar]) used importance sampling for calculating bootstrap distributions of one-dimensional statistics. Realizing that their methods can not be extended easily to multi-dimensional statistics, Fuh and Hu (2004 Fuh , C. D. , Hu , I. ( 2004 ). Efficient importance sampling for events of moderate deviations with applications . Biometrika 91 : 471 – 490 .[Crossref], [Web of Science ®] , [Google Scholar]) proposed an exponential tilting formula for statistics of multi-dimension, which is optimal in the sense that the asymptotic variance is minimized for estimating tail probabilities of asymptotically normal statistics. For one-dimensional statistics, Hu and Su (2008 Hu , J. , Su , Z. ( 2008 ). Adaptive resampling algorithms for estimating bootstrap distributions . Journal of Statistical Planning and Inference 138 ( 6 ): 1763 – 1777 .[Crossref], [Web of Science ®] , [Google Scholar]) proposed a multi-step variance minimization approach that can be viewed as a generalization of the two-step variance minimization approach proposed by Do and Hall (1991 Do , K. A. , Hall , P. ( 1991 ). On importance sampling for the bootstrap . Biometrika 78 : 161 – 167 .[Crossref], [Web of Science ®] , [Google Scholar]). In this article, we generalize the approach of Hu and Su (2008 Hu , J. , Su , Z. ( 2008 ). Adaptive resampling algorithms for estimating bootstrap distributions . Journal of Statistical Planning and Inference 138 ( 6 ): 1763 – 1777 .[Crossref], [Web of Science ®] , [Google Scholar]) to multi-dimensional statistics, which applies to general statistics and does not resort to asymptotics. Empirical results on a real survival data set show that the proposed algorithm provides significant computational efficiency gains.     

A Simple Method for Constructing Multidimensional Distributions of Correlated Categorical Data

We suggest a method for constructing a multidimensional distribution of correlated categorical data with fixed marginal distributions and specified degrees of association based on the log-linear models. A convex combination approach by Lee (1997 Lee , A. J. ( 1997 ). Some simple methods for generating correlated categorical variates . Computational Statistics & Data Analysis 26 : 133 – 148 .[Crossref], [Web of Science ®] , [Google Scholar]) is applied to get a joint distribution with fixed Pearson chi-square coefficient. By using the suggested method, we can generate three-dimensional distributions which have a fixed association among three variables. Therefore, the suggested method could be extended to higher dimensions.