Cálculo manual PACF

Estoy tratando de replicar el cálculo que hacen SAS y SPSS para la función de autocorrelación parcial (PACF). En SAS se produce a través de Proc Arima. Los valores PACF son los coeficientes de una autorregresión de la serie de interés sobre los valores rezagados de la serie. Mi variable de interés...