Abstract
We analyze predictions of future recruitment to a multicenter clinical trial based on a maximum-likelihood fitting of a commonly used hierarchical Poisson-gamma model for recruitments at individual centers. We consider the asymptotic accuracy of quantile predictions in the limit as the number of recruitment centers grows large and find that, in an important sense, the accuracy of the quantiles does not improve as the number of centers increases. When predicting the number of further recruits in an additional time period, the accuracy degrades as the ratio of the additional time to the census time increases, whereas when predicting the amount of additional time to recruit a further patients, the accuracy degrades as the ratio of to the number recruited up to the census period increases. Our analysis suggests an improved quantile predictor. Simulation studies verify that the predicted pattern holds for typical recruitment scenarios in clinical trials and verify the much improved coverage properties of prediction intervals obtained from our quantile predictor. In the process of extending the applicability of our methodology, we show that in terms of the accuracy of all integer moments it is always better to approximate the sum of independent gamma random variables by a single gamma random variable matched on the first two moments than by the moment-matched Gaussian available from the central limit theorem.
Highlights
1 Introduction randomised controlled trials represent the gold standard for evaluating the safety and efficacy of a new healthcare intervention or treatment [Akobeng, 2005]
A potential recruit is enroled at a given centre; after a lag, the potential recruit is screened for suitability, and if suitable that patient is randomised onto a particular treatment
Methods for predicting future recruitment usually model the probability of screening success separately, so we focus on the second stage of the process
Summary
Let C be the number of centres, and for c = 1, . C, let tc and Nc represent the time for which centre c was open before the census time and number recruited in centre c during the time tc. When all centres are open for the same time we denote that time by t. For Objective One, let t+ be the additional time ahead at which predictions will be made, and let Nc+ be the number recruited in centre c in that time, with N+ =. The negative binomial distribution of the number of successes until there are a failures when the probability of success is p is denoted NB(a, p). We use the notation →p and ⇒ to indicate convergence in probability and in distribution, respectively, and Φ to indicate the cumulative distribution function of a N(0, 1) random variable
Published Version (Free)
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have