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

without tree

Browse Source
This commit is contained in:
Hollos Roland 2021-06-30 08:49:56 +02:00
parent c257e965bd
commit a051b465aa
1 changed files with 10 additions and 10 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

20
RBBGCMuso/R/multiSite.R
View File

@ -216,17 +216,17 @@ multiSiteCalib <- function(measurements,
results <- (rbind(res0,resultsSans0)) results <- (rbind(res0,resultsSans0))
write.csv(results,"result.csv") write.csv(results,"result.csv")
calibrationPar <- future::value(fut[[1]], stdout = FALSE, signal=FALSE)[["calibrationPar"]] 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 <- numeric(0)
notForTree <- c(notForTree,which(sapply(seq_along(calibrationPar),function(i){sd(results[,i])==0}))) # notForTree <- c(notForTree,which(sapply(seq_along(calibrationPar),function(i){sd(results[,i])==0})))
treeData <- results[,-notForTree] # treeData <- results[,-notForTree]
treeData["failType"] <- as.factor(results$failType) # treeData["failType"] <- as.factor(results$failType)
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) # rpart.plot(rp)
dev.off() # dev.off()
} # }
origModOut <- future::value(fut[[1]], stdout = FALSE, signal=FALSE)[["origModOut"]] origModOut <- future::value(fut[[1]], stdout = FALSE, signal=FALSE)[["origModOut"]]
# Just single objective version TODO:Multiobjective # Just single objective version TODO:Multiobjective
results <- results[results[,"Const"] == 1,] results <- results[results[,"Const"] == 1,]