From a051b465aa5a2a684dbb7dba230c7fc7d798d1fc Mon Sep 17 00:00:00 2001 From: Hollos Roland Date: Wed, 30 Jun 2021 08:49:56 +0200 Subject: [PATCH] without tree --- RBBGCMuso/R/multiSite.R | 20 ++++++++++---------- 1 file changed, 10 insertions(+), 10 deletions(-) diff --git a/RBBGCMuso/R/multiSite.R b/RBBGCMuso/R/multiSite.R index 2bbe2bd..1bfd904 100644 --- a/RBBGCMuso/R/multiSite.R +++ b/RBBGCMuso/R/multiSite.R @@ -216,17 +216,17 @@ multiSiteCalib <- function(measurements, results <- (rbind(res0,resultsSans0)) write.csv(results,"result.csv") calibrationPar <- future::value(fut[[1]], stdout = FALSE, signal=FALSE)[["calibrationPar"]] - notForTree <- c(seq(from = (length(calibrationPar)+1), length.out=3)) + # notForTree <- c(seq(from = (length(calibrationPar)+1), length.out=3)) # notForTree <- numeric(0) - notForTree <- c(notForTree,which(sapply(seq_along(calibrationPar),function(i){sd(results[,i])==0}))) - treeData <- results[,-notForTree] - treeData["failType"] <- as.factor(results$failType) - if(ncol(treeData) > 4){ - rp <- rpart(failType ~ .,data=treeData,control=treeControl) - svg("treeplot.svg") - rpart.plot(rp) - dev.off() - } + # notForTree <- c(notForTree,which(sapply(seq_along(calibrationPar),function(i){sd(results[,i])==0}))) + # treeData <- results[,-notForTree] + # treeData["failType"] <- as.factor(results$failType) + # if(ncol(treeData) > 4){ + # rp <- rpart(failType ~ .,data=treeData,control=treeControl) + # svg("treeplot.svg") + # rpart.plot(rp) + # dev.off() + # } origModOut <- future::value(fut[[1]], stdout = FALSE, signal=FALSE)[["origModOut"]] # Just single objective version TODO:Multiobjective results <- results[results[,"Const"] == 1,]