Define the sufficient data sums, all over ,
Multiplying the likelihood function by the independent standard-normal priors and collecting the quadratic terms gives
Let and . This precision matrix is , where , so it is positive definite even for a short or singular design. Completing the square proves
The fully normalized posterior density is
This is Gaussian conjugacy for an initialized AR(2) regression.
Completing each one-dimensional square, or using the supplied conditional-normal identity, gives
When , all regressor sums vanish and the posterior remains the independent standard-normal prior. No stationary-parameter restriction is imposed: the specified prior is on all of , and the finite initialized chain is defined for all .