Some Results on Fractional Factorial Split-Plot Designs with Multi-Level Factors

Fractional factorial split-plot (FFSP) designs have received much attention in recent years. In this article, the matrix representation for FFSP designs with multi-level factors is first developed, which is an extension of the one proposed by Bingham and Sitter (1999b Bingham , D. , Sitter , R. R. ( 1999b ). Some theoretical results for fractional factorial split-plot designs . Ann. Statist. 27 : 1240 – 1255 . [Google Scholar]) for the two-level case. Based on this representation, periodicity results of maximum resolution and minimum aberration for such designs are derived. Differences between FFSP designs with multi-level factors and those with two-level factors are highlighted.     

Optimum Mixture Design Using Deficiency Criterion

In a mixture experiment the measured response is assumed to depend only on the relative proportion of ingredients or components present in the mixture. Scheffe (1958 Scheffe , H. ( 1958 ). Experiments with mixtures . J. R. Statist. Soc. B 20 : 344 – 360 . [Google Scholar], 1963 Scheffe , H. ( 1963 ). Simplex – centroid design for experiments with mixtures . J. R. Statist. Soc. B 25 : 235 – 263 . [Google Scholar]) first systematically considered this problem and introduced different models and designs suitable in such situations. Optimum designs for the estimation of parameters of different mixture models are available in the literature. The problem of estimating the optimum proportion of mixture components is of great practical importance. Pal and Mandal (2006 Pal , M. , Mandal , N. K. ( 2006 ). Optimum designs for optimum mixtures . Statist. Probab. Lett. 76 ( 13 ): 1369 – 1379 .[Crossref], [Web of Science ®] , [Google Scholar], 2007 Pal , M. , Mandal , N. K. (2007). Optimum mixture design via equivalence theory. communicated. JCISS 32:107–126. [Google Scholar]) attempted to find a solution to this problem by adopting a pseudo-Bayesian approach and using the trace criterion. Subsequently, Pal and Mandal (2008 Pal , M. , Mandal , N. K. ( 2008 ). Minimax designs for optimum mixtures . To appear in Statist. Probab. Lett.  [Google Scholar]) solved the problem using minimax criterion. In this article, the deficiency criterion due to Chatterjee and Mandal (1981 Chatterjee , S. K. , Mandal , N. K. ( 1981 ). Response surface designs for estimating the optimal point . Calcutta Statist. Assoc. Bull. . 30 : 145 – 169 .[Crossref] , [Google Scholar]) has been used as a measure for comparing the performance of competing designs.     

Computer-Generated Change-Over Designs for Correlated Observations

Change-over designs with independently distributed errors in the model have been studied extensively in the literature. Martin and Eccleston (2001 Martin , R. J. , Eccleston , J. A. ( 2001 ). Optimal and near optimal designs for dependent observations . Statist. Applic. 3 : 101 – 116 . [Google Scholar]) gave an algorithm for the generation of efficient change-over designs when the errors are correlated. This article proposes an algorithm for the generation of efficient change-over designs for estimation of direct effects of treatments in the presence of first-order residual effects in the model and when the errors are correlated.     

Selection of Binary Variables and Classification by Boosting

We adopt boosting for classification and selection of high-dimensional binary variables for which classical methods based on normality and non singular sample dispersion are inapplicable. Boosting seems particularly well suited for binary variables. We present three methods of which two combine boosting with the relatively classical variable selection methods developed in Wilbur et al. (2002 Wilbur , J. D. , Ghosh , J. K. , Nakatsu , C. H. , Brouder , S. M. , Doerge , R. W. ( 2002 ). Variable selection in high-dimensional multivariate binary data with application to the analysis of microbial community DNA fingerprints . Biometrics 58 : 378 – 386 . [Google Scholar]). Our primary interest is variable selection in classification with small misclassification error being used as validation of proposed method for variable selection. Two of the new methods perform uniformly better than Wilbur et al. (2002 Wilbur , J. D. , Ghosh , J. K. , Nakatsu , C. H. , Brouder , S. M. , Doerge , R. W. ( 2002 ). Variable selection in high-dimensional multivariate binary data with application to the analysis of microbial community DNA fingerprints . Biometrics 58 : 378 – 386 . [Google Scholar]) in one set of simulated and three real life examples.     

Reverse Foldovers for Blocked 2 m−k Fractional Factorial Designs

Two-level regular fractional factorial designs are often used in industry as screening designs to help identify early on in an experimental process those experimental or system variables which have significant effects on the process being studied. When the experimental material to be used in the experiment is heterogenous or the experiment must be performed over several well-defined time periods, blocking is often used as a means to improve experimental efficiency by removing the possible effects of heterogenous experimental material or possible time period effects. In a recent article, Li and Jacroux (2007 Li , F. , Jacroux , M. (2007). Optimal foldover plans for blocked 2 ^m?k fractional factorial designs. J. Statsist. Plann. Infer 137:2434–2452. [Google Scholar]) suggested a strategy for constructing optimal follow-up designs for blocked fractional factorial designs using the well-known foldover technique in conjunction with several optimality criteria. In this article, we consider the reverse foldover problem for blocked fractional factorial designs. In particular, given a 2^(m+p)?(p+k) blocked fractional factorial design D, we derive simple sufficient conditions which can be used to determine if there exists a 2^{(m+p?1)?(p?1+k+1)} initial fractional factorial design d which yields D as a foldover combined design as well how to generate all such d. Such information is useful in developing an overall experimental strategy in situations where an experimenter wants an overall blocked fractional factorial design with “desirable” properties but also wants the option of analyzing the observed data at the halfway mark to determine if the significant experimental variables are obvious (and the experiment can be terminated) or if a different path of experimentation should be taken from that initially planned.     

Optimum Designs for Estimation of Optimum Amount and Proportions in a Mixture-Amount Model

This article studies a mixture-amount model, which is quadratic both in the proportions of mixing components and the amount of mixture. Using the pseudo-Bayesian approach of Pal and Mandal (2006 Pal , M. , Mandal , N. K. ( 2006 ). Optimum designs for optimum mixtures . Statist. Probab. Lett. 76 : 1369 – 1379 .[Crossref], [Web of Science ®] , [Google Scholar]), it attempts to find the A-optimal design for the estimation of the optimum mixing proportions and the optimum amount.     

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.     

Comparison of Frailty Models for Acute Leukemia Data under Gompertz Baseline Distribution

In this article, we consider two different shared frailty regression models under the assumption of Gompertz as baseline distribution. Mostly assumption of gamma distribution is considered for frailty distribution. To compare the results with gamma frailty model, we consider the inverse Gaussian shared frailty model also. We compare these two models to a real life bivariate survival data set of acute leukemia remission times (Freireich et al., 1963 Freireich, E.J., Gehan, E., Frei, E., Schroeder, L.R., Wolman, I.J., Anbari, R., Burgert, E.O., Mills, S.D., Pinkel, D., Selawry, O.S., Moon, J.H., Gendel, B.R., Spurr, C.L., Storrs, R., Haurani, F., Hoogstraten, B., Lee, S. (1963). The effect of 6-mercaptopurine on the duration of steroid-induced remissions in acute leukemia: a model for evaluation of other potentially useful therapy. Blood 21:699–716.[Web of Science ®] , [Google Scholar]). Analysis is performed using Markov Chain Monte Carlo methods. Model comparison is made using Bayesian model selection criterion and a well-fitted model is suggested for the acute leukemia data.     

A Generalized Correlation Coefficient for the General Two-Treatment Crossover Design

This article extends the correlation methodology developed by Chinchilli et al. (2005 Chinchilli , V. M. , Phillips , B. R. , Mauger , D. T. , Szefler , S. J. ( 2005 ). A general class of correlation coefficients for the 2 × 2 crossover design . Biometr. J. 47 : 1 – 10 . [Google Scholar]) for the 2 × 2 crossover design to more complex crossover designs for clinical trials. We describe how the methodology can be adapted to a general type of two-treatment crossover design which includes either at least two sequences or at least two treatment periods or both. We then derive the asymptotic theory for the corresponding correlation statistics, investigate the statistical accuracy of the estimators via bootstrap analyses, and demonstrate their use with two real data examples.     

Construction and Nonparametric Testing of Orthogonal Arrays through Algebraic Strata and Inequivalent Permutation Matrices

In this article we use a new methodology, based on algebraic strata, to generate the class of all the orthogonal arrays of given size and strength. From this class we extract all the non isomorphic orthogonal arrays. Then, using all these non isomorphic orthogonal arrays, we suggest a method based on the inequivalent matrices permutations testing procedures Basso et al. (2004 Basso , D. , Evangelaras , H. , Koukouvinos , C. , Salmaso , L. ( 2004 ). Nonparametric testing for main effects on inequivalent designs. Proc. 7th Int. Workshop Model-Oriented Design Anal. Heeze, Netherlands, June 14–18 . [Google Scholar]) in order to obtain separate permutation tests for the effects in unreplicated mixed level fractional factorial designs. In order to validate the proposed method we perform a Monte Carlo simulation study and find out that the permutation tests appear to be a valid solution for testing effects, in particular when the usual normality assumptions cannot be justified.     

Optimal block designs for double cross experiments

ABSTRACT

This paper presents systematic methods of construction of optimal block designs for a double cross experiments for both even and odd values of “p” parental lines. The both even and odd values of designs are derived by using initial block of unreduced balanced incomplete block designs and initial block of row–column designs given by Bose et al. (1953 Bose, R.C., Shrikhande, S.S., Bhattacharya, K.N. (1953). On the construction of group divisible incomplete block design. Ann. Math. Stat. 24:167–195.[Crossref] , [Google Scholar]) and Gupta and Choi (1998 Gupta, S., Choi, K.C. (1998). Optimal row-column designs for complete diallel crosses. Commun. Stat. Theory Methods 27(11):2827–2835.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]), respectively. In this attempt we have found some new universally optimal block designs for double cross experiments.     

Marginal Correlation from Logit- and Probit-Beta-Normal Models for Hierarchical Binary Data

In hierarchical data settings, be it of a longitudinal, spatial, multi-level, clustered, or otherwise repeated nature, often the association between repeated measurements attracts at least part of the scientific interest. Quantifying the association frequently takes the form of a correlation function, including but not limited to intraclass correlation. Vangeneugden et al. (2010 Vangeneugden, T., Molenberghs, G., Laenen, A., Geys, H., Beunckens, C., Sotto, C. (2010). Marginal correlation in longitudinal binary data based on generalized linear mixed models. Communi. Stati. Theory &; Methods. 39:3542–3557. [Google Scholar]) derived approximate correlation functions for longitudinal sequences of general data type, Gaussian and non-Gaussian, based on generalized linear mixed-effects models. Here, we consider the extended model family proposed by Molenberghs et al. (2010 Molenberghs, G., Verbeke, G., Demétrio, C., Vieira, A. (2010). A family of generalized linear models for repeated measures with normal and conjugate random effects. Stat. Sci. 25:325–347.[Crossref], [Web of Science ®] , [Google Scholar]). This family flexibly accommodates data hierarchies, intra-sequence correlation, and overdispersion. The family allows for closed-form means, variance functions, and correlation function, for a variety of outcome types and link functions. Unfortunately, for binary data with logit link, closed forms cannot be obtained. This is in contrast with the probit link, for which such closed forms can be derived. It is therefore that we concentrate on the probit case. It is of interest, not only in its own right, but also as an instrument to approximate the logit case, thanks to the well-known probit-logit ‘conversion.’ Next to the general situation, some important special cases such as exchangeable clustered outcomes receive attention because they produce insightful expressions. The closed-form expressions are contrasted with the generic approximate expressions of Vangeneugden et al. (2010 Vangeneugden, T., Molenberghs, G., Laenen, A., Geys, H., Beunckens, C., Sotto, C. (2010). Marginal correlation in longitudinal binary data based on generalized linear mixed models. Communi. Stati. Theory &; Methods. 39:3542–3557. [Google Scholar]) and with approximations derived for the so-called logistic-beta-normal combined model. A simulation study explores performance of the method proposed. Data from a schizophrenia trial are analyzed and correlation functions derived.     

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.     

Marginal correlation from an extended random-effects model for repeated and overdispersed counts

Vangeneugden et al. [15 Vangeneugden, T., Molenberghs, G., Laenen, A., Geys, H., Beunckens, C. and Sotto, C. 2007. Marginal correlation in longitudinal binary data based on generalized linear mixed models, Tech. Rep., Hasselt University. submitted for publication [Google Scholar]] derived approximate correlation functions for longitudinal sequences of general data type, Gaussian and non-Gaussian, based on generalized linear mixed-effects models (GLMM). Their focus was on binary sequences, as well as on a combination of binary and Gaussian sequences. Here, we focus on the specific case of repeated count data, important in two respects. First, we employ the model proposed by Molenberghs et al. [13 Molenberghs, G., Verbeke, G. and Demétrio, C. G.B. 2007. An extended random-effects approach to modeling repeated, overdispersed count data. Lifetime Data Anal., 13: 513–531. [Crossref], [PubMed], [Web of Science ®] [Google Scholar]], which generalizes at the same time the Poisson-normal GLMM and the conventional overdispersion models, in particular the negative-binomial model. The model flexibly accommodates data hierarchies, intra-sequence correlation, and overdispersion. Second, means, variances, and joint probabilities can be expressed in closed form, allowing for exact intra-sequence correlation expressions. Next to the general situation, some important special cases such as exchangeable clustered outcomes are considered, producing insightful expressions. The closed-form expressions are contrasted with the generic approximate expressions of Vangeneugden et al. [15 Vangeneugden, T., Molenberghs, G., Laenen, A., Geys, H., Beunckens, C. and Sotto, C. 2007. Marginal correlation in longitudinal binary data based on generalized linear mixed models, Tech. Rep., Hasselt University. submitted for publication [Google Scholar]]. Data from an epileptic-seizures trial are analyzed and correlation functions derived. It is shown that the proposed extension strongly outperforms the classical GLMM.     

On Universal Optimality of Circular Weakly Neighbor Balanced Designs Under an Interference Model

In this article, the optimal design problem in a fixed effects interference model with left-neighbor effects is studied. It is known (Druilhet, 1999 Druilhet , P. ( 1999 ). Optimality of circular neighbor balanced designs . J. Statist. Plann. Infer. 81 : 141 – 152 .[Crossref], [Web of Science ®] , [Google Scholar]) that circular neighbor balanced designs (CNBDs) are universally optimal in such a model. We prove the universal optimality of circular weakly neighbor balanced designs (CWNBDs), which require a smaller number of blocks than CNBDs. CWNBDs with the number of blocks smaller than the number of treatments belong to the class of partially neighbor balanced designs (PNBDs) defined by Wilkinson et al. (1983 Wilkinson , G. N. , Eckert , R. , Hancock , T. W. , Mayo , O. ( 1983 ). Nearest neighbour (NN) analysis of field experiments . J. Roy. Statist. Soc. Ser. B 45 : 151 – 210 . [Google Scholar]). We give a construction method for some CWNBDs, with examples.     

A Method for Analyzing Supersaturated Designs with a Block Orthogonal Structure

Supersaturated designs is a large class of factorial designs which can be used for screening out the important factors from a large set of potentially active variables. The huge advantage of these designs is that they reduce the experimental cost drastically, but their critical disadvantage is the confounding involved in the statistical analysis. In this article, we propose a method for analyzing data using a specific type of supersaturated designs. This method heavily uses the special block orthogonal structure of the supersaturated designs given by Tang and Wu (1997 Tang , B. , Wu , C. F. J. ( 1997 ). A method for constructing supersaturated designs and its Es ²-optimality . Canadian J. Statist. 25 : 191 – 201 .[Crossref], [Web of Science ®] , [Google Scholar]). Also, we compare our method with several known statistical analysis methods by using some of the existing supersaturated designs. The comparison is performed by some simulating experiments and the Type I and Type II error rates are calculated. The results are presented in tables and the discussion to follow.     

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.     

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.     

Whole-Plot Exchange Algorithms for Constructing D-Optimal Multistratum Designs

Multistratum experiments contain several different sizes of experimental units. Examples include split-plot, strip-plot designs, and randomized block designs. We propose a strategy for constructing a D-optimal multistratum design by improving a randomly generated design through a sequence of whole-plot exchanges. This approach preserves the design structure and simplifies updates to the information and is applicable to any multistratum design where the largest-sized experimental unit is either a whole plot or a block. Two whole-plot exchange algorithms inspired by the point-exchange strategies of Fedorov (1972 Fedorov , V. V. ( 1972 ). Theory of Optimal Experiments . New York : Academic Press . [Google Scholar]) and Wynn (1972 Wynn , H. P. ( 1972 ). Results in the theory and construction of D-optimum experimental designs . Journal of the Royal Statistical Society Series B Statistics Methodology 34 : 133 – 147 . [Google Scholar]) are described. The application of the algorithms to several design problems is discussed.     

Population Attributable Fraction as a Measure of Heritability in Dichotomous Twin Data

In many genetic analyses of dichotomous twin data, odds ratios have been used to test hypotheses on heritability and shared common environment effects of a given disease (Lichtenstein et al., 2000 Lichtenstein , P. , Holm , N. , Verkasalo , P. , Iliadou , A. , Kaprio , J. , Koskenvuo , M. , Pukkala , E. , Skytthe , A. , Hemminki , K. ( 2000 ). Environmental and heritable factors in the causation of cancer . New England Journal of Medicine 343 : 78 – 85 .[Crossref], [Web of Science ®] , [Google Scholar]; Ahlbom et al., 1997 Ahlbom , A. , Lichtenstein , P. , Malmström , H. , Feychting , M. , Hemminki , K. , Pedersen , N. L. ( 1997 ). Cancer in twins: genetic and nongenetic familial risk factors . Journal of the National Cancer Institute 89 : 287 – 93 . [Google Scholar]; Ramakrishnan et al., 1992 Ramakrishnan , V. , Goldberg , J. , Henderson , W. , Elsen , S. , True , W. , Lyons , M. , Tsuang , M. ( 1992 ). Elementary methods for the analysis of dichotomous outcomes in unselected samples of twins . Genetic Epidemiology 9 : 273 – 287 . [Google Scholar], 4). However, estimates of these two effects have not been dealt with in the literature. In epidemiology, the attributable fraction (AF), a function of the odds ratio and the prevalence of the risk factor has been used to describe the contribution of a risk factor to a disease in a given population (Leviton, 1973 Leviton , A. ( 1973 ). Definitions of attributable risk . American Journal of Epidemiology 98 : 231 . [Google Scholar]). In this article, we adapt the AF to quantify the heritability and the shared common environment. Twin data on cancer, gallstone disease and phobia are used to illustrate the applicability of the AF estimate as a measure of heritability.