NOTE THIS FUNCTION IS IN TESTING AT PRESENT
exceedProbsAggregated(threshold, lg = NULL, lastonly = TRUE)vector of threshold levels for the indicator function
an object of class aggregatedPredict
logical, whether to only compute the exceedances for the last time point. default is TRUE
a function of Y that computes the indicator function I(exp(Y)>threshold) evaluated for each cell of a matrix Y, but with values aggregated to regions If several tresholds are specified an array is returned with the [,,i]th slice equal to I(exp(Y)>threshold[i])
This function computes regional exceedance probabilities after MCMC has finished, it requires the information to have been dumped to disk, and to have been computed using the function lgcpPredictAggregated $$P[\exp(Y_{t_1:t_2})>k],$$ that is the probability that the relative risk exceeds threshold \(k\). Note that it is possible to pass vectors of tresholds to the function, and the exceedance probabilities will be computed for each of these.