1

这里有关于如何为二进制数据运行 pc 算法的文档 ( https://www.rdocumentation.org/packages/pcalg/versions/2.7-3/topics/pc )。

library(pcalg)
##################################################
## Using discrete data
##################################################
## Load data
data(gmD)
V <- colnames(gmD$x)
## define sufficient statistics
suffStat <- list(dm = gmD$x, nlev = c(3,2,3,4,2), adaptDF = FALSE)
## estimate CPDAG
pc.D <- pc(suffStat,
           ## independence test: G^2 statistic
           indepTest = disCItest, alpha = 0.01, labels = V, verbose = TRUE)
if (require(Rgraphviz)) {
  ## show estimated CPDAG
  par(mfrow = c(1,2))
  plot(pc.D, main = "Estimated CPDAG")
  plot(gmD$g, main = "True DAG")
}

我想知道是否无论如何我可以使用该pc.D对象在ggdag（https://cran.r-project.org/web/packages/ggdag/vignettes/intro-to-ggdag.html）中可视化它。

有办法吗？

4

0 回答 0