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

图片

这张才对

第二张是dx.y算出了的数值图

看样子dx.y不是一个连续的函数啊，这个用integrate来积分可能会有问题。

哦，原来这样，那我怎么办呢

我把超过阈值的数据y3放进dx.yh<-function(x){
                a=0
                j=1
                for(i in 1:length(x))
  if(x<u){a[j]=x
                        j=j+1}
                a}
y3<-h(y1)
dx.y<-function(x){dcopula(rcop.bb7.fit,gpd.2p(y3,gpd.y1.fit2),gpd.2p(y2.VaR,gpd.y2.fit2))*dx(x)}
plot(dx.y(y3),type="l")得到的图形是这样的

我应该怎么办，运行下去的时候程序出现
> f <- function(x)
{
        out <- 0 * x
        for(i in 1:length(x))
                out[i] <- integrate(dx.y, lower = -90, upper = x[i])$integral -
                        0.05
        out
}
> CoVaR = uniroot(f, c(-71, 0))$root

Problem in dcopula(rcop.bb7.fit, gpd.2p(y3, gpd.y1...: length of longer operand (71) should be a multiple of length of shorter (15)
Use traceback() to see the call stack

楼上的还在吗

看样子问题在你调用dcopula的时候两个向量的长度不一致，单独运行一下这个函数试试。

这两个变量需要怎么设置呢