nlv2000
高手看过来！请教一个SAS互谱分析的问题！ 我在参考陈磊发表在统计研究2001年第九期的文章《我国宏观经济指标周期波动相关性的互谱分析》中看到互谱分析结果中（表1）
周期长度(月) 相干谱平方 相位谱(弧度) 时差(月)
∞ 0·99 0 0
216 0·88 0·39 13·5
108 0·98 0·43 7·3
72 0·85 0·51 5·8
54 0·52 0·88 7·6
43·2 0·72 0·35 2·4
36 0·93 0·35 2·0
30·8 0·70 0·72 3·5
27 0·69 0·33 1·4
24 0·82 0 0
21·6 0·70 -0·3 -0·9
19·6 0·79 -0·10 -0·4
18 0·82 0·22 0·6
16·6 0·85 0·58 1·5
15·4 0·52 0·18 0·4
14·4 0·12 -1·13 -2·6
13·5 0·40 1·3 2·8
12·7 0·83 1·42 2·9
12 0·60 2·19 4·2
注:负号表示领先于工业总产值,正号则反之,0表示两者同步,以下同。
我不知道这个时差如何计算得出来。
我用SAS做时，振幅、相干谱和相位谱（弧度）等都可以直接得到。请问一下这个时差指的具体含义是什么。因为一直没有看到定义。
还有就是选用Tukey-Hanning窗时，如何把定义截断点M（比如=5）值写进程序中！我参考了SAS的说明中好像没有这方面的具体说明。我用SPSS做是可以做的，但是横坐标值太大，结果曲线都变成直线了，看不出结果来。
谢谢各位了！
nlv2000
我编的程序如下：
title "Cross Spectral Analysis of dlnrfe &dlnrgdp";
data gdpe;
input quarter dlnrfe dlnrgdp @;
datalines;
。。。
;
symbol1 i=splines v=dot;
proc spectra data=gdpe out=g cross coef a k p ph s;
var dlnrfe dlnrgdp;
weights tukey ;
run;
proc contents data=g position;
run;
symbol1 i=splines v=dot;
proc gplot data=g;
plot a_01_02 * freq;
plot k_01_02 * freq;
plot ph_01_02 * freq;
run;
proc gplot data=g;
plot a_01_02 * period;
plot k_01_02 * period;
plot ph_01_02 * period;
where period < 28;
run;
各位看看有什么问题吗？
