Evaluation of the type I error rate when using parametric bootstrap analysis of a cluster randomized controlled trial with binary outcomes and a small number of clusters
Lilian Golzarri-Arroyo,Stephanie L Dickinson,Yasaman Jamshidi-Naeini,Roger S Zoh,Andrew W Brown,Arthur H Owora,Peng Li,J Michael Oakes,David B Allison,Stephanie L. Dickinson,Roger S. Zoh,Andrew W. Brown,Arthur H. Owora,J. Michael Oakes,David B. Allison
DOI: https://doi.org/10.1016/j.cmpb.2022.106654
IF: 6.1
2022-03-01
Computer Methods and Programs in Biomedicine
Abstract:BACKGROUND: Cluster randomized controlled trials (cRCTs) are increasingly used but must be analyzed carefully. We conducted a simulation study to evaluate the validity of a parametric bootstrap (PB) approach with respect to the empirical type I error rate for a cRCT with binary outcomes and a small number of clusters.METHODS: We simulated a case study with a binary (0/1) outcome, four clusters, and 100 subjects per cluster. To compare the validity of the test with respect to error rate, we simulated the same experiment with K=10, 20, and 30 clusters, each with 2,000 simulated datasets. To test the null hypothesis, we used a generalized linear mixed model including a random intercept for clusters and obtained p-values based on likelihood ratio tests (LRTs) using the parametric bootstrap method as implemented in the R package "pbkrtest".RESULTS: The PB test produced error rates of 9.1%, 5.5%, 4.9%, and 5.0% on average across all ICC values for K=4, K=10, K=20, and K=30, respectively. The error rates were higher, ranging from 9.1% to 36.5% for K=4, in the models with singular fits (i.e., ignoring clustering) because the ICC was estimated to be zero.CONCLUSION: Using the parametric bootstrap for cRCTs with a small number of clusters results in inflated error rates and is not valid.
engineering, biomedical,computer science, interdisciplinary applications,medical informatics, theory & methods