wavePurifyVx(SpatialVx)
wavePurifyVx()所属R语言包:SpatialVx
Apply traditional forecast verification after wavelet denoising
采用传统的预测验证后,小波消噪
译者:生物统计家园网 机器人LoveR
描述----------Description----------
Apply traditional forecast verification after wavelet denoising ala Briggs and Levine (1997).
应用传统的的预报检验后,小波消噪阿拉·布里格斯和Levine(1997)。
用法----------Usage----------
wavePurifyVx(x, y, object = NULL, climate=NULL, which.stats = c("bias", "ts", "ets", "pod", "far", "f", "hk", "mse"), thresholds = NULL, return.fields = FALSE, verbose = FALSE, ...)
## S3 method for class 'wavePurifyVx'
plot(x, ...)
## S3 method for class 'wavePurifyVx'
summary(object, ...)
参数----------Arguments----------
参数:x,y
For wavePurifyVx, m X n matrices giving the verification and forecast fields, resp. If either is not supplied, then object must be supplied, and the fields given by object are used instead (even if one of x or y are supplied. If one of these is missing, and object not supplied, the function will error out. For plot.wavePurifyVx, list object as output from wavePurifyVx.
对于wavePurifyVx,m×n的矩阵验证和预测领域分别。要么是不提供的,然后object必须提供,并领域的给出的object来代替(即使x或y提供。这些失踪,对于object,列表对象,输出plot.wavePurifyVx。wavePurifyVx不提供的,该功能会出错。
参数:object
For wavePurifyVx, a list object as returned by hoods2dPrep. Not used if both x and y are given. For summary.wavePurifyVx, list object as returned by wavePurifyVx.
wavePurifyVx,返回一个列表对象hoods2dPrep。未使用,如果没有x和y是给定的。对于summary.wavePurifyVx,列表对象返回的wavePurifyVx。
参数:climate
m X n matrix defining a climatology field. If not NULL, then the anamoly correlation coefficient will be applied to the wavelet denoised fields.
M×N矩阵定义一个气候学领域。如果不是NULL,则的anamoly相关系数将施加到小波去噪后字段。
参数:which.stats
character describing which traditional verification statistics to calculate on the wavelet denoised fields. This is the argument passed to the argument of the same name in vxstats.
性格描述,传统的检验统计计算的小波消噪领域。这是参数传递的参数中的同名vxstats。
参数:thresholds
If object is used, then this argument is ignored. Otherwise, it is used to define events for all of the verification statistics except MSE. However, if supplied or other statistics are to be computed, then MSE will be calculated for the fields at values >= thresholds. If only MSE is to be computed, and thresholds is NULL, then no thresholding is applied. If NULL, and any of the statistics besides MSE are to be calculated, then default values of the 0, 0.1, 0.25, 0.33, 0.5, 0.66, 0.75, 0.9 and 0.95 quantiles (for each field, so that the thresholds may differ between fields) are used. The same holds for anamoly correlation coefficient.
如果object使用,则此参数将被忽略。否则,它被用于定义事件可用于所有的检验统计除了MSE。但是,如果提供的或其他的统计数据计算,那么MSE将计算出的值字段> =thresholds。如果只MSE的计算,thresholds是NULL,则没有阈值。如果为NULL,除了MSE的任何的统计计为,那么默认值0,0.1,0.25,0.33,0.5,0.66,0.75,0.9和0.95分位数(为每个字段,从而使阈值可显着差异字段)使用。这同样适用于anamoly相关系数。
参数:return.fields
logical, should the denoised fields be returned (e.g., for subsequent plotting)?
逻辑,应去噪领域退还(例如,为后续的绘制)?
参数:verbose
logical, should progress information (including total run time) be printed to the screen?
逻辑,应打印到屏幕上的信息(包括总运行时间)的进展吗?
参数:...
For wavePurifyVx, optional additional arguments to denoise.dwt.2d (or denoise.modwt.2d) from package waveslim. Note that if the argument J is not passed, then it will be determined as J=log2(min(m,n)). If the fields are dyadic, then the usual DWT is used, otherwise the maximal overlap DWT is used instead. For the plot and summary method functions, these are not used.
wavePurifyVx,可选的附加参数denoise.dwt.2d(或denoise.modwt.2d)从包装waveslim的。请注意,如果参数J不通过,那么将被定为J = LOG2分钟(M,N)。如果信息双值的,那么通常的载重吨被使用,,否则最大重叠载重吨代替。对于图和汇总方法的功能,这些都没有用。
Details
详细信息----------Details----------
If the fields are dyadic, then the denoise.dwt.2d function from package waveslim is applied to each field before calculating the chosen verification statistics. Otherwise denoise.modwt.2d from the same package is used. The result is that high-frequency fluctuations in the two fields are removed before calculating verification statistics so that the resulting statistics are less susceptible to small-scale errors (see Briggs and Levine, 1997). See Percival and Guttorp (1994) and Lindsay et al. (1996) for more on this type of wavelet analysis including maximal overlap DWT.
如果该字段是双值的,然后denoise.dwt.2d从包装waveslim功能被应用于各个领域,然后计算所选择的检验统计。否则denoise.modwt.2d使用相同的包。结果是,在这两个领域的高频率波动中删除,然后再计算检验统计,但所得的统计数字是不容易受到小规模的错误(见布里格斯和Levine,1997年)。 ,请参阅波斯富街和Guttorp的(1994)和林赛等。 (1996)更多关于这种类型的小波分析,包括最大重叠离散小波变换。
值----------Value----------
A list object of class "wavePurifyVx" is returned with possible components (depending on what is supplied in the arguments, etc.):
一个List对象返回的类“wavePurifyVx”可能的组成部分(根据提供的参数等):
参数:X2,Y2
m X n matrices of the denoised verification and forecast fields, resp. (only if 'return.fields' is TRUE).
m×n的矩阵的的去噪验证和预报场分别。 (只有当return.fieldsTRUE)。
参数:thresholds
q X 2 matrix of thresholds applied to the forecast (first column) and verification (second column) fields, resp. If 'climate' is not NULL, then the same thresholds for the forecast field are applied to the climatology.
Q X 2的阈值应用于预测的矩阵(第一列)和验证(第二列)字段,分别。如果“气候”不为NULL,那么相同的阈值的预测领域应用的气候。
参数:qs
If object and thresholds are NULL, and statistics other than MSE or ACC are desired, then this will be created along with the thresholds, and is just a character version of the trhesholds.
如果对象和阈值是NULL,和所需的统计MSE或ACC以外的,那么这将一起创建和设置的阈值,而仅仅是一个字符的trhesholds版本。
参数:args
list object containing all the optional arguments passed into ..., and the value of J used (e.g., even if not passed into ...).
列表对象,其中包含的所有可选参数传递到...,和使用的J值(例如,即使不会进入...)。
参数:bias,ts,ets,pod,far,f,hk,mse,acc
numeric vectors of length q (i.e., the number of thresholds) giving the associated verification statistics.
为长度q的数值向量(即,阈值的数目),给出了相关的检验统计。
(作者)----------Author(s)----------
Eric Gilleland
参考文献----------References----------
IEEE Transactions on Geoscience and Remote Sensing, 34 (3), 771–787.
and P. Kumar, Eds., New York: Academic, pp. 325–343.
参见----------See Also----------
denoise.dwt.2d, denoise.modwt.2d, waverify2d, mowaverify2d, waveIS, hoods2dPrep, vxstats
denoise.dwt.2d,denoise.modwt.2d,waverify2d,mowaverify2d,waveIS,hoods2dPrep,vxstats
实例----------Examples----------
data(UKobs6)
data(UKfcst6)
data(UKloc)
hold <- surrogater2d(UKobs6, n=1, maxiter=50, verbose=TRUE)
hold <- matrix(hold, 256, 256)
UKobj <- hoods2dPrep("UKfcst6", "UKobs6", thresholds=c(0.1,2,5,10), loc=UKloc, units="mm/h")
look <- wavePurifyVx(object=UKobj, climate=hold, return.fields=TRUE, verbose=TRUE)
plot(look)
summary(look)
## Not run: [#不运行:]
data(pert004)
data(pert000)
hold <- hoods2dPrep( "pert004", "pert000", thresholds=c(1,10,50), loc=ICPg240Locs, units="mm/h")
look <- wavePurifyVx(object=hold, return.fields=TRUE, verbose=TRUE)
summary( look)
## End(Not run)[#(不执行)]
转载请注明:出自 生物统计家园网(http://www.biostatistic.net)。
注:
注1:为了方便大家学习,本文档为生物统计家园网机器人LoveR翻译而成,仅供个人R语言学习参考使用,生物统计家园保留版权。
注2:由于是机器人自动翻译,难免有不准确之处,使用时仔细对照中、英文内容进行反复理解,可以帮助R语言的学习。
注3:如遇到不准确之处,请在本贴的后面进行回帖,我们会逐渐进行修订。
|