STA457 Time Series Analysis Assignment
R Language Assignment Question
Statistical properties of (moving-average) rule returns
1. Introduction
Technical indicator is widely used to generate trading signals by practitioners to make trading decisions. The usual rule is to trade with the trend. In this case, the trader initiates a position early in the trend and maintains that position as long as the trend continues.
In this assignment, you are asked to study the statistical properties of returns for applying the oldest and most widely used method in technical indicatorsโmoving averages.1
The structure of this paper is given as follows. Section 2 defines the trading rule (or strategy). In Section 3 and 4, we formulate the trading return based on a given trading rule and state the corresponding statistical properties, respectively. The questions for you to answer are listed in Section 5. Finally, references and appendix are given in Section 6 and 7, respectively.
2. Trading rule
Suppose that at each time ๐ก, market participants predict the direction of the trend of asset prices using a price-based forecast ๐น , where ๐น is a function of past asset prices
##
๐น =๐(๐,โฆ,๐ ,โฆ). # # #*+,-
1 The simplest rule of this family is the single moving average which says when the rate penetrates from below (above) a moving average of a given length, a buy (sell) signal is generated.
2 The above predictor is then converted to buy and sell trading signals ๐ต#: buy (+1) and sell (-1)
using, i.e.
{"{"}โ๐๐๐๐โ โ ๐ต#=โ1 โ ๐น#<0 โ๐ต๐ข๐ฆโ โ ๐ต#=+1 โ ๐น#>0
Note that the signal of a trading rule is completely defined by one of the inequalities giving a sell or buy order (if the position is not short, it is long).
For example, consider a trading rule based on the moving average of order five rule (๐ = 5). In this case, ๐ is given by
โC ๐ ๐(๐,โฆ,๐ )=๐โ BDE #*B.
# #*+,- #
5 In this case, we buy the asset (๐ต# = +1) at time ๐ก + 1when
โC ๐
๐น >0โบ๐ > BDE #*B;
and sell the asset (๐ต# = โ1) when
โC ๐
๐น <0โบ๐ < BDE #*B.
The Figure below illustrates the dynamics of the above 5-periods moving average methodโ when the rate penetrates from below (above) the moving average of order five, a buy (sell) signal is generated.
##
5
##
5
2.45
24
2.35
Price
2.25
2.2
Signal
006
006 _
004 _ IN
9 13 17 21 25 29 33 37 II 45 Day
a
5 9 13 17 21 25 29 33 37 41 45 Day
Return oscillator
i A
iโ\ 002 โ i 1
i\r,
it\ Logretum t
.402
-004
-036
โข
Settโ, โ,s1
\ 13:.C,. 1 01/
N,โ0
saxโ
i 3.11
I
9 13 17 21 25
33 37 41 45
3
:Simple moving average method 5-days moving average
For your assignment, we consider ๐น based on a moving-average technical indicator. In #
= ๐๐(๐ โ๐
where ๐
# ##*- K
Figure 3.1: The simple moving average method.
general, for a given moving-average indicator, ๐น may be expressed as (a function of log #
a. price series, b. signal time series, c return oscillator
+*M
returns):
), ๐ฟ and ๐ are defined by a given trading rule (See Appendix for more details). For this assignment, we assume ๐ฟ = 0.
๐น=๐ฟ+I๐๐ , (1)
#
KDE
K #*K
48
3. Rule returns
For the period [๐ก โ 1, ๐ก), a trader following a given technical rule establishes a position (long or short) at time ๐ก โ 1, ๐ต#*-. The returns at time t made by applying such a decision rule is called โruled returnsโ and denoted as ๐ # . Their value can be expressed as
๐ # =๐ต#*-๐# โR๐ # =โ๐# ๐๐ ๐ต#*- =โ1T, ๐ #=+๐# ๐๐ ๐ต#*-=+1
where ๐ = ๐๐(๐ โ๐ ) denote the logarithm return over this period (assume no dividend # ##*-
payout during period ๐ก).
Remark: ๐
# is unconditional and unrealized returns. By unrealized we mean that rule returns
are recorded every day even if the position is neither closed nor reversed, but simply carries on. Remark: We may define the realized returns as
W
๐ U # = I ๐ # , V , VD-
where ๐ท represents the stochastic duration of the position which will last ๐ days if {'{๐ท=๐}'}โ{'{๐ต#*- โ ๐ต#,๐ต# =๐ต#,- =โฏ=๐ต#,W*-,๐ต#*W,- โ ๐ต#,W }'}.
4
5
4. Statistical properties of rule returns
Under the assumption that ๐# follows a stationary Gaussian process, several statistical properties of rule returns can be derived:
1. Unconditional expected return:
2 ๐M ๐
๐ธ(๐ )=^ ๐ โ ๐๐๐๐(๐,๐น )โ ๐๐ฅ๐iโ k l+๐m1โ2๐ทoโ kpq, (2) #๐b##*-2๐M ๐
kk
where ๐ท(h) = โซx tโ{'2๐v*- ๐๐ฅ๐{โ๐ฅM/2}'} ๐๐ฅ, ๐ = ๐ธ(๐ ), ๐ = ๐ฃ๐๐(๐ ), ๐ = ๐ธ(๐น ), and
*y
2. Unconditional variance:
๐ฃ๐๐(๐ ) = ๐ธ(๐M) โ ๐ธ(๐ )M = ๐M + ๐M โ ๐ธ(๐ )M. ##bb#
Additionally, Kedem (1986) shows that the expected zero crossing rate for a stationary process as the expected zero-crossing rate for a discrete-time, zero-mean, stationary Gaussian sequence ๐# is given by
1 ๐๐๐ *- {'๐(1)'}, ๐
where ๐(1) denotes the autocorrelation function of {'{๐#}'} at lag one. Using the same assumption, we can show that ๐น is stationary. Using this result, we may approximate the
expected length of the holding period2 for a given trading rule as
2 Intuitively, the longer holding period, the larger the expected return on a trading rule.
๐M = ๐ฃ๐๐(๐น ). k#
k#b#k#
#
>๐ป= ๐ . (3) ๐๐๐ *- ๐k(1)
5. Questions
1. Derive the variance of the predictor ๐น given in Equation (1). #
Hint: ๐M = ๐ฃ๐๐(โ+*M ๐ ๐ ). k BDE B #*B
2. Derive the expectation of the predictor ๐น . #
Hint: ๐ = ๐ธ(โ+*M ๐ ๐ ). k BDE B #*B
3. Derive the autocorrelation function at lag one for the predictor. Hint: ๐ (1) = ๐๐๐๐(๐น , ๐น ).
6
k # #*-
4. Write a R function to calculate the expectation of the rule return for a given double MA trading rule (See Appendix) and the expected length of the holding period.
5. Use a R function to download daily, weekly S&P500 index from Oct/01/2009 to Sep/30/2018 from yahoo finance
Hint: adjusted Close and R quantmod library.
6. Write a R function to choose the optimal daily and weekly double MA trading rules (that maximize the expected rule returns) for S&P500 index.
7. Write a R function to calculate the in-sample trading statistics (cumulative return and holding time) of your choice and compare them with your theoretical results.
8. (Optional) Run and back-test your daily trading rule using six months of rolling window. Show the empirical trading statistics and show the difference between the theoretical results.
Hint: Given asset price time series and a pair of integers, ๐ and ๐ (function arguments),
your function calculates the expected rule return ๐ธ(๐ #) and the expected length of holding
periods ๐ป
Hint: Find the ๐ and ๐ pair that has the highest ๐ธ(๐ #). For simplicity, let the maximum
values of ๐ be 250 and 52 for daily and weekly data, respectively
Hint: Use the ratio of the cumulative return over the number of trading periods as the
estimate ๐ธ(๐ #).
6. Reference
- Acar, E. (1993). Economic evaluation of financial forecasting. (Unpublished Doctoral thesis, City University London.)
- Acar E. (200?), โAdvanced trading ruleโ, Second edition. (Chapter 4. Expected returns of directional forecasters).
- Kedem (1986), โSpectral analysis and discrimination by zero-crossingsโ, Proceedings of IEEE, Vol 74, No. 11, page 1477-1493.
7. Appendix
134 Advanced Trading Rules
Table 1: Return/Price signals equivalence
Table 4.1
Return/price signals equivalence
Return sell signals
dj ๔ฐ๔ฐm๔ฐj๔ฐ๔ฐm๔ฐj๔ฐ1๔ฐ 2
see generalization
dj๔ฐ1
Xm๔ฐ2
j๔ฐ0
dj ๔ฐ๔ฐm๔ฐr๔ฐ๔ฐj๔ฐ 1๔ฐ
for0๔ฐ j๔ฐ r๔ฐ1 dj ๔ฐr๔ฐm๔ฐj๔ฐ1๔ฐ
forr๔ฐ j๔ฐ m๔ฐ1 m๔ฐ2
๔ฐ X dj Xt๔ฐj ` 0Ywith:
j๔ฐ0 m๔ฐ1 m๔ฐ1 dj ๔ฐ ๔ฐ X ai and ๔ฐ X aj
i๔ฐj๔ฐ 1 j๔ฐ0
Parameter(s) |
Price sell signals |
Pt ` Xm๔ฐ1 j๔ฐ0 ajPt๔ฐj | |
m๔ฐ2 |
a j ๔ฐ m1 |
m๔ฐ2 |
aj๔ฐ m๔ฐj ๔ฐm๔ฐm ๔ฐ 1๔ฐ๔ฐa2 |
1 b a b 0 Ym ๔ฐ 2 |
aj ๔ฐa๔ฐ1๔ฐa๔ฐj |
m๔ฐ2 |
aj ๔ฐ1forj๔ฐm๔ฐ1Y aj ๔ฐ0forj6๔ฐm๔ฐ1 |
r๔ฐ1 m๔ฐ1 XbjPt๔ฐj `XajPt๔ฐj j๔ฐ0 j๔ฐ0 | |
mbr๔ฐ1 |
bj ๔ฐ1rY aj ๔ฐm1 |
Xm๔ฐ1 ajPt๔ฐj `0 j๔ฐ0 |
Rule
Simple order
Simple MA
Weighted MA
Exponential MA
Momentum
Double orders
Double MA
Generalization
Xm๔ฐ2 j๔ฐ0
djXt๔ฐj `0 dj ๔ฐ๔ฐm๔ฐj๔ฐ1๔ฐ
djXt๔ฐj `0
for the simple and weighted moving average rule.4 As can be seen from Tables 4.2 and 4.3, signals are di๔ฐerent in less than 3 per cent of all cases for simulated series. The largest deviation comes from โ ๔ฐ 0.001 and โ ๔ฐ 0.03, for m ๔ฐ 200. Even then, returns signals are identical to price signals in at least 97 per cent of all cases. This case represents an upper bound in terms of both volatility and average returns over ten years for ยฎnancial series (Taylor, 1986: Tables 3 and 4, pp. 33 and 34). On the basis of the empirical results presented in Tables 4.2 and 4.3, one can safely conclude that return signals lead to essentially the same investment strategies as price signals for values of m as large as 200.