This paper shows the validity of the (arbitrary) higher order Edgeworth expansion for the distribution of estimators of the coefficient parameter θ∈(−1, 1) in the AR(1) process {Xt} with a possibly nonzero initial value X0=x. The stationary case of X0∼N(0, 1/(1−θ2)) is also treated.