Cálculo manual do PACF

Estou tentando replicar o cálculo que o SAS e o SPSS fazem para a função de autocorrelação parcial (PACF). No SAS, é produzido através do Proc Arima. Os valores PACF são os coeficientes de uma auto-regressão da série de interesse em valores defasados ​​da série. Minha variável de interesse é...