GEVparameters_from_models {SpatialModelsZAMG} | R Documentation |
this function calculates the GEV parameters according to the linear models
specified by coefficients
GEVparameters_from_models(covariables, coefficients)
covariables |
a named matrix or vector with the covariables. each row corresponds to one location, columns are the corresponding covariables (names of the coefficients) from the linear models; the intercept is added seperately and thus doesn't have to be included |
coefficients |
a list with the named coefficients of the linear models including |
GEVparam |
a matrix or vector with the GEV parameters. each row corresponds to one location, columns are loc, scale and shape |
GEVparameters_map
, returnlevels
, cond_returnlevels
, returnlevel_plot
,
cond_returnlevel_plot
# get covariables lon = get(data(lon.at)) lon = as.vector(t(lon)) lat = get(data(lat.at)) lat = as.vector(t(lat)) alt = get(data(alt.at)) alt = as.vector(t(alt)) x = get(data(sample_grid_data)) mdday = x$mdday sd_mmax = x$mmsd swe_mmax = x$mmswe # take only locations of the (Austrian) domain lon = lon[which(!is.na(mdday))] lat = lat[which(!is.na(mdday))] alt = alt[which(!is.na(mdday))] mdday = mdday[which(!is.na(mdday))] sd_mmax = sd_mmax[which(!is.na(sd_mmax))] swe_mmax = swe_mmax[which(!is.na(swe_mmax))] # define matrix 'covariables' covariables = cbind("lon" = lon, "lat" = lat, "alt" = alt, "mdday" = mdday, "sd_mmax" = sd_mmax, "swe_mmax" = swe_mmax) # load function output from optimizer_biv_hr_model data("optim_hr") # define coefficients sd_coeff = optim_hr$coefficients$sd_coeff swe_coeff = optim_hr$coefficients$swe_coeff cor_coeff = optim_hr$coefficients$cor_coeff # Calculate the GEV parameters from linear models sd_GEVparam = GEVparameters_from_models(covariables = covariables, coefficients = sd_coeff) swe_GEVparam = GEVparameters_from_models(covariables = covariables, coefficients = swe_coeff)