fixing sensi bug

RBBGCMuso/R/calibration.R

@@ -109,7 +109,7 @@ optiMuso <- function(measuredData, parameters = NULL, startDate = NULL,
 
     modellOut <- numeric(iterations + 1) # single variable solution
     rmse <- numeric(iterations + 1)
-    origModellOut <- calibMuso(settings=settings,silent=TRUE, skipSpinup = skipSpinup)
+    origModellOut <- calibMuso(settings=settings,silent=TRUE, skipSpinup = skipSpinup,postProcString=postProcString)
 
 
     write.csv(x=origModellOut, file=paste0(pretag,1,".csv"))

RBBGCMuso/R/musoMonte.R

@@ -32,6 +32,7 @@ musoMonte <- function(settings=NULL,
                      keepEpc = FALSE,
                      constrains = NULL,
                      skipZero = TRUE,
+                     postProcString=NULL,
                      ...){
 
 
@@ -86,6 +87,9 @@ musoMonte <- function(settings=NULL,
         outVarNames <- sapply(outVars, musoMapping)
     }
     
+    if(!is.null(postProcString)){
+        outVarNames <- c(outVarNames,gsub("\\s","",unlist(strsplit(procString,"<-"))[1]))
+    }
 
     parameterNames <- gsub("([\\s]|\\-epc)","",parameters[,1],perl=TRUE)
     # settings$calibrationPar <- A[,1] #:LATER:
@@ -145,7 +149,7 @@ musoMonte <- function(settings=NULL,
                              skipSpinup = skipSpinup,
                              keepEpc = keepEpc,
                              debugging = debugging,
-                             outVars = outVars), error = function (e) NA)
+                             outVars = outVars,postProcString=postProcString), error = function (e) NA)
             
             if(length(dim(tmp))>=1){
                 for(j in 1:numVars){

RBBGCMuso/R/musoSensi.R

@@ -31,6 +31,7 @@ musoSensi <- function(monteCarloFile = NULL,
                      plotTitle = "Sensitivity",
                      skipSpinup = TRUE,
                      skipZero = TRUE,
+                     postProcString=NULL,
                      dpi=300){
 
     if(is.null(parameters)){
@@ -48,6 +49,7 @@ musoSensi <- function(monteCarloFile = NULL,
         # browser()
         npar <- ncol(M)-1
         M <- M[which(!is.na(M[,ncol(M)])),]
+        M <- M[-1,]
         y <- M[,(npar+1)]
         colnames(M) <- gsub("\\.epc","-epc",colnames(M))
         M <- M[,colnames(M) %in% parameters[,1]]
@@ -93,7 +95,8 @@ musoSensi <- function(monteCarloFile = NULL,
                       fun = fun,
                       varIndex = varIndex,
                       skipSpinup = skipSpinup,
-                      skipZero=skipZero
+                      skipZero=skipZero,
+                      postProcString=postProcString
                       )
         M <- cbind(seq_along(M[,1]),M)
         yInd <-  grep("mod.", colnames(M))[varIndex]

RBBGCMuso/R/postProcString.R

@@ -0,0 +1,11 @@
+postProcMuso  <- function(modelData, procString){
+    modelDat <- modelData[,-(1:4)]
+    cNames <- colnames(modelData)
+    tocalc <- gsub("(@)(\\d)","modelDat[,\\2]",procString)
+    newVarName <- gsub("\\s","",unlist(strsplit(procString,"<-"))[1])
+    assign(newVarName,eval(parse(text = unlist(strsplit(tocalc,"<-"))[2])))
+    modelData <- cbind.data.frame(modelData,eval(parse(text = newVarName)))
+    colnames(modelData) <- c(cNames,newVarName)
+    modelData
+}
+