R积分出现"non-finite function value"错误

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

用R计算函数积分时，出现了“non-finite function value”的错误。
具体代码如下：

1. TT <- c(0.5120,0.5990,0.6860,0.7300,0.8170,0.9040,0.9920,1.0790,1.1660,1.2530)
   
2. 
   
3. C1_B <- function(T){integrate(function(tau)f(tau),lower=0.01*T,upper=0.99*T)$value}
   
4. 
   
5. f <- function(tau) {integrate(function(tau1)sqrt(1/(tau-tau1)),lower=0.01*tau,upper=0.99*tau)$value}
   
6. 
   
7. C <- matrix(rep(0,length(TT)*2),nrow=length(TT),ncol=2)
   
8. for (j in 1:length(TT)){
   
9.     T=TT[j]
   
10.     C[j,] <- c(T,C1_B(T))
    
11. }

复制代码

请教大家该如何解决上述代码里的问题呢？先谢谢了！

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：function finite value nite alue function

f <- function(tau) {integrate(function(tau1)sqrt(1/(tau-tau1)),lower=0.01*tau,upper=0.99*tau)$value}
##这条命令错了，没有规定取值范围，根号下面不能为负数，也就是说tau必须大于tau1

取值可以定为c(-Inf,tau)

机械唯物主义123 发表于 2016-12-24 17:47
f

谢谢你的回复。但你的回复可能不对。
我已经设置tau1的范围在0.01*tau，0.99*tau
那么tau-tau1就在0.99*tau，0.01*tau之间，而tau>0的。

抱歉，我没看出来，不过问题肯定出在根号上，大概是下面那样

Error in integrate(function(tau1) sqrt(1/(tau - tau1)), lower = 0.01 *  :
  non-finite function value
此外: Warning message:
In sqrt(1/(tau - tau1)) : 产生了NaNs