Scatter
Posted: Fri Jul 16, 2010 6:05 am
Code: Select all
open data sum.txt
calendar(q) 1947:1
data(format=free,org=columns) 1947:1 2009:1 lgdp
set lgdp / = 100*lgdp
compute nobs=2009:1-1949:1+1
Because I need to make some changes to the graphic option, I follow the default setting in @spectrum but I did not use the .src directly
Code: Select all
@bkfilter(nodiff) lgdp / bkcycle
set frequencies 1 256 = (t-1.0)/256
frequency 2 512
rtoc 1949:1 2009:1 1
# bkcycled
# 1
fft 1
cmult(scale=1.0/(2.0*%pi*nobs)) 1 1
window(type=flat) 1 / 2
ctor 1 256
# 2
# bk_spect
Code: Select all
scatter(hshade=||0.05,0.257||) 1
# frequencies bk_spect 1 256