ST 520 Statistical Principles of Clinical Trials - NCSU Statistics ...
ST 520 Statistical Principles of Clinical Trials - NCSU Statistics ...
ST 520 Statistical Principles of Clinical Trials - NCSU Statistics ...
You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
CHAPTER 9 <strong>ST</strong> <strong>520</strong>, A. TSIATIS and D. Zhang<br />
is the value <strong>of</strong> the non-centrality parameter necessary so that a non-central chi-square distributed<br />
random variable with K − 1 degrees <strong>of</strong> freedom and non-centrality parameter φ 2 (α, β, K − 1)<br />
will exceed the value χ 2 α;K−1 with probability (1 − β). Tables <strong>of</strong> φ2 (α, β, K − 1) for α = .05 were<br />
provided in chapter 7.<br />
For example, if we take K = 3, then in order to ensure that we have at least 90% power to detect<br />
a hazard ratio between any two treatments that may exceed 1.5, using a logrank test at the .05<br />
level <strong>of</strong> significance, we would need the total number <strong>of</strong> deaths to exceed<br />
d =<br />
2 × 3 × 12.654<br />
= 462.<br />
{log(1.5)} 2<br />
We can contrast this to a two-sample comparison which needs 256 events. As in the two-sample<br />
problem, the computations during the design stage will involve the best guesses for the accrual<br />
rate, accrual period, follow-up period, and underlying treatment-specific survival distributions<br />
which can be translated to the desired number <strong>of</strong> failures. Thus we can experiment with different<br />
values <strong>of</strong><br />
• accrual rate a(u)<br />
• underlying treatment-specific failure time distributions Fj(t) = P(T ≤ t|A = j) = 1 −<br />
Sj(t), j = 1, . . .,K under the alternative hypothesis <strong>of</strong> interest (we may take these at the<br />
least favorable configuration)<br />
• the accrual period A<br />
• the length <strong>of</strong> study L<br />
so that<br />
K�<br />
j=1<br />
dj{a(·), Fj(·), A, L} = 2Kφ2 (α, β, K − 1)<br />
γ2 ,<br />
A<br />
where dj{a(·), Fj(·), A, L} denotes the expected number <strong>of</strong> deaths in treatment group j as a<br />
function <strong>of</strong> a(·), Fj(·), A, L, computed using equation (9.5).<br />
CALGB 8541 Example<br />
We now return to the data from CALGB 8541 which compared three treatments in a randomized<br />
study <strong>of</strong> node positive stage II breast cancer patients. The three treatments were<br />
PAGE 161