Description Usage Arguments Details Value Author(s) References See Also Examples

`nSurvival()`

is used to calculate the sample size for a clinical trial with a time-to-event endpoint. The Lachin and Foulkes (1986) method is used.
`nEvents`

uses the Schoenfeld (1981) approximation to provide sample size and power in terms of the underlying hazard ratio and the number of events observed in a survival analysis.
The functions `hrz2n()`

, `hrn2z()`

and `zn2hr()`

also use the Schoenfeld approximation to provide simple translations between hazard ratios, z-values and the number of events in an analysis; input variables can be given as vectors.

1 2 3 4 5 6 7 8 9 10 | ```
nSurvival(lambda1=1/12, lambda2=1/24, Ts=24, Tr=12, eta = 0, ratio = 1,
alpha = 0.025, beta = 0.10, sided = 1, approx = FALSE,
type = c("rr", "rd"), entry = c("unif", "expo"), gamma = NA)
## S3 method for class 'nSurvival'
print(x,...)
nEvents(hr = .6, alpha = .025, beta = .1, ratio = 1, sided = 1,
hr0 = 1, n = 0, tbl = FALSE)
hrn2z(hr, n, ratio=1, hr0=1, hr1=.7)
hrz2n(hr, z, ratio=1, hr0=1)
zn2hr(z, n, ratio=1, hr0=1, hr1=.7)
``` |

`lambda1, lambda2` |
event hazard rate for placebo and treatment group respectively. |

`eta` |
equal dropout hazard rate for both groups. |

`ratio` |
randomization ratio between placebo and treatment group. Default is balanced design, i.e., randomization ratio is 1. |

`Ts` |
maximum study duration. |

`Tr` |
accrual (recruitment) duration. |

`alpha` |
type I error rate. Default is 0.025 since 1-sided testing is default. |

`beta` |
type II error rate. Default is 0.10 (90% power). Not needed for |

`sided` |
one or two-sided test? Default is one-sided test. |

`approx` |
logical. If |

`type` |
type of sample size calculation: risk ratio (“rr”) or risk difference (“rd”). |

`entry` |
patient entry type: uniform entry ( |

`gamma` |
rate parameter for exponential entry. |

`x` |
An object of class "nSurvival" returned by |

`hr` |
Hazard ratio. For |

`hr0` |
Hazard ratio under the null hypothesis (>0, for |

`hr1` |
Hazard ratio under the alternate hypothesis for |

`n` |
Number of events. For |

`tbl` |
Indicator of whether or not scalar (vector) or tabular output is desired for |

`z` |
A z-statistic. |

`...` |
Allows additional arguments for |

`nSurvival()`

produces an object of class "nSurvival" with
the number of subjects and events for a set of
pre-specified trial parameters, such as accrual duration and follow-up
period. The calculation is based on Lachin and Foulkes (1986) method and can
be used for risk ratio or risk difference. The function also consider
non-uniform (exponential) entry as well as uniform entry.

If the logical `approx`

is `TRUE`

, the variance under alternative
hypothesis is used to replace the variance under null hypothesis.
For non-uniform entry, a non-zero value of `gamma`

for exponential entry
must be supplied. For positive `gamma`

, the entry distribution is
convex, whereas for negative `gamma`

, the entry distribution is concave.

`nEvents()`

uses the Schoenfeld (1981) method to approximate the number of events `n`

(given `beta`

) or the power (given `n`

). Arguments may be vectors or scalars, but any vectors must have the same length.

The functions `hrz2n`

, `hrn2z`

and `zn2hr`

also all apply the Schoenfeld approximation for proportional hazards modelling.
This approximation is based on the asymptotic normal distribtuion of the logrank statistic as well as related statistics are asymptotically normal.
Let *λ* denote the underlying hazard ratio (`lambda1/lambda2`

in terms of the arguments to `nSurvival`

). Further, let *n* denote the number of events observed when computing the statistic of interest and *r* the ratio of the sample size in an experimental group relative to a control. The estimated natural logarithm of the hazard ratio from a proportional hazards ratio is approximately normal with a mean of *log{λ}* and variance *(1+r)^2/nr*.
Let *z* denote a logrank statistic (or a Wald statistic or score statistic from a proportional hazards regression model).
The same asymptotic theory implies *z* is asymptotically equivalent to a normalized estimate of the hazard ratio *λ* and thus *z* is asymptotically normal with variance 1 and mean

*\frac{log{λ}r}{(1+r)^2}.*

Plugging the estimated hazard ratio into the above equation allows approximating any one of the following based on the other two: the estimate hazard ratio, the number of events and the z-statistic. That is,

*\hat{λ}= \exp(z(1+r)/√{rn})*

*z=log(\hat{λ})√{nr}/(1+r)*

*n= (z(1+r)/log(\hat{λ}))^2/r.*

`hrz2n()`

translates an observed interim hazard ratio and interim z-value into the number of events required for the Z-value and hazard ratio to correspond to each other. `hrn2z()`

translates a hazard ratio and number of events into an approximate corresponding Z-value. `zn2hr()`

translates a Z-value and number of events into an approximate corresponding hazard ratio. Each of these functions has a default assumption of an underlying hazard ratio of 1 which can be changed using the argument `hr0`

. `hrn2z()`

and `zn2hr()`

also have an argument `hr1`

which is only used to compute the sign of the computed Z-value in the case of `hrn2z()`

and whether or not a z-value > 0 corresponds to a hazard ratio > or < the null hazard ratio `hr0`

.

`nSurvival`

produces a list with the following component returned:

`type` |
As input. |

`entry` |
As input. |

`n` |
Sample size required (computed). |

`nEvents` |
Number of events required (computed). |

`lambda1` |
As input. |

`lambda2` |
As input. |

`eta` |
As input. |

`ratio` |
As input. |

`gamma` |
As input. |

`alpha` |
As input. |

`beta` |
As input. |

`sided` |
As input. |

`Ts` |
As input. |

`Tr` |
As input. |

`nEvents`

produces a scalar or vector of sample sizes (or powers) when `tbl=FALSE`

or, when `tbl=TRUE`

a data frame of values with the following columns:

`hr` |
As input. |

`n` |
If |

`alpha` |
As input. |

`beta` |
If |

`Power` |
One minus the output |

`delta` |
Standardized effect size represented by input difference between null and alternative hypothesis hazard ratios. |

`ratio` |
Ratio of experimental to control sample size where 'experimental' is the same as the group with hazard represented in the numerator of the hazard ratio. |

`se` |
Estimated standard error for the observed log(hazard ratio) with the given sample size. |

`hrz2n`

outputs a number of events required to approximately have the input hazard ratio, z-statistic and sample size correspond.
`hrn2z`

outputs an approximate z-statistic corresponding to an input hazard ratio and number of events.
`zn2hr`

outputs an approximate hazard ratio corresponding to an input z-statistic and number of events.

Shanhong Guan [email protected], Keaven Anderson [email protected]

Lachin JM and Foulkes MA (1986),
Evaluation of Sample Size and Power for Analyses of Survival
with Allowance for Nonuniform Patient Entry, Losses to Follow-Up,
Noncompliance, and Stratification. *Biometrics*, 42, 507-519.

Schoenfeld D (1981),
The Asymptotic Properties of Nonparametric Tests for Comparing Survival Distributions. *Biometrika*,
68, 316-319.

gsDesign package overview, Plots for group sequential designs, gsDesign, gsHR

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 | ```
# consider a trial with
# 2 year maximum follow-up
# 6 month uniform enrollment
# Treatment/placebo hazards = 0.1/0.2 per 1 person-year
# drop out hazard 0.1 per 1 person-year
# alpha = 0.025 (1-sided)
# power = 0.9 (default beta=.1)
ss <- nSurvival(lambda1=.2 , lambda2=.1, eta = .1, Ts = 2, Tr = .5,
sided=1, alpha=.025)
# group sequential translation with default bounds
# note that delta1 is log hazard ratio; used later in gsBoundSummary summary
x<-gsDesign(k = 5, test.type = 2, n.fix=ss$nEvents, nFixSurv=ss$n,
delta1=log(ss$lambda2/ss$lambda1))
# boundary plot
plot(x)
# effect size plot
plot(x, plottype = "hr")
# total sample size
x$nSurv
# number of events at analyses
x$n.I
# print the design
x
# overall design summary
cat(summary(x))
# tabular summary of bounds
gsBoundSummary(x,deltaname="HR",Nname="Events",logdelta=TRUE)
# approximate number of events required using Schoenfeld's method
# for 2 different hazard ratios
nEvents(hr=c(.5, .6), tbl=TRUE)
# vector output
nEvents(hr=c(.5, .6))
# approximate power using Schoenfeld's method
# given 2 sample sizes and hr=.6
nEvents(hr=.6, n=c(50, 100), tbl=TRUE)
# vector output
nEvents(hr=.6, n=c(50, 100))
# approximate hazard ratio corresponding to 100 events and z-statistic of 2
zn2hr(n=100,z=2)
# same when hr0 is 1.1
zn2hr(n=100,z=2,hr0=1.1)
# same when hr0 is .9 and hr1 is greater than hr0
zn2hr(n=100,z=2,hr0=.9,hr1=1)
# approximate number of events corresponding to z-statistic of 2 and
# estimated hazard ratio of .5 (or 2)
hrz2n(hr=.5,z=2)
hrz2n(hr=2,z=2)
# approximate z statistic corresponding to 75 events
# and estimated hazard ratio of .6 (or 1/.6)
# assuming 2-to-1 randomization of experimental to control
hrn2z(hr=.6,n=75,ratio=2)
hrn2z(hr=1/.6,n=75,ratio=2)
``` |

gsDesign documentation built on May 31, 2017, 2:15 a.m.

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.