summary.bcdagCE {BCDAG}

bcdagCE object summary

Description

This method produces summaries of the input and output of the function get_causaleffect().

Usage

## S3 method for class 'bcdagCE'
summary(object, ...)

Arguments

object

a bcdagCE object for which a summary is desired

...

additional arguments affecting the summary produced

Value

A printed message listing the inputs given to learn_DAG() and get_causaleffect() and summary statistics of the posterior distribution.

Examples

q = 8
w = 0.2
set.seed(123)
DAG = rDAG(q = q, w = w)
outDL = rDAGWishart(n = 1, DAG = DAG, a = q, U = diag(1, q))
L = outDL$L; D = outDL$D
Sigma = solve(t(L))%*%D%*%solve(L)
n = 200
# Generate observations from a Gaussian DAG-model
X = mvtnorm::rmvnorm(n = n, sigma = Sigma)
# Run the MCMC (set S = 5000 and burn = 1000 for better results)
out_mcmc = learn_DAG(S = 500, burn = 100, a = q, U = diag(1,q)/n, data = X, w = w,
                     fast = TRUE, save.memory = FALSE, verbose = FALSE)
out_ce <- get_causaleffect(out_mcmc, targets = c(4,6), response = 1)
# summary(out_ce)