Splus 中的integrate 怎么用 急需 赏金10元

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

有程序如下：
> f <- function(x)
{
        out <- 0 * x
        for(i in 1:length(x))
                out <- integrate(dx.y, lower = -90, upper = x)$integral -0.05
        out
}
> CoVaR = uniroot(f, c(-72, 0))$root
结果：
Problem in qf15(f, lower, upper, aux = optargs(list..: f is not vectorized - length(f(x))!=length(x) Use traceback() to see the call stack请问是哪里出错了呢？急

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：integrate Rate PLUS Inte int aux function 程序

dx.y是什么？

#dx.y表示f(x|y)是x关于y=us.VaR的密度函数
dx.y<-function(x)
dcopula(cop.bb7.fit,gpd.2p(x,foreign.gpd),gpd.2p(us.VaR,us.gpd))*dx(x)

qoiqpwqr 发表于 2012-3-12 01:08
dx.y是什么？

#参数设定
xi<-foreign.gpd$lower.par.ests[1]
beta<-foreign.gpd$lower.par.ests[2]
u<-foreign.gpd$lower.thresh
#x代表foreign，dx表示x的密度函数
dx<-function(x) (k/n)*(1/beta)*(1-xi*(x-u)/beta)^(-1-1/xi)
#dx.y表示f(x|y)是x关于y=us.VaR的密度函数
dx.y<-function(x)
dcopula(cop.bb7.fit,gpd.2p(x,foreign.gpd),gpd.2p(us.VaR,us.gpd))*dx(x)
#求解显著水平5%下的foreign的CoVaR
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
#显著水平从0.1-0.01变化时，%CoVaR值
pCoVaR=0
fq<-function(x,q){
out <- 0*x
for(i in 1:length(x))
out<-integrate(dx.y,lower=-101,upper=x)$integral-q
out
}
for(i in 1:10){
j=i/100
uVaR=gpd.2q(j,us.gpd)
fVaR=gpd.2q(j,foreign.gpd)
xi<-foreign.gpd$lower.par.ests[1]
beta<-foreign.gpd$lower.par.ests[2]
u<-foreign.gpd$lower.thresh
#x代表foreign，dx表示x的密度函数
dx<-function(x) (k/n)*(1/beta)*(1-xi*(x-u)/beta)^(-1-1/xi)
#dx.y表示f(x|y)是x关于y=us.VaR的密度函数
dx.y<-function(x)
dcopula(cop.bb7.fit,gpd.2p(x,foreign.gpd),gpd.2p(uVaR,us.gpd))*dx(x)
CoVaR=uniroot(fq,c(-100,0),q=j)$root
pCoVaR=(CoVaR-fVaR)/fVaR
}

能加你qq联系嘛？

我的邮箱是zhuangdedong0925@163.com

qoiqpwqr 发表于 2012-3-12 01:08
dx.y是什么？
#参数设定
xi<-foreign.gpd$lower.par.ests[1]
beta<-foreign.gpd$lower.par.ests[2]
u<-foreign.gpd$lower.thresh
#x代表foreign，dx表示x的密度函数
dx<-function(x) (k/n)*(1/beta)*(1-xi*(x-u)/beta)^(-1-1/xi)
#dx.y表示f(x|y)是x关于y=us.VaR的密度函数
dx.y<-function(x)
dcopula(cop.bb7.fit,gpd.2p(x,foreign.gpd),gpd.2p(us.VaR,us.gpd))*dx(x)
#求解显著水平5%下的foreign的CoVaR
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
#显著水平从0.1-0.01变化时，%CoVaR值
pCoVaR=0
fq<-function(x,q){
out <- 0*x
for(i in 1:length(x))
out<-integrate(dx.y,lower=-101,upper=x)$integral-q
out
}
for(i in 1:10){
j=i/100
uVaR=gpd.2q(j,us.gpd)
fVaR=gpd.2q(j,foreign.gpd)
xi<-foreign.gpd$lower.par.ests[1]
beta<-foreign.gpd$lower.par.ests[2]
u<-foreign.gpd$lower.thresh
#x代表foreign，dx表示x的密度函数
dx<-function(x) (k/n)*(1/beta)*(1-xi*(x-u)/beta)^(-1-1/xi)
#dx.y表示f(x|y)是x关于y=us.VaR的密度函数
dx.y<-function(x)
dcopula(cop.bb7.fit,gpd.2p(x,foreign.gpd),gpd.2p(uVaR,us.gpd))*dx(x)
CoVaR=uniroot(fq,c(-100,0),q=j)$root
pCoVaR=(CoVaR-fVaR)/fVaR
}

把f函数里面的$integral改成$value试试看

在R里面是这样的。

qoiqpwqr 发表于 2012-3-12 08:15
把f函数里面的$integral改成$value试试看

在R里面是这样的。

还是不行，在splus中应该写成$integral

qoiqpwqr 发表于 2012-3-12 08:15
把f函数里面的$integral改成$value试试看

在R里面是这样的。

结果提示是：

> f(y1)
Error: Expressions nested beyond limit (256)
        only 29 of 88 frames dumped