| exceedProbsAggregated {lgcp} | R Documentation |
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 |
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.
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])