rmdiff | R Documentation |
Obtains the estimate of restricted mean survival time difference between two treatment groups.
rmdiff(
data,
rep = "",
stratum = "",
treat = "treat",
time = "time",
event = "event",
milestone = NA_real_,
rmstDiffH0 = 0,
conflev = 0.95,
biascorrection = 0L
)
data |
The input data frame that contains the following variables:
|
rep |
The name of the replication variable in the input data. |
stratum |
The name of the stratum variable in the input data. |
treat |
The name of the treatment variable in the input data. |
time |
The name of the time variable in the input data. |
event |
The name of the event variable in the input data. |
milestone |
The milestone time at which to calculate the restricted mean survival time. |
rmstDiffH0 |
The difference in restricted mean survival times under the null hypothesis. Defaults to 0 for superiority test. |
conflev |
The level of the two-sided confidence interval for the difference in restricted mean survival times. Defaults to 0.95. |
biascorrection |
Whether to apply bias correction for the variance estimate of individual restricted mean survival times. Defaults to no bias correction. |
A data frame with the following variables:
rep
: The replication number.
milestone
: The milestone time relative to randomization.
rmstDiffH0
: The difference in restricted mean survival times
under the null hypothesis.
rmst1
: The estimated restricted mean survival time for
the treatment group.
rmst2
: The estimated restricted mean survival time for
the control group.
rmstDiff
: The estimated difference in restricted mean
survival times.
vrmst1
: The variance for rmst1.
vrmst2
: The variance for rmst2.
vrmstDiff
: The variance for rmstDiff.
rmstDiffZ
: The Z-statistic value.
rmstDiffPValue
: The one-sided p-value.
lower
: The lower bound of confidence interval.
upper
: The upper bound of confidence interval.
conflev
: The level of confidence interval.
biascorrection
: Whether to apply bias correction for the
variance estimate of individual restricted mean survival times.
Kaifeng Lu, kaifenglu@gmail.com
df <- rmdiff(data = rawdata, rep = "iterationNumber",
stratum = "stratum", treat = "treatmentGroup",
time = "timeUnderObservation", event = "event",
milestone = 12)
head(df)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.