Numpy净现值计算

Question

我正在用numpy和我自己的代码计算NPV，结果是不同的。我一定是在哪里犯了个错误。有指针吗？

// Solution 1
r = .06
flows = {0:1200, 3:-450, 6:-450, 15:-450}
print  sum([C/(1+r)**i for i,C in flows.iteritems()])
// => 317

// Solution using numpy's npv function
flows = zeros(16)
flows[0] = 1200
flows[3] = -450
flows[6] = -450
flows[15]= -450
print np.npv(r, flows)
// => 299

Answer 1

看起来(尽管它说的是在医生里)，np.npv开始用t= 1，而不是t=0求和：

In [56]: r = 0.06

In [57]: R = r+1

In [58]: (1200/R**0 - 450/R**3 - 450/R**6 - 450/R**15)
Out[58]: 317.16980210661666

In [59]: (1200/R**0 - 450/R**3 - 450/R**6 - 450/R**15)/R
Out[59]: 299.21679444020435

In [64]: np.npv(r, flows)*(1+r)
Out[64]: 317.16980210661683

的确，np.npv是这样定义

def npv(rate, values):
    values = np.asarray(values)
    return (values / (1+rate)**np.arange(1,len(values)+1)).sum(axis=0)

Answer 2

这是在numpy 1.8中修复的

修正到financial.npv npv函数有一个错误。与文档声明的相反，它从索引1到M，而不是从0到M1。修复更改了返回的值。mirr函数称为npv函数，但它解决了这个问题，因此这也是固定的，并且mirr函数的返回值保持不变。

NumPy 1.8.0发行说明

如果我在np.npv中使用numpy 1.8，我会得到：

import numpy as np

r = .06

flows = np.zeros(16)
flows[0] = 1200
flows[3] = -450
flows[6] = -450
flows[15]= -450

result = np.npv(r, flows)

结果：

>>> result
317.16980210661666