# NOT RUN {
## Data
data(ecoli)
data <- ecoli[,-c(1,9)] ## remove sequece.name and class
X <- TrainingandTestData(data, percentage_test = 0.95)
Xtraining <- X$Training
Xtest <- X$Test
## DAG
dag <- LearningHC(data)
plot(dag)
## Means and desviations
colnames(data)
m <- sapply(data, mean)
m <- m[-which(is.na(m))]
names(m)
d <- sapply(data, sd)
d <- d[-which(is.na(m))]
names(d)
## Prior Dataset
n <- 5600
priorData <- generateNormalPriorData(dag, data = Xtraining, size = n, means = m)
summary(priorData)
ncol(priorData)
nrow(priorData)
class(priorData)
# }
Run the code above in your browser using DataLab