library(flowReMix)
ncpus <- 15
filenames <- as.list(dir(path = 'results', pattern="rv144_2_*"))
filenames <- lapply(filenames, function(x) paste0('results/', x))
post <- list()
random <- list()
assign <-list()
for(i in 1:length(filenames)) {
load(file = filenames[[i]])
post[[i]] <- fit$posteriors[, -1]
random[[i]] <- fit$randomEffectSamp
assign[[i]] <- fit$assignmentList
}
post <- Reduce("+", post) / length(filenames)
assign <- do.call("c", assign)
random <- do.call("c", random)
fit$posteriors[, -1] <- post
fit$randomEffectSamp <- random
fit$assignmentList <- assign
stability <- stabilityGraph(fit, type = "ising", cv = FALSE, reps = 500,
cpus = ncpus, AND = TRUE, gamma = 0.25)
save(stability, file = "results/rv144AggStab2.Robj")
stability <- stabilityGraph(fit, type = "randomEffects", cv = FALSE, reps = 500,
cpus = ncpus, AND = TRUE, gamma = 0.25)
save(stability, file = "results/rv144AggStab2.Robj")
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.