profile2d(SpatialExtremes)
profile2d()所属R语言包:SpatialExtremes
Method for profiling (in 2d) fitted max-stable objects
分析方法(在2D)合身的最大的稳定的物体
译者:生物统计家园网 机器人LoveR
描述----------Description----------
Computes profile surfaces for fitted max-stable models.
计算型材表面为合身的最大稳定模型。
用法----------Usage----------
## S3 method for class 'maxstab'
profile2d(fitted, params, ranges, n = 10, plot = TRUE,
...)
参数----------Arguments----------
参数:fitted
An object of class “maxstab”. Most often, it will be the output of the function fitmaxstab.
的对象类“maxstab”。大多数情况下,这将是功能fitmaxstab的输出。
参数:params
A character vector giving the two model parameters that are to be profiled.
给两个模型参数是要分析一个字符向量。
参数:ranges
A matrix corresponding to the ranges for the profiled model parameters that must be explored. Each row corresponds to one model parameter range.
对应的矩阵的分析模型参数,必须探索的范围。每一行对应于一个模型参数范围。
参数:n
Integer. The number of profiled model parameter that must be considered.
整数。 ,异形模型参数的数目是必须考虑的。
参数:plot
Logical. If TRUE (default), the profile surface is plotted.
逻辑。如果TRUE(默认),在型材表面绘制。
参数:...
Extra options that must be passed to the plot function.
额外的选项,必须通过plot功能。
值----------Value----------
A list with two arguments: coord and llik. coord is a matrix representing the grid where the profiled model parameters are fixed. llik the corresponding pairwise log-likelihood.
列表有两个参数:coord和llik。 coord是一个矩阵的网格,其中的异形模型的参数是固定的。 llik相应的成对数似然。
警告----------Warnings----------
This function can be really time consuming!
此功能可真是费时!
(作者)----------Author(s)----------
Mathieu Ribatet
实例----------Examples----------
## Not run: [#不运行:]
##Define the coordinates of each location[#定义的每个位置的坐标]
n.site <- 30
locations <- matrix(rnorm(2*n.site, sd = sqrt(.2)), ncol = 2)
colnames(locations) <- c("lon", "lat")
##Simulate a max-stable process - with unit Frechet margins[#模拟一个最大稳定的过程 - 与单位的Frechet空间]
data <- rmaxstab(30, locations, cov.mod = "whitmat", nugget = 0, range = 30,
smooth = 0.5)
##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 + (lat^2)
shape.form <- shape ~ 1
##Fit a max-stable process[#适合一个最大稳定的过程]
## 1- using the Schlather representation[#1 - 使用Schlather表示]
fitted <- fitmaxstab(data, locations, "whitmat", loc.form, scale.form,
shape.form)
##Plot the profile pairwise log-likelihood for the smooth parameter[#图的档案成对的顺利参数的对数似然]
ranges <- rbind(c(9,11), c(.3, .8))
profile2d(fitted, c("range", "smooth"), ranges = ranges)
## End(Not run)[#(不执行)]
转载请注明:出自 生物统计家园网(http://www.biostatistic.net)。
注:
注1:为了方便大家学习,本文档为生物统计家园网机器人LoveR翻译而成,仅供个人R语言学习参考使用,生物统计家园保留版权。
注2:由于是机器人自动翻译,难免有不准确之处,使用时仔细对照中、英文内容进行反复理解,可以帮助R语言的学习。
注3:如遇到不准确之处,请在本贴的后面进行回帖,我们会逐渐进行修订。
|