拓端tecdat|R语言动量交易策略分析调整后的数据

原文链接：http://tecdat.cn/?p=18686

 

用于动量策略中所谓的动量（Momentum），是指某一对象所具有的一种倾向于保持其原有属性或特征的性质，也可以简单理解成一种惰性（Inertia）。股票的动量，简单地说就是涨的还会接着涨，跌的还会接着跌；过去涨得越猛，未来涨的也就越猛；过去跌得越狠，未来也会跌的越狠。

下面，本文将尝试将动量策略应用于调整后的数据（历史价格，股息和拆分数据，以预测公司的未来或获得市场见解）：

  1.  
    #*****************************************************************
  2.  
    # 加载一天结束时的历史数据
  3.  
    #*****************************************************************
  4.  
     
  5.  
    # 股票代码
  6.  
    tickers = '
  7.  
    AGG
  8.  
    DBC
  9.  
    EEM
  10.  
    EFA
  11.  
    '
  12.  
     
  13.  
     
  14.  
     
  15.  
    getSymbols.extra(tickers, src = 'yahoo', from = '1970-01-01', env = data, set.symbolnames = T, auto.assign = T)
  16.  
     
  17.  
    #*****************************************************************
  18.  
    # 构建另一个不含股利的价格调整后的回测环境
  19.  
     
  20.  
     
  21.  
    #*****************************************************************
  22.  
     
  23.  
    for(i in data$symbolnames) data.price[[i]] = adjustOHLC(data[[i]], symbol.name=i, adjust='split', use.Adjusted=F)
  24.  
     
  25.  
     
  26.  
    #*****************************************************************
  27.  
    # 调整价格
  28.  
    #*****************************************************************
  29.  
    for(i in data$symbolnames) data[[i]] = adjustOHLC(data[[i]], use.Adjusted=T)
  30.  
    prep(data, align='keep.all')

我们构建了两个环境：

  • 包含拆分和股息调整后的价格
  • 仅包含调整后的价格
  1.  
     
  2.  
     
  3.  
     
  4.  
    #*****************************************************************
  5.  
    # 排名60天变化率
  6.  
     
  7.  
    #******************************************************************
  8.  
    return = prices / mlag(prices,60) - 1
  9.  
    position.score = iif(return < 0, NA, return)
  10.  
    data$weight[period.ends0,] = ntop(position.score[period.ends1,], 1)
  11.  
     
  12.  
     
  13.  
     
  14.  
    #*****************************************************************
  15.  
    # 60/120天的排名变化率
  16.  
    #******************************************************************
  17.  
    return = prices / mlag(prices,60) - 1 + prices / mlag(prices,120) - 1
  18.  
    position.score = iif(return < 0, NA, return)
  19.  
     
  20.  
     
  21.  
    data$weight[period.ends0,] = ntop(position.score[period.ends1, , 1)
  22.  
     
  23.  
     
  24.  
    #*****************************************************************
  25.  
    # 创建报表
  26.  
    #*****************************************************************
  27.  
    plot(models)

 

 

  1.  
     
  2.  
    mom60 mom120
  3.  
    Period Jan2002 - Mar2015 Jan2002 - Mar2015
  4.  
    Cagr 18.05 15.47
  5.  
    Sharpe 0.89 0.76
  6.  
    DVR 0.85 0.7
  7.  
    R2 0.95 0.92
  8.  
    Volatility 21.26 22.21
  9.  
    MaxDD -33.49 -51.75
  10.  
    Exposure 94.36 94.36
  11.  
    Win.Percent 63.76 62.42
  12.  
    Avg.Trade 1.69 1.49
  13.  
    Profit.Factor 2.24 2
  14.  
    Num.Trades 149 149

 

  1.  
    models$mom60 weight entry.date exit.date nhold entry.price exit.price return
  2.  
    SPY 100 2013-06-28 2013-07-31 33 155.04 163.06 5.17
  3.  
    XLV 100 2013-07-31 2013-08-30 30 49.88 48.12 -3.53
  4.  
    XLV 100 2013-08-30 2013-09-30 31 48.12 49.66 3.20
  5.  
    FXI 100 2013-09-30 2013-10-31 31 35.92 36.40 1.34
  6.  
    EEM 100 2013-10-31 2013-11-29 29 41.16 41.05 -0.27
  7.  
    XLV 100 2013-11-29 2013-12-31 32 54.24 54.64 0.75
  8.  
    SPY 100 2013-12-31 2014-01-31 31 180.35 173.99 -3.53
  9.  
    XLV 100 2014-01-31 2014-02-28 28 55.16 58.59 6.22
  10.  
    IYR 100 2014-02-28 2014-03-31 31 65.72 65.81 0.14
  11.  
    IYR 100 2014-03-31 2014-04-30 30 65.81 67.81 3.04
  12.  
    EEM 100 2014-04-30 2014-05-30 30 40.42 41.62 2.97
  13.  
    EEM 100 2014-05-30 2014-06-30 31 41.62 42.62 2.40
  14.  
    IYR 100 2014-06-30 2014-07-31 31 70.41 70.33 -0.11
  15.  
    FXI 100 2014-07-31 2014-09-30 61 39.96 37.80 -5.41
  16.  
    UUP 100 2014-09-30 2014-10-31 31 22.87 23.09 0.96
  17.  
    XLV 100 2014-10-31 2014-11-28 28 67.02 69.35 3.48
  18.  
    XLV 100 2014-11-28 2014-12-31 33 69.35 68.38 -1.40
  19.  
    IYR 100 2014-12-31 2015-01-30 30 76.84 81.23 5.71
  20.  
    IYR 100 2015-01-30 2015-02-27 28 81.23 79.12 -2.60
  21.  
    FXI 100 2015-02-27 2015-03-30 31 43.76 44.74 2.24

接下来，我们考虑以下3种设置：

  • “已调整”-拆分和股息调整后的价格数据。信号和损益仅从调整后的数据时间序列中得出。
  • “实际”-仅使用拆分后的调整数据，尚未对股息进行调整。信号和损益仅从唯一的拆分调整后的数据时间序列中得出。
  • “混合”-使用“实际”时间序列生成信号，使用“调整”时间序列进行损益计算。
  1.  
    #*****************************************************************
  2.  
    # 辅助函数
  3.  
    #******************************************************************
  4.  
     
  5.  
     
  6.  
     
  7.  
    return = prices / mlag(prices,60) - 1
  8.  
    position.score = iif(return < 0, NA, return)
  9.  
     
  10.  
     
  11.  
    return = prices / mlag(prices,60) - 1 + prices / mlag(prices,120) - 1
  12.  
    position.score = iif(return < 0, NA, return)
  13.  
     
  14.  
     
  15.  
    models[[paste0('mom120',name)]] =run(data, clean.signal=F, commission = commission, trade.summary=T, silent=T)
  16.  
     
  17.  
    models
  18.  
     
  19.  
    #*****************************************************************
  20.  
    # 安装
  21.  
    #******************************************************************
  22.  
     
  23.  
    all.models = c(all.models, strategy(prices.adj, data, 'Adjusted'))
  24.  
     
  25.  
     
  26.  
    #*****************************************************************
  27.  
    # 创建报表
  28.  
    #*****************************************************************
  29.  
     
  30.  
    plot(models

 

print(plotbt.strateg(models, make.plot=F,

 

  1.  
     
  2.  
    mom60Adjusted mom60Actual mom60Hybrid
  3.  
    Period Jan2002 - Mar2015 Jan2002 - Mar2015 Jan2002 - Mar2015
  4.  
    Cagr 18.05 15.02 17.22
  5.  
    Sharpe 0.89 0.76 0.85
  6.  
    DVR 0.85 0.73 0.82
  7.  
    R2 0.95 0.96 0.96
  8.  
    Volatility 21.26 21.33 21.33
  9.  
    MaxDD -33.49 -35.99 -33.74
  10.  
    Exposure 94.36 92.44 92.44
  11.  
    Win.Percent 63.76 62.33 63.01
  12.  
    Avg.Trade 1.69 1.49 1.67
  13.  
    Profit.Factor 2.24 1.98 2.15
  14.  
    Num.Trades 149 146 146
  1.  
     
  2.  
     
  3.  
    plotbt(models, plotX = T, log = 'y'

 

print(plotstrategy)
  1.  
     
  2.  
    mom120Adjusted mom120Actual mom120Hybrid
  3.  
    Period Jan2002 - Mar2015 Jan2002 - Mar2015 Jan2002 - Mar2015
  4.  
    Cagr 15.47 13.61 15.66
  5.  
    Sharpe 0.76 0.69 0.77
  6.  
    DVR 0.7 0.64 0.73
  7.  
    R2 0.92 0.92 0.94
  8.  
    Volatility 22.21 22.08 22.07
  9.  
    MaxDD -51.75 -49.93 -47.15
  10.  
    Exposure 94.36 92.47 92.47
  11.  
    Win.Percent 62.42 60.96 61.64
  12.  
    Avg.Trade 1.49 1.37 1.54
  13.  
    Profit.Factor 2 1.9 2.02
  14.  
    Num.Trades 149 146 146

经调整后的数据表现优于实际数据和混合数据。

最受欢迎的见解

1.用机器学习识别不断变化的股市状况—隐马尔科夫模型(HMM)的应用

2.R语言GARCH-DCC模型和DCC（MVT）建模估计

3.R语言实现 Copula 算法建模依赖性案例分析报告

4.R语言COPULAS和金融时间序列数据VaR分析

5.R语言多元COPULA GARCH 模型时间序列预测

6.用R语言实现神经网络预测股票实例

7.r语言预测波动率的实现：ARCH模型与HAR-RV模型

8.R语言如何做马尔科夫转换模型markov switching model

9.matlab使用Copula仿真优化市场风险

 

posted @ 2021-01-02 21:10  拓端tecdat  阅读(144)  评论(0)  编辑  收藏  举报