This function calculates the posterior predictive density at points x based on a matrix of posterior density parameters
.gev.postpred(x, posterior, Nyr = 100, type = c("density", "quantile"))a vector of values for the posterior predictive density or quantile at x, depending on the value of type
n vector of points
n by 3 matrix of posterior samples
number of years to extrapolate
string indicating whether to return the posterior density or the quantile.