Backtesting Part 4: random strategies
[This article was first published on Modern Toolmaking, and kindly contributed to R-bloggers]. (You can report issue about the content on this page here)
Want to share your content on R-bloggers? click here if you have a blog, or here if you don't.
Note: This post is NOT financial advice! This is just a fun way to explore some of the capabilities R has for importing and manipulating data. Want to share your content on R-bloggers? click here if you have a blog, or here if you don't.
In part 2, we found that our 200-day high, hold 100 days strategy yielded average annual returns on the S&P 500 index of 7% in a backtest going back to 1950. That’s pretty good, but it’s entirely possible that we got lucky and this was due to chance. One way to test this is to compare our strategy to a strategy that randomly chooses long or short each day. I tested just such a strategy on the S&P 500 index, again going back to 1950. I also adjusted the returns series for splits and dividends, and assumed trade costs were 0.5%. These are the exact same conditions I tested the “high-and-hold” strategy under.
Here is a histogram of average annual returns for these random strategies:
As you can see these strategies perform pretty dismally: the 0.5% trading cost quickly eats up all the available capital.
It would probably be a good idea to model a set dollar amount, and have trades costs a fixed amount, something like $7. It might also be a good idea to make the random strategies auto-correlated, which would probably better reflect actual investor behavior. The code for this test is below the fold, if anyone wants to modify it and make these improvements.
This file contains bidirectional Unicode text that may be interpreted or compiled differently than what appears below. To review, open the file in an editor that reveals hidden Unicode characters.
Learn more about bidirectional Unicode characters
#Setup | |
rm(list = ls(all = TRUE)) | |
set.seed(1) | |
nsims <- 1000 | |
library(quantmod) | |
library(PerformanceAnalytics) | |
#Load Data | |
getSymbols('^GSPC',from='1900-01-01') | |
spyReturns <- dailyReturn(GSPC$GSPC.Adjusted, type = "arithmetic") | |
#Make 1000 random position vectors | |
randomPositions <- sample(c(0,1),nsims*length(spyReturns),TRUE) | |
randomPositions <- matrix(randomPositions,ncol=nsims) | |
#Determine trades and calculate trade costs | |
trades <- apply(randomPositions,2,Lag,1) != randomPositions | |
trades[1,] <- trades[1,]==1 | |
trades <- ifelse(trades,0.001,0) #Assume the commision is 0.5% | |
#Determine dailty returns of random strats | |
randomStrats <- randomPositions*matrix(rep(spyReturns,nsims),ncol=nsims) | |
randomStrats <- randomStrats-trades | |
randomStrats.return <- apply(randomStrats,2,Return.annualized,scale=252) | |
hist(randomStrats.return, freq=FALSE) | |
#Plot some of the strategies | |
library(PerformanceAnalytics) | |
charts.PerformanceSummary(na.omit(cbind(spyReturns,randomStrats)[,1:15]),colorset=redfocus,) | |
#Same plot, removes buy&hold | |
#charts.PerformanceSummary(na.omit(cbind(spyReturns,randomStrats)[,2:15]),colorset=redfocus,) |
To leave a comment for the author, please follow the link and comment on their blog: Modern Toolmaking.
R-bloggers.com offers daily e-mail updates about R news and tutorials about learning R and many other topics. Click here if you're looking to post or find an R/data-science job.
Want to share your content on R-bloggers? click here if you have a blog, or here if you don't.