基于信号的特定日期股票投资组合的同等加权重新配置

时间:2014-12-09 09:12:17

标签: r quantmod stock portfolio algorithmic-trading

我想在特定日期重新分配策略组合:

require(PerformanceAnalytics)
require(TTR)
require(quantmod)

获取资产价格并获取每日离散收益

tickers = c("ABI.BR","AI.PA","AIR.PA","ALV.DE","ASML.AS")
getSymbols(tickers, from="2012-01-01", to="2013-12-01")

close.prices = do.call(merge, lapply(tickers, function(x) Cl(get(x))))
colnames(close.prices) = c("Anheuser-Busch InBev",
                            "L'Air Liquide","AIRBUS GROUP","Allianz","ASML HLDG")

assets.ret = ROC(close.prices,type="discrete")[-1]

现在我通过将RSI函数应用于每个资产来获取RSI信号

rsi.fct = function(x) RSI(x, n=20, maType = SMA) 
rsi     = xts(apply(close.prices, 2, rsi.fct), 
              order.by=index(rsi.fct(close.prices[,1]) ) )

> tail(rsi)
           Anheuser-Busch InBev L'Air Liquide AIRBUS GROUP  Allianz ASML HLDG
2013-11-22             51.15171      49.36494     60.25836 61.07143  46.84159
2013-11-25             54.95495      50.82237     63.54717 61.07143  49.63168
2013-11-26             49.65470      52.55102     58.29563 58.18182  48.59023
2013-11-27             54.60575      61.81980     57.94677 62.05674  52.11640
2013-11-28             46.52778      60.76994     57.85061 63.35616  45.70000
2013-11-29             50.99905      61.90476     56.09756 65.49296  48.82479

策略如下:当RSI <&lt;时,我购买资产。 30当RSI> = 30

时不要购买
ret.mat.rsi = lag(ifelse (rsi < 30, 1, 0))*assets.ret

现在这是我遇到问题的部分。 ret.mat.rsi的返回值是每日返回值。 假设我想在每月的第一天查看rsi矩阵,例如

> rsi[110]
           Anheuser-Busch InBev L'Air Liquide AIRBUS GROUP  Allianz ASML HLDG
2012-06-01             39.66126       31.1599     30.39443 17.17647  43.85172

由于RSI低于30,我想购买同样加权到我投资组合中的前4个资产 并在本月剩余时间内保持头寸不变(无论RSI信号如何),直到下个月的第一天:

> rsi[131]
           Anheuser-Busch InBev L'Air Liquide AIRBUS GROUP  Allianz ASML HLDG
2012-07-02             84.69529      73.87205     66.25561 74.52642  71.65021

我选择不购买任何资产。

现在的问题是如何优雅地编写投资组合的自动重新分配 在特定日期,即每月初(也可能是每周或每三周)。投资组合回报应仅包括在重新分配日期满足指标条件(此处为RSI <30)的那些资产。

1 个答案:

答案 0 :(得分:3)

我将如何编写您的示例:

require(quantmod)
tickers <- c("ABI.BR","AI.PA","AIR.PA","ALV.DE","ASML.AS")
myEnv <- new.env()
getSymbols(tickers, from="2012-01-01", to="2013-12-01", env=myEnv)

close.prices <- do.call(merge, eapply(myEnv, Cl))
close.prices <- close.prices[,pmatch(tickers,colnames(close.prices))]
colnames(close.prices) <- c("Anheuser-Busch InBev",
  "L'Air Liquide","AIRBUS GROUP","Allianz","ASML HLDG")

assets.ret <- ROC(close.prices,type="discrete")[-1]

rsi.fct <- function(x) RSI(x, n=20, maType = SMA) 
rsi <- xts(apply(close.prices, 2, rsi.fct), index(close.prices))

现在,要回答您的问题,请使用GSee's startpoints function获取每个月的第一个RSI值。 startpoints允许您选择任意数周的周,月,季等作为重新平衡期。

startpoints <- function (x, on = "months", k = 1) {
  head(endpoints(x, on, k) + 1, -1)
}
# get the signal at the beginning of each month
rsi.signal <- lag(ifelse(rsi < 30, 1, 0))[startpoints(rsi),]
# rsi.signal is monthly; we need a daily series where each day has the
# value from the first day of the month, so we merge with an empty xts
# object that has the daily index and use na.locf to fill the gaps
rsi.signal <- merge(rsi.signal, xts(,index(rsi)), fill=na.locf)
# now calculate returns
rsi.ret <- rsi.signal * assets.ret