- Research article
- Open Access
- Open Peer Review

# Comparing methods to estimate treatment effects on a continuous outcome in multicentre randomized controlled trials: A simulation study

- Rong Chu
^{1, 2}Email author, - Lehana Thabane
^{1, 2, 4}, - Jinhui Ma
^{1, 2}, - Anne Holbrook
^{1, 3, 4}, - Eleanor Pullenayegum
^{1, 2, 4}and - Philip James Devereaux
^{1}

**11**:21

https://doi.org/10.1186/1471-2288-11-21

© Chu et al; licensee BioMed Central Ltd. 2011

**Received:**26 May 2010**Accepted:**21 February 2011**Published:**21 February 2011

## Abstract

### Background

Multicentre randomized controlled trials (RCTs) routinely use randomization and analysis stratified by centre to control for differences between centres and to improve precision. No consensus has been reached on how to best analyze correlated continuous outcomes in such settings. Our objective was to investigate the properties of commonly used statistical models at various levels of clustering in the context of multicentre RCTs.

### Methods

Assuming no treatment by centre interaction, we compared six methods (ignoring centre effects, including centres as fixed effects, including centres as random effects, generalized estimating equation (GEE), and fixed- and random-effects centre-level analysis) to analyze continuous outcomes in multicentre RCTs using simulations over a wide spectrum of intraclass correlation (ICC) values, and varying numbers of centres and centre size. The performance of models was evaluated in terms of bias, precision, mean squared error of the point estimator of treatment effect, empirical coverage of the 95% confidence interval, and statistical power of the procedure.

### Results

While all methods yielded unbiased estimates of treatment effect, ignoring centres led to inflation of standard error and loss of statistical power when within centre correlation was present. Mixed-effects model was most efficient and attained nominal coverage of 95% and 90% power in almost all scenarios. Fixed-effects model was less precise when the number of centres was large and treatment allocation was subject to chance imbalance within centre. GEE approach underestimated standard error of the treatment effect when the number of centres was small. The two centre-level models led to more variable point estimates and relatively low interval coverage or statistical power depending on whether or not heterogeneity of treatment contrasts was considered in the analysis.

### Conclusions

All six models produced unbiased estimates of treatment effect in the context of multicentre trials. Adjusting for centre as a random intercept led to the most efficient treatment effect estimation across all simulations under the normality assumption, when there was no treatment by centre interaction.

## Keywords

- Mean Square Error
- Centre Size
- Generalize Estimate Equation
- Balance Design
- Empirical Coverage

## Background

A multicentre randomized control trial (RCT) is an experimental study "conducted according to a single protocol but at more than one site and, therefore, carried out by more than one investigator"[1]. Multicentre RCTs are usually carried out for two main reasons. First, they provide a feasible way to accrue sufficient participants to achieve reasonable statistical power to detect the effect of an experimental treatment compared with some control treatment. Second, by enrolling participants of more diverse demographics from a broader spectrum of geographical locations and various clinical settings, multicentre RCTs increase generalizability of the experimental treatment for future use [1].

Randomization is the most important feature of RCTs, for on average it balances known and unknown baseline prognostic factors between treatment groups, in addition to minimizing selection bias. Nevertheless, randomization does not guarantee complete balance of participant characteristics especially when the sample size is moderate or small. Stratification is a useful technique to guard against potential bias introduced by imbalance in key prognostic factors. In multicentre RCTs, investigators often use a stratified randomization design to achieve balance over key differences in study population (e.g. environmental, socio-economic or demographical factors) and management team (e.g. patient administration and management) at centre level to improve precision of statistical analysis [2]. Regulatory agencies recommend that stratification variables in design should usually be accounted for in analysis, unless the potential value of adjustment is questionable (e.g. very few subjects per centre) [1].

The current study was motivated by the COMPETE II trial which was designed to determine if an integrated computerized decision support system shared by primary care providers and patients could improve management of diabetes [3]. A total number of 511 patients were recruited from 46 family physician practices. Individual patients were randomized to one of the two intervention groups stratified by physician practice using permuted blocks of size 6.The number of patients treated by one physician varied from 1 to 26 (interquartiles = 7.25, 11, 15; mean = 11; standard deviation [SD] = 6). The primary outcome was a continuous variable representing the change of a 10-point process composite score based on eight diabetes-related component variables from baseline to a mean of 5.9 months' follow-up. A positive change indicated a favourable result. During the study, the possibility of clustering within physician practice and its consequence on statistical analysis was a concern to the investigators. The phenomenon of clustering emerges when outcomes observed from patients managed by the same centre, practice or physician are more similar than outcomes from different centres, practices or physicians. Clustering often arises in situations where patients are selective about which centre they belong to, patients in a centre or practice are managed according to the same clinical care paths, or patients influence each other in the same cluster [4]. Intraclass (or intracentre) correlation (ICC) is often used to quantify the average correlation between any two outcomes within the same cluster [5]. It is a number between zero and one. A large value indicates that within-cluster observations are similar relative to observations from other clusters and each observation within cluster contains less unique information. This implies that the independence assumption which many standard statistical models are based on is violated. An ICC of zero indicates that individual observations within the same clusters are uncorrelated and different clusters on average have similar observations.

Through a literature review, we identified six statistical methods that were sometimes employed to analyze continuous outcomes in multicentre RCTs: A. simple linear regression (two sample t-test), B. fixed-effects regression, C. mixed-effects regression, D. generalized estimating equations (GEE), E-1. fixed-effects centre-level analysis, and E-2. random-effects centre-level analysis. The first four methods use patient as unit of analysis, yet address centre effects differently [6–8]. Simple linear regression completely ignores centre effects that are likely to arise from two sources: (1) possible differences in environmental, socio-economic or treatment factors between centres, and (2) potential correlation among patients within centres. Although stratified randomization attempts to minimize the impact of centre on standard error of the treatment effect by ensuring that the treated and control groups are largely balanced with respect to centre, failure to control for stratification in analysis will likely inflate variance of the effect estimate. The fixed-effects model treats each participating centre as a fixed intercept to control for possible population or environmental differences among centres. This model assumes that study subjects from the same centre have independent outcomes, i.e. the intraclass correlation is fixed at zero. The mixed-effects model incorporates dependence of outcomes within a centre and treats centres as random intercepts. Proposed by Liang and Zeger [9], the generalized estimating equation (GEE) model extends generalized linear regression with continuous, categorical or count outcomes to correlated observations within cluster. Under a commonly used and perhaps oversimplified assumption, that the degree of similarity between any two outcomes from a centre is equal, an exchangeable correlation structure can be used to assess treatment effect in Model C and D. Though the within- and between-centre variances (${\sigma}_{e}^{2}$ and ${\sigma}_{b}^{2}$) are estimated differently in these two models. Method E-1 and E-2 are routinely employed to combine information from different studies in meta-analysis [10]. One can also apply them to aggregate treatment effects over multiple centres [11–13]. The overall effect is obtained as the average within-centre effect differences over centre, using inverse-variance weighting.

To date, only a few studies have been carried out to compare the performance of statistical models in analyzing multicentre RCTs using Monte Carlo simulation [6, 7, 14], whereas many studies assessed the impact of ICC in cluster randomization trials. Moerbeek et al [6] compared the simple linear regression model, fixed-effects regression and fixed-effects centre-level analysis with equal centre size. Pickering et al [7] examined the bias, precision and power of three methods: simple regression, fixed-effects and mixed-effects regression assuming block randomization of size 2 or 4 on a continuous outcome. In the presence of imbalance and non-orthogonality, they found ignoring centres or incorporating them as random-effects led to greater precision and smaller type II error compared with treating centres as fixed effects. Performance of the GEE approach and centre-level methods were not investigated in that work. Jones et al [14] compared the fixed-effects and random-effects regression models to a two-step Frequentist procedure as well as a Bayesian model, in the presence of treatment by centre interaction, and recommended fixed-effects weighted method for future analysis of multicentre trials. The investigation was further expanded to assessing correlated survival outcomes from large multicentre cancer trials. A series of random-effects approaches were proposed to account for centre or treatment by centre heterogeneity in proportional hazards models [15, 16].

A lack of definitive evidence on which models perform the best in various situations led to this comprehensive simulation study to examine the performance of all six commonly used models with continuous outcomes. The objective was to assess their comparative performance in terms of bias, precision (simulation standard deviation (SD) and average estimated SE), and mean squared error (MSE) of the point estimator of the treatment effect, empirical coverage of the 95% confidence interval (CI) and the empirical statistical power, over a wide spectrum of ICC value and centre size. We did not consider treatment by centre interaction this study, partly because clinicians and trialists have been making efforts to standardize the conduct and management of multicentre trials via, for instance, uniform patient selection criteria, staff training, and trial monitoring and auditing to reduce heterogeneity of treatment effects among centres. Furthermore it is uncommon to find clinical trials designed with sufficient power to detect treatment by covariate interactions.

In this paper, we survey six methods to investigate the effect of a treatment in multicentre RCTs in detail. We outline the design and analysis of an extensive simulation study, and report how model performance varies with ICC, centre size and the number of centres. We also present the estimated effect of the computer-aid decision support system on management of diabetes using different methods.

## Methods

### Approaches assessing treatment effects

We investigated six statistical approaches to evaluating effect of an experimental treatment on a continuous outcome compared with the control, for multicentre RCTs. Assuming baseline prognostic characteristics are approximately balanced between the treatment and control groups via randomization, we do not consider covariates other than centre effects in the models. The first four approaches use individual patient as unit of analysis, while centre is the unit of analysis in the last two approaches.

#### Simple linear regression (Model A)

where *Y*
_{
ij
} is the outcome of the i-th patient in the j-th centre, *X*
_{
ij
} stands for the treatment assignment (*X*
_{
ij
} = 1 for the treatment, *X*
_{
ij
} = 0 for the control), and *e*
_{
ij
} is the random error assumed to follow a normal distribution with mean 0 and variance ${\sigma}_{e}^{2}$. The intercept, *β*
_{0} , represents the mean outcome for the control group in all participating centres, and the slope *β*
_{1} represents effect of the treatment on the mean outcome.

#### Fixed-effects regression (Model B)

*β*

_{0j }) as a fixed effect by restricting the scope of statistical inference to the sample of participating centres in a RCT. Interpretation for

*β*

_{1}remains the same as in Model A. Model A and B were fitted using the linear model procedure 'lm( )' in

**R**.

#### Mixed-effects regression (Model C)

*β*

_{0j }=

*β*

_{0}+

*b*

_{0j }follows a normal distribution N(

*β*

_{0}, ${\sigma}_{b}^{2}$), and is thus random effect. In Equation 3,

*b*

_{0j }is the random deviation from the mean intercept

*β*

_{0}, specific for each centre.

Similar to the previous models, the within-centre variability is reflected by ${\sigma}_{e}^{2}$. The variability of outcome between-centre is captured by ${\sigma}_{b}^{2}$ in Model C. The variance and covariance of outcomes in the same or different centres can be expressed as: $\text{Var}({Y}_{ij})={\sigma}_{b}^{2}+{\sigma}_{e}^{2}$, $\text{Cov}({Y}_{ij},{Y}_{i\text{'}j})={\sigma}_{b}^{2}$, $\text{Cov}({Y}_{ij},{Y}_{i\text{'}j\text{'}})=0$. The intraclass correlation that measures the correlation among outcomes within centre is given by $\frac{{\sigma}_{b}^{2}}{{\sigma}_{e}^{2}+{\sigma}_{b}^{2}}$, assumed equal across all centres. We fitted Model C in **R** via linear mixed-effects procedure 'lme( )' using the restricted maximum likelihood (REML) method [17, 18].

#### Generalized estimating equations (Model D)

The GEE method has gained increasing popularity among health science researchers for its availability in most statistical software. As opposed to the mixed-effects method that estimates treatment difference between arms and individual centre effects, the GEE approach models the marginal population-average treatment effects in two steps: 1) it fits a naïve linear regression assuming independence between observations within and across centres, and 2) it estimates parameters of the working correlation matrix using residuals in the naïve model and refit regression model to adjust standard error and confidence interval for within-centre dependence [19]. As a result, the estimated impact of treatment on the outcome in GEE model reflects the "combined" within- and between-centre relationship. GEE employs quasi-likelihood to estimate regression coefficients iteratively, and a working correlation needs to be supplied to approximate the within centre correlation. When the working correlation is mis-specified, the sandwich-based covariance estimator will lead to a robust yet less efficient estimate of treatment effect in GEE model [9]. Recently, statisticians found that variance of the estimated treatment effect could be underestimated when the number of centres was small [20]. We therefore assessed the efficiency of GEE models using procedure 'gee( )' in library(gee) in **R**. As in the mixed-effects model, an exchangeable correlation structure was assumed in fitting Model D.

#### Centre-level fixed-effects model (Model E - 1)

*n*

_{ j }=

*n*

_{ tj }+

*n*

_{ cj }, given equal numbers of patients in two arms. Here

*n*

_{ tj }and

*n*

_{ cj }represent the number of patients in the treatment and control group, respectively, in the j-th centre. This form of the weighted analysis (without adjustment for covariates) was discussed extensively by many researchers [21–23]. We implemented Models E - 1 using the fixed-effects method for meta-analysis provided by the 'metacont( )' procedure in

**R**.

#### Centre-level random-effects model (Model E - 2)

A random-effects approach is used to aggregate mean effect differences not only across all participating centres but also across a population of centres represented by the sample. This model factors heterogeneity of treatment effect among centres (i.e. random treatment by centre interaction) into its weighting scheme and captures within- and between-centre variation of the outcome. One should not confuse this method with the mixed-effects model using patient-level data (Model C). For Model E-2, the underlying true treatment effects are not a fixed single value for all centres, rather they are considered random effects, normally distributed around a mean treatment effect with between-centre variation. Model C, on the other hand, treats centres as random intercepts and postulates the same treatment effect across all centres. Model E-2 does not serve as a fair comparator to the alternatives listed here which assume no treatment by centre interaction. Preliminary investigation suggested E-2 could outperform E-1 in some situations; we therefore included E-2 in the study to advance understanding of these models. Details of centre level models are provided in Figure 1. Model E - 2 was carried out using DerSimonian-Laird random-effects [24] method using the 'metacont( )' procedure in **R**. The confidence interval for Model E - 2 was constructed based on the within- and between-centre variances.

### Study data simulation

We used Monte Carlo simulation to assess performance of statistical models to analyze parallel group multicentre RCTs with a continuous outcome. We simulated outcome, *Y*, using the mixed-effects linear regression model (Model C): *Y*
_{
ij
} = *β*
_{0} + *b*
_{0j
}+ *β*
_{1}
*X*
_{
ij
} + *e*
_{
ij
} for the i-th patient in the j-th centre, where *X*
_{
ij
} (= 0, 1) is the dummy variable for treatment allocation (i = 1...*m*
_{
j
} , j = 1... J). We generated random error, *e*, from $N(0,{\sigma}_{e}^{2}=1)$. We set the true treatment effect (*β*
_{1}) to be 0.5 residual standard deviation (*σ*
_{
e
} ), an effect size suggested by the COMPETE II trial. This corresponds to a medium effect size according to Cohen's criterion [25]. To simulate centre effects, we employed the relationship between ICC and ${\sigma}_{b}^{2}:\text{ICC}=\frac{{\sigma}_{b}^{2}}{{\sigma}_{e}^{2}+{\sigma}_{b}^{2}}$: To fully study the behaviour of candidate models at various ICC levels, we considered the following values of ICC for completeness: 0.00, 0.01, 0.05, 0.10, 0.15, 0.20, 0.25, 0.30, 0.35, 0.40, 0.45, 0.50 and 0.75. This in turn set the corresponding ${\sigma}_{b}^{2}$ values to be 0, 1/99, 1/19, 1/9, 3/17, 1/4, 1/3, 3/7, 7/13, 2/3, 9/11, 1 and 3. However, we focused interpretation of the results on lower values of ICC as they were more likely to occur in practice [26–28].

Catalogue of simulation designs

Design | Scenario | Number of centres | Centre size | ICC |
---|---|---|---|---|

Balance | 1 | 45 | 4 | |

2 | 18 | 10 | ||

3 | 6 | 30 | ||

Chance imbalance | 4 | 45 | 4 | |

5 | 18 | 10 | 0 - 0.75 | |

6 | 6 | 30 | ||

Blocking (size = 2) | 7 | 17 | 1, 1, 4, 5, 5, 5, 8, 8, 10, 10, 10 | |

Blocking (size = 4) | 8 | 17 | 10, 15, 15, 20, 25, 28 | |

Blocking (size = 6) | 9 | 46 | Same as Table 2 | 0.125 |

*β*

_{0}= 1.34,

*β*

_{1}= 1.26, ${\sigma}_{b}^{2}=1$, ${\sigma}_{e}^{2}=7$, ICC = 0.125.

Centre composition of the COMPETE II trial

Number of patients per centre | Number of Centres |
---|---|

1 | 3 |

2 | 0 |

3 | 1 |

4 | 4 |

5 | 1 |

6 | 1 |

7 | 2 |

8 | 3 |

9 | 4 |

10 | 3 |

11 | 5 |

12 | 3 |

13 | 2 |

14 | 0 |

15 | 3 |

16 | 3 |

17 | 0 |

18 | 2 |

19 | 2 |

20 | 1 |

21 | 0 |

22 | 1 |

23 | 1 |

24 | 0 |

25 | 1 |

Estimates of intervention effects in COMPETE II trial

Model | Estimate of intervention effect | SE | 95% CI | Variance component |
---|---|---|---|---|

A: Simple linear regression | 1.270 | 0.246 | (0.787, 1.753) | ${\sigma}_{e}^{2}=7.712$ |

B: Fixed-effects regression | 1.291 | 0.231 | (0.836, 1.745) | ${\sigma}_{e}^{2}=6.682$ |

C: Mixed-effects regression | 1.263 | 0.230 | (0.811, 1.714) | ${\sigma}_{e}^{2}=6.678$ ${\sigma}_{b}^{2}=1.069$ |

D: GEE | 1.263 | 0.193 | (0.884, 1.641) | |

E - 1: centre-level Fixed-effects model | 1.397 | 0.219 | (0.967, 1.826) | |

E - 2: centre-level Random-effects model | 1.397 | 0.219 | (0.967, 1.826) |

We generated 1000 simulations for each of the 13 ICC values under each of the first eight scenarios and 1000 simulations for the specified ICC value for the ninth scenario. Separate sets of centre effects were simulated for each scenario and each simulation 1-1000. We chose to simulate 1000 replicates so that the simulation standard deviation for the empirical power at a nominal level of 90% in the absence of clustering was controlled at 1%. This also ensured that standard deviations of the coverage of the confidence interval and the empirical power not exceed 1.6%.

### Comparison of analytic models

We applied six statistical models to each simulated dataset. For each model, we calculated the bias, simulation standard deviation (SD), average of estimated standard error (SE) and mean squared error (MSE) of the point estimator of treatment effect (i.e. *β*
_{1}), empirical coverage of the 95% confidence interval around *β*
_{1} and the empirical statistical power. We constructed confidence intervals based on t-test for Models A - C, and Wald interval based on normal approximation for Models D and E. We estimated bias as the difference between the average estimate of *β*
_{1} over 1000 simulated datasets and the true effect. The simulation or empirical SD was calculated as the standard deviation of the estimated *β*
_{1}s across simulations, indicating precision of the estimator. We also obtain average of the estimated SEs from 1000 simulations to assess accuracy of variance estimator from each simulation dataset. The overall error rate of the point estimator was captured by the estimated MSE, enumerated by the average squared difference between the estimated *β*
_{1} and true value across the 1000 datasets. Furthermore, we reported performance of the interval estimators in each model. The empirical coverage was estimated as the proportion of 95% confidence intervals that covered the true *β*
_{1}, and the empirical power was the proportion of confidence intervals that rejected a false null hypothesis, i.e. zero lies outside CI. All datasets were simulated and analyzed in **R** version 2.4.1[29].

## Results

### Analysis of COMPETE II trial data

We applied all six models to the COMPETE II data and reported results in Table 3. Approximately equal numbers of patients were randomized to the intervention and control groups within each family doctor, leading to 253 and 258 patients in the intervention and control group, respectively. Among 46 family physicians, 11 physicians (24%) treated equal numbers of patients in two arms, 24 physicians (52%) treated one more patient in the intervention or control arm, 10 physicians (22%) managed 2 more patients in either arm, and one physician (2%) managed 3 more patients in one arm compared with the other.

All baseline characteristics were roughly balanced between arms [3]. The analyses using patient-level data produced similar estimates for *β*
_{1} and the effect size was around 0.5 times the corresponding residual standard deviation. The standard error of the estimated *β*
_{1} reduced from 0.25 (Model A) to 0.23 (Models B, C) then 0.19 (Model D) when centre effects were adjusted, leading to narrower CIs around estimated *β*
_{1} in Models B - D. The intraclass correlation was estimated 0.138 in Model C and 0.124 in Model D. The two centre-level analyses returned slightly larger estimates of *β*
_{1} than those from the individual patient-level models. In fact the minimal variance between physicians indicated no noticeable heterogeneity between physicians (τ^{2} = 0, *I*
^{2} = 0), resulting in same estimates from E-1 and E-2. Zero was not contained in the 95% confidence intervals, therefore all models led to the conclusion that the experimental intervention significantly improved patient management over usual care based on the change of composite process score.

### Balanced design with equal centre size

#### Properties of point estimates

*β*

_{1}were unbiased in all six models for all ICC values. Upon review, it was surprising that the point estimates in Model A, ignoring stratification and clustering, were invariant of ICC, and that the same estimates were returned by four patient-level models for each simulation. In fact, when treatments are allocated in same proportion in all centres, centre has no association with the treatment allocation, hence adjusting for centre effect or not has little impact on point estimate of the treatment - response relationship given a continuous response variable. For this reason, different ways to incorporate between-centre information (Models B-D) led to same estimates of treatment contrast in a balanced design. Same point estimates led to same empirical SD and overall error rate (measured by MSE) of the estimator in Models A - D regardless of ICC. Across different ICC values and scenarios 1 - 3, Models B and C yielded accurate estimates of the standard error of ${\widehat{\beta}}_{1}$ that approximated the empirical SD and the true standard deviation, 0.149, calculated using the best linear unbiased estimator of the simulation model, i.e. Model C [18]. From Table 4 we found that the standard error of ${\widehat{\beta}}_{1}$, in Model A increased with ICC in each scenario, deviating from the corresponding empirical SD. The standard error could be slightly underestimated in Model D when the number of centres was small (Table 4 scenario 2 and 3 comparing empirical SD and average SE). This agreed with previous work concerning small sample properties of the GEE model [20].

Properties of point estimates of the treatment effect from Models A - E in scenarios 1 to 3

ICC = 0.01 | ICC = 0.05 | ICC = 0.20 | |||||||
---|---|---|---|---|---|---|---|---|---|

Model | Mean effect (SD) | Ave. SE | MSE | Mean effect (SD) | Ave. SE | MSE | Mean effect (SD) | Ave. SE | MSE |

Scenario 1 - balanced design, 45 centres each with 4 subjects | |||||||||

A | 0.496 (0.148) | 0.149 | 0.022 | 0.499 (0.146) | 0.152 | 0.021 | 0.497 (0.151) | 0.167 | 0.023 |

B | 0.496 (0.148) | 0.148 | 0.022 | 0.499 (0.146) | 0.148 | 0.021 | 0.497 (0.151) | 0.149 | 0.023 |

C | 0.496 (0.148) | 0.147 | 0.022 | 0.499 (0.146) | 0.148 | 0.021 | 0.497 (0.151) | 0.149 | 0.023 |

D | 0.496 (0.148) | 0.146 | 0.022 | 0.499 (0.146) | 0.146 | 0.021 | 0.497 (0.151) | 0.147 | 0.023 |

E-1 | 0.496 (0.494) | 0.066 | 0.244 | 0.491 (0.454) | 0.066 | 0.206 | 0.506 (0.447) | 0.065 | 0.200 |

E-2 | 0.499 (0.163) | 0.172 | 0.027 | 0.497 (0.166) | 0.170 | 0.027 | 0.494 (0.162) | 0.170 | 0.026 |

Scenario 2 - balanced design, 18 centres each with 10 subjects | |||||||||

A | 0.490 (0.149) | 0.150 | 0.022 | 0.504 (0.155) | 0.152 | 0.024 | 0.498 (0.145) | 0.166 | 0.021 |

B | 0.490 (0.149) | 0.149 | 0.022 | 0.504 (0.155) | 0.149 | 0.024 | 0.498 (0.145) | 0.149 | 0.021 |

C | 0.490 (0.149) | 0.148 | 0.022 | 0.504 (0.155) | 0.148 | 0.024 | 0.498 (0.145) | 0.149 | 0.021 |

D | 0.490 (0.149) | 0.142 | 0.022 | 0.504 (0.155) | 0.143 | 0.024 | 0.498 (0.145) | 0.142 | 0.021 |

E-1 | 0.490 (0.178) | 0.130 | 0.032 | 0.501 (0.180) | 0.130 | 0.032 | 0.498 (0.171) | 0.130 | 0.029 |

E-2 | 0.492 (0.164) | 0.154 | 0.027 | 0.503 (0.165) | 0.155 | 0.027 | 0.498 (0.158) | 0.153 | 0.025 |

Scenario 3 - balanced design, 6 centres each with 30 subjects | |||||||||

A | 0.496 (0.149) | 0.149 | 0.022 | 0.492 (0.149) | 0.152 | 0.022 | 0.504 (0.151) | 0.164 | 0.023 |

B | 0.496 (0.149) | 0.149 | 0.022 | 0.492 (0.149) | 0.149 | 0.022 | 0.504 (0.151) | 0.149 | 0.023 |

C | 0.496 (0.149) | 0.149 | 0.022 | 0.492 (0.149) | 0.149 | 0.022 | 0.504 (0.151) | 0.149 | 0.023 |

D | 0.496 (0.149) | 0.130 | 0.022 | 0.492 (0.149) | 0.130 | 0.022 | 0.504 (0.151) | 0.149 | 0.023 |

E-1 | 0.497 (0.153) | 0.144 | 0.023 | 0.491 (0.154) | 0.144 | 0.024 | 0.508 (0.156) | 0.144 | 0.024 |

E-2 | 0.497 (0.151) | 0.163 | 0.023 | 0.491 (0.151) | 0.163 | 0.023 | 0.507 (0.153) | 0.161 | 0.023 |

The centre-level analyses produced larger empirical SE and MSE for ${\widehat{\beta}}_{1}$ compared with the patient-level analyses given small or moderate centre sizes (Table 4). The difference reduced as centre size increased. When only a few patients were enrolled per centre, the fixed-effects centre-level point estimator in Model E - 1 had large sampling variation that was severely underestimated at all ICC values. The random-effects model (E - 2) based on DerSimonian-Laird method on the other hand seemed to yield valid SE for ${\widehat{\beta}}_{1}$ that was on average greater than SEs from the patient-level models. The average estimate of SE for ${\widehat{\beta}}_{1}$ over all simulations in Model E - 2 was always larger than estimates of SE in Models B and C, followed by the SE estimated in Model E - 1 across different combinations of centre size and number of centres. In this study, although datasets were generated so that the treatment effects were homogeneous among centres (i.e. no treatment by centre interaction), random-effects analysis using centre-level data outperformed the fixed-effects analysis when the centre size was small, for Model E - 2 took into account the observed "heterogeneity" due to imprecise estimation of the centre mean difference and the associated standard error.

#### Properties of interval estimates

Coverage of the 95% interval estimate of the treatment effect and statistical power of Models A - E in scenarios 1 to 3

ICC = 0.01 | ICC = 0.05 | ICC = 0.20 | ||||
---|---|---|---|---|---|---|

Model | Cover. of CI | Power | Cover. of CI | Power | Cover. of CI | Power |

Scenario 1 - balanced design, 45 centres each with 4 subjects | ||||||

A | 0.952 | 0.901 | 0.953 | 0.912 | 0.973 | 0.862 |

B | 0.947 | 0.905 | 0.945 | 0.924 | 0.951 | 0.899 |

C | 0.947 | 0.907 | 0.944 | 0.924 | 0.951 | 0.899 |

D | 0.941 | 0.911 | 0.936 | 0.931 | 0.933 | 0.902 |

E-1 | 0.286 | 0.902 | 0.294 | 0.920 | 0.320 | 0.912 |

E-2 | 0.933 | 0.810 | 0.921 | 0.818 | 0.938 | 0.821 |

Scenario 2 - balanced design, 18 centres each with 10 subjects | ||||||

A | 0.955 | 0.899 | 0.941 | 0.903 | 0.973 | 0.881 |

B | 0.954 | 0.906 | 0.935 | 0.906 | 0.954 | 0.916 |

C | 0.951 | 0.908 | 0.935 | 0.906 | 0.954 | 0.916 |

D | 0.929 | 0.909 | 0.902 | 0.919 | 0.940 | 0.924 |

E-1 | 0.845 | 0.904 | 0.835 | 0.917 | 0.857 | 0.924 |

E-2 | 0.921 | 0.868 | 0.905 | 0.886 | 0.938 | 0.875 |

Scenario 3 - balanced design, 6 centres each with 30 subjects | ||||||

A | 0.953 | 0.905 | 0.949 | 0.901 | 0.966 | 0.888 |

B | 0.948 | 0.907 | 0.947 | 0.906 | 0.952 | 0.918 |

C | 0.948 | 0.907 | 0.946 | 0.906 | 0.952 | 0.918 |

D | 0.860 | 0.915 | 0.854 | 0.931 | 0.867 | 0.929 |

E-1 | 0.939 | 0.918 | 0.931 | 0.910 | 0.926 | 0.927 |

E-2 | 0.952 | 0.867 | 0.949 | 0.846 | 0.953 | 0.880 |

Overall, Models B and C had very close performance that outweighed other models in balanced design. Models C and D converged to a solution in all simulations.

### Design with equal centre size and chance imbalance

#### Properties of point estimates

Properties of point estimates of the treatment effect from Models A - E in scenarios 4 to 6

ICC = 0.01 | ICC = 0.05 | ICC = 0.20 | |||||||
---|---|---|---|---|---|---|---|---|---|

Model | Mean effect (SD) | Ave. SE | MSE | Mean effect (SD) | Ave. SE | MSE | Mean effect (SD) | Ave. SE | MSE |

Scenario 4 - chance imbalance, 45 centres each with 4 subjects | |||||||||

A | 0.502 (0.146) | 0.150 | 0.021 | 0.511 (0.154) | 0.154 | 0.024 | 0.494 (0.168) | 0.166 | 0.028 |

B | 0.506 (0.162) | 0.172 | 0.026 | 0.510 (0.174) | 0.172 | 0.030 | 0.496 (0.180) | 0.172 | 0.032 |

C | 0.502 (0.145) | 0.149 | 0.021 | 0.511 (0.155) | 0.152 | 0.024 | 0.496 (0.165) | 0.159 | 0.027 |

D | 0.501 (0.146) | 0.146 | 0.021 | 0.511 (0.155) | 0.149 | 0.024 | 0.496 (0.165) | 0.155 | 0.027 |

E-1 | 0.492 (0.525) | 0.122 | 0.275 | 0.504 (0.544) | 0.126 | 0.296 | 0.481 (0.482) | 0.127 | 0.232 |

E-2 | 0.506 (0.274) | 0.265 | 0.075 | 0.515 (0.284) | 0.269 | 0.081 | 0.490 (0.285) | 0.260 | 0.081 |

Scenario 5 - chance imbalance, 18 centres each with 10 subjects | |||||||||

A | 0.495 (0.148) | 0.150 | 0.022 | 0.498 (0.150) | 0.153 | 0.023 | 0.497 (0.169) | 0.166 | 0.028 |

B | 0.494 (0.156) | 0.157 | 0.024 | 0.498 (0.152) | 0.157 | 0.023 | 0.500 (0.161) | 0.157 | 0.026 |

C | 0.495 (0.148) | 0.150 | 0.022 | 0.498 (0.149) | 0.151 | 0.022 | 0.499 (0.159) | 0.153 | 0.025 |

D | 0.494 (0.148) | 0.142 | 0.022 | 0.498 (0.150) | 0.144 | 0.022 | 0.499 (0.159) | 0.148 | 0.025 |

E-1 | 0.488 (0.206) | 0.130 | 0.042 | 0.498 (0.199) | 0.130 | 0.039 | 0.503 (0.204) | 0.130 | 0.042 |

E-2 | 0.490 (0.177) | 0.163 | 0.031 | 0.501 (0.172) | 0.162 | 0.030 | 0.501 (0.178) | 0.164 | 0.032 |

Scenario 6 - chance imbalance, 6 centres each with 30 subjects | |||||||||

A | 0.499 (0.153) | 0.149 | 0.023 | 0.502 (0.150) | 0.153 | 0.022 | 0.510 (0.165) | 0.164 | 0.027 |

B | 0.499 (0.155) | 0.150 | 0.024 | 0.501 (0.150) | 0.151 | 0.022 | 0.507 (0.151) | 0.152 | 0.023 |

C | 0.499 (0.153) | 0.149 | 0.024 | 0.503 (0.149) | 0.150 | 0.022 | 0.507 (0.151) | 0.151 | 0.023 |

D | 0.498 (0.154) | 0.129 | 0.024 | 0.503 (0.150) | 0.131 | 0.022 | 0.508 (0.151) | 0.134 | 0.023 |

E-1 | 0.498 (0.159) | 0.146 | 0.025 | 0.502 (0.156) | 0.146 | 0.024 | 0.507 (0.157) | 0.146 | 0.025 |

E-2 | 0.498 (0.157) | 0.165 | 0.025 | 0.502 (0.153) | 0.164 | 0.023 | 0.507 (0.154) | 0.167 | 0.024 |

Coverage of the 95% interval estimate of the treatment effect and statistical power of Models A - E in scenarios 4 to 6

ICC = 0.01 | ICC = 0.05 | ICC = 0.20 | ||||
---|---|---|---|---|---|---|

Model | Cover. of CI | Power | Cover. of CI | Power | Cover. of CI | Power |

Scenario 4 - chance imbalance, 45 centres each with 4 subjects | ||||||

A | 0.954 | 0.910 | 0.949 | 0.918 | 0.942 | 0.845 |

B | 0.966 | 0.846 | 0.946 | 0.822 | 0.931 | 0.810 |

C | 0.954 | 0.912 | 0.945 | 0.917 | 0.934 | 0.870 |

D | 0.948 | 0.924 | 0.934 | 0.926 | 0.934 | 0.878 |

E-1 | 0.411 | 0.782 | 0.417 | 0.793 | 0.424 | 0.745 |

E-2 | 0.897 | 0.468 | 0.900 | 0.501 | 0.887 | 0.468 |

Scenario 5 - chance imbalance, 18 centres each with 10 subjects | ||||||

A | 0.954 | 0.898 | 0.949 | 0.900 | 0.942 | 0.843 |

B | 0.946 | 0.874 | 0.959 | 0.891 | 0.937 | 0.891 |

C | 0.952 | 0.898 | 0.951 | 0.904 | 0.939 | 0.895 |

D | 0.922 | 0.916 | 0.932 | 0.911 | 0.910 | 0.900 |

E-1 | 0.776 | 0.868 | 0.794 | 0.890 | 0.791 | 0.895 |

E-2 | 0.905 | 0.810 | 0.918 | 0.834 | 0.918 | 0.839 |

Scenario 6 - chance imbalance, 6 centres each with 30 subjects | ||||||

A | 0.950 | 0.897 | 0.953 | 0.905 | 0.961 | 0.860 |

B | 0.949 | 0.892 | 0.954 | 0.907 | 0.961 | 0.916 |

C | 0.950 | 0.897 | 0.952 | 0.905 | 0.959 | 0.910 |

D | 0.856 | 0.911 | 0.879 | 0.918 | 0.874 | 0.908 |

E-1 | 0.922 | 0.904 | 0.931 | 0.921 | 0.931 | 0.913 |

E-2 | 0.944 | 0.831 | 0.951 | 0.867 | 0.955 | 0.857 |

#### Properties of interval estimates

Similar results were observed relative to the balanced design. Patient - level models A - C guaranteed nominal coverage of confidence intervals at different ICC values, whereas the other models were likely to produce lower coverage under certain conditions. Among all models, Models C and D achieved the best empirical power that was closest to the nominal value of 91.8% across different centre sizes. When centre size was small and the number of centre was large (scenario 4), power for Models C and D also decreased with ICC, a pattern that was less obvious in scenarios 5 and 6. Models C and D achieved convergence in analyzing all simulated datasets.

### Design with unequal centre sizes and chance imbalance

Properties of point and 95% interval estimates calculated from Models A - E based on 1000 simulated datasets in scenario 9 - unbalanced, 46 centres, same centre composition as the COMPETE II trial.

SCENARIO 9 | |||||
---|---|---|---|---|---|

Model | Mean effect (SD) | Ave. SE | MSE | Cover. of CI | Power |

ICC = 0.125 | |||||

A | 1.254 (0.236) | 0.249 | 0.056 | 0.965 | 0.999 |

B | 1.253 (0.240) | 0.236 | 0.058 | 0.952 | 1 |

C | 1.253 (0.237) | 0.235 | 0.056 | 0.949 | 0.999 |

D | 1.253 (0.237) | 0.230 | 0.056 | 0.944 | 0.999 |

E-1 | 1.256 (0.405) | 0.207 | 0.165 | 0.787 | 0.991 |

E-2 | 1.257 (0.270) | 0.261 | 0.073 | 0.935 | 0.995 |

## Discussion

In this paper, we investigated six modelling strategies in a Frequentist framework to study the effect of an experimental treatment compared to the control treatment in the context of multicentre RCTs with a continuous outcome. We focused on three designs with equal or varying centre sizes and a treatment allocation ratio of 1:1 in the absence of treatment by centre interaction. Results of this simulation study showed that, when the proportion of patients allocated to the experimental treatment was the same in each centre or subject to chance imbalance only, models using patient-level and centre-level data yielded unbiased point estimates of treatment effect across a wide spectrum of ICC values. Ignoring stratification by centre or within-centre correlation did not bias the estimated treatment effects even when ICC was large. In fact, Parzen et al showed that mathematically the usual two-sample t-test, naively assuming independent observations of the response within centre was asymptotically unbiased in this context [30].

The simulation study also indicated that these models produced different standard errors of ${\widehat{\beta}}_{1}$, and the properties of interval estimates were affected by several factors: whether and how centre effects were incorporated in analysis, the combination of centre size and number of participating centres, and the level of non-orthogonality of the observed data. Treating centre as a random intercept resulted in the most precise estimate, and nominal values of coverage and power were attained in all circumstances. The fixed-effects model had extremely similar performance compared with the mixed-effects model in balanced design, but was slightly less efficient when the number of centres was large (J > 20) in an unbalanced design. Pickering and Weatherall observed the same pattern in their simulation study comparing three patient-level models with small ICC values [7]. The GEE model using information sandwich covariance method tended to underestimate the standard error across centre effects when the sample of centres was small, a property noticed by researchers [20, 31]. This resulted in higher statistical power. That is, the treatment effect estimate was more likely to be significant with a smaller standard error, but was associated with a lower coverage of the conference interval. Marray et al suggested that at least 40 centres should be used to ensure reliable estimate of standard error in the context of cluster randomized trials [32]. Our simulation results suggested that such cut value was also applicable to multicentre RCTs. Failure to control for centre effects in any form resulted in inflation of standard error, falsely high interval coverage and sizable drop of power, as ICC increased. Parzen et al quantified the impact of correlation among observations within centre on the variance of ${\widehat{\beta}}_{1}$ in Model A as 1/(1-ICC) [30]. Alternatively, one may consider a variant of robust variance estimation or a GEE model with an independent working correlation to control for the impact of ICC on variance estimation using t-test. Centre-level models generally produced larger standard errors, lower coverage or power than the patient-level models. Centre-level random-effects model incorporated variability of the treatment effect over centres, and was not a fair comparator to other models. Interestingly, this model seemed to fare better than the centre-level fixed-effects model in terms of precision and coverage even though the simulated datasets contained no treatment by centre interaction. Despite that the random-effects centre-level model may be a reasonable alternative for patient-level models when the number of patients per centre is large (≥30), centre-level models cannot adjust for patient-level covariates, a potential fatal drawback in the presence of patient prognostic imbalance.

Statisticians have different viewpoints on treating centre effects and treatment by centre interaction as fixed or random effects when analyzing multicentre RCTs [12, 13, 21, 33]. Our simulation results demonstrated the advantage of treating centres as random intercepts in the absence of treatment by centre interaction. When many centres enrol a few patients and allocation is unbalanced, the random intercept models can give more precise estimates of the treatment effect than the fixed intercept models, because they recover inter-centre information in unbalanced situations. For instance, in a multicentre RCT consisting of 45 centres each recruiting 4 patients, the empirical variance of the estimator of the treatment effect resulting from the fixed-effects model was 24.8% and 26.0% greater than that from the random-effects model when the ICC was 0.01 and 0.05, respectively. In the sentence alluded to, we need to compare the empirical variance of 0.162^{2} with the value of 0.145^{2} for ICC = 0.01, and 0.174^{2} to 0.155^{2} for ICC = 0.05 (Table 6 scenario 4). We therefore take the same position as Grizzle [33] and Agresti and Hartzel [12] that, "Although the clinics are not randomly chosen, the assumption of random clinic effect will result in tests and confidence intervals that better capture the variability inherent in the system more realistically than clinical effects are considered fixed".

Our results have some implications for the design of multicentre RCTs in the absence of treatment by centre interaction. First, regardless of the pre-determined allocation ratio, permutated block randomization (of relatively small block sizes) should be used to maintain approximate balance or orthogonality (i.e. same treatment allocation proportion across centres [7]) between treatments and centres, so that their individual effects can be evaluated independently. Variable block sizes can be used to strengthen allocation concealment. Second, for a given sample size, the number of patients randomized in majority of centres should be sufficiently large to ensure reliable estimate of within-centre variation. Third, it is essential for investigators to obtain a rough estimate of ICC for within-centre responses, through literature review or a pilot study. To reach nominal power of 80% or 90% (in the absence of clustering), centre effects should be taken into consideration in sample size assessment. When centre effects are included without treatment by centre interaction, the analysis becomes more powerful than a two-sample t-test. One method to assess sample size is to start with a two sample t-test for continuous outcomes (ignoring centre effect) then multiple the original estimated error variance by an variation inflation factor of 1/(1-ICC). This factor would have the effect of increasing the required sample size. Ignoring centre effects results in the larger sample size in the absence of interaction. Sample size determined using information sandwich covariance of GEE model could lead to slight loss of power, when the number of centres is small (≥40) and no proper adjustment is done. Lastly, there is no particular reason to require equal numbers of patients being enrolled in all participating centres and this is seldom the case in practice. Throughout the simulations, we observed similar results for studies of equal and varying centre sizes. In the study, we considered three scenarios representing the particular centre composition of the COMPETE II trial. For discussion on potential impact of enrolment patterns on the point and interval estimates of treatment effect, readers can refer to the publications on random enrolment verse determined enrolment, and relative efficiency between equal and unequal cluster sizes in the reference list [34, 35].

The current ICH E9 guideline recommends that researchers investigate treatment effect using a model that allows for centre differences in the absence of treatment by centre interaction [1]. However, it is implausible or impractical to include centre effects in statistical modelling or stratify randomization by centre, when it is anticipated from the start that trials may have very few subjects per centre. As it is acknowledged in the document, these recommendations are based on fixed-effects models. Mixed-effects models on the other hand may also be used to explore the centre and centre by interaction effects, especially when the number of centres is large [1]. Our simulation results indicated that when a considerable number of centres contains only a few patients, adjusting for centre as a fixed effect may lead to reduced precision (depending on distribution of patients between arms) compared with the naïve unadjusted analysis. Our work complements the ICH E9 guideline, by studying the impact of intraclass correlation on the assessment of treatment effects - a challenge that is seldom discussed, although routinely faced by investigators in reality. Our investigation suggests that, (1) ignoring centre effects completely may cause substantial overestimation of the standard error, faulty increase of coverage of the confidence interval and reduction of power; and (2) mixed-effects models and GEE models, if employed appropriately, can produce accurate and precise effect estimates, regardless of the degree of clustering. We recommend consider these methods in developing future guidelines.

When the number of patients per centre is very small, it is not practical to include centre as a fixed effect to analyze patient-level data, as centre effects cannot be reliably estimated, and precision of the treatment effect will be compromised. In fact for extremely small centres, all patients may be allocated to the same treatment group, and such centres will be ignored by the fixed-effects model [36–39]. The alternatives include collapsing all centres to perform a two-sample t-test, collapsing smaller centres to create an artificial centre and treating it as a fixed effect, and exploring other models discussed above. The mixed-effects model utilizes small centres more efficiently by "borrowing" information from larger centres. The GEE approach models the average treatment difference across all centres and adjusts for centre effects through a uniform correlation structure. This is an intuitively more efficient model which unfortunately does not always converge when the number of patients per centre was highly variable (simulation scenarios 7 and 8). In the current study, non-convergence problems were more likely to arise for very small or large ICC values (less than 0.1 or greater than 0.4 for block size 2 or 4) due to non-positive definite working correlation matrices, and the frequency could be as big as 10% after 2000 iterations. Conversely, convergence problems did not occur for the mixed-effects models in any scenarios. Our results show that analysis of trials consisting of very small centres (i.e. those containing less than 2 patients per arm) using centre-level models may not be an optimal strategy, because the within-centre standard deviation of treatment difference cannot be estimated for such centres, and consequently these very small centres are excluded from the analysis.

Results of two large empirical studies and one systematic review of cluster RCTs in primary care clinics suggested that most ICC values on physical, functional and social measures were less than 0.10 [26–28]. The estimated ICC in the COMPETE II trial using GEE and linear mixed-effects model, on the other hand, was 0.124 and 0.138, respectively. We chose to include rare yet possible large ICC values (0-0.75) in this simulation to examine the overall trend of model performance by ICC, and for the purpose of completeness and generalizability. Readers should anticipate the ICC values likely to emerge from their studies when interpreting these results. Throughout the work, we quantified correlation among subjects within centre using ICC, the most commonly used concept to assess clustering in biomedical literature. As indicated in previous sections, ICC reflects the interplay of two variance components in multicentre data: the between-centre variance and within-centre variance. These variance components are relatively easy to interpret for analysis of continuous outcomes using linear models. For analysis of binary or time-to-event data from multicentre trials using generalized mixed and frailty models, interpretation of centre heterogeneity can present challenges because random effects are linked to the outcome via nonlinear functions [40]. Reparameterization of the probability density function may be used to assess the impact of within- and between-centre variance. Interested readers can refer to Duchateau and Janssen [40] for more details.

A major limitation of the study is that it did not address model performance when the treatment by centre interaction exists. The interactions may be due to different patient populations or variable standard of care. Interested readers may read Moerbeek et al [6] for formulas of variance of ${\widehat{\beta}}_{1}$ in different models and Jones et al [14] for simulation results. Future studies addressing interaction effects in multicentre RCTs are needed. Datasets in the current paper were generated based on a moderate treatment effect reflected by the standardized mean difference between the treatment and control group. More or less prominent treatment effects are also likely to occur in clinical studies and similar findings are expected. The current study investigated on continuous outcomes in two groups from a Frequentist perspective. The models discussed above can be naturally extended to compare three or more treatments. Agresti and Hartzel [12] surveyed different methods to evaluate treatments for binary outcomes in multicentre RCTs. Non-parametric approaches and Bayesian methods are also available to obtain treatment contrast. Interested readers can refer to Aitkin [41], Gould [11], Smith et al [42], Legrand et al [16], and Louis [43], to name a few.

## Conclusions

We used simulations to investigate the performance of six statistical approaches that have been advocated to analyze continuous outcomes in multicentre RCTs. Our simulation study showed that all six models produced unbiased estimates of treatment effect in individual patient randomization multicentre trials. Adjusting for centre as random effects resulted in more efficient effect estimates in all scenarios over a wide spectrum of ICC values and various centre compositions. Fixed-effects model performed comparably to the mixed-effects model under most circumstances but lost efficiency when many centres contained a relatively small number of patients. The GEE model underestimated standard error of the effect estimates when a small number of centres were involved, and did not always converge when the centre size was variable for very large or small ICC values. Two-sample t-test severely overestimated standard error given moderate to large ICC values. The relative efficiencyof statistical modelling of treatment contrasts was also affected by ICC, distribution of patient enrolment, centre size and the number of centres.

## Declarations

### Acknowledgements

We thank the reviewers and Dr. Werner Vach for their constructive comments that have helped us to improve this work. We thank Ms. Ruthanne Cameron for proofreading the manuscript. Dr. Thabane is a clinical trials mentor for the Canadian Institutes of Health Research. Rong Chu was supported in part by funding from the CANNeCTIN (Canadian Network and Centre for Trials Internationally) training award, Father Sean O'Sullivan Research Center (FSORC) Studentship award and the Canadian Institute of Health Research (CIHR) doctoral research award.

## Authors’ Affiliations

## References

- International Conference on Harmonisation E9 Expert Working Group: ICH Harmonised Tripartite Guideline. Statistical principles for clinical trials. >Stat Med. 1999, 18 (15): 1905-1942.Google Scholar
- Lachin JM: Biostatistical methods: the assessment of relative risks:. 2000, New York: John Wiley and Sons, 1View ArticleGoogle Scholar
- Holbrook A, Thabane L, Keshavjee K, Dolovich L, Bernstein B, Chan D, Troyan S, Foster G, Gerstein H, COMPETE II Investigators: Individualized electronic decision support and reminders to improve diabetes care in the community: COMPETE II randomized trial. CMAJ. 2009, 181 (1-2): 37-44.View ArticlePubMedPubMed CentralGoogle Scholar
- Donner A, Klar N: Design and analysis of cluster randomization trials in health research:. 2000, London: Arnold, 1Google Scholar
- Kerry SM, Bland JM: The intracluster correlation coefficient in cluster randomisation. BMJ. 1998, 316 (7142): 1455-View ArticlePubMedPubMed CentralGoogle Scholar
- Moerbeek M, van Breukelen GJ, Berger MP: A comparison between traditional methods and multilevel regression for the analysis of multicenter intervention studies. J Clin Epidemiol. 2003, 56 (4): 341-350. 10.1016/S0895-4356(03)00007-6.View ArticlePubMedGoogle Scholar
- Pickering RM, Weatherall M: The analysis of continuous outcomes in multi-centre trials with small centre sizes. Stat Med. 2007, 26 (30): 5445-5456. 10.1002/sim.3068.View ArticlePubMedGoogle Scholar
- Worthington H: Methods for pooling results from multi-center studies. J Dent Res. 2004, 83 (Spec No C): C119-21.View ArticlePubMedGoogle Scholar
- Liang KY, Zeger SL: Longitudinal data analysis using generalized linear models. Biometrika. 1986, 73: 13-22. 10.1093/biomet/73.1.13.View ArticleGoogle Scholar
- Whitehead A: Meta-analysis of controlled clinical trials:. 2002, Chichester: John Wiley and Sons, 1View ArticleGoogle Scholar
- Gould AL: Multi-centre trial analysis revisited. Stat Med. 1998, 17 (15-16): 1779-97. 10.1002/(SICI)1097-0258(19980815/30)17:15/16<1779::AID-SIM979>3.0.CO;2-7. discussion 1799-800View ArticlePubMedGoogle Scholar
- Agresti A, Hartzel J: Strategies for comparing treatments on a binary response with multi-centre data. Stat Med. 2000, 19 (8): 1115-1139. 10.1002/(SICI)1097-0258(20000430)19:8<1115::AID-SIM408>3.0.CO;2-X.View ArticlePubMedGoogle Scholar
- Fleiss JL: Analysis of data from multiclinic trials. Control Clin Trials. 1986, 7 (4): 267-275. 10.1016/0197-2456(86)90034-6.View ArticlePubMedGoogle Scholar
- Jones B, Teather D, Wang J, Lewis JA: A comparison of various estimators of a treatment difference for a multi-centre clinical trial. Stat Med. 1998, 17 (15-16): 1767-77. 10.1002/(SICI)1097-0258(19980815/30)17:15/16<1767::AID-SIM978>3.0.CO;2-H. discussion 1799-800View ArticlePubMedGoogle Scholar
- Glidden DV, Vittinghoff E: Modelling clustered survival data from multicentre clinical trials. Stat Med. 2004, 23 (3): 369-388. 10.1002/sim.1599.View ArticlePubMedGoogle Scholar
- Legrand C, Ducrocq V, Janssen P, Sylvester R, Duchateau L: A Bayesian approach to jointly estimate centre and treatment by centre heterogeneity in a proportional hazards model. Stat Med. 2005, 24 (24): 3789-3804. 10.1002/sim.2475.View ArticlePubMedGoogle Scholar
- Brown HK, Kempton RA: The application of REML in clinical trials. Stat Med. 1994, 13 (16): 1601-1617. 10.1002/sim.4780131602.View ArticlePubMedGoogle Scholar
- McLean RA, Sanders WL: Approximating the degrees of freedom for SE's in mixed linear models. 1988, Proceedings of the Statistical Computing Section of the American Statistical Association. New Orleans, LouisianaGoogle Scholar
- Twisk J: Applied longitudinal data analysis for epidemiology: a practical guide. 2003, Cambridge: Cambridge University PressGoogle Scholar
- Ukoumunne OC, Carlin JB, Gulliford MC: A simulation study of odds ratio estimation for binary outcomes from cluster randomized trials. Stat Med. 2007, 26 (18): 3415-3428. 10.1002/sim.2769.View ArticlePubMedGoogle Scholar
- Senn S: Some controversies in planning and analysing multi-centre trials. Stat Med. 1998, 17 (15-16): 1753-65. 10.1002/(SICI)1097-0258(19980815/30)17:15/16<1753::AID-SIM977>3.0.CO;2-X. discussion 1799-800View ArticlePubMedGoogle Scholar
- Lin Z: An issue of statistical analysis in controlled multi-centre studies: how shall we weight the centres?. Stat Med. 1999, 18 (4): 365-373. 10.1002/(SICI)1097-0258(19990228)18:4<365::AID-SIM46>3.0.CO;2-2.View ArticlePubMedGoogle Scholar
- Kallen A: Treatment-by-centre interaction: what is the issue. Drug Info J. 1997, 31: 927-936.Google Scholar
- DerSimonian R, Laird N: Meta-analysis in clinical trials. Control Clin Trials. 1986, 7 (3): 177-188. 10.1016/0197-2456(86)90046-2.View ArticlePubMedGoogle Scholar
- Cohen J: Statistical Power Analysis for the Behavioral Sciences:. 1988, Hillsdale, NJ: Lawrence Erlbaum Associates, 2Google Scholar
- Adams G, Gulliford MC, Ukoumunne OC, Eldridge S, Chinn S, Campbell MJ: Patterns of intra-cluster correlation from primary care research to inform study design and analysis. J Clin Epidemiol. 2004, 57 (8): 785-794. 10.1016/j.jclinepi.2003.12.013.View ArticlePubMedGoogle Scholar
- Parker DR, Evangelou E, Eaton CB: Intraclass correlation coefficients for cluster randomized trials in primary care: the cholesterol education and research trial (CEART). Contemp Clin Trials. 2005, 26 (2): 260-267. 10.1016/j.cct.2005.01.002.View ArticlePubMedGoogle Scholar
- Smeeth L, Ng ES: Intraclass correlation coefficients for cluster randomized trials in primary care: data from the MRC Trial of the Assessment and Management of Older People in the Community. Control Clin Trials. 2002, 23 (4): 409-421. 10.1016/S0197-2456(02)00208-8.View ArticlePubMedGoogle Scholar
- R Core Development Team: R: A language and environment for statistical computing:. 2005, Vienna: R Foundation for Statistical Computing, 1Google Scholar
- Parzen M, Lipsitz S, Dear K: Does clustering affect the usual test statistics of no treatment effect in a randomized clinical trial?. Biomtrc J. 1998, 40: 385-402. 10.1002/(SICI)1521-4036(199808)40:4<385::AID-BIMJ385>3.0.CO;2-#.View ArticleGoogle Scholar
- Mancl LA, DeRouen TA: A covariance estimator for GEE with improved small-sample properties. Biometrics. 2001, 57 (1): 126-134. 10.1111/j.0006-341X.2001.00126.x.View ArticlePubMedGoogle Scholar
- Murray DM, Varnell SP, Blitstein JL: Design and analysis of group-randomized trials: a review of recent methodological developments. Am J Public Health. 2004, 94 (3): 423-432. 10.2105/AJPH.94.3.423.View ArticlePubMedPubMed CentralGoogle Scholar
- Grizzle JE: Letter to the editor. Control Clin Trials. 1987, 8 (4): 392-393. 10.1016/0197-2456(87)90158-9.View ArticlePubMedGoogle Scholar
- van Breukelen GJ, Candel MJ, Berger MP: Relative efficiency of unequal cluster sizes for variance component estimation in cluster randomized and multicentre trials. Stat Methods Med Res. 2008, 17 (4): 439-458. 10.1177/0962280206079018.View ArticlePubMedGoogle Scholar
- Fedorov V, Jones B: The design of multicentre trials. Stat Methods Med Res. 2005, 14 (3): 205-248. 10.1191/0962280205sm399oa.View ArticlePubMedGoogle Scholar
- Localio AR, Berlin JA, Ten Have TR, Kimmel SE: Adjustments for center in multicenter studies: an overview. Ann Intern Med. 2001, 135 (2): 112-123.View ArticlePubMedGoogle Scholar
- Breslow NE, Day NE: Statistical Methods in Cancer Research. Volume I: The Analysis of Case-Control Studies. 1980, Lyon: International Agency for Research on CancerGoogle Scholar
- Cox DR, Hinkley DV: Theoretical Statistics:. 1974, London: Chapman & HallView ArticleGoogle Scholar
- Graubard BI, Korn EL: Regression analysis with clustered data. Stat Med. 1994, 13 (5-7): 509-522. 10.1002/sim.4780130514.View ArticlePubMedGoogle Scholar
- Duchateau L, Janssen P: Understanding heterogeneity in generalized mixed and frailty models. The American Statistician. 2005, 59 (2): 143-146. 10.1198/000313005X43236.View ArticleGoogle Scholar
- Aitkin M: A general maximum likelihood analysis of variance components in generalized linear models. Biometrics. 1999, 55 (1): 117-128. 10.1111/j.0006-341X.1999.00117.x.View ArticlePubMedGoogle Scholar
- Smith TC, Spiegelhalter DJ, Thomas A: Bayesian approaches to random-effects meta-analysis: a comparative study. Stat Med. 1995, 14 (24): 2685-2699. 10.1002/sim.4780142408.View ArticlePubMedGoogle Scholar
- Louis TA: Using empirical Bayes methods in biopharmaceutical research. Stat Med. 1991, 10 (6): 811-27. 10.1002/sim.4780100604. discussion 828-9View ArticlePubMedGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2288/11/21/prepub

### Pre-publication history

## Copyright

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.