Estoy usandobnlearn ypcalg R paquetes para obtener el mapa de causalidad de los conjuntos de datos. Hay un algoritmo independiente del orden que dice ser independiente de las variables que se proporcionan como entrada. Cuando cambio el orden de las variables, las direcciones de las puntas de flecha están cambiando. A continuación se muestra el código que estoy 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)

El código anterior me da el siguiente resultado:

Ahora ejecuto el mismo código con diferente orden de datos.

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 puede ver allí, las flechas han cambiado la dirección de v6 y v7. ¿Me estoy perdiendo de algo? Nota: Soy consciente de que el esqueleto no ha cambiado (gráfico sin flechas).