tests/test_causalEffect.R

In pcalg: Methods for Graphical Models and Causal Inference

library(pcalg)
suppressWarnings(RNGversion("3.5.0"))
set.seed(34)
g <- randomDAG(5,0.8)
trMat <- matrix(0, 5,5)
trMat[1,5] <- 0.305
trMat[1,4] <- 0.863
trMat[1,2] <- 0.354
trMat[2,4] <- 0.392
trMat[2,5] <- 0.495
trMat[3,4] <- 0.278


## eff 1->5: 0.305 + 0.354*0.495
trEff <- 0.305 + 0.354*0.495
estEff <- causalEffect(g, 5,1)

if (!(round(trEff,3) == round(estEff, 3))) {
  stop("Test in wgtMatrix: Weights have wrong value!")
}