mcX2TestB50k: Wrapper around the bnlearn mc-x2 test (B=50k)

In ddgraph: Distinguish direct and indirect interactions with Graphical Modelling

Description

Version of mcX2Test() with 50000 Monte Carlo replicates.

Usage

mcX2TestB50k(x, y, S, suffStat)

Arguments

x	the index of the first variable
y	the index of the second variable
S	the conditioning set
suffStat	the sufficient statistics to do the test, in this case a list of one element: dm where the values matrix is stored

Value

p value of the test

Examples

suffStat <- list(dm = cbind("a"=c(0,1,0,0,1,0), "b"=c(1,0,0,0,1,0), "c"=c(0,0,0,1,1,1)))
# test if a is independent of b
mcX2TestB50k(1, 2, NULL, suffStat)
# test if a is independent of b conditioned on c
mcX2TestB50k(1, 2, 3, suffStat)

ddgraph documentation built on Nov. 17, 2017, 10:50 a.m.