LDTestimate: Function for calculating LTD estimates In DTR: Estimation and Comparison of Dynamic Treatment Regimes

Description

The function computes the survival estimates and estimated standard errors for dynamic treatment regimes (DTRs) at the observed event times as proposed in Lunceford, Davidian and Tsiatis (2002) Equation (3) and Equation (10).

Usage

 `1` ```LDTestimate(data, L = .Machine\$double.xmax) ```

Arguments

 `data` a data frame (X, R, Z, U, delta) representing the data from a two-stage randomization design 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 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 event time if delta=1, and U is censoring time if delta=0 delta: censoring indicator, delta=1 for event, and delta=0 for censored `L` restricted survival time. Default is .Machine\$double.xmax, which is the largest double value of R. Set L to a numeric number smaller than the maximum follow-up time if restricted follow-up time up to L is considered.

Details

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 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.

Value

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

References

Lunceford JK, Davidian M, Tsiatis AA: Estimation of survival distributions of treatment policies in two-stage randomization designs in clinical trials. Biometrics 58:48-57, 2002

`simLDTdata`, `DTR.object`, `print.DTR`, `summary.DTR`, `print.summary.DTR`, `plot.DTR`
 ```1 2 3 4 5``` ```## Not run: data("LDTdata") est <- LDTestimate(data=LDTdata) est ## End(Not run) ```