Which DickeyFuller test should I apply to this time series?by pythonscript Tags: apply, dickeyfuller, series, test, time 

#1
Nov2912, 06:48 PM

P: 7

I have a time series of climate data that I'm testing for stationarity. Based on previous research, I expect the model underlying the data to have an intercept term, a positive linear time trend, and some normally distributed error term. In other words, I expect the underlying model to look something like this:
y_{t} = a_{0} + a_{1}t + β y_{t1} + u_{t} $ where u_{t} is normally distributed. Since I'm assuming the underlying model has both an intercept and a linear time trend, I tested for a unit root with equation #3 of the simple DickeyFuller test, as shown: ∇y_{t} = α_{0}+α_{1}t+δ y_{t1}+u_{t} This test returns a critical value that would lead me to reject the null hypothesis and conclude that the underlying model is nonstationary. However, I question if I'm applying this correctly, since even though the underlying model is assumed to have an intercept and a time trend, this does not imply that the first difference ∇y_{t} will as well. Quite the opposite, in fact, if my math is correct. Calculating the first difference based on the equation of the assumed underlying model gives: ∇y_{t} = y_{t}  y_{t1} = [a_{0} + a_{1} + β y_{t1} + u_{t}]  [a_{0} + a_{1}(t1) + β y_{t2} + u_{t1}] ∇y_{t} = [a_{0}  a_{0}] + [a_{1}t  a_{1}(t1)] + β[y_{t1}  y_{t2}] + [u_{t}  u_{t1}] ∇y_{t} = a_{1} + β * ∇y_{t1} + u_{t}  u_{t1}$ Therefore, the first difference ∇y_{t} appears to only have an intercept, not a time trend. Because the underlying model has an intercept and a time trend, should I use the DickeyFuller test that includes an intercept and time trend when it tests for a unit root, or should I use the DickeyFuller test that only includes an intercept because the first difference of the original time series only has an intercept? 



#2
Nov2912, 11:14 PM

Sci Advisor
P: 3,173

You are correct that the difference operation a linear trend term (k)(t) produces a constant, not another linear trend term.
It's interesting to try to read that Wikipedia article you linked. It's about some hypothesis tests, but it never manages to state the null hypotheses and the test statistics. Apparently you understand the test statistics  which is more than I know. However, I did find this page: http://stats.stackexchange.com/quest...eyfullertest where the reply mentions that using the test with term of the form (k)(t) implies we are investigating a model with a quadratic trend. So my guess is that you don't use that form to test your null hypothesis. 



#3
Nov2912, 11:25 PM

P: 7

As far as I know, the null hypothesis is that a unit root exists and that the process is therefore stationary. I had found that thread you linked, and I'm the user who started the short but ongoing discussion in the comments. The linked post (listed in the righthand column) is also mine, and is quite similar to my post here.
http://stats.stackexchange.com/q/44647 I don't understand why using a test with a trend term a_{1}t would be used with a quadratic term, however. Can you shed any light on that? 



#4
Nov3012, 01:14 AM

Sci Advisor
P: 3,173

Which DickeyFuller test should I apply to this time series?
The only light I can shed is that if [itex] f(t) = kt^2 [/itex] then
[itex] \triangle f(t) = f(t+1)  f(t) = k(t+1)^2  kt^2 = k( t^2 + 2t + 1)  kt^2 = 2kt + k [/itex] So the [itex] \triangle [/itex] of a model with a [itex] t^2 [/itex] term has a term linear in [itex] t [/itex]. 


Register to reply 
Related Discussions  
Geometric Series Test Rather Than Integral Test  Calculus & Beyond Homework  4  
Integral test and ratio test on haromonic series.  Calculus  1  
How to test if 2 timeseries variables are cointegrated?  Set Theory, Logic, Probability, Statistics  0  
Does the discriminant test apply to this conic?  Calculus & Beyond Homework  21  
Infinite Series Test (Ratio Test get 1)  Calculus & Beyond Homework  13 