exceedProbsAggregated: exceedProbsAggregated function
Description
NOTE THIS FUNCTION IS IN TESTING AT PRESENT
Usage
exceedProbsAggregated(threshold, lg = NULL, lastonly = TRUE)
Arguments
threshold
vector of threshold levels for the indicator function
lg
an object of class aggregatedPredict
lastonly
logical, whether to only compute the exceedances for the last time point. default is TRUE
Value
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])
Details
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.