Two nuisance fits and two fluctuations, one for the outcome and one for the mediator. Targeting both buys the double robustness: a mis-specified mediator model can be rescued by a correct outcome model and the other way round, which a product-of-coefficients estimator cannot do.
Examples
TmlMd(Y = c(1, 2, 3, 4, 5, 6, 7, 8), X = c(1, 2, 3, 4, 5, 6, 7, 8), M = 5L,
Cc = c(1, 2, 3, 4, 5, 6, 7, 8))
#> $estimate
#> [1] -7.112651e+15
#>
#> $nie
#> [1] NA
#>
#> $total
#> [1] NA
#>
#> $se
#> [1] 219.6523
#>
#> $eps
#> [1] -0.01354023
#>
#> $n
#> [1] 8
#>
#> $method
#> [1] "TMLE for natural direct and indirect effects"
#>
#> attr(,"class")
#> [1] "morie_rich_result" "list"
