全部版块 我的主页
› 论坛 › 数据科学与人工智能 › 数据分析与数据科学 › R语言论坛
在线等：请问数值积分出错是什么原因啊？
ssjwzhang
1843 7
2012-10-10
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


请问是怎么回事？怎么解决?
二维码

扫码加我 拉你入群

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

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

全部回复
qoiqpwqr
2012-10-10 19:40:32
可能是你的积分限里面的无穷的问题，换个数试试看。
二维码

扫码加我 拉你入群

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

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

ssjwzhang
2012-10-11 15:34:48
把无穷换成1000还是同样问题，我是想算0到无穷的积分的
二维码

扫码加我 拉你入群

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

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

qoiqpwqr
2012-10-11 19:39:28
var_y不能是一个向量，必须是一个数
二维码

扫码加我 拉你入群

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

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

ssjwzhang
2012-10-12 07:56:55
发现程序粘贴的时候少了下表，程序其实是下面的：
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
        }
二维码

扫码加我 拉你入群

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

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

ssjwzhang
2012-10-12 07:58:56
不好意思啊，版主。其实我是做模拟的，有时候报错，所以想搞明白。附件是完整程序。
macro.doc
大小:(28 KB)

 马上下载





二维码

扫码加我 拉你入群

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

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

点击查看更多内容…
相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

微信QQ空间QQ微博
扫码加好友，拉您进群
各岗位、行业、专业交流群