Genetic.optim.IWKMEIV: The optimal treatment regime based on the (S)IWKMEIV...

In otrKM: Optimal Treatment Regimes in Survival Contexts with Kaplan-Meier-Like Estimators

The optimal treatment regime based on the (S)IWKMEIV estimator.

Description

Given a predetermined t0, estimate the optimal treatment regime by maximizing t0-year survival probability based on the (S)IWKMEIV estimator.

Usage

Genetic.optim.IWKMEIV(datalist, ps, t0, smooth = TRUE)

Arguments

datalist	A list used to calculate the (S)IWKMEIV estimator including instrument named z, treatment named a, observed time named obs.t, censoring indicator (0, censored) named delta, and baseline covariates used to assign treatment named l. Notice that all the data in the datalist should be ordered by observed time.
ps	A list including the probability of receiving instrument given baseline covariates named fzl and the difference between fal1 and fal0 named deltal, where fal0 denotes the probability of receiving treatment given baseline covariates and instrument equaling 0, and fal1 denotes the probability of receiving treatment given baseline covariates and instrument equaling 1. Fps.IWKMEIV can produce ps by positing logistic model.
t0	A predetermined time.
smooth	A logic variable indicating wether a smoothed version should be used.

Details

More details can be found in references.

Value

A numeric vector in which the last number is the estimated optimal t0-year survival probability and others are the estimated parameter of the optimal treatment regime.

References

Xia, J., Zhan, Z., Zhang, J. (2022) Estimating optimal treatment regime in survival contexts using an instrumental variable. Under Review.

Examples

# load data
data(simulation)
simulation=simulation[order(simulation$Survival),]

# convert the data into a datalist
datalist=list(z=simulation$Instrument,a=simulation$Treatment,
              obs.t=simulation$Survival,delta=simulation$Status,
              l=cbind(simulation$Covariate1,simulation$Covariate2))

# predetermined t0
t0=1

# calculate ps and prep
ps=Fps.IWKMEIV(datalist, t0)

Genetic.optim.IWKMEIV(datalist, ps, t0, smooth=TRUE)

otrKM documentation built on July 9, 2023, 5:35 p.m.