cluster/rv144/fixResults.R

cpus <- 1
args <- commandArgs(TRUE)
eval(parse(text=args[[1]]))
setting <- as.numeric(setting)

library(flowReMix)

files <- as.list(dir(path = 'results', pattern="rv144_*"))
dirfiles <- lapply(files, function(x) paste0('results/', x))
file <- dirfiles[[setting]]
load(file = file)

stability <- stabilityGraph(fit, type = "ising", cv = FALSE, reps = 200,
                            cpus = 1, AND = TRUE, gamma = 0.25)
stabfile <- paste("results2/stab_", files[[setting]], sep ="")
save(stability, file = stabfile)

fit$assignmentList <- NULL
fit$randomEffectSamp <- NULL
savefile <- paste("results2/", files[[setting]], sep ="")
save(fit, file = savefile)
RGLab/flowReMix documentation built on May 8, 2019, 5:55 a.m.