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.