PACF ручной расчет

Я пытаюсь повторить расчет, который SAS и SPSS делают для функции частичной автокорреляции (PACF). В SAS производится через Proc Arima. Значения PACF - это коэффициенты авторегрессии интересующего ряда по запаздывающим значениям ряда. Моя переменная интереса - продажи, поэтому я вычисляю lag1, lag2...