变量顺序更改时因果关系图更改

Question

我正在使用bnlearn和pcalg R包从数据集中获取因果关系图。有一种与顺序无关的算法，该算法声称与变量无关，作为输入给出。当我更改变量的顺序时，箭头的方向也在改变。下面是我正在使用的代码：

data("gmG")
set.seed(101)
cols = sample(ncol(gmG8$x))
suffStat <- list(C = cor(gmG8$x[,cols]), n = nrow(gmG8$x))
pc.gmG <- pc(suffStat, indepTest = gaussCItest,
               labels = colnames(gmG8$x)[cols], alpha = 0.01)
plot(pc.gmG)

以上代码为我提供了以下输出：

现在，我以不同的数据顺序运行相同的代码。

data("gmG")
set.seed(102)
cols = sample(ncol(gmG8$x))
suffStat <- list(C = cor(gmG8$x[,cols]), n = nrow(gmG8$x))
pc.gmG <- pc(suffStat, indepTest = gaussCItest,
               labels = colnames(gmG8$x)[cols], alpha = 0.01)
plot(pc.gmG)

正如人们所看到的，箭头已经改变了v6和v7的方向。 我在这里想念什么吗？ 注意：我知道骨架没有改变（图形没有箭头）。

Answer 1

我找到了答案。 通过使用conservative将maj.rule或TRUE设置为solve.confl = TRUE，可以使因果关系图完全与顺序无关。

pc.gmG <- pc(suffStat, indepTest = gaussCItest,skel.method = 'stable',
                  conservative = TRUE,solve.confl = TRUE,
               labels = colnames(gmG8$x)[cols], alpha = 0.01)