Bayesian Analysis of Two-Level Fractional Factorial Experiments with Non-Normal Responses

An intractable issue on screening experiments is to identify significant effects and to select the best model when the number of factors is large, especially for fractional factorial experiments with non-normal responses. In such cases, a three-stage Bayesian approach based on generalized linear models (GLMs) is proposed to identify which effects should be included in the target model where the principles of effect sparsity, hierarchy, and heredity are successfully considered. Three simulation experiments with non-normal responses which follow Poisson, binomial, and gamma distributions are presented to illustrate the performance of the proposed approach.     

Three Parallel Flats Designs for Two-level Factorial Experiments

This paper investigates the properties of the class of three parallel flats designs for two-level factorial experiments. It shows that the designs constructed from this class of designs can have a very simple correlation structure. The correlation of any pair of best linear unbiased estimators of factorial effects is 0, ⅓ or ¼. Furthermore, the designs obtained also have high D-efficiency. Finally, a class of designs is generated with run-size N = 12 to illustrate the use of the theorem.     

Goodness-of-fit Tests for GEE with Correlated Binary Data

The marginal logistic regression, in combination with GEE, is an increasingly important method in dealing with correlated binary data. As for independent binary data, when the number of possible combinations of the covariate values in a logistic regression model is much larger than the sample size, such as when the logistic model contains at least one continuous covariate, many existing chi-square goodness-of-fit tests either are not applicable or have some serious drawbacks. In this paper two residual based normal goodness-of-fit test statistics are proposed: the Pearson chi-square and an unweighted sum of residual squares. Easy-to-calculate approximations to the mean and variance of either statistic are also given. Their performance, in terms of both size and power, was satisfactory in our simulation studies. For illustration we apply them to a real data set.     

Determining the effects of intraclass correlation on factorial experiments

This paper considers the effects of the presence of eqmi-correlation between observations in a k-way factorial experiment, A technique to provide unbiased F-tests is also constructed.     

A Note on Tests for Zero-Inflation in Correlated Count Data

Zero-inflated Poisson mixed regression models are popular approaches to analyze clustered count data with excess zeros. Prior to application of these models, it is essential to examine the necessity of the adjustment for zero outcomes. The existing literature, however, has focused only on score tests for testing the suitability of zero-inflated models for correlated count data. In view of the observed bias and non-optimal size of score tests, it deserves further investigation of other alternative ways for the test. This article aims to explore the use of the null Wald and likelihood ratio tests for zero-inflation in correlated count data. Our simulation study shows that both the null Wald and likelihood ratio tests outperform the score test of Xiang et al. (2006 Xiang , L. , Lee , A. H. , Yau , K. K. W. , McLachlan , G. J. ( 2006 ). A score test for zero-inflation in correlated count data . Statistics in Medicine 25 : 1660 – 1671 . [Google Scholar]) in terms of statistical power, regardless of the computational convenience of the score test. A bootstrap null Wald statistic is also proposed, which results in improved performance in terms of the size and power of the test.     

Joint GEEs for multivariate correlated data with incomplete binary outcomes

This study considers a fully-parametric but uncongenial multiple imputation (MI) inference to jointly analyze incomplete binary response variables observed in a correlated data settings. Multiple imputation model is specified as a fully-parametric model based on a multivariate extension of mixed-effects models. Dichotomized imputed datasets are then analyzed using joint GEE models where covariates are associated with the marginal mean of responses with response-specific regression coefficients and a Kronecker product is accommodated for cluster-specific correlation structure for a given response variable and correlation structure between multiple response variables. The validity of the proposed MI-based JGEE (MI-JGEE) approach is assessed through a Monte Carlo simulation study under different scenarios. The simulation results, which are evaluated in terms of bias, mean-squared error, and coverage rate, show that MI-JGEE has promising inferential properties even when the underlying multiple imputation is misspecified. Finally, Adolescent Alcohol Prevention Trial data are used for illustration.     

Rank Tests for Clustered Survival Data

In a clinical trial, we may randomize subjects (called clusters) to different treatments (called groups), and make observations from multiple sites (called units) of each subject. In this case, the observations within each subject could be dependent, whereas those from different subjects are independent. If the outcome of interest is the time to an event, we may use the standard rank tests proposed for independent survival data, such as the logrank and Wilcoxon tests, to test the equality of marginal survival distributions, but their standard error should be modified to accommodate the possible intracluster correlation. In this paper we propose a method of calculating the standard error of the rank tests for two-sample clustered survival data. The method is naturally extended to that for K-sample tests under dependence.     

Simultaneous Confidence Intervals for the Population Cell Means,for Two-by-Two Factorial Data,that Utilize Uncertain Prior Information

Consider a two-by-two factorial experiment with more than one replicate. Suppose that we have uncertain prior information that the two-factor interaction is zero. We describe new simultaneous frequentist confidence intervals for the four population cell means, with simultaneous confidence coefficient 1 ? α, that utilize this prior information in the following sense. These simultaneous confidence intervals define a cube with expected volume that (a) is relatively small when the two-factor interaction is zero and (b) has maximum value that is not too large. Also, these intervals coincide with the standard simultaneous confidence intervals obtained by Tukey’s method, with simultaneous confidence coefficient 1 ? α, when the data strongly contradict the prior information that the two-factor interaction is zero. We illustrate the application of these new simultaneous confidence intervals to a real data set.     

Optimum designs for estimation of regression parameters in a balanced treatment incomplete block design set-up

The use of covariates in block designs is necessary when the experimental errors cannot be controlled using only the qualitative factors. The choice of values of the covariates for a given set-up attaining minimum variance for estimation of the regression parameters has attracted attention in recent times. In this paper, optimum covariate designs (OCD) have been considered for the set-up of the balanced treatment incomplete block (BTIB) designs, which form an important class of test-control designs. It is seen that the OCDs depend much on the methods of construction of the basic BTIB designs. The series of BTIB designs considered in this paper are mainly those as described by Bechhofer and Tamhane (1981) and Das et al. (2005). Different combinatorial arrangements and tools such as Hadamard matrices and different kinds of products of matrices viz Khatri-Rao product and Kronecker product have been conveniently used to construct OCDs with as many covariates as possible.     

MultiOrd: An R Package for Generating Correlated Ordinal Data

In this article, operational details of an R package MultiOrd that is designed for the generation of correlated ordinal data are described, and examples of some important functions are given. The package provides a valuable and needed tool that has been lacking for generating multivariate ordinal data.     

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.     

Estimating and Testing a Structured Covariance Matrix for Three-Level Multivariate Data

This article considers an approach to estimating and testing a new Kronecker product covariance structure for three-level (multiple time points (p), multiple sites (u), and multiple response variables (q)) multivariate data. Testing of such covariance structure is potentially important for high dimensional multi-level multivariate data. The hypothesis testing procedure developed in this article can not only test the hypothesis for three-level multivariate data, but also can test many different hypotheses, such as blocked compound symmetry, for two-level multivariate data as special cases. The tests are implemented with two real data sets.     

Variance least squares estimators for multivariate linear mixed model

Non-iterative, distribution-free, and unbiased estimators of variance components by least squares method are derived for multivariate linear mixed model. A general inter-cluster variance matrix, a same-member only general inter-response variance matrix, and an uncorrelated intra-cluster error structure for each response are assumed. Projection method is suggested when unbiased estimators of variance components are not nonnegative definite matrices. A simulation study is conducted to investigate the properties of the proposed estimators in terms of bias and mean square error with comparison to the Gaussian (restricted) maximum likelihood estimators. The proposed estimators are illustrated by an application of gene expression familial study.     

A basic lemma and the analysis of block and Kronecker product designs

In this paper the analysis of the class of block designs whose C matrix can be expressed in terms of the Kronecker product of some elementary matrices is considered. The analysis utilizes a basic result concerning the spectral decomposition of the Kronecker product of symmetric matrices in terms of the spectral decomposition of the component matrices involved in the Kronecker product. The property (A) of Kurkjian and Zelen (1963) is generalised and the analysis of generalised property (A) designs is given. It is proved that a design is balanced factorially if and only if it is a generalised property (A) design. A method of analysis of Kronecker product block designs whose component designs are equi-replicate and proper is also suggested.