R语言POT超阈值模型和极值理论EVT分析

简介: R语言POT超阈值模型和极值理论EVT分析

本文依靠EVT对任何连续分布的尾部建模。尾部建模,尤其是POT建模,对于许多金融和环境应用至关重要。

POT模型其主要动机是为高洪水流量的概率模型提供实用工具。但是,EVT的优势在于结果不取决于要建模的过程。因此,人们可以使用POT来分析降水,洪水,金融时间序列,地震等。

 

特征

POT软件包可以执行单变量和双变量极值分析;一阶马尔可夫链也可以考虑。例如,目前使用18个 估算器拟合(单变量)GPD  。这些估算器依靠三种不同的技术:

  • 极大似然:MLE,LME,MPLE
  • 动量法:MOM,PWM,MED
  • 距离最小化:MDPD和MGF估计器。

与单变量情况相反,没有用于对超过阈值的双变量超出进行建模的有限参数化。POT允许对双变量GPD进行6种参数化:对数模型,负对数模型和混合模型-以及它们各自的不对称版本。

最后,可以使用二元GPD拟合一阶马尔可夫链,以实现两个连续观测值的联合分布。

在本节中,我们明确介绍了软件包中一些最有用的功能。 但是,对于完整的描述,用户可能希望查看软件包的小插图和软件包的html帮助。


GPD 计算:


模拟来自GPD(0,1,0.2)的样本:
x <- rgpd(100, 0, 1, 0.2)
##评估x = 3时的密度和不超过的概率:
dgpd(3, 0, 1, 0.2); pgpd(3, 0, 1, 0.2)


#计算非超出概率为0.95的分位数:


qgpd(0.95, 0, 1, 0.2)


y <- rbvgpd(100, mo
##评估不超过(5,14)的可能性


pbvgpd(c(3,15), mode
GPD 拟合


##最大似然估计(阈值= 0):
mle <- fgpd(x, 0)
##最大似然估计(阈值= 0):
pwu <- fgpd(x, 0, "pwmu")
##最大拟合优度估算器:
adr <- fgpd(x, 0, "mgf"
##指定已知参数:
fgpd(x, 0, "mple",
##指定数值优化的起始值:
fgpd(x, 0, "mdpd", start =
##拟合具有逻辑依存关系的双变量GPD:
log <- fitbv

绘图用于单变量和双变量情况的通用函数:


plot(mle); plot(log)

返回等级图:

 

概率图和QQ图


qq(mle)

绘制密度

绘制Pickands的依赖函数:

 

光谱密度图:



##对数似然(分位数):confint(mle, prob = 0.95)



##对数似然(参数):
confint(mle, "shape")

 

相关文章
|
7天前
|
资源调度 前端开发 数据可视化
R语言GARCH模型对股市sp500收益率bootstrap、滚动估计预测VaR、拟合诊断和蒙特卡罗模拟可视化
R语言GARCH模型对股市sp500收益率bootstrap、滚动估计预测VaR、拟合诊断和蒙特卡罗模拟可视化
27 2
|
11天前
POT超阈值模型和极值理论EVT分析
POT超阈值模型和极值理论EVT分析
|
3天前
MATLAB用GARCH-EVT-Copula极值理论模型VaR预测分析股票投资组合
MATLAB用GARCH-EVT-Copula极值理论模型VaR预测分析股票投资组合
|
4天前
|
数据可视化 数据建模
【视频】R语言极值理论EVT:基于GPD模型的火灾损失分布分析|数据分享(下)
【视频】R语言极值理论EVT:基于GPD模型的火灾损失分布分析|数据分享
|
4天前
|
NoSQL 数据处理
【视频】R语言极值理论EVT:基于GPD模型的火灾损失分布分析|数据分享(上)
【视频】R语言极值理论EVT:基于GPD模型的火灾损失分布分析|数据分享
11 0
|
4天前
|
数据挖掘 C语言
极值理论 EVT、POT超阈值、GARCH 模型分析股票指数VaR、条件CVaR:多元化投资组合预测风险测度分析
极值理论 EVT、POT超阈值、GARCH 模型分析股票指数VaR、条件CVaR:多元化投资组合预测风险测度分析
11 1
|
6天前
|
前端开发
数据分享|R语言零膨胀泊松回归ZERO-INFLATED POISSON(ZIP)模型分析露营钓鱼数据实例估计IRR和OR
数据分享|R语言零膨胀泊松回归ZERO-INFLATED POISSON(ZIP)模型分析露营钓鱼数据实例估计IRR和OR
37 4
|
6天前
|
机器学习/深度学习 算法 vr&ar
PYTHON用时变马尔可夫区制转换(MARKOV REGIME SWITCHING)自回归模型分析经济时间序列
PYTHON用时变马尔可夫区制转换(MARKOV REGIME SWITCHING)自回归模型分析经济时间序列
16 4
|
7天前
|
算法 数据可视化
R语言马尔可夫MCMC中的METROPOLIS HASTINGS,MH算法抽样(采样)法可视化实例
R语言马尔可夫MCMC中的METROPOLIS HASTINGS,MH算法抽样(采样)法可视化实例
26 1
R语言马尔可夫MCMC中的METROPOLIS HASTINGS,MH算法抽样(采样)法可视化实例
|
12天前
|
算法 Linux Python
R语言BUGS序列蒙特卡罗SMC、马尔可夫转换随机波动率SV模型、粒子滤波、Metropolis Hasting采样时间序列分析
R语言BUGS序列蒙特卡罗SMC、马尔可夫转换随机波动率SV模型、粒子滤波、Metropolis Hasting采样时间序列分析
40 8

热门文章

最新文章

http://www.vxiaotou.com