qqgev(SpatialExtremes)
qqgev()所属R语言包:SpatialExtremes
QQ-plot for the GEV parameters
QQ-GEV参数积为
译者:生物统计家园网 机器人LoveR
描述----------Description----------
This function compares the GEV parameters estimated separately for each location to the ones estimated from a fitted spatial model.
此功能比较GEV分别为每个位置估计从合身的空间模型的参数估计。
用法----------Usage----------
qqgev(fitted, xlab, ylab, ...)
参数----------Arguments----------
参数:fitted
An object of class maxstab or spatgev. Most often, this will be the output of fitmaxstab, fitcovmat, fitcovariance or fitspatgev.
对象的类maxstab或spatgev。大多数情况下,这将是输出fitmaxstab,fitcovmat,fitcovariance或fitspatgev。
参数:xlab,ylab
The x and y-axis labels. May be missing.
的x和y轴的标签。可能会丢失。
参数:...
Optional arguments to be passed to the plot function.
到传递给plot功能的可选参数。
值----------Value----------
A QQ-plot.
一个QQ的图。
(作者)----------Author(s)----------
Mathieu Ribatet
参考文献----------References----------
Fields. Extremes 5:1, 33–44.
multivariate and spatial extreme values: Properties and inference. Biometrika 90(1):139–156.
extremes. Unpublished manuscript.
参见----------See Also----------
qqextcoeff
qqextcoeff
实例----------Examples----------
##Define the coordinate of each location[#定义的每个位置的坐标]
n.site <- 30
locations <- matrix(runif(2*n.site, 0, 10), ncol = 2)
colnames(locations) <- c("lon", "lat")
##Simulate a max-stable process - with unit Frechet margins[#模拟一个最大稳定的过程 - 与单位的Frechet空间]
data <- rmaxstab(50, locations, cov.mod = "gauss", cov11 = 100, cov12 =
25, cov22 = 220)
##Now define the spatial model for the GEV parameters[#GEV参数定义的空间模型]
param.loc <- -10 + 2 * locations[,2]
param.scale <- 5 + 2 * locations[,1] + locations[,2]^2
param.shape <- rep(0.2, n.site)
##Transform the unit Frechet margins to GEV [#变换单位的Frechet空间GEV]
for (i in 1:n.site)
data[,i] <- frech2gev(data[,i], param.loc[i], param.scale[i], param.shape[i])
##Define a model for the GEV margins to be fitted[#定义一个模型,的GEV利润的安装]
##shape ~ 1 stands for the GEV shape parameter is constant[#的形状~1的看台为GEV形状参数是常数]
##over the region[#以上的区域]
loc.form <- loc ~ lat
scale.form <- scale ~ lon + I(lat^2)
shape.form <- shape ~ 1
fitted <- fitspatgev(data, locations, loc.form = loc.form, scale.form =
scale.form, shape.form = shape.form)
qqgev(fitted)
转载请注明:出自 生物统计家园网(http://www.biostatistic.net)。
注:
注1:为了方便大家学习,本文档为生物统计家园网机器人LoveR翻译而成,仅供个人R语言学习参考使用,生物统计家园保留版权。
注2:由于是机器人自动翻译,难免有不准确之处,使用时仔细对照中、英文内容进行反复理解,可以帮助R语言的学习。
注3:如遇到不准确之处,请在本贴的后面进行回帖,我们会逐渐进行修订。
|