## Not run:
# library(pcalg)
# set.seed(4)
# n <- 300
# data <- NULL
# x1 <- 2*(runif(n)-0.5)
# x2 <- x1 + runif(n)-0.5
# x3 <- x1^2 + 0.6*runif(n)
# x4 <- rnorm(n)
# x5 <- x3 + x4^2 + 2*runif(n)
# x6 <- 10*(runif(n)-0.5)
# x7 <- x6^2 + 5*runif(n)
# x8 <- 2*x7^2 + 1.5*rnorm(n)
# x9 <- x7 + 4*runif(n)
# data <- cbind(x1,x2,x3,x4,x5,x6,x7,x8,x9)
# true <- matrix(0,9,9)
# true[c(1),c(2,3)]<-true[c(3,4),5]<-true[c(6),c(7)]<-true[c(7),c(8)]<-true[7,9]<-1
#
# pc <- pc(suffStat = list(C = cor(data), n = 9),
# indepTest = gaussCItest,
# alpha = 0.9,
# labels = colnames(data),
# u2pd = "relaxed",
# skel.method = "stable",
# verbose = TRUE)
# kpc1 <- kpc(suffStat = list(data=data, ic.method="dcc.perm"),
# indepTest = kernelCItest,
# alpha = 0.1,
# labels = colnames(data),
# u2pd = "relaxed",
# skel.method = "stable",
# verbose = TRUE)
# kpc2 <- kpc(suffStat = list(data=data, ic.method="hsic.gamma"),
# indepTest = kernelCItest,
# alpha = 0.1,
# labels = colnames(data),
# u2pd = "relaxed",
# skel.method = "stable",
# verbose = TRUE)
# kpc3 <- kpc(suffStat = list(data=data, ic.method="hsic.perm"),
# indepTest = kernelCItest,
# alpha = 0.1,
# labels = colnames(data),
# u2pd = "relaxed",
# skel.method = "stable",
# verbose = TRUE)
# kpc4 <- kpc(suffStat = list(data=data, ic.method="hsic.clust"),
# indepTest = kernelCItest,
# alpha = 0.1,
# labels = colnames(data),
# u2pd = "relaxed",
# skel.method = "stable",
# verbose = TRUE)
#
# if (require(Rgraphviz)) {
# par(mfrow=c(2,3))
# plot(pc,main="pc")
# plot(kpc1,main="dpc.perm")
# plot(kpc2,main="kpc.gamma")
# plot(kpc3,main="kpc.perm")
# plot(kpc4,main="kpc.clust")
# plot(as(true,"graphNEL"),main="True DAG")
# }
# ## End(Not run)
Run the code above in your browser using DataLab