Bayesian Sample size Determination for Longitudinal Studies with Continuous Response using Marginal Models

Author
Taban Baghfalaki
Tarbiat Modares University
10.29252/mmr.4.1.1
Abstract
Introduction

Longitudinal study designs are common in a lot of scientific researches, especially in medical, social and economic sciences. The reason is that longitudinal studies allow researchers to measure changes of each individual over time and often have higher statistical power than cross-sectional studies. Choosing an appropriate sample size is a crucial step in a successful study. A study with insufficient sample size may have a small statistical power to detect significant effects and may lead to incorrect answers to many important research questions. On the other hand, a study with more than what a sample size should be wastes the resources and budget. In longitudinal studies, because of the complexity of the design of experiment, including the selection of the number of individuals and the number of repeated measurements, the sample size determination is less studied. This paper uses a simulation-based method to determine the sample size from a Bayesian perspective. For this purpose several Bayesian criteria of sample size determination for a longitudinal study using marginal model are used. Most of the methods of determining the sample size are based on creation of one hypothesis. In this paper, in addition to using this method, we also present a method to determine the sample size for multiple hypothesis testing. Using several examples the proposed Bayesian methods are discussed.

Material and method

Based on the Bayesian sample size determination approach proposed by Wang and Gelfand (2002), which is a simulation-based approach, two sets of priors are considered. The first set is called the “fitting” or “analysis” priors which are used for analyzing data. The other set of priors is called “design” or “sampling” sets which drawing upon expertise, we may speculate upon a variety of informative scenarios regarding the unknown parameters and capture each with a suitable sampling prior. In the first step of the simulation-based approach, one generates the parameters from the design priors then a sample from the generated parameters is obtained. Note that the frequentist/classical viewpoint of the sample size determination approach require a point estimate of the variance and the smallest meaningful difference (Diggle et al., 2002). This information is elicited based on pilot data or the opinion of experts. The use of design priors in the proposed approach of Wang and Gelfand (2002) replaces this part of the classical sample size determination approach. After generating data using the design priors in this stage, they are analyzed by the fitting priors. In this paper, non-informative priors are used as fitting priors, as we assume there is no prior information or any expert opinion to be used to construct a prior distribution. If there is any prior information to be used to elicit an appropriate prior distribution, one may use it in the step of fitting prior and for sure, if the specification of the prior is correct, will achieve a better result. Also, in this paper, four criteria are used for determining sample size in longitudinal studies: Bayesian power criterion (BPC), average length criterion (ALC), average posterior variance criterion (APVC) and average converge criterion (ACC).

The simulation-based approach of Wang and Gelfand (2002) can be summarized as follows:


Specify the value of the effect of an important regression coefficient which is the interest parameter and specify the design and fitting priors.
For each sample size, the following steps are repeated M times:

Generate values of the unknown parameters from their design priors.
Simulate the values of covariates from continuous or discrete distributions and the response variable from its distribution.
Analyze the generated data set of step (ii) using the fitting priors.
Calculate BPC.


Fit a curve or surface through Bayesian power values and find an adequate sample size for a desired power using interpolation. In this paper, the curve is fitted using a polynomial regression.

In step 2(iv), one can calculate ACC, APVC, or ALC and determine the sample size based on these criteria.

References


Diggle, P., Heagerty, P., Liang, K. and Zeger, S. (2002). Analysis of Longitudinal Data. Oxford:Oxford University Press.
Wang, F., and Gelfand, A. E. (2002). A simulation-based approach to Bayesian sample size determination for performance under a given model and for separating models. Statistical Science, 193-208../files/site1/files/41/1Extended_Abstract.pdf
Keywords
Bayesian analysis
Longitudinal data
marginal model
posterior distribution
Bayesian power criterion
1. Muller K. E., Lavange L. M., Ramey S. L., Ramey C. T. "Power calculations for general linear multivariate models including repeated measures applications", Journal of the American Statistical Association, 87 (420) (1992) 1209-1226. 2. Muller K. E., Stewart P. W. "Linear model theory: univariate, multivariate, and mixed models", John Wiley & Sons, Canada (2006). 3. Guo Y., Logan H. L., Glueck D. H., Muller K. E. "Selecting a sample size for studies with repeated measures", BMC medical research methodology, 13(1) (2013) 1. 4. Gurka M. J., Edwards L. J., Muller K. E. "Avoiding bias in mixed model inference for fixed effects", Statistics in medicine, 30 (22) (2011) 2696-2707. 5. Cheng J., Edwards L. J., Maldonado Molina, M. M., Komro K. A., Muller K. E. "Real longitudinal data analysis for real people: building a good enough mixed model", Statistics in medicine, 29 (4) (2010) 504-520. 6. Muller K. E., Edwards L. J., Simpson S. L., Taylor D. J. "Statistical tests with accurate size and power for balanced linear mixed models", Statistics in Medicine, 26 (19) (2007) 3639-3660. 7. Johnson J. L., Muller K. E., Slaughter J. C., Gurka M. J., Gribbin M. J., Simpson S. L. "POWERLIB: SAS/IML software for computing power in multivariate linear models", Journal of statistical software, 30 (5) (2009). 8. Liang K. Y., Zeger S. L. "Longitudinal data analysis using generalized linear models", Biometrika, 73 (1) (1986) 13-22. 9. Liu G., Liang K. Y., "Sample size calculations for studies with correlated observations", Biometrics (1997) 937-947. 10. Overall J. E., Doyle S. R. "Estimating sample sizes for repeated measurement designs", Controlled clinical trials, 15 (2) (1994) 100-123. 11. Hedeker D., Gibbons R. D., Waternaux C. "Sample size estimation for longitudinal designs with attrition: comparing time-related contrasts between two groups", Journal of Educational and Behavioral Statistics, 24 (1) (1999) 70-93. 12. Comulada W. S., Weiss R. E. "Sample size and power calculations for correlations between bivariate longitudinal data", Statistics in medicine, 29 (27) (2010) 2811-2824. 13. Hintze J. L. "Power Analysis and Sample Size (PASS) for Windows User’s Guide", Kaysville, Utah, NCSS (2000). 14. Adcock C. J. "Sample size determination: a review", Journal of the Royal Statistical Society: Series D (The Statistician), 46(2) (1997) 261-283. 15. Ahn C., Heo M., Zhang S., "Sample size calculations for clustered and longitudinal outcomes in clinical research", CRC Press. Boca Raton (2014). 16. Joseph L., Belisle P., "Bayesian sample size determination for normal means and differences between normal means", Journal of the Royal Statistical Society: Series D (The Statistician), 46(2) (1997) 209-226. 17. Rahme E., Joseph L., Gyorkos T. W., "Bayesian sample size determination for estimating binomial parameters from data subject to misclassification", Journal of the Royal Statistical Society: Series C (Applied Statistics), 49 (1) (2000) 119-128. 18. Wang F., Gelfand A. E., "A simulation-based approach to Bayesian sample size determination for performance under a given model and for separating models", Statistical Science (2002) 193-208. 19. Brutti P., De Santis F., Gubbiotti S., "Robust Bayesian sample size determination in clinical trials", Statistics in Medicine, 27 (13) (2008) 2290-2306. 20. Pham-Gia T., Turkkan N., "Sample size determination in Bayesian analysis", The Statistician (1992) 389-397. 21. Joseph L., Wolfson D. B., Berger R., "Some comments on Bayesian sample size determination", The Statistician (1995) 167-171. 22. Pezeshk H., "Bayesian techniques for sample size determination in clinical trials: a short review", Statistical Methods in Medical Research, 12 (6) (2003) 489-504. 23. Christensen R., Johnson W., Branscum A., Hanson T., "Bayesian Ideas and Data Analysis", CHAPMAN & HALL/CRC, Boca Raton (2012). 24. Joseph L., Berger R., Belisle P., "Bayesian and mixed Bayesian/likelihood criteria for sample size determination", Statistics in medicine, 16 (7) (1997) 769-781. 25. Gelman A., Rubin D. B., "Inference from iterative simulation using multiple sequences. Statistical Science, 7 (1992) 457-511. 26. Diggle P. J., Heagerty P., Liang K. Y., Zeger S. L., "Analysis of longitudinal data", Oxford University Press (2002).