Confidence regions for repeated measures ANOVA power curves based on estimated covariance

BackgroundUsing covariance or mean estimates from previous data introduces randomness into each power value in a power curve. Creating confidence intervals about the power estimates improves study planning by allowing scientists to account for the uncertainty in the power estimates. Driving examples arise in many imaging applications.MethodsWe use both analytical and Monte Carlo simulation methods. Our analytical derivations apply to power for tests with the univariate approach to repeated measures (UNIREP). Approximate confidence intervals and regions for power based on an estimated covariance matrix and fixed means are described. Extensive simulations are used to examine the properties of the approximations.ResultsClosed-form expressions are given for approximate power and confidence intervals and regions. Monte Carlo simulations support the accuracy of the approximations for practical ranges of sample size, rank of the design matrix, error degrees of freedom, and the amount of deviation from sphericity. The new methods provide accurate coverage probabilities for all four UNIREP tests, even for small sample sizes. Accuracy is higher for higher power values than for lower power values, making the methods especially useful in practical research conditions. The new techniques allow the plotting of power confidence regions around an estimated power curve, an approach that has been well received by researchers. Free software makes the new methods readily available.ConclusionsThe new techniques allow a convenient way to account for the uncertainty of using an estimated covariance matrix in choosing a sample size for a repeated measures ANOVA design. Medical imaging and many other types of healthcare research often use repeated measures ANOVA.

[1]  M. Dagenais Comment on Zeger and Brookmeyer , 1991 .

[2]  R. H. Browne On the use of a pilot sample for sample size determination. , 1995, Statistics in medicine.

[3]  M. Gribbin,et al.  Better power methods for the univariate approach to repeated measures , 2007 .

[4]  B. Lecoutre A Correction for the ε̃ Approximate Test in Repeated Measures Designs With Two or More Independent Groups , 1991 .

[5]  Richard F. Gunst,et al.  Regression and ANOVA: An Integrated Approach Using SAS Software , 2003, Technometrics.

[6]  Douglas J. Taylor,et al.  BIAS IN LINEAR MODEL POWER AND SAMPLE SIZE CALCULATION DUE TO ESTIMATING NONCENTRALITY. , 1996, Communications in statistics: theory and methods.

[7]  S. Geisser,et al.  On methods in the analysis of profile data , 1959 .

[8]  Keith E. Muller,et al.  Linear Model Theory: Univariate, Multivariate, and Mixed Models , 2006 .

[9]  K. Muller,et al.  Statistical tests with accurate size and power for balanced linear mixed models , 2007, Statistics in medicine.

[10]  G. Box Some Theorems on Quadratic Forms Applied in the Study of Analysis of Variance Problems, I. Effect of Inequality of Variance in the One-Way Classification , 1954 .

[11]  S. Maxwell The persistence of underpowered studies in psychological research: causes, consequences, and remedies. , 2004, Psychological methods.

[12]  R. Davies The distribution of a linear combination of 2 random variables , 1980 .

[13]  J. Cavanaugh Biostatistics , 2005, Definitions.

[14]  Keith E Muller,et al.  Properties of internal pilots with the univariate approach to repeated measures , 2003, Statistics in medicine.

[15]  H. Kraemer,et al.  Caution regarding the use of pilot studies to guide power calculations for study proposals. , 2006, Archives of general psychiatry.

[16]  H. Huynh,et al.  Estimation of the Box Correction for Degrees of Freedom from Sample Data in Randomized Block and Split-Plot Designs , 1976 .

[17]  Keith E Muller,et al.  Computing Confidence Bounds for Power and Sample Size of the General Linear Univariate Model. , 1995, The American statistician.

[18]  G. Box Some Theorems on Quadratic Forms Applied in the Study of Analysis of Variance Problems, II. Effects of Inequality of Variance and of Correlation Between Errors in the Two-Way Classification , 1954 .

[19]  D. F. Morrison,et al.  Multivariate Statistical Methods , 1968 .

[20]  Keith E Muller,et al.  Analytic, Computational, and Approximate Forms for Ratios of Noncentral and Central Gaussian Quadratic Forms , 2006, Journal of computational and graphical statistics : a joint publication of American Statistical Association, Institute of Mathematical Statistics, Interface Foundation of North America.

[21]  S. Geisser,et al.  An Extension of Box's Results on the Use of the $F$ Distribution in Multivariate Analysis , 1958 .

[22]  Keith E. Muller,et al.  Approximate Power for Repeated-Measures ANOVA Lacking Sphericity , 1989 .

[23]  Sean L Simpson,et al.  POWERLIB: SAS/IML Software for Computing Power in Multivariate Linear Models. , 2009, Journal of statistical software.

[24]  Keith E. Muller,et al.  BIAS IN LINEAR MODEL POWER AND SAMPLE SIZE DUE TO ESTIMATING VARIANCE. , 1997, Communications in statistics: theory and methods.

[25]  Yueh-Yun Chi,et al.  Global hypothesis testing for high‐dimensional repeated measures outcomes , 2012, Statistics in medicine.

[26]  J. Wishart THE GENERALISED PRODUCT MOMENT DISTRIBUTION IN SAMPLES FROM A NORMAL MULTIVARIATE POPULATION , 1928 .