PAMfmado {ExtremalDep} | R Documentation |
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, |
threshold |
Threshold corresponding to the quantile level for the resampling. The resulting quantile is printed (when |
max.min |
A threshold to remove very small values. For example, some raingauges cannot go below 2 mm. By default, |
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)