WRSEestimate {DTR}R Documentation

Function for calculating WRSE estimates


The function computes the weighted risk set estimates (WRSE) of the survival functions and their estimated standard errors for dynamic treatment regimes (DTRs) at observed event times as proposed in Guo and Tsiatis (2005) Equation (3) and Equation (16).





a data frame (X, TR, R, Z, U, delta) representing the data from a two-stage randomized designs with therapies A1 and A2 available at the first stage, and B1 and B2 available at the second stage.
X: treatment arm, X=0 for A1, and X=1 for A2
TR: the time to response for responders (R=1)
R: response status, R=1 for responders, and R=0 for non-responders
Z: second-stage indicator, Z=0 if assigned to B1, and Z=1 if assigned to B2
U: observed survival time, U is death time if delta=1, and U is censoring time if delta=0
delta: censoring indicator, delta=1 for died, and delta=0 for censored


In sequentially randomized designs, there could be more than two therapies available at each stage. For simplicity, and to maintain similarity to the most common sequentially randomized clinical trials, a simple two-stage randomization design allowing two treatment options at each stage is used in the current version of the package. In detail, patients are initially randomized to either A1 or A2 at the first stage. Based on their response status, they are then randomized to either B1 or B2 at the second stage. Therefore, there are a total of four DTRs: A1B1, A1B2, A2B1, and A2B2.


The function returns an object of class DTR. See DTR.object for details.


Guo X, Tsiatis AA: A weighted risk set estimator for survival distributions in two-stage randomization designs with censored survival data. Int. J. Biostatistics 1:1-15, 2005

See Also

simWRSEdata, DTR.object, print.DTR, summary.DTR, print.summary.DTR, plot.DTR


## Not run: 
est <- WRSEestimate(data=WRSEdata)
## End(Not run)

[Package DTR version 1.7 Index]