R/cTMed-delta-med-std-dot.R

In cTMed: Continuous-Time Mediation

.DeltaMedStd <- function(delta_t,
                         phi,
                         sigma,
                         vcov_theta,
                         from,
                         to,
                         med) {
  constructor <- function(delta_t,
                          from,
                          to,
                          med) {
    function(x) {
      .MedStdVec(
        v = x,
        delta_t = delta_t,
        from = from,
        to = to,
        med = med
      )
    }
  }
  func <- constructor(
    delta_t = delta_t,
    from = from,
    to = to,
    med = med
  )
  v <- c(
    .Vec(phi),
    .Vech(sigma)
  )
  jacobian <- numDeriv::jacobian(
    func = func,
    x = v
  )
  list(
    delta_t = delta_t,
    jacobian = jacobian,
    est = func(x = v),
    vcov = jacobian %*% vcov_theta %*% t(jacobian)
  )
}