在线等：请问数值积分出错是什么原因啊？

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

integ_fw<-0
        for (i in 1:3){
               
                fw<-function(t) c*lambda^r/gamma(r)*(c*t+k-1)^(r-1)*exp(-lambda*(c*t+k-1))*(1/(var_y+t))
        
        integ_fw<-integrate(fw,0,Inf)$value
        }

上面一段程序中，fw函数表达式中除t外均是已知常数，var_y是一个已知的3维向量，错误提示是：
Error in integ_fw<-integrate(fw,0,Inf)$value :
  non-finite function value


请问是怎么回事？怎么解决?

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：是什么原因 数值积分 在线等 Inte int 原因 积分 function 表达式 Error

可能是你的积分限里面的无穷的问题，换个数试试看。

把无穷换成1000还是同样问题，我是想算0到无穷的积分的

var_y不能是一个向量，必须是一个数

发现程序粘贴的时候少了下表，程序其实是下面的：
integ_fw<-0
        for (i in 1:3){
               
                fw<-function(t) c*lambda^r/gamma(r)*(c*t+k-1)^(r-1)*exp(-lambda*(c*t+k-1))*(1/(var_y[i]+t))
        
        integ_fw[i]<-integrate(fw,0,Inf)$value
        }

不好意思啊，版主。其实我是做模拟的，有时候报错，所以想搞明白。附件是完整程序。 macro.doc (28 KB)

2012-10-12 08:28:19 上传

你的函数是不是要返回integ_fw的值？

我试了一下，模拟了很多次都没有任何问题啊。

版主对的，我也不太明白，这个程序是从原来的大程序中抽出来改了一下的，可能原来的程序中某个小问题没发现的。现在这样确实可以啊，谢谢！