楼主: alxzjiayou
2268 1

[问答] 利用EVIM软件包在MATLAB中估计GPD参数出现的问题 [推广有奖]

  • 15关注
  • 0粉丝

已卖:2份资源

博士生

29%

还不是VIP/贵宾

-

威望
0
论坛币
16 个
通用积分
8.9000
学术水平
0 点
热心指数
0 点
信用等级
0 点
经验
1133 点
帖子
8
精华
0
在线时间
525 小时
注册时间
2013-7-3
最后登录
2023-1-6

楼主
alxzjiayou 发表于 2017-3-9 12:45:54 |AI写论文
10论坛币
对于POT模型,我用峰度法选取阈值,然后利用MATLAB中的EVIM软件包估计参数(样本数853),所有文件运行也是在EVIM软件包目录下完成的,出现如下问题
(1)收益率数据乘以100后,即阈值也被放大100倍,运行报错
data=xlsread('data1.xlsx');
out=gpd(data,3.944405445840538,[])
运行结果:
错误使用 fminusub (line 16)
Objective function is undefined at initial point. Fminunc cannot continue.

出错 fminunc (line 439)
   [x,FVAL,GRAD,HESSIAN,EXITFLAG,OUTPUT] = fminusub(funfcn,x, ...

出错 gpd (line 67)
[res.par_ests,res.funval,res.terminated,res.details] = fminunc('negloglikgpd',res.par_ests,opts,excess);

出错 rm1 (line 4)
out=gpd(data,3.944405445840538,[])

(2)收益率序列不乘以100,运行结果如下:
data=xlsread('data1.xlsx');
out=gpd(data,0.039444054458405,[])

运行结果

out =

         par_ests: [-1.248790482490188 0.082325514512838]
           funval: -1.182592767021184e+02
       terminated: 2
          details: [1x1 struct]
           varcov: [2x2 double]
          par_ses: [NaN NaN]
        threshold: 0.039444054458405
             data: [41x1 double]
    p_less_thresh: 0.951934349355217


为何显示: par_ses: [NaN NaN]
本人不才,刚刚学习极值理论,很多东西还不太懂,请各位帮帮忙,非常感谢!

关键词:excess 软件包 收益率 模型 样本

沙发
rajeshy 发表于 2020-12-4 09:48:47
老铁,这个包能提供一下么?我没有找到呢。谢谢啦

您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注cda
拉您进交流群
GMT+8, 2025-12-26 11:05