如何在 R 中选择峰值超过阈值方法的阈值

Question

我有一个财务数据的数据集，并已将其转换为日志返回，因此所有数据都在-1和1之间。

我正在尝试使用以下代码：

getSymbols("^IXIC",from="2016-11-11",to="2018-11-08")
index<-as.vector(IXIC$IXIC.Close)
logret <- diff(log(index))

使用包绘制不同的阈值

library("ismev")
library("evir")
gpd.fitrange(logreturns, umin=0,umax=1)

但它给了我一个错误，我不知道如何编码如何选择合适的阈值。如何选择合适的阈值？

Answer 1

我更喜欢POT包中的tcplot函数。它绘制一个阈值选择图，即它估计阈值不同值的比例和形状参数。 然后，您可以使用这些图获得良好的阈值：良好的阈值的特征是具有低方差的稳定估计值。

POT::tcplot(logret, which = 1, nt = 25, u.range = range(logret)* c(1, 0.99), ask = F)
POT::tcplot(logret, which = 2, nt = 25, u.range = range(logret)* c(1, 0.99), ask = F)

第一个函数绘制相对于刻度参数的阈值选择图，第二个函数绘制相对于形状参数的阈值选择图。nt指定图中的点数，u.range可以根据您的数据指定。