PAMfmado: Clustering of maxima

In ExtremalDep: Extremal Dependence Models

Clustering of maxima

Description

Clustering of times series of maxima based on the pam package tailored for the F-madogram distance

Usage

PAMfmado(x, K, J = 0, threshold = 0.99, max.min = 0)

Arguments

x	x a matrix of maxima. For example, number of stations = ncol(x) and time series length = nrow(x) for weekly maxima of precipitation.
K	number of clusters
J	number of resampling for which the stations are randomly moved to break the dependence. By default, J=0 means no resampling.
threshold	Threshold corresponding to the quantile level for the resampling. The resulting quantile is printed (when J does not take value 0).
max.min	A threshold to remove very small values. For example, some raingauges cannot go below 2 mm. By default, max.min=0.

Value

an object of class "pam" representing the clustering. See ?pam.object for details.

Author(s)

Philippe Naveau

References

Bernard E., Naveau P., Vrac M. and Mestre O. (2013). Clustering of maxima: Spatial dependencies among heavy rainfall in France. Journal of Climate 26, 7929–7937.

Naveau, P., A. Guillou, D. Cooley, and J. Diebolt (2009). Modeling pairwise dependence of maxima in space. Biometrika 96(1).

Cooley, D., P. Naveau, and P. Poncet (2006). Variograms for spatial max-stable random fields. In: Bertail, P., Soulier, P., Doukhan, P. (eds) Dependence in Probability and Statistics. Lecture Notes in Statistics, vol 187. Springer, New York, NY .

Reynolds, A., Richards, G., de la Iglesia, B. and Rayward-Smith, V. (1992). Clustering rules: A comparison of partitioning and hierarchical clustering algorithms. Journal of Mathematical Modelling and Algorithms 5, 475-504.

See Also

See the function as pam in the package cluster

Examples

data(PrecipFrance)
attach(PrecipFrance)
PAMmado <- PAMfmado(precip,7) 

ExtremalDep documentation built on April 4, 2025, 1:44 a.m.