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

warning fixed in GLUE results, more messages for the users

Browse Source
This commit is contained in:
Roland Hollós 2022-11-22 10:18:50 +01:00
parent 6e75619255
commit 92cbd78854
1 changed files with 4 additions and 1 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

5
RBBGCMuso/R/calibration.R
View File

@ -270,10 +270,13 @@ musoGlue <- function(presCalFile, w, delta = 0.17, settings=setupMuso(), paramet
par(pari) par(pari)
dev.off() dev.off()
maxParValues <- preservedCalibtop5[which.max(preservedCalibtop5$combined),] maxParValues <- unlist(preservedCalibtop5[which.max(preservedCalibtop5$combined),])[1:length(paramIndex)]
maxParIndexes <- paramIndex maxParIndexes <- paramIndex
maxLikelihoodParameters <- data.frame(parameter_index=maxParIndexes,parameter_value=maxParValues)
write.csv(cbind.data.frame(calibrationPar=maxParValues,parameters=maxParIndexes),"maxLikelihood.csv") write.csv(cbind.data.frame(calibrationPar=maxParValues,parameters=maxParIndexes),"maxLikelihood.csv")
cat("\n\n- A file containing the parameters with the maximum likelihood (maxlikelihood_parameters.csv) has been created.\n")
write.csv(optRanges,"optRanges.csv") write.csv(optRanges,"optRanges.csv")
cat("- GLUE interval values have been written into optRanges.csv\n")
# browser() # browser()
# There are some serious problems with this implementation. The uncertainity bouns are not for the parameters, but for the output values. The median is pointwise median for all simulation. # There are some serious problems with this implementation. The uncertainity bouns are not for the parameters, but for the output values. The median is pointwise median for all simulation.
# And the 95 and 5 percentile also. # And the 95 and 5 percentile also.