ipc: Estimate an IDAG

In dhessing/ipcalg: Causal discovery for IDAGs

Description

Estimate an IDAG for the causal effect of x on y from a fitted PC function.

Usage

ipc(pc.fit, data, x, y, alpha = 0.01)

Arguments

pc.fit	fitted PC object
data	dataset
x	intervention
y	outcome
alpha	significance level for regressions

Value

fitted IPC object

Examples

library(tibble)
library(pcalg)
n <- 1000
data <- tibble::tibble(
  A = rnorm(n),
  B = rnorm(n),
  Y = A + B + A * B + rnorm(n)
)

pc.fit <- pc(
  suffStat = list(C = cor(data), n = nrow(data)),
  indepTest = gaussCItest,
  alpha = 0.01,
  labels = colnames(data)
)

ipc.fit <- ipc(pc.fit, data, x = 'A', y = 'Y', alpha = 0.01)

dhessing/ipcalg documentation built on Dec. 19, 2021, 11:07 p.m.