Estou usandobnlearn epcalgacotes @ R para obter o mapa de causalidade dos conjuntos de dados. Existe um algoritmo independente da ordem que afirma ser independente das variáveis é dado como entrada. Quando mudo a ordem das variáveis, as direções das pontas das setas estão mudando. Abaixo está o código que estou usando:

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)

código @Above fornece a seguinte saída:

gora, eu corro o mesmo código com uma ordem diferente de dado

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)

Como se pode ver lá, as setas mudaram a direção das v6 e v7. Estou faltando alguma coisa aqui? Nota: Estou ciente de que o esqueleto não foi alterado (gráfico sem setas