用蒙特卡罗模拟计算净现值

Question

我想知道是否有人能提供我的链接或想法，我如何可以计算随机净现值后，蒙特卡罗模拟，我如何能够计算概率的npv>0？我们首先用所有的假设计算确定性净现值，然后取一些重要的参数，在这里我可以分配不确定性，然后给它们分配均匀分布(runif)。但是，正向净现值的概率似乎是0/1，介于两者之间，我计算正向净现值的概率有什么问题吗?或者我是如何计算npv_veci的。

..。

a<- runif(100,10,20)
b<- runif(100,20,30)
npv_vec <- rep(NA,ndraw)
profit_vec <- rep(NA,ndraw)


for(i in 1:ndraw) {
  npv_vec[i] <- NPV_fun(a_vec[i],b_vec[i])
  profit_vec[i] <- ifelse(npv_vec[i]>0,1,0)
}


# calculate the probability of positive npv
pb_profit <- mean(profit_vec)
pb_profit

..。

Answer 1

在一枚硬币上，硬币的正面或反面都会出现。这并不意味着头部的概率要么是0，要么是1。要估计这个概率，你必须对硬币翻转进行多次试验，并确定掷硬币的比例。

类似地，NPV>0的概率为0或1，两者之间没有任何可能性。与硬币翻转一样，您可以根据多次试验确定概率，并计算有NPV>0的试验的比例。