From 92cbd78854cb972902880823f2268238ec832cfd Mon Sep 17 00:00:00 2001 From: =?UTF-8?q?Roland=20Holl=C3=B3s?= Date: Tue, 22 Nov 2022 10:18:50 +0100 Subject: [PATCH] warning fixed in GLUE results, more messages for the users --- RBBGCMuso/R/calibration.R | 5 ++++- 1 file changed, 4 insertions(+), 1 deletion(-) diff --git a/RBBGCMuso/R/calibration.R b/RBBGCMuso/R/calibration.R index d8f76fd..c3c8e45 100644 --- a/RBBGCMuso/R/calibration.R +++ b/RBBGCMuso/R/calibration.R @@ -270,10 +270,13 @@ musoGlue <- function(presCalFile, w, delta = 0.17, settings=setupMuso(), paramet par(pari) dev.off() - maxParValues <- preservedCalibtop5[which.max(preservedCalibtop5$combined),] + maxParValues <- unlist(preservedCalibtop5[which.max(preservedCalibtop5$combined),])[1:length(paramIndex)] maxParIndexes <- paramIndex + maxLikelihoodParameters <- data.frame(parameter_index=maxParIndexes,parameter_value=maxParValues) 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") + cat("- GLUE interval values have been written into optRanges.csv\n") # 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. # And the 95 and 5 percentile also.