b34sexec matrix; call load(data_acf); call print(data_acf); ar=array(:.7,-.2); n=3000; start=array(:.1 .1); testar1a=genarma(ar,ma,1.0,start,.1,n,10); ar=array(:-.7,.2); testar1b=genarma(ar,ma,1.0,start,.1,n,10); acf1a=acf(testar1a,20,se,pacf1a); acf1b=acf(testar1b,20,se,pacf1b); call data_acf(testar1a,'model 1',300); call data_acf(testar1b,'model 2',300); call autobj(testar1a :print :autobuild); call data_acf(%res,'res 1',300); call spectral(%res,sinx,cosx,px,sx1a,freq:1 2 3 2 1); call graph(freq,sx1a:heading 'Spectrum of res1' :plottype xyplot); call autobj(testar1b :print :autobuild); call data_acf(%res,'res 2',300); call spectral(%res,sinx,cosx,px,sx1a,freq:1 2 3 2 1); call graph(freq,sx1a:heading 'Spectrum of res2' :plottype xyplot); call spectral(testar1a,sinx,cosx,px,sx1a,freq:1 2 3 2 1); call spectral(testar1b,sinx,cosx,px,sx1b,freq:1 2 3 2 1); call graph(freq,sx1a:heading 'Spectrum of ar(.9)' :plottype xyplot); call graph(freq,sx1b:heading 'Spectrum of ar(-.9)' :plottype xyplot); call tabulate(acf1a,acf1b,pacf1a,pacf1b); b34srun;