Создать функцию в R, чтобы найти ставку купона для облигации

Question

Пока что я написал следующую функцию:

c_rate <- function (bond_value, par, ttm, y) {

   t <- seq(1, ttm, 1)

  pv_factor <- 1 / (1 + y)^t

 cr <- (bond_value - par / (1+y)^t) / (par*sum(pv_factor))

 cr
 }

однако это дает несколько результатов.

Как я могу обновить функцию, чтобы получить только один конечный индекс только?

Answer 1

Я думаю, что вы хотите поднять (1+y)^t до степени ttm вместо:

cr <- (bond_value - par / (1+y)^ttm) / (par*sum(pv_factor))