The ctmle from jucheng1992

ctmle's People

Contributors

Stargazers

Watchers

ctmle's Issues

what is wrong with V=10?

Doesn't work for V=10

set.seed(10) 
x1= rbinom(N, 1, .3)
  x2= rbinom(N, 1, .5)
  x3= rbinom(N, 1, .7)
  x4 = rbinom(N, 1, .5)
  W = data.frame(x1,x2,x3,x4)
  head(W)
  A = rbinom(N, 1, plogis(.4+.3*x1-.4*x2+.5*x3-2*x4))
  Y = rbinom(N, 1, plogis(.4+.3*x1-.4*x2 + A))
  
  ctmleDiscrete(Y, A, W, Wg = W, Q = NULL, preOrder = FALSE, order = NULL,
                patience = FALSE, Qbounds = NULL, cvQinit = FALSE, Qform = NULL,
                SL.library = c("SL.glm"), alpha = 0.995, family = "binomial", gbound = 0.025,
                like_type = "RSS", fluctuation = "logistic", verbose = FALSE,
                detailed = FALSE, PEN = FALSE, V = 10, folds = NULL,
                stopFactor = 10^6)

yields
Error in if (nterms1[v] <= totalEst) { :
missing value where TRUE/FALSE needed

However, this works for V=5

x1= rbinom(N, 1, .3)
 x2= rbinom(N, 1, .5)
 x3= rbinom(N, 1, .7)
 x4 = rbinom(N, 1, .5)
 W = data.frame(x1,x2,x3,x4)
 head(W)
 A = rbinom(N, 1, plogis(.4+.3*x1-.4*x2+.5*x3-2*x4))
 Y = rbinom(N, 1, plogis(.4+.3*x1-.4*x2 + A))
 
 ctmleDiscrete(Y, A, W, Wg = W, Q = NULL, preOrder = FALSE, order = NULL,
               patience = FALSE, Qbounds = NULL, cvQinit = FALSE, Qform = NULL,
               SL.library = c("SL.glm"), alpha = 0.995, family = "binomial", gbound = 0.025,
               like_type = "RSS", fluctuation = "logistic", verbose = FALSE,
               detailed = FALSE, PEN = FALSE, V = 5, folds = NULL,
               stopFactor = 10^6)

Recommend Projects