基于动态无风险率的Tidyquant滚动夏普比率实现问题

我明白你现在的痛点：想用tq_mutate配合修改后的PortfolioAnalytics包里的rollapply，让无风险率的滚动窗口（比如120个月）和投资组合收益的窗口完全同步，但SharpeRatio.annualized默认的Rf参数只能传单一数值或简单公式，没法直接传入滚动的时间序列。下面是具体的解决方案：

核心思路 #

问题的关键在于，要让rollapply的每个滚动窗口同时包含投资组合收益和对应的无风险率数据，然后在自定义函数里对每个窗口内的两组数据进行匹配计算，而不是分开传递收益和无风险率。

具体实现步骤 #

假设你的数据框df包含三列：date（日期列）、portfolio_return（投资组合月度收益）、rf_monthly（月度无风险率），且数据已经按日期排序。

1. 加载必要包

library(tidyquant)
library(PortfolioAnalytics) # 确保是你修改后的版本

2. 使用tq_mutate+rollapply实现同步滚动计算

df_rolling_sharpe <- df %>%
  tq_mutate(
    x = ., # 传入整个数据框，而非单独的收益列
    FUN = rollapply,
    width = 120, # 设置滚动窗口大小（120个月）
    align = "right", # 窗口对齐方式，默认右对齐（取最近的120期）
    # 自定义窗口处理函数，确保收益和无风险率窗口完全同步
    FUN = function(window_df) {
      # 提取当前窗口内的收益和无风险率序列
      r_port <- window_df$portfolio_return
      r_rf <- window_df$rf_monthly
      
      # 调用你修改后的SharpeRatio.annualized函数
      # 这里假设你已经修改函数，使其支持传入与R等长的Rf向量
      SharpeRatio.annualized(R = r_port, Rf = r_rf)
    },
    by.column = FALSE, # 必须设为FALSE，因为我们传入的是整个数据框
    col_rename = "rolling_sharpe_120m" # 给计算结果列命名
  )

关键注意点 #

• 函数修改验证：确保你修改后的SharpeRatio.annualized能够接受与R长度一致的Rf向量。原函数可能只支持单一数值，你需要调整逻辑，先计算每个时期的超额收益（R - Rf），再进行年化夏普比率的计算。
• 窗口数据完整性：如果是月度数据，120个月的窗口意味着前119行的计算结果会是NA（因为窗口不足），这是正常现象。
• 数据排序：必须确保数据按日期升序排列，否则滚动窗口的时间范围会混乱。

内容的提问来源于stack exchange，提问作者Drew Thornton