Alan Kimber

Excess burden of constipation in Parkinson's disease: A pilot study

An analysis was undertaken of clinic‐based questionnaires that asked people with Parkinsons disease and a control group of older people without a known neurological condition about their experiences of constipation. People with Parkinsons disease report higher constipation on a validated objective measure, the Rome criterion (59% vs. 20.9%); a behavioral indicator, laxative‐taking (38.4% vs. 14.2%); and subjective self‐report of being always or often concerned by it (33.4% vs. 6.1%). Many people with Parkinsons disease experience constipation problems but they may not bring these to the attention of their healthcare providers. More research is required to understand the causes and management options.

A Statistical Analysis of Batting in Cricket

The batting average is ubiquitous in cricket. In this paper we show that the traditional batting average depends on an unrealistic parametric assumption. We propose a nonparametric approach based on runs scored for assessing batting performance. The methods have been applied to a large sample of players at various levels of cricket, examples of which are featured in this paper. The statistical methodology employed is akin to that used in reliability and survival analysis.

Bayesian analysis of an inverse Gaussian correlated frailty model

In survival analysis frailty is often used to model heterogeneity between individuals or correlation within clusters. Typically frailty is taken to be a continuous random effect, yielding a continuous mixture distribution for survival times. A Bayesian analysis of a correlated frailty model is discussed in the context of inverse Gaussian frailty. An MCMC approach is adopted and the deviance information criterion is used to compare models. As an illustration of the approach a bivariate data set of corneal graft survival times is analysed.

Meta‐analysis of clinical trials with rare events

Meta-analysis of rare event studies has recently become a subject of controversy and debate. We will argue and demonstrate in this paper that the occurrence of zero events in clinical trials or cohort studies, even if zeros occur in both arms (the case of a double-zero trial), is less problematic, at least from a statistical perspective, if the available statistical tools are applied in the appropriate way. In particular, it is neither necessary nor advisable to exclude studies with zero events from the meta-analysis. In terms of statistical tools, we will focus here on Mantel-Haenszel techniques, mixed Poisson regression and related regression models.

Proportional hazards models with discrete frailty

We extend proportional hazards frailty models for lifetime data to allow a negative binomial, Poisson, Geometric or other discrete distribution of the frailty variable. This might represent, for example, the unknown number of flaws in an item under test. Zero frailty corresponds to a limited failure model containing a proportion of units that never fail (long-term survivors). Ways of modifying the model to avoid this are discussed. The models are illustrated on a previously published set of data on failures of printed circuit boards and on new data on breaking strengths of samples of cord.

Detection of frailty in Weibull lifetime data using Outlier tests

Heterogeneity in lifetime data may be modelled by multiplying an individuals hazard by an unobserved frailty. We test for the presence of frailty of this kind in univariate and bivariate data with Weibull distributed lifetimes, using statistics based on the ordered Cox–Snell residuals from the null model of no frailty. The form of the statistics is suggested by outlier testing in the gamma distribution. We find through simulation that the sum of the k largest or k smallest order statistics, for suitably chosen k, provides a powerful test when the frailty distribution is assumed to be gamma or positive stable, respectively. We provide recommended values of k for sample sizes up to 100 and simple formulae for estimated critical values for tests at the 5% level.

A Parametric Model for the Interval Censored Survival Times of Acacia Mangium Plantation in a Spacing Trial

Abstract Survival times for the Acacia mangium plantation in the Segaliud Lokan Project, Sabah, East Malaysia were analysed based on 20 permanent sample plots (PSPs) established in 1988 as a spacing experiment. The PSPs were established following a complete randomized block design with five levels of spacing randomly assigned to units within four blocks at different sites. The survival times of trees in years are of interest. Since the inventories were only conducted annually, the actual survival time for each tree was not observed. Hence, the data set comprises censored survival times. Initial analysis of the survival of the Acacia mangium plantation suggested there is block by spacing interaction; a Weibull model gives a reasonable fit to the replicate survival times within each PSP; but a standard Weibull regression model is inappropriate because the shape parameter differs between PSPs. In this paper we investigate the form of the non-constant Weibull shape parameter. Parsimonious models for the Weibull survival times have been derived using maximum likelihood methods. The factor selection for the parameters is based on a backward elimination procedure. The models are compared using likelihood ratio statistics. The results suggest that both Weibull parameters depend on spacing and block.

A comparison of different ways of including baseline counts in negative binomial models for data from falls prevention trials

A common design for a falls prevention trial is to assess falling at baseline, randomize participants into an intervention or control group, and ask them to record the number of falls they experience during a follow-up period of time. This paper addresses how best to include the baseline count in the analysis of the follow-up count of falls in negative binomial (NB) regression. We examine the performance of various approaches in simulated datasets where both counts are generated from a mixed Poisson distribution with shared random subject effect. Including the baseline count after log-transformation as a regressor in NB regression (NB-logged) or as an offset (NB-offset) resulted in greater power than including the untransformed baseline count (NB-unlogged). Cook and Weis conditional negative binomial (CNB) model replicates the underlying process generating the data. In our motivating dataset, a statistically significant intervention effect resulted from the NB-logged, NB-offset, and CNB models, but not from NB-unlogged, and large, outlying baseline counts were overly influential in NB-unlogged but not in NB-logged. We conclude that there is little to lose by including the log-transformed baseline count in standard NB regression compared to CNB for moderate to larger sized datasets.

Model robust designs for survival trials

The exponential-based proportional hazards model is often assumed in time-to-event experiments but may only approximately hold. Deviations in different neighbourhoods of this model are considered that include other widely used parametric proportional hazards models and the data are assumed to be subject to censoring. Minimax designs are then found explicitly, based on criteria corresponding to classical c- and D-optimality. Analytical characterisations of optimal designs are provided which, unlike optimal designs for related problems in the literature, have finite support and thus avoid the issues of implementing a density-based design in practice. Finally, the proposed designs are compared with the balanced design that is traditionally used in practice, and recommendations for practitioners are given.

A covariance-based test for shared frailty in multivariate lifetime data

We decompose the score statistic for testing for shared finite variance frailty in multivariate lifetime data into marginal and covariance-based terms. The null properties of the covariance-based statistic are derived in the context of parametric lifetime models. Its non-null properties are estimated using simulation and compared with those of the score test and two likelihood ratio tests when the underlying lifetime distribution is Weibull. Some examples are used to illustrate the covariance-based test. A case is made for using the covariance-based statistic as a simple diagnostic procedure for shared frailty in a parametric exploratory analysis of multivariate lifetime data and a link to the bivariate Clayton–Oakes copula model is shown.