workfile corr u 1 500
series result
for !i=1 to 500
smpl 1 100 series x=nrnd series y=nrnd
series xx series yy
scalar sum1=0 scalar sum2=0
for !counter=1 to 100 sum1=sum1+x(!counter) sum2=sum2+y(!counter) xx(!counter)=sum1 yy(!counter)=sum2
next
scalar r=@cor(xx,yy) result(!i)=r
next
result.hist | 定义一个非时间序列(u)工作文件,corr,容量为500。 定义一个空序列result,用来存储相关系数的计算结果。 !i为控制变量,通过一个for循环语句使计算进行500次。 把样本范围设置成100。 生成两个互不相关的白噪声序列x、y,样本容量100。
定义两个空的序列xx和yy,样本容量也是100。
定义两个标量sum1和sum2,初始值为0。
!counter为控制变量,在这个for循环中,分别对序列x和y进行一次累加生成两个一阶单整的序列,将结果分别放到序列xx和yy中。
累加一次。
计算序列xx和yy的相关系数,并将结果放到标量r中。
将相关系数计算结果放到序列result中,在这个for循环中,这个操作要进行500次。
显示序列result的直方图以及有关统计量。 |