Sample size determination in clinical trials with multiple co‐primary endpoints including mixed continuous and binary variables

In the field of pharmaceutical drug development, there have been extensive discussions on the establishment of statistically significant results that demonstrate the efficacy of a new treatment with multiple co-primary endpoints. When designing a clinical trial with such multiple co-primary endpoints, it is critical to determine the appropriate sample size for indicating the statistical significance of all the co-primary endpoints with preserving the desired overall power because the type II error rate increases with the number of co-primary endpoints. We consider overall power functions and sample size determinations with multiple co-primary endpoints that consist of mixed continuous and binary variables, and provide numerical examples to illustrate the behavior of the overall power functions and sample sizes. In formulating the problem, we assume that response variables follow a multivariate normal distribution, where binary variables are observed in a dichotomized normal distribution with a certain point of dichotomy. Numerical examples show that the sample size decreases as the correlation increases when the individual powers of each endpoint are approximately and mutually equal.

[1]  Toshimitsu Hamasaki,et al.  A convenient formula for sample size calculations in clinical trials with multiple co‐primary continuous endpoints , 2012, Pharmaceutical statistics.

[2]  G Molenberghs,et al.  Evaluation of surrogate endpoints in randomized experiments with mixed discrete and continuous outcomes , 2001, Statistics in medicine.

[3]  Robert F. Tate,et al.  Correlation Between a Discrete and a Continuous Variable. Point-Biserial Correlation , 1954 .

[4]  Bartlett Ms The use of transformations. , 1947 .

[5]  R. Prentice,et al.  Correlated binary regression with covariates specific to each binary observation. , 1988, Biometrics.

[6]  A. R. de Leon Pairwise likelihood approach to grouped continuous model and its extension , 2005 .

[7]  Toshimitsu Hamasaki,et al.  Sample size determination in clinical trials with multiple co‐primary binary endpoints , 2010, Statistics in medicine.

[8]  Feng Gao,et al.  Power and sample size for clinical trials when efficacy is required in multiple endpoints: application to an Alzheimer's treatment trial , 2005, Clinical trials.

[9]  Sik-Yum Lee,et al.  Maximum likelihood estimation of multivariate polyserial and polychoric correlation coefficients , 1987 .

[10]  J. Dale Global cross-ratio models for bivariate, discrete, ordered responses. , 1986, Biometrics.

[11]  R. Muirhead,et al.  Multiple Co-primary Endpoints: Medical and Statistical Solutions: A Report from the Multiple Endpoints Expert Team of the Pharmaceutical Research and Manufacturers of America , 2007 .

[12]  Hamid Pezeshk,et al.  Sample Size Determination in Clinical Trials , 1999 .

[13]  F. Yates Contingency Tables Involving Small Numbers and the χ2 Test , 1934 .

[14]  M. Piedmonte,et al.  A Method for Generating High-Dimensional Multivariate Binary Variates , 1991 .

[15]  R. Tate Applications of Correlation Models for Biserial Data , 1955 .

[16]  Toshimitsu Hamasaki,et al.  Sample Size Determination in Superiority Clinical Trials with Multiple Co-Primary Correlated Endpoints , 2011, Journal of biopharmaceutical statistics.

[17]  Isao Yoshimura,et al.  Power and Sample Size Calculations in Clinical Trials with Multiple Primary Variables , 2006 .

[18]  Karl Pearson,et al.  ON A NEW METHOD OF DETERMINING CORRELATION BETWEEN A MEASURED CHARACTER A, AND A CHARACTER B, OF WHICH ONLY THE PERCENTAGE OF CASES WHEREIN B EXCEEDS (OR FALLS SHORT OF) A GIVEN INTENSITY IS RECORDED FOR EACH GRADE OF A , 1909 .

[19]  Sue-Jane Wang,et al.  Some Controversial Multiple Testing Problems in Regulatory Applications , 2009, Journal of biopharmaceutical statistics.

[20]  Joseph Lev,et al.  The Point Biserial Coefficient of Correlation , 1949 .

[21]  James X. Song Sample size for simultaneous testing of rate differences in non-inferiority trials with multiple endpoints , 2009, Comput. Stat. Data Anal..

[22]  W. Laskey,et al.  Nephrotoxicity of iodixanol versus ioversol in patients with chronic kidney disease: the Visipaque Angiography/Interventions with Laboratory Outcomes in Renal Insufficiency (VALOR) Trial. , 2008, American heart journal.

[23]  F. Breedveld,et al.  The PREMIER study: A multicenter, randomized, double-blind clinical trial of combination therapy with adalimumab plus methotrexate versus methotrexate alone or adalimumab alone in patients with early, aggressive rheumatoid arthritis who had not had previous methotrexate treatment. , 2006, Arthritis and rheumatism.

[24]  R. Plackett A Class of Bivariate Distributions , 1965 .

[25]  A. Genz Numerical Computation of Multivariate Normal Probabilities , 1992 .

[26]  D. E. Walters In Defence of the Arc Sine Approximation , 1979 .

[27]  Morris L. Eaton,et al.  On a multiple endpoints testing problem , 2007 .

[28]  Robert F. Tate,et al.  THE THEORY OF CORRELATION BETWEEN TWO CONTINUOUS VARIABLES WHEN ONE IS DICHOTOMIZED , 1955 .