pacf(x, lag.max, plot, na.action, ...)

有两个办法：
1，要求x与y的偏相关系数，如果Z是第三个相关变量或者多个变量的统称。
只要拿x和y分别对z求线性回归，然后得到相应的残差序列，然后求两个残差序列
的相关系数，这就是偏相关系数。
这个方法要用到:lm(),resid(model),cor(x1,x2)。

2，第二种办法比较简单，就是求出所有变量的相关系数矩阵，然后求这个矩阵的逆矩阵。里面对应
的元素就是偏相关系数。这个只要用到：cor(data.frame),solve(matrix);