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

怎么办呢？有哪位高手帮我解决吧，非常感谢

dedongzhuang 发表于 2012-3-12 10:18
结果提示是：

> f(y1)

就是说你的f函数根本就算不出来一个数值，是吧？那问题就是出在dx.y上了。先检查一个dx.y函数，比如先画个图看看。

qoiqpwqr 发表于 2012-3-12 10:55
就是说你的f函数根本就算不出来一个数值，是吧？那问题就是出在dx.y上了。先检查一个dx.y函数，比如先画个 ...

xi<-gpd.y1.fit2$lower.par.ests[1]
beta<-gpd.y1.fit2$lower.par.ests[2]
u<-gpd.y1.fit2$lower.thresh

dx<-function(x){(k/n)*(1-xi*(x-u)/beta)^(-1/xi)}

dx.y<-function(x){dcopula(rcop.bb7.fit,gpd.2p(y1,gpd.y1.fit2),gpd.2p(y2.VaR,gpd.y2.fit2))/dx(x)}
plot(dx.y(y1),type="l")

dedongzhuang 发表于 2012-3-12 11:35
xi

你的dx.y的图能画出来吗？

dedongzhuang 发表于 2012-3-12 11:35
xi

结果出来的密度函数图没有什么异样

qoiqpwqr 发表于 2012-3-12 11:37
你的dx.y的图能画出来吗？

可以画出来，怎么可以发给你看

你可以把整个Splus程序传上来让大家看一看，包括数据

数据怎么能上传啊

前面已经给出了整个程序了

dedongzhuang 发表于 2012-3-12 11:43
可以画出来，怎么可以发给你看

点高级模式，然后上传图片