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

back to the original

Browse Source
This commit is contained in:
hollorol 2019-02-17 20:23:38 +01:00
parent 2318804ea0
commit 4ab2b60f49
1 changed files with 2 additions and 2 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

4
RBBGCMuso/R/musoRand.R
View File

@ -171,10 +171,10 @@ musoRand <- function(parameters, constrains = NULL, iterations=3000){
h <- c(Gh0$h,h) h <- c(Gh0$h,h)
E <- do.call(rbind,lapply(Ef,function(x){x$E})) E <- do.call(rbind,lapply(Ef,function(x){x$E}))
f <- do.call(c,lapply(Ef,function(x){x$f})) f <- do.call(c,lapply(Ef,function(x){x$f}))
randVal <- suppressWarnings(limSolve::xsample(G=G,H=h,E=E,F=f,iter = iterations, type = "rda"))$X randVal <- suppressWarnings(limSolve::xsample(G=G,H=h,E=E,F=f,iter = iterations))$X
} else{ } else{
Gh0<-genMat0(dependences) Gh0<-genMat0(dependences)
randVal <- suppressWarnings(xsample(G=Gh0$G,H=Gh0$h, iter = iterations, type = "rda"))$X randVal <- suppressWarnings(xsample(G=Gh0$G,H=Gh0$h, iter = iterations))$X
} }
results <- list(INDEX =dependences$INDEX, randVal=randVal) results <- list(INDEX =dependences$INDEX, randVal=randVal)