library(huge)
library(xtable)
devtools::install_github("mattdelhey/bmar")
library(bmar)
## Parameters
f.out.graphs <- "simulation-graphs.pdf"
reps <- 2
## Simulations
#sim.plot(sim.objects, f.out.graphs)
## sim.result[method, error, model, replicate]
sim.results <- simplify2array(lapply(1:reps, function(i) {
sim.random <- huge.generator(n = 200, d = 300, v = 0.3, u = 0.1, graph = "random", prob = 0.006)
sim.cluster <- huge.generator(n = 200, d = 300, v = 0.3, u = 0.1, graph = "cluster")
sim.band <- huge.generator(n = 200, d = 300, g = 2, graph = "band")
sim.sclfree <- huge.generator(n = 200, d = 300, graph = "scale-free")
sim.objects <- list(sim.random, sim.cluster, sim.band, sim.sclfree)
sim.summarize(sim.objects, nlambda = 1000, K = 5, stars.thresh = 0.05,
stars.subsample.ratio = NULL, rep.num = 10,
l = 10, niter = 10000, burnin = 2000)
}))
## Mean, sd
means <- apply(sim.results, c(1, 2, 3), mean, na.rm = TRUE)
apply(means[, c(-3, -4), ], 3, xtable)
sds <- apply(sim.results, c(1, 2, 3), sd, na.rm = TRUE)
apply(sds[, c(-3, -4), ], 3, xtable)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.