R/delta_method_function.R

In gpapadog/Interference: Causal inference with interference for realistic regimes

#' @export
delta_method <- function(x, vec = c(1, - 1))  {
  
  if (length(vec) != dim(x)[1] | length(vec) != dim(x)[2]) {
    stop('Wrong dimensions.')
  }
  
  vec <- matrix(vec, ncol = 1)
  return(t(vec) %*% x %*% vec)
}