Learn R Programming
DynTxRegime (version 3.01)

.newOptimalSeq: Value Search Method for Optimal Treatment Regime.

Description

Value search methods to estimate the optimal treatment. Method not exported.

Usage

.newOptimalSeq(moPropen, moMain, moCont, regimes, fSet, ...)
# S4 method for modelObj,NULL,NULL,Regime,NULL
.newOptimalSeq(moPropen, moMain, moCont, regimes, fSet, data, response, txName, 
               suppress, argsList)
# S4 method for modelObj,NULL,NULL,Regime,function
.newOptimalSeq(moPropen, moMain, moCont, regimes, fSet, data, response, txName, 
               suppress, argsList)
# S4 method for ModelObj_SubsetList,NULL,NULL,Regime,function
.newOptimalSeq(moPropen, moMain, moCont, regimes, fSet, data, response, txName, 
               suppress, argsList)
# S4 method for modelObj,modelObj,modelObj,Regime,NULL
.newOptimalSeq(moPropen, moMain, moCont, regimes, fSet, iter, data, response, 
               txName, suppress, argsList)
# S4 method for modelObj,modelObj,NULL,Regime,NULL
.newOptimalSeq(moPropen, moMain, moCont, regimes, fSet, iter, data, response,  
               txName, suppress, argsList)
# S4 method for modelObj,NULL,modelObj,Regime,NULL
.newOptimalSeq(moPropen, moMain, moCont, regimes, fSet, iter, data, response,  
               txName, suppress, argsList)
# S4 method for modelObj,modelObj,modelObj,Regime,function
.newOptimalSeq(moPropen, moMain, moCont, regimes, fSet, iter, data, response,  
               txName, suppress, argsList)
# S4 method for modelObj,modelObj,NULL,Regime,function
.newOptimalSeq(moPropen, moMain, moCont, regimes, fSet, iter, data, response,  
               txName, suppress, argsList)
# S4 method for modelObj,NULL,modelObj,Regime,function
.newOptimalSeq(moPropen, moMain, moCont, regimes, fSet, iter, data, response,  
               txName, suppress, argsList)
# S4 method for ModelObj_SubsetList,modelObj,modelObj,Regime,function
.newOptimalSeq(moPropen, moMain, moCont, regimes, fSet, iter, data, response,  
               txName, suppress, argsList)
# S4 method for ModelObj_SubsetList,modelObj,NULL,Regime,function
.newOptimalSeq(moPropen, moMain, moCont, regimes, fSet, iter, data, response,  
               txName, suppress, argsList)
# S4 method for ModelObj_SubsetList,NULL,modelObj,Regime,function
.newOptimalSeq(moPropen, moMain, moCont, regimes, fSet, iter, data, response,  
               txName, suppress, argsList)
# S4 method for ModelObj_SubsetList,ModelObj_SubsetList,ModelObj_SubsetList,Regime,function
.newOptimalSeq(moPropen, moMain, moCont, regimes, fSet, iter, data, response,  
               txName, suppress, argsList)
# S4 method for ModelObj_SubsetList,ModelObj_SubsetList,NULL,Regime,function
.newOptimalSeq(moPropen, moMain, moCont, regimes, fSet, iter, data, response,  
               txName, suppress, argsList)
# S4 method for ModelObj_SubsetList,NULL,ModelObj_SubsetList,Regime,function
.newOptimalSeq(moPropen, moMain, moCont, regimes, fSet, iter, data, response,  
               txName, suppress, argsList)
# S4 method for modelObj,ModelObj_SubsetList,ModelObj_SubsetList,Regime,function
.newOptimalSeq(moPropen, moMain, moCont, regimes, fSet, iter, data, response,  
               txName, suppress, argsList)
# S4 method for modelObj,ModelObj_SubsetList,NULL,Regime,function
.newOptimalSeq(moPropen, moMain, moCont, regimes, fSet, iter, data, response,  
               txName, suppress, argsList)
# S4 method for modelObj,NULL,ModelObj_SubsetList,Regime,function
.newOptimalSeq(moPropen, moMain, moCont, regimes, fSet, iter, data, response,  
               txName, suppress, argsList)
# S4 method for ModelObj_DecisionPointList,NULL,NULL,Regime_DecisionPointList,list
.newOptimalSeq(moPropen, moMain, moCont, regimes, fSet, data, response, txName, 
               suppress, argsList)
# S4 method for ModelObj_DecisionPointList,NULL,NULL,Regime_DecisionPointList,NULL
.newOptimalSeq(moPropen, moMain, moCont, regimes, fSet, data, response, txName, 
               suppress, argsList)
# S4 method for ModelObj_SubsetList_DecisionPointList,NULL,NULL,Regime_DecisionPointList,list
.newOptimalSeq(moPropen, moMain, moCont, regimes, fSet, data, response, 
               txName, suppress, argsList)
# S4 method for ModelObj_DecisionPointList,ModelObj_DecisionPointList,ModelObj_DecisionPointList,Regime_DecisionPointList,NULL
.newOptimalSeq(moPropen, moMain, moCont, regimes, fSet, iter, data, response,  
               txName, suppress, argsList)
# S4 method for ModelObj_DecisionPointList,ModelObj_DecisionPointList,NULL,Regime_DecisionPointList,NULL
.newOptimalSeq(moPropen, moMain, moCont, regimes, fSet, iter, data, response,  
               txName, suppress, argsList)
# S4 method for ModelObj_DecisionPointList,NULL,ModelObj_DecisionPointList,Regime_DecisionPointList,NULL
.newOptimalSeq(moPropen, moMain, moCont, regimes, fSet, iter, data, response,  
               txName, suppress, argsList)
# S4 method for ModelObj_DecisionPointList,ModelObj_DecisionPointList,ModelObj_DecisionPointList,Regime_DecisionPointList,list
.newOptimalSeq(moPropen, moMain, moCont, regimes, fSet, iter, data, response,  
               txName, suppress, argsList)
# S4 method for ModelObj_DecisionPointList,ModelObj_DecisionPointList,NULL,Regime_DecisionPointList,list
.newOptimalSeq(moPropen, moMain, moCont, regimes, fSet, iter, data, response,  
               txName, suppress, argsList)
# S4 method for ModelObj_DecisionPointList,NULL,ModelObj_DecisionPointList,Regime_DecisionPointList,list
.newOptimalSeq(moPropen, moMain, moCont, regimes, fSet, iter, data, response,  
               txName, suppress, argsList)
# S4 method for ModelObj_SubsetList_DecisionPointList,ModelObj_DecisionPointList,ModelObj_DecisionPointList,Regime_DecisionPointList,list
.newOptimalSeq(moPropen, moMain, moCont, regimes, fSet, iter, data, response,  
               txName, suppress, argsList)
# S4 method for ModelObj_SubsetList_DecisionPointList,ModelObj_DecisionPointList,NULL,Regime_DecisionPointList,list
.newOptimalSeq(moPropen, moMain, moCont, regimes, fSet, iter, data, response,  
               txName, suppress, argsList)
# S4 method for ModelObj_SubsetList_DecisionPointList,NULL,ModelObj_DecisionPointList,Regime_DecisionPointList,list
.newOptimalSeq(moPropen, moMain, moCont, regimes, fSet, iter, data, response,  
               txName, suppress, argsList)
# S4 method for ModelObj_SubsetList_DecisionPointList,ModelObj_SubsetList_DecisionPointList,ModelObj_SubsetList_DecisionPointList,Regime_DecisionPointList,list
.newOptimalSeq(moPropen, moMain, moCont, regimes, fSet, iter, data, response,  
               txName, suppress, argsList)
# S4 method for ModelObj_SubsetList_DecisionPointList,ModelObj_SubsetList_DecisionPointList,NULL,Regime_DecisionPointList,list
.newOptimalSeq(moPropen, moMain, moCont, regimes, fSet, iter, data, response,  
               txName, suppress, argsList)
# S4 method for ModelObj_SubsetList_DecisionPointList,NULL,ModelObj_SubsetList_DecisionPointList,Regime_DecisionPointList,list
.newOptimalSeq(moPropen, moMain, moCont, regimes, fSet, iter, data, response,  
               txName, suppress, argsList)
# S4 method for ModelObj_DecisionPointList,ModelObj_SubsetList_DecisionPointList,ModelObj_SubsetList_DecisionPointList,Regime_DecisionPointList,list
.newOptimalSeq(moPropen, moMain, moCont, regimes, fSet, iter, data, response,  
               txName, suppress, argsList)
# S4 method for ModelObj_DecisionPointList,ModelObj_SubsetList_DecisionPointList,NULL,Regime_DecisionPointList,list
.newOptimalSeq(moPropen, moMain, moCont, regimes, fSet, iter, data, response,  
               txName, suppress, argsList)
# S4 method for ModelObj_DecisionPointList,NULL,ModelObj_SubsetList_DecisionPointList,Regime_DecisionPointList,list
.newOptimalSeq(moPropen, moMain, moCont, regimes, fSet, iter, data, response,  
               txName, suppress, argsList)

Arguments

moPropen
"modelObj" for propensity score regression.
moMain
"modelObj" for main effects of outcome regression.
moCont
"modelObj" for contrasts of outcome regression.
regimes
"function" defining class of regimes
fSet
"function" defining feasible treatment rules
...
Used to pass arguments that are required but the class of which does not determine the method that is selected.
data
"data.frame" of covariates and treatment history.
response
Outcome of interest.
txName
"character" name of column header in data containing treatment variable.
iter
Maximum number of iterations for iterative algorithm.
suppress
T/F indicating if screen prints are generated.
argsList
"list" of additional arguments to be passed to rgenoud.

Examples

Run this code

data(bmiData)

y <- -(bmiData$month12BMI - bmiData$baselineBMI) / bmiData$baselineBMI * 100

bmiData$A2 <- as.factor(bmiData$A2)

moPropen <- buildModelObj(model = ~1, 
                          solver.method = 'glm',
                          solver.args = list("family" = "binomial"),
                          predict.args = list("type" = "response"))

moMain <- buildModelObj(model = ~parentBMI+baselineBMI+month4BMI, 
                        solver.method = 'lm')
moCont <- buildModelObj(model = ~parentBMI+baselineBMI+month4BMI, 
                        solver.method = 'lm')

regime <- function(a,data){
            vec <- character(nrow(data))
            vec[data$month4BMI < a] <- "MR"
            vec[data$month4BMI >= a] <- "CD"
            return( vec )
          }

regime <- DynTxRegime:::.newRegime(object = regime)

regList <- list(Domains = matrix(c(20,40),ncol=2), 
                pop.size = 500, 
                starting.values = 30)

## Not run: ------------------------------------
# obj <- DynTxRegime:::.newOptimalSeq(moPropen = moPropen,
#                                     moMain = NULL, 
#                                     moCont = NULL, 
#                                     regime = regime,
#                                     fSet = NULL,
#                                     data = bmiData,
#                                     response = y,
#                                     txName = 'A2', 
#                                     suppress = TRUE,
#                                     argsList = regList)
# 
# is(obj)
# 
# coef(obj)
# DTRstep(obj)
# estimator(obj)
# fitObject(obj)
# optTx(obj)
# optTx(obj,bmiData)
# outcome(obj)
# plot(obj)
# propen(obj)
# show(obj)
# summary(obj)
## ---------------------------------------------

Run the code above in your browser using DataLab