Exploratory Data Analysis

Time Series Plots, Autocorrelation, and Stationarity Tests

Published

December 3, 2023

Modified

September 15, 2026


Purpose

This section examines each series before any model is fitted. For each one it shows four things: the series over time, its autocorrelation function (ACF), its partial autocorrelation function (PACF), and the ACF of its first difference. It then applies the augmented Dickey-Fuller test to the level and to the first difference.

The question these tools answer is whether a series is stationary. A stationary series has a constant mean and variance and an autocorrelation that depends only on the lag, so its past behaviour is informative about its future. A series with a unit root does not revert to a mean and must be differenced before it can be modelled with ARMA methods. The answer determines the order of differencing used in the sections that follow.

The series fall into two groups by their relation to the study’s question. The three spreads are Treasury yields net of the Federal Funds Rate, so their stationarity is a statement about the policy rate: a spread that reverts to a mean is a market rate that stays a fixed distance from the Federal Reserve’s, and one with a unit root is not. The remaining series are the targets the policy rate is expected to influence, and the order of differencing found here is the order used when the policy rate is added to them in the sections that follow.

Twelve series are examined: the three short-term spreads, the three long-maturity yields, and six economic series. The 5-year spread is used only in the ARIMAX section and consumer sentiment only in the deep learning section, so they are examined there. Each series is read from its own source file, described on the data source page, and analysed over its full available sample.

Stationarity

The augmented Dickey-Fuller test (Dickey and Fuller 1979) tests the null hypothesis that a series has a unit root against the alternative that it is stationary around a linear trend. The implementation is adf.test from the tseries package (Trapletti and Hornik 2026), which fits a regression with a constant and a linear trend and sets the lag order to the integer part of \((n-1)^{1/3}\). Its p-values are interpolated from the tables in Banerjee et al. (1993) and are reported only within the range 0.01 to 0.99, so a value shown as “< 0.01” is at the lower bound. A p-value below 0.05 rejects the unit root at the 5% level. The table reports the test statistic and p-value for the level of each series and for its first difference.

Series Frequency Sample n ADF Level p ADF First Difference p Verdict
T3MFF Daily Jan 1982 to Oct 2023 10,442 -8.28 < 0.01 -30.71 < 0.01 I(0)
T6MFF Daily Jan 1982 to Oct 2023 10,442 -7.60 < 0.01 -30.45 < 0.01 I(0)
T1YFF Daily Jan 1962 to Oct 2023 15,423 -8.96 < 0.01 -26.60 < 0.01 I(0)
DGS10 Daily Jan 1962 to Oct 2023 15,425 -2.36 0.426 -22.85 < 0.01 I(1)
DGS20 Daily Oct 1993 to Oct 2023 7,507 -1.95 0.599 -17.60 < 0.01 I(1)
DGS30 Daily Feb 1977 to Oct 2023 11,655 -3.20 0.088 -20.60 < 0.01 I(1)
MORTGAGE30US Weekly Apr 1971 to Sep 2023 2,740 -2.18 0.501 -14.02 < 0.01 I(1)
FEDFUNDS Monthly Jul 1954 to Aug 2023 830 -3.18 0.091 -8.41 < 0.01 I(1)
CPIAUCSL Monthly Jan 1947 to Aug 2023 920 -0.82 0.960 -5.43 < 0.01 I(1)
UNRATE Monthly Jan 1948 to Aug 2023 908 -3.92 0.013 -9.72 < 0.01 I(0)
GDPC1 Quarterly 1947 Q1 to 2023 Q2 306 -1.05 0.929 -7.62 < 0.01 I(1)
DTWEXBGS Daily Jan 2006 to Nov 2023 4,488 -2.71 0.276 -15.41 < 0.01 I(1)

The results divide the series into two groups. The three spreads reject the unit root in levels, with test statistics between −7.6 and −9.0. They are stationary and can be modelled without differencing. This is expected: a spread between two rates that move together is bounded. The ARMA section models quarterly averages of the 3-month and 6-month spreads.

Eight of the remaining nine series do not reject the unit root in levels at the 5% level, and every series rejects it after a single difference, with statistics between −5.4 and −22.8 among those eight. They are integrated of order one and enter the models in first differences. Three results deserve comment. The 30-year yield and the Federal Funds Rate have p-values of 0.088 and 0.091, so they reject at the 10% level but not at 5%; both are treated as I(1). The unemployment rate rejects at the 5% level, with p = 0.013, on the full sample from 1948. Unemployment is bounded and reverts toward a long-run level, so the result is plausible, but the series is highly persistent, with an autocorrelation above 0.5 out to 16 months, and it is differenced in the sections that use it. CPI rejects the unit root after one difference, but the differenced series is itself strongly autocorrelated, with a lag-one coefficient of 0.59, which is why the ARMA section fits the logarithm of the index with a second difference.

Every level series apart from the spreads has a first-order autocorrelation above 0.96. Their ACF plots below decay so slowly that the bars remain near one across the entire plotted range, which is the visual signature of a unit root. The ACF of the first difference, in the lower-right panel of each figure, is the plot that carries information about the model order.

Treasury Series

The 3-month spread runs daily from January 1982 to October 2023, 10,442 observations. It ranges from −10.32 points on 30 December 1986 to 1.98 points on 24 October 1984, and its mean is −0.21. The low is a single-day year-end spike in the effective Federal Funds Rate, discussed on the visualization page. The ADF statistic in levels is −8.28, which rejects the unit root. The ACF starts at 0.83 and falls below 0.5 within 44 trading days, about two months, so shocks to the spread fade within a quarter. The PACF has a large first spike and small values after it. The first difference has a lag-one autocorrelation of −0.20, the negative value that results from differencing a stationary series.

The 6-month spread covers the same sample. Its range is −10.24 to 2.37 points, with both extremes on the same days as the 3-month spread, and its mean is −0.03. The ADF statistic is −7.60 and the unit root is rejected. The ACF starts at 0.83 and falls below 0.5 within 35 trading days. The lag-one autocorrelation of the first difference is −0.20. The 6-month spread behaves as the 3-month spread does, with slightly faster decay.

The 1-year spread has the longest sample of the three, daily from January 1962 to October 2023, 15,423 observations. It ranges from −10.11 points, again on 30 December 1986, to 3.61 points on 8 December 1965, and its mean is 0.04. The ADF statistic in levels is −8.96 and the unit root is rejected. Its autocorrelation decays most slowly of the three: the ACF starts at 0.89 and does not fall below 0.5 until lag 64, about three months. The longer horizon leaves more room for expected rate changes to accumulate, so the spread moves in longer swings. The first difference has a lag-one autocorrelation of −0.27.

The 10-year yield runs daily from January 1962 to October 2023, 15,425 observations. It peaked at 15.84% on 30 September 1981 and reached its low of 0.52% on 4 August 2020. The ADF statistic in levels is −2.36 with p = 0.43, so the unit root is not rejected. The ACF at lag one is 0.9997 and does not fall below 0.5 until lag 1,887, about seven and a half years. After one difference the statistic is −22.85 and the unit root is rejected. The differenced series has a lag-one autocorrelation of 0.06, close to white noise, which is what a random walk implies for daily yield changes.

The 20-year series has a gap. The Treasury discontinued the 20-year constant-maturity yield at the end of 1986 and reinstated it on 1 October 1993, so the level plot is blank for those years. Joining the two segments end to end would produce a false jump and a compressed time axis, so the autocorrelation plots and the test use only the continuous segment from October 1993, 7,507 observations. Over that period the yield ranges from 8.30% on 4 November 1994 to 0.87% on 9 March 2020. The ADF statistic is −1.95 with p = 0.60 in levels and −17.60 after one difference. The lag-one autocorrelation of the first difference is −0.01.

The 30-year yield runs from February 1977 to October 2023, 11,655 observations, from a peak of 15.21% on 26 October 1981 to a low of 0.99% on 9 March 2020. The ADF statistic in levels is −3.20 with p = 0.088, which rejects at 10% but not at 5%. Given a lag-one autocorrelation of 0.9997 and an ACF that does not fall below 0.5 until lag 1,685, the series is treated as I(1). After one difference the statistic is −20.60 and the lag-one autocorrelation is 0.03.

Economic Series

The 30-year mortgage rate is weekly, from April 1971 to September 2023, 2,740 observations. It peaked at 18.63% on 9 October 1981 and reached its low of 2.65% on 7 January 2021. The ADF statistic in levels is −2.18 with p = 0.50. The ACF does not fall below 0.5 until lag 364, seven years. After one difference the statistic is −14.02. Unlike the daily yields, the weekly changes are positively autocorrelated, with a lag-one coefficient of 0.27: a move in the mortgage rate in one week tends to be followed by a move in the same direction the next, because the rate is a weekly survey average that adjusts to the bond market with a lag.

The effective Federal Funds Rate is monthly, from July 1954 to August 2023, 830 observations. It ranges from 0.05% in April 2020 to 19.10% in June 1981. The ADF statistic in levels is −3.18 with p = 0.091, which rejects at 10% but not at 5%, and the series is treated as I(1). The ACF at lag one is 0.99 and does not fall below 0.5 until lag 72, six years. After one difference the statistic is −8.41. The monthly changes have a lag-one autocorrelation of 0.39, the highest among the rate series, because the Fed moves the rate in sequences of steps in the same direction.

The Consumer Price Index is monthly, from January 1947 to August 2023, 920 observations, rising from 21.48 to 306.27. The series is seasonally adjusted at source. The ADF statistic in levels is −0.82 with p = 0.96, the weakest rejection in the table, as expected for a series that has risen in almost every month of the sample. After one difference the statistic is −5.43 and the unit root is rejected, but the differenced series has a lag-one autocorrelation of 0.59 and an ACF that decays slowly, since the monthly change in the index is inflation, and inflation is itself persistent. The ARMA section uses a second difference for this reason.

The unemployment rate is monthly, from January 1948 to August 2023, 908 observations, seasonally adjusted at source. It ranges from 2.5% in May 1953 to 14.7% in April 2020. The ADF statistic in levels is −3.92 with p = 0.013, so the unit root is rejected at the 5% level on the full sample. The ACF at lag one is 0.97 and falls below 0.5 at lag 17, faster than any other level series, which is consistent with a bounded series that returns toward a long-run rate. After one difference the statistic is −9.72 and the lag-one autocorrelation is 0.04.

Real GDP is quarterly, from the first quarter of 1947 to the second quarter of 2023, 306 observations, seasonally adjusted at source. It ranges from $2,172 billion in the third quarter of 1947 to $22,225 billion in the second quarter of 2023, in chained 2017 $. The ADF statistic in levels is −1.05 with p = 0.93. After one difference the statistic is −7.62 and the unit root is rejected. The lag-one autocorrelation of the differenced series is −0.12, so quarterly changes in the level are close to uncorrelated once the trend is removed. The ARMA section models the logarithm of the series.

The nominal broad dollar index is daily, from January 2006 to November 2023, 4,488 observations. It ranges from 85.47 on 26 July 2011 to 128.32 on 26 September 2022. The ADF statistic in levels is −2.71 with p = 0.28, so the unit root is not rejected. The ACF does not fall below 0.5 until lag 801, about three years. After one difference the statistic is −15.41 and the lag-one autocorrelation is 0.04. The financial time series section models the volatility of these daily changes.

References

Banerjee, Anindya, Juan J. Dolado, John W. Galbraith, and David F. Hendry. 1993. Co-Integration, Error Correction, and the Econometric Analysis of Non-Stationary Data. Oxford University Press.
Dickey, David A., and Wayne A. Fuller. 1979. “Distribution of the Estimators for Autoregressive Time Series with a Unit Root.” Journal of the American Statistical Association 74 (366): 427–31. https://doi.org/10.1080/01621459.1979.10482531.
Trapletti, Adrian, and Kurt Hornik. 2026. Tseries: Time Series Analysis and Computational Finance. https://doi.org/10.32614/CRAN.package.tseries.