A Generalized Logrank-type Test for Comparison of Treatment Regimes in Sequential Multiple Assignment Randomized Trials

Anastasios A. Tsiatis,Marie Davidian

2024-03-25

Abstract:The sequential multiple assignment randomized trial (SMART) is the ideal study design for the evaluation of multistage treatment regimes, which comprise sequential decision rules that recommend treatments for a patient at each of a series of decision points based on their evolving characteristics. A common goal is to compare the set of so-called embedded regimes represented in the design on the basis of a primary outcome of interest. In the study of chronic diseases and disorders, this outcome is often a time to an event, and a goal is to compare the distributions of the time-to-event outcome associated with each regime in the set. We present a general statistical framework in which we develop a logrank-type test for comparison of the survival distributions associated with regimes within a specified set based on the data from a SMART with an arbitrary number of stages that allows incorporation of covariate information to enhance efficiency and can also be used with data from an observational study. The framework provides clarification of the assumptions required to yield a principled test procedure, and the proposed test subsumes or offers an improved alternative to existing methods. We demonstrate performance of the methods in a suite of simulation studies. The methods are applied to a SMART in patients with acute promyelocytic leukemia.

Methodology

What problem does this paper attempt to address?

The problem that this paper attempts to solve is to compare the effects of different treatment regimes in Sequential, Multiple Assignment, Randomized Trials (SMART). Specifically, the author focuses on comparing the survival distributions of these treatment regimes based on time - to - event (such as survival time) outcomes. The paper proposes a generalized Logrank - type test method for comparing the effects of multiple treatment regimes in SMART design, especially when the data may contain censoring information. In addition, this method can also incorporate covariate information to improve the test efficiency and is applicable to observational study data.

### Core problems of the paper
1. **Comparing the effects of treatment regimes**: How to effectively compare the time - to - event outcomes of different treatment regimes in SMART, especially in the presence of censored data.
2. **Improving test efficiency**: How to improve the statistical power of the test by incorporating covariate information.
3. **Universality of the method**: The proposed test method is applicable not only to SMART but also to observational study data.

### Solutions
- **Generalized Logrank - type test**: The author has developed a new statistical framework and proposed a generalized Logrank - type test method for comparing the survival distributions of different treatment regimes. This method can handle any number of decision points and can improve the test efficiency in the presence of covariates.
- **Inverse probability weighting**: Censored data are processed by inverse probability weighting (IPW) technique to ensure the validity of the test.
- **Covariate adjustment**: Baseline and intermediate covariate information are combined to improve the statistical power of the test.

### Application examples
- **Acute promyelocytic leukemia (APL) study**: The proposed method in the paper is applied to analyze a SMART study (Study C9710) on APL patients and compare the effects of four embedded treatment regimes.

### Formula explanations
- **Survival function and hazard rate**:
- \( T^*(d) \) represents the potential survival time if an individual follows treatment regime \( d \).
- \( \lambda(u, d) \) represents the hazard rate of an individual following treatment regime \( d \) at time \( u \).
- \( \Lambda(u, d)=\int_{0}^{u} \lambda(s, d) \, ds \) represents the cumulative hazard function.

- **Counting process and risk set**:
- \( N^*(u, d)=I\{T^*(d) \leq u\} \) represents the event - occurrence counting process up to time \( u \).
- \( Y^*(u, d)=I\{T^*(d) \geq u\} \) represents the individuals still in the risk set up to time \( u \).

- **Inverse probability weighting**:
- \( C(u, d) \) and \( \pi(u, d) \) respectively represent the indicator function and inverse probability weight of an individual receiving treatment consistent with treatment regime \( d \).
- \( \Omega(u, d)=\frac{C(u, d) I(U \geq u)}{\pi(u, d) K_c(u | H(u))} \) represents the weight after inverse probability weighting.

- **Test statistic**:
- \( T_j=\sum_{i = 1}^{n} \int_{0}^{\infty} \Omega_i(u, d_j)\{dN_i(u)-\hat{\Lambda}_0(u) dY_i(u)\} \) represents the score vector of the \( j \) - th treatment regime.
- \( Z = n^{-1}T^T\hat{\Sigma}^{-1}T \) represents the final test statistic, where \( \hat{\Sigma} \) is the estimated value of the asymptotic covariance matrix of \( T \).

Through these methods, the paper provides a systematic and comprehensive framework for comparing the effects of different treatment regimes in SMART, especially in the presence of censored data.