全部版块 我的主页
论坛 计量经济学与统计论坛 五区 计量经济学与统计软件 经管代码库
1277 1
2017-06-17
利用Copula-EVT方法求解CoVaR,利用S-plus软件求解积分公式见图片

代码如下:
dx<-function(x)   (k/n)*(1/beta)*(1-xi*(x-u)/beta)^(-1-1/xi)
dx.y<-function(x)
     dcopula(cop.normal.fit,gpd.2p(x,rm.gpd),gpd.2p(rn.VaR,rn.gpd))*dx(x)
f<-function(x)
{
        out<-0*x
        for(i in 1:length(x))
                out<-integrate(dx.y,lower=-101,upper=x)$integral-0.05
out
}
CoVaR=uniroot(f,c(-100,0))$root

请问为何积分下限为负无穷,这里写为-101,这是参考别人论文附上的程序。

附件列表
QQ图片20170617205523.png

原图尺寸 5.06 KB

QQ图片20170617205523.png

二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

全部回复
2017-9-11 15:04:05
您好,请问S-PLUS的软件包可否共享?不胜感激
二维码

扫码加我 拉你入群

请注明:姓名-公司-职位

以便审核进群资格,未注明则拒绝

相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群