ssizeCT {powerSurvEpi} | R Documentation |
Sample Size Calculation in the Analysis of Survival Data for Clinical Trials
Description
Sample size calculation for the Comparison of Survival Curves Between Two Groups under the Cox Proportional-Hazards Model for clinical trials. Some parameters will be estimated based on a pilot data set.
Usage
ssizeCT(formula,
dat,
power,
k,
RR,
alpha = 0.05)
Arguments
formula |
A formula object, e.g. |
dat |
a data frame representing the pilot data set and containing at least 3 columns: (1) survival/censoring time; (2) censoring indicator;
(3) group indicator which is a factor object in R and takes only two possible values ( |
power |
numeric. power to detect the magnitude of the hazard ratio as small as that specified by |
k |
numeric. ratio of participants in group E (experimental group) compared to group C (control group). |
RR |
numeric. postulated hazard ratio. |
alpha |
numeric. type I error rate. |
Details
This is an implementation of the sample size calculation method described in Section 14.12 (page 807) of Rosner (2006). The method was proposed by Freedman (1982).
The movitation of this function is that some times we do not have information about or
and
available, but we have a pilot data set that can be used to estimate
and
hence
, where
is the expected total number of events over both groups,
and
are numbers of participants in group E (experimental group) and group C (control group), respectively.
is the probability of failure in group E (experimental group) over the maximum time period of the study (t years).
is the probability of failure in group C (control group) over the maximum time period of the study
(t years).
Suppose we want to compare the survival curves between an experimental group () and
a control group (
) in a clinical trial with a maximum follow-up of
years.
The Cox proportional hazards regression model is assumed to have the form:
Let be the number of participants in the
group
and
be the number of participants in the
group.
We wish to test the hypothesis
versus
not equal to 1,
where
underlying hazard ratio
for the
group versus the
group. Let
be the postulated hazard ratio,
be the significance level. Assume that the test is a two-sided test.
If the ratio of participants in group
E compared to group C
, then the number of participants needed in each group to
achieve a power of
is
where
and
is the
-th percentile of
the standard normal distribution
.
and
can be calculated from the following formulaes:
where ,
,
,
,
. And
is the probability of failure at time
i
among participants in the
control group, given that a participant has survived to time and is not censored at time
,
i.e., the approximate hazard time
in the control group,
;
is the probability of failure at time
i
among participants in the
experimental group, given that a participant has survived to time and is not censored at time
,
i.e., the approximate hazard time
in the experimental group,
;
is the prbability that a participant is censored at time
given that he was
followed up to time
and has not failed,
, which is assumed the same in each group.
Value
mat.lambda |
a matrix with 9 columns and |
mat.event |
a matrix with 5 columns and |
pC |
estimated probability of failure in group C (control group) over the maximum time period of the study (t years). |
pE |
estimated probability of failure in group E (experimental group) over the maximum time period of the study (t years). |
ssize |
a two-element vector. The first element is |
Note
(1) The estimates of . That is,
RRlambda
is not directly estimated based on data
from the experimental group;
(2) The sample size formula assumes that the central-limit theorem is valid and hence is appropriate for large samples.
(3) and
will be rounded up to integers.
References
Freedman, L.S. (1982). Tables of the number of patients required in clinical trials using the log-rank test. Statistics in Medicine. 1: 121-129
Rosner B. (2006). Fundamentals of Biostatistics. (6-th edition). Thomson Brooks/Cole.
See Also
Examples
# Example 14.42 in Rosner B. Fundamentals of Biostatistics.
# (6-th edition). (2006) page 809
library(survival)
data(Oph)
res <- ssizeCT(formula = Surv(times, status) ~ group,
dat = Oph,
power = 0.8,
k = 1,
RR = 0.7,
alpha = 0.05)
# Table 14.24 on page 809 of Rosner (2006)
print(round(res$mat.lambda, 4))
# Table 14.12 on page 787 of Rosner (2006)
print(round(res$mat.event, 4))
# the sample size
print(res$ssize)