tests/testthat/test_basic.R

In podTockom/tstmle01: Targeted Maximum Likelihood Estimation of Causal Effects of Interventions on a Single Binary Time Series

library(Hmisc)
library(doParallel)  
library(parallel)
library(foreach)
library(tstmle01)

set.seed(2)
data("ts_samp_data_ex")
data=ts_samp_data_ex

# Get initial estimates of the process:
fit <- initEst(data, t=5, freqW = 2, freqA = 2, freqY = 2)

# Get estimates:
est1 <- mainTMLE(fit, t = 5, Anode = 3, intervention = 1, B = 100, N = 100, MC = 100, maxIter = 20)
est2 <- mainTMLE(fit, t = 5, Anode = 3, B = 100, N = 100, MC = 50, maxIter = 20)