hollorol/RBBGCMuso
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

constraints results table: constRes.csv

Browse Source
This commit is contained in:
Hollos Roland 2021-07-06 19:01:40 +02:00
parent 2ec8e1fccc
commit 53e1c0ec4f
1 changed files with 13 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

15
RBBGCMuso/R/multiSite.R
View File

@ -46,15 +46,19 @@ compose <- function(expr){
eval(parse(text=finalExpression),envir=penv) eval(parse(text=finalExpression),envir=penv)
} }
compoVect <- function(mod, constrTable){ compoVect <- function(mod, constrTable, fileToWrite = "const_results.data"){
with(as.data.frame(mod), { with(as.data.frame(mod), {
nexpr <- nrow(constrTable) nexpr <- nrow(constrTable)
filtered <- numeric(nexpr) filtered <- numeric(nexpr)
vali <- numeric(nexpr)
for(i in 1:nexpr){ for(i in 1:nexpr){
val <- compose(constrTable[i,1]) val <- compose(constrTable[i,1])
filtered[i] <- (val <= constrTable[i,3]) && filtered[i] <- (val <= constrTable[i,3]) &&
(val >= constrTable[i,2]) (val >= constrTable[i,2])
vali[i] <- val
} }
write(paste(vali,collapse=","), fileToWrite, append=TRUE)
filtered filtered
}) })
} }
@ -134,6 +138,7 @@ multiSiteCalib <- function(measurements,
} else { } else {
print("copy skipped") print("copy skipped")
file.remove(file.path(list.dirs("tmp",recursive=FALSE),"progress.txt")) file.remove(file.path(list.dirs("tmp",recursive=FALSE),"progress.txt"))
file.remove(file.path(list.dirs("tmp", recursive=FALSE), "const_results.data"))
} }
# ____ _ _ _ # ____ _ _ _
@ -207,6 +212,9 @@ multiSiteCalib <- function(measurements,
# | |__| (_) | | | | | | |_) | | | | | __/ # | |__| (_) | | | | | | |_) | | | | | __/
# \____\___/|_| |_| |_|_.__/|_|_| |_|\___| # \____\___/|_| |_| |_|_.__/|_|_| |_|\___|
resultFiles <- list.files(pattern="preservedCalib.*csv$",recursive=TRUE) resultFiles <- list.files(pattern="preservedCalib.*csv$",recursive=TRUE)
constRes <- file.path(list.dirs("tmp", recursive=FALSE), "const_results.data")
constRes <- lapply(constRes, function(f){read.csv(f, stringsAsFactors=FALSE, header=FALSE)})
write.csv(constRes, "constRes.csv")
res0 <- read.csv(grep("thread_1/",resultFiles, value=TRUE),stringsAsFactors=FALSE) res0 <- read.csv(grep("thread_1/",resultFiles, value=TRUE),stringsAsFactors=FALSE)
resultFilesSans0 <- grep("thread_1/", resultFiles, value=TRUE, invert=TRUE) resultFilesSans0 <- grep("thread_1/", resultFiles, value=TRUE, invert=TRUE)
# results <- do.call(rbind,lapply(resultFilesSans0, function(f){read.csv(f, stringsAsFactors=FALSE)})) # results <- do.call(rbind,lapply(resultFilesSans0, function(f){read.csv(f, stringsAsFactors=FALSE)}))
@ -223,7 +231,10 @@ multiSiteCalib <- function(measurements,
if(ncol(treeData) > 4){ if(ncol(treeData) > 4){
rp <- rpart(failType ~ .,data=treeData,control=treeControl) rp <- rpart(failType ~ .,data=treeData,control=treeControl)
svg("treeplot.svg") svg("treeplot.svg")
rpart.plot(rp) tryCatch(rpart.plot(rp), error = function(e){
print(e)
})
dev.off() dev.off()
} }
origModOut <- future::value(fut[[1]], stdout = FALSE, signal=FALSE)[["origModOut"]] origModOut <- future::value(fut[[1]], stdout = FALSE, signal=FALSE)[["origModOut"]]