在网上搜索到相关的R脚本，见下面的地址http://www.stat.pitt.edu/stoffer/tsa2/Examples.htm

由于该脚本存在一个问题，即滞后数和对应的数值错位的问题，我做了小小修改，在plot中采用LAG+1，使得lag值和统计值在图中对应。
现将我修改后的代码文件附在后面。