1

我正在尝试使用 R（类似于名为 XIRR 的 Excel 函数）计算非定期现金流的 IRR（内部收益率）。唯一的问题是它不起作用。我当然可以使用一些帮助来弄清楚出了什么问题..

R>cashflow <- c(275,275,275,275,275,275,275,275,275,275,275,275,416.833333333333,416.833333333333,416.833333333333,416.833333333333,416.833333333333,416.833333333333,416.833333333333,416.833333333333,416.833333333333,416.833333333333,416.833333333333,416.833333333333,452.666666666667,452.666666666667,452.666666666667,452.666666666667,452.666666666667,452.666666666667,452.666666666667,452.666666666667,452.666666666667,452.666666666667,452.666666666667,452.666666666667,534.583333333333,534.583333333333,534.583333333333,534.583333333333,534.583333333333,534.583333333333,534.583333333333,534.583333333333,534.583333333333,534.583333333333,534.583333333333,534.583333333333,-9246.969046)
R>dates <- c("1997-01-31","1997-02-28","1997-03-31","1997-04-30","1997-05-31","1997-06-30","1997-07-31","1997-08-31","1997-09-30","1997-10-31","1997-11-30","1997-12-31","1998-01-31","1998-02-28","1998-03-31","1998-04-30","1998-05-31","1998-06-30","1998-07-31","1998-08-31","1998-09-30","1998-10-31","1998-11-30","1998-12-31","1999-01-31","1999-02-28","1999-03-31","1999-04-30","1999-05-31","1999-06-30","1999-07-31","1999-08-31","1999-09-30","1999-10-31","1999-11-30","1999-12-31","2000-01-31","2000-02-29","2000-03-31","2000-04-30","2000-05-31","2000-06-30","2000-07-31","2000-08-31","2000-09-30","2000-10-31","2000-11-30","2000-12-31","2001-01-31")

R>irr <- xirr(cashflow, dates)

它抛出错误：

Error in while (pchg >= tol & abs(fval) > tol & iter <= itmax) { : missing value where TRUE/FALSE needed

这是功能：

xirr <- function(cf, dates) {

 # Secant method.
secant <- function(par, fn, tol=1.e-07, itmax = 100, trace=FALSE, ...) { 
    # par = a starting vector with 2 starting values 
    # fn = a function whose first argument is the variable of interest 
    if (length(par) != 2) stop("You must specify a starting parameter vector of length 2") 
    p.2 <- par[1] 
    p.1 <- par[2] 
    f <- rep(NA, length(par)) 
    f[1] <- fn(p.1, ...) 
    f[2] <- fn(p.2, ...) 
    iter <- 1 
    pchg <- abs(p.2 - p.1) 
    fval <- f[2] 
    if (trace) cat("par: ", par, "fval: ", f, "\n") 
    while (pchg >= tol & abs(fval) > tol & iter <= itmax) { 
        p.new <- p.2 - (p.2 - p.1) * f[2] / (f[2] - f[1]) 
        pchg <- abs(p.new - p.2) 
        fval <- fn(p.new, ...) 
        p.1 <- p.2 
        p.2 <- p.new 
        f[1] <- f[2] 
        f[2] <- fval 
        iter <- iter + 1 
        if (trace) cat("par: ", p.new, "fval: ", fval, "\n") 
    } 
    list(par = p.new, value = fval, iter=iter) 
} 

# Net present value.
npv <- function(irr, cashflow, times) sum(cashflow / (1 + irr)^times) 

times <- as.numeric(difftime(dates, dates[1], units="days"))/365.24 

r <- secant(par=c(0,0.1), fn=npv, cashflow=cf, times=times) 

return(r$par)
}
4

3 回答 3

1

你有一个数学问题，而不是编程问题；

接近 -100% 的回报率时，净现值接近无穷大，符号为最后现金流量，因此无法计算内部回报率。

http://en.wikipedia.org/wiki/Internal_rate_of_return

于 2013-01-15T10:52:10.577 回答
0

首先让我感谢你发表这篇文章，因为我一直在寻找这样的功能很长一段时间。您的功能仍然无法控制NaN. [fval]您应该将while循环中的表达式更改为

fval <- ifelse(is.na(fn(p.new, ...)),1,fn(p.new, ...))

希望能帮助到你。

于 2015-02-18T18:25:18.247 回答
0

包中有一个xirr功能tvm。

library(tvm)

xirr_value <- xirr(cashflow, dates)

希望这可以帮助。

于 2019-10-05T11:20:43.167 回答