Binary Responses and the Three-period Cross-over

Problems with carry-over effects in the simple two-period cross-over have lead to interest in more complex cross-over designs. A method for analysing the optimum two-treatment three-period design with binary response variables is given by making a simple extension to Gart's logistic model. The method gives independent tests for, and estimates of the difference in treatment and first-order carry-over effects. An example of the analysis is given, using the loglinear models facility in GLIM.     

Modeling continuous auxiliary covariate data in generalized linear mixed models using the kernel smoother

Auxiliary covariate data are often collected in biomedical studies when the primary exposure variable is only assessed on a subset of the study subjects. In this study, we investigate a semiparametric‐estimated likelihood estimation for the generalized linear mixed models (GLMM) in the presence of a continuous auxiliary variable. We use a kernel smoother to handle continuous auxiliary data. The method can be used to deal with missing or mismeasured covariate data problems in a variety of applications when an auxiliary variable is available and cluster sizes are not too small. Simulation study results show that the proposed method performs better than that which ignores the random effects in GLMM and that which only uses data in the validation data set. We illustrate the proposed method with a real data set from a recent environmental epidemiology study on the maternal serum 1,1‐dichloro‐2,2‐bis(p‐chlorophenyl) ethylene level in relationship to preterm births.     

Score Tests for Exploring Complex Models: Application to HIV Dynamics Models

In biostatistics, more and more complex models are being developed. This is particularly the case in system biology. Fitting complex models can be very time‐consuming, since many models often have to be explored. Among the possibilities are the introduction of explanatory variables and the determination of random effects. The particularity of this use of the score test is that the null hypothesis is not itself very simple; typically, some random effects may be present under the null hypothesis. Moreover, the information matrix cannot be computed, but only an approximation based on the score. This article examines this situation with the specific example of HIV dynamics models. We examine the score test statistics for testing the effect of explanatory variables and the variance of random effect in this complex situation. We study type I errors and the statistical powers of this score test statistics and we apply the score test approach to a real data set of HIV‐infected patients.     

Modeling heterogeneity in nest survival data

Current statistical methods for estimating nest survival rates assume that nests are identical in their propensity to succeed. However, there are several biological reasons to question this assumption. For example, experience of the nest builder, number of nest helpers, genetic fitness of individuals, and site effects may contribute to an inherent disparity between nests with respect to their daily mortality rates. Ignoring such heterogeneity can lead to incorrect survival estimates. Our results show that constant survival models can seriously underestimate overall survival in the presence of heterogeneity. This paper presents a flexible random-effects approach to model heterogeneous nest survival data. We illustrate our methods through data on redwing blackbirds.     

Memory in coal tits: an alternative model

Jolliffe and Jolliffe (1997, Biometrics 53, 1136-1142) proposed various models for data from an experiment on memory in coal tits. This article describes an alternative model, which fits equally well and which may be simpler to interpret.     

Empirical Bayes estimation of random effects parameters in mixed effects logistic regression models

We extend an approach for estimating random effects parameters under a random intercept and slope logistic regression model to include standard errors, thereby including confidence intervals. The procedure entails numerical integration to yield posterior empirical Bayes (EB) estimates of random effects parameters and their corresponding posterior standard errors. We incorporate an adjustment of the standard error due to Kass and Steffey (KS; 1989, Journal of the American Statistical Association 84, 717-726) to account for the variability in estimating the variance component of the random effects distribution. In assessing health care providers with respect to adult pneumonia mortality, comparisons are made with the penalized quasi-likelihood (PQL) approximation approach of Breslow and Clayton (1993, Journal of the American Statistical Association 88, 9-25) and a Bayesian approach. To make comparisons with an EB method previously reported in the literature, we apply these approaches to crossover trials data previously analyzed with the estimating equations EB approach of Waclawiw and Liang (1994, Statistics in Medicine 13, 541-551). We also perform simulations to compare the proposed KS and PQL approaches. These two approaches lead to EB estimates of random effects parameters with similar asymptotic bias. However, for many clusters with small cluster size, the proposed KS approach does better than the PQL procedures in terms of coverage of nominal 95% confidence intervals for random effects estimates. For large cluster sizes and a few clusters, the PQL approach performs better than the KS adjustment. These simulation results agree somewhat with those of the data analyses.     

A mixed model for two-state Markov processes under panel observation

Many chronic medical conditions can be meaningfully characterized in terms of a two-state stochastic process. Here we consider the problem in which subjects make transitions among two such states in continuous time but are only observed at discrete, irregularly spaced time points that are possibly unique to each subject. Data arising from such an observation scheme are called panel data, and methods for related analyses are typically based on Markov assumptions. The purpose of this article is to present a conditionally Markov model that accommodates subject-to-subject variation in the model parameters by the introduction of random effects. We focus on a particular random effects formulation that generates a closed-form expression for the marginal likelihood. The methodology is illustrated by application to a data set from a parasitic field infection survey.     

Genetic diversity (RAPD-PCR) of lactobacilli isolated from "Almagro" eggplant fermentations from two seasons

Genetic diversity of 323 strains of lactobacilli isolated from an Almagro eggplant manufacturing plant was analyzed by using random amplified polymorphic DNA (RAPD). Thirty-four distinct RAPD patterns were obtained with 95% of isolates grouped into 18 main clusters. Genetic diversity was higher in brines from season II, Lactobacillus plantarum and Lactobacillus fermentum/cellobiosus being the species with the greatest number of genotypes. A single L. fermentum/cellobiosus genotype comprised isolates from both seasons and could be considered endemic to that factory.     

Analyzing a randomized trial on breast self-examination with noncompliance and missing outcomes

Recently, instrumental variables methods have been used to address non-compliance in randomized experiments. Complicating such analyses is often the presence of missing data. The standard model for missing data, missing at random (MAR), has some unattractive features in this context. In this paper we compare MAR-based estimates of the complier average causal effect (CACE) with an estimator based on an alternative, nonignorable model for the missing data process, developed by Frangakis and Rubin (1999, Biometrika, 86, 365-379). We also introduce a new missing data model that, like the Frangakis-Rubin model, is specially suited for models with instrumental variables, but makes different substantive assumptions. We analyze these issues in the context of a randomized trial of breast self-examination (BSE). In the study two methods of teaching BSE, consisting of either mailed information about BSE (the standard treatment) or the attendance of a course involving theoretical and practical sessions (the new treatment), were compared with the aim of assessing whether teaching programs could increase BSE practice and improve examination skills. The study was affected by the two sources of bias mentioned above: only 55% of women assigned to receive the new treatment complied with their assignment and 35% of the women did not respond to the post-test questionnaire. Comparing the causal estimand of the new treatment using the MAR, Frangakis-Rubin, and our new approach, the results suggest that for these data the MAR assumption appears least plausible, and that the new model appears most plausible among the three choices.     

Cell balance equation for chemotactic bacteria with a biphasic tumbling frequency

Alts three-dimensional cell balance equation characterizing the chemotactic bacteria was analyzed under the presence of one-dimensional spatial chemoattractant gradients. Our work differs from that of others who have developed rather general models for chemotaxis in the use of a non-smooth anisotropic tumbling frequency function that responds biphasically to the combined temporal and spatial chemoattractant gradients. General three-dimensional expressions for the bacterial transport parameters were derived for chemotactic bacteria, followed by a perturbation analysis under the planar geometry. The bacterial random motility and chemotaxis were summarized by a motility tensor and a chemotactic velocity vector, respectively. The consequence of invoking the diffusion-approximation assumption and using intrinsic one-dimensional models with modified cellular swimming speeds was investigated by numerical simulations. Characterizing the bacterial random orientation after tumbles by a turn angle probability distribution function, we found that only the first-order angular moment of this turn angle probability distribution is important in influencing the bacterial long-term transport. Mathematics Subject Classification (2000):60G05, 60J60, 82A70