The Dickey-Fuller test aims to test the presence of unit roots in a process. If I've the following process:
$$y_t = a_1y_{t-1} +\epsilon_t$$
it tests $H_0\colon \ a_1=1$.
If I specify $\theta =a_1-1$ I can rewrite the process as:
$$(a_1-1)y_{t-1} +\epsilon_t = \theta y_{t-1} + \epsilon_t$$
and test $H_0\colon \ \theta=0$.
My question is: why does the distribution of the $t$-statistic have to be computed using Monte Carlo method?