使用循环为R中的多个债券定价

library(quantmod)

bonds <- data.frame(
  code = c("BOND1", "BOND2", "BOND3"),
  face_value = c(1000, 2000, 1500),
  maturity = as.Date(c("2022-12-31", "2023-06-30", "2024-05-31")),
  coupon_frequency = c(2, 4, 1),
  coupon_dates = list(as.Date(c("2022-06-30", "2022-12-31")), 
                      as.Date(c("2023-03-31", "2023-06-30", "2023-09-30", "2023-12-31")), 
                      as.Date(c("2024-03-31", "2024-05-31"))),
  yield = c(0.05, 0.06, 0.04)
)

calculate_bond_price <- function(face_value, maturity, coupon_frequency, coupon_dates, yield) {
  cashflows <- c(rep(face_value * coupon_frequency, length(coupon_dates) - 1), face_value + face_value * coupon_frequency)
  pv <- sum(cashflows / ((1 + yield / coupon_frequency) ^ (seq_along(coupon_dates) - 1)))
  return(pv)
}

bond_prices <- list()
for (i in seq_len(nrow(bonds))) {
  bond_price <- calculate_bond_price(bonds$face_value[i], bonds$maturity[i], bonds$coupon_frequency[i], bonds$coupon_dates[[i]], bonds$yield[i])
  bond_prices[[i]] <- bond_price
}

for (i in seq_len(nrow(bonds))) {
  cat("债券", bonds$code[i], "的定价为", bond_prices[[i]], "\n")
}