RDP 9601: Why Does the Australian Dollar Move so Closely with the Terms of Trade? Appendix C: Monte Carlo Simulations
May 1996
 Download the Paper 204KB
This appendix outlines the Monte Carlo simulations used to generate results reported in Tables 5a, 5b and 6 in the text. Simulations are conducted assuming four different data generating processes (DGPs), which we examine in turn.
C.1 The Real Exchange Rate
To begin, we test the null hypotheses that our models of the real exchange rate have no explanatory power, that is, that = 0 in each of the equations:
To test these null hypotheses, we assume a data generating process (DGP_{1}) with the log terms of trade, tot_{t}, following the preferred specification identified in Section 2.2, and embodying the hypothesis that the real exchange rate is unforecastable:
Let and the estimates be . Both equations are estimated by OLS, with equation (C2) estimated over the full terms of trade sample, 1969:Q3 to 1994:Q2, and (C3) estimated over the postfloat period, 1984:Q1 to 1994:Q2. The offdiagonal elements of V are derived from the correlation coefficient, _{1,2}, of the estimated errors, _{1,t} and _{2,t} over the postfloat period 1984:Q1 to 1994:Q2. The estimates are:
and
The Monte Carlo distributions are generated by running 5,000 trials with each trial (i = 1, …, 5,000) proceeding as follows.
 Draw a vector sequence of observations from a bivariate normal distribution with mean 0, covariance matrix , n = 58 and T = 44. (n is the length of the prefloat terms of trade sample; T is the postfloat period, 1984:Q1 to 1994:Q4.)

Generate sequences of observations ,
according to ,
using (tot_{5},tot_{4},tot_{3},tot_{2},tot_{1}) to start the autoregression.  For time t, (t = 71 corresponds to 1987:Q1; t = 100 to 1994:Q2) use the sequence of synthetic terms of trade data, to estimate the terms of trade model (C2) and thereby generate a vector of parameter estimates , and hence derive , k = 0,…, 8, t = 71,…, 100.
 For time t, , use the sequences of synthetic data, and to estimate a preferred real exchange rate model. For model A, this is equation (6) in the text, . For model B, it is one of the three specifications in equation (8), depending on the time, t. Use this preferred real exchange rate model, together with estimates of , derived above, to generate , x = A,B, k = 1,…, 8, t = 71,…, 100.
 Estimate the regressions , t = 71,…,100, by OLS for all values k = 1,…, 8, and for the two models and generate the OLS tstatistics, , x = A, β.^{[10]}
These 5,000 observations of t_{OLS} and t_{OLS} form the Monte Carlo distributions under the null hypothesis that the models of the real exchange rate have no explanatory power. Based on these Monte Carlo distributions, the results of the hypothesis tests, H_{01}:β = 0, for k = 1,…, 8, are shown in column (x) in Table 5a for model A and in Table 5b for model B.
We turn now to the second data generating process (DGP_{2}) which assumes that the real exchange rate model, estimated over the postfloat period 1984:Q1 to 1994:Q2, is the true model. For this datagenerating process, we again assume that the terms of trade follow model (C2) with parameter values (C4). The real exchange rate model (C3) is, however, replaced by:
with parameter values, estimated by OLS over the postfloat period 1984:Q1 to 1994:Q2, given by:
The Monte Carlo distributions under this data generating process are again generated by running 5,000 trials with each trial (i = 1, …, 5,000) proceeding as before, but with the following modifications. In step 2, the sequences of observations , are now derived from and .
In step 5, after running the regressions , x = A,B, k = 1,…,8 we now collect the coefficients, .
These 5,000 observations of for the two models, x = A,B and for each value of k form the Monte Carlo distributions under the assumption that the real exchange rate is described by model (C6) with parameter values (C7). They are used to test the null hypotheses, H_{02}:β = 1, for k = 1,…, 8, with the results shown in column (xi) of Tables 5a and 5b.
C.2 Excess Returns
The third data generating process, DGP_{3}, combines the terms of trade and real exchange rate models defined by equations (C2) and (C6) with parameter values (C4) and (C7) – used for DGP_{2} – with models for the Australianforeign quarterly inflation differential, , and the oneyear (fourquarter) Australianforeign nominal interest differential, . is defined by = − , where = − ) and is the tradeweighted foreign quarterly inflation rate, =, and is defined analogously.
Simple bivariate autoregressive timeseries models are fitted for and by OLS using data from 1987:Q1 to 1994:Q4 and allowing up to six lags for each variable. Eliminating insignificant lags leads to these models:
where ε3_{t} and ε4_{t} are uncorrelated i.i.d. N(0,1) random variables.
Assuming this data generating process, DGP_{3}, Monte Carlo distributions are again generated by running 5,000 trials with each trial i proceeding as for DGP_{2}, but with the following modifications. As well as generating sequences of observations , , we now also generate sequences of observations and using (C8) and (C9).
As for DGP_{2}, for time t, , we use the sequences of synthetic data, and to estimate real exchange rate model B (which involves estimating one of the three specifications in equation (8), depending on the time, t). With this real exchange rate model, together with estimates of E_{t}(), we generate , t = 71,…, 100.
At time t, the expected and actual excess returns to holding a oneyear Australian dollar bond rather than a tradeweighted basket of foreign bonds, are and which are equations (12) and (10) in the text. By construction, the expected real interest differential using backwardlooking inflationary expectations is , while the actual real interest differential is . Thus, for trial i, the expected and actual excess returns are (approximately) and .
For trial i, we use the values of E_{t}() and , to run the regression = α^{i} + β^{i} E_{t}() + ε_{t}. The 5,000 values of β^{i} form the Monte Carlo distribution with which the hypotheses, H_{01}:β = 0 and H_{02}:β = 1 are tested. The results are shown in the first row of Table 6.
The fourth data generating process, DGP_{4}, combines the terms of trade and real exchange rate models defined by equations (C2) and (C6) with parameter values (C4) and (C7), with models for the Australianforeign quarterly inflation differential, , and the twoyear (eightquarter) Australianforeign nominal interest differential, .
As before, simple bivariate autoregressive timeseries models are fitted for and by OLS using data from 1987:Q1 to 1994:Q4 and allowing up to six lags for each variable. Eliminating insignificant lags leads to these models:
where, again, ε_{3t} and ε_{4t} are uncorrelated i.i.d. N(0,1) random variables.
The Monte Carlo distributions are derived as for DGP_{3}, with the following modifications. The expected and actual excess returns to holding a twoyear Australian dollar bond are and which are equations (14) and (13) in the text. The expected real interest differential using backwardlooking inflationary expectations is , while the actual real interest differential is . Thus, for trial i, the expected and actual excess returns are (approximately) and .
For trial i, we use the values of E_{t} () and , to run the regression . Again, the 5,000 values of β^{i} form the Monte Carlo distribution with which the hypotheses, H_{01}:β = 0 and H_{02}:β = 1 are tested. The results are shown in the second row of Table 6.
Footnote
As explained earlier, the actual exchange rate data ends in 1994:Q4, and hence the sample size is 30 for k = 1,2, but progressively less for forecasts further ahead (for k = 8, it is only 24). This pattern is replicated for the synthetic data. [10]