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

fixing musoSensi bug

Browse Source
This commit is contained in:
Roland Hollós 2019-11-18 13:20:37 +01:00
parent 4cc5e5b371
commit d90a29f017
1 changed files with 3 additions and 0 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

3
RBBGCMuso/R/musoSensi.R
View File

@ -44,9 +44,11 @@ musoSensi <- function(monteCarloFile = NULL,
}} }}
doSensi <- function(M){ doSensi <- function(M){
browser()
npar <- ncol(M)-1 npar <- ncol(M)-1
M <- M[which(!is.na(M[,ncol(M)])),] M <- M[which(!is.na(M[,ncol(M)])),]
y <- M[,(npar+1)] y <- M[,(npar+1)]
colnames(M) <- gsub("\\.epc","-epc",colnames(M))
M <- M[,colnames(M) %in% parameters[,1]] M <- M[,colnames(M) %in% parameters[,1]]
npar <- ncol(M) npar <- ncol(M)
M <- apply(M[,1:npar],2,function(x){x-mean(x)}) M <- apply(M[,1:npar],2,function(x){x-mean(x)})
@ -103,6 +105,7 @@ musoSensi <- function(monteCarloFile = NULL,
yInd <- grep("mod.", colnames(M))[varIndex] yInd <- grep("mod.", colnames(M))[varIndex]
parNames <- grep("mod.",colnames(M), invert=TRUE, value = TRUE) parNames <- grep("mod.",colnames(M), invert=TRUE, value = TRUE)
M <- M[,c(grep("mod.", colnames(M),invert=TRUE),yInd)] M <- M[,c(grep("mod.", colnames(M),invert=TRUE),yInd)]
# browser()
return(doSensi(M)) return(doSensi(M))
} }
} }