powered by
Conditional Gaussian likelihood (good enough for sampling steps)
llike_pacf(data, psi, sigma2)