library(mgcv) library(dHSIC) source("utils.R") XX <- read.csv('./RestrictedSCMsData1.txt', sep = "\t") XX <- XX[1:1000,] Y <- XX[,2] X <- XX[,1] plot(X, Y, pch = 19, cex = .8) load(file = "./RestrictedSCMsData2.RData") # loads dat load(file = "./RestrictedSCMsallDagsWith4Nodes.RData") # loads all DAGs as allDags load(file = "./RestrictedSCMsallFullDagsWith4Nodes.RData") # loads all fully connected DAGs as allFullDags