计量上机指南实验模拟.pdfVIP

上机实验14：Eviews中的模拟

一、随机数的生成

（1）生成服从某种分布的随机数序列

【例】生成标准正态分布、指数分布、poisson分布、t分布的随机数序列

EViews程序如下：

生成服从某种分布的随机数序列（T=1000），并存入random1文件。

workfilerandom1u11000

seriesY1=@rnorm

seriesY2=nrnd

seriesY3=@rexp(3)

seriesY4=@rpoisson(3)

seriesY5=@rtdist(3)

（3）生成服从某种分布随机数序列基础上求分位数值，求频数分布表。

计算分位数的EViews命令为：

scalarq=@quantile(y,τ,s)，

其中y表示求分位数所用的序列；τ表示要计算的分位数；s取1~6，依次

表示6种计算分位数序数方法，上命令的含义是计算第τ分位数所对应的分

位数值，存入标量q中。

【例】以生成的Y1为例，求0.2分位数和中位数值。

在空白处键入命令：

scalarq1=@quantile(Y1,.2,1)

scalarq2=@quantile(Y1,.5,1)

注意：q1和q2以标量形式存入工作文件。双击之，EViews6可以在窗口

下方显示该值。EViews7则用矩阵窗口的形式显示。

（4）查看随机数列的频数、频率（百分数）、累积频数、频率（百分数）直方图

【例】以Y1为例，查看频数、频率（百分数）、累积频数、频率（百分数）直方图。

激活序列窗口，点击View选one-waytabulation功能，点击OK键。

上述合并在一起，变成如下：

生成服从某种分布的随机数序列（T=1000），并存入random1文件。

workfilerandom1u11000

seriesY1=@rnorm

seriesY2=nrnd

seriesY3=@rexp(3)

seriesY4=@rpoisson(3)

seriesY5=@rtdist(3)

seriesZ=nrnd

seriesXZ*2+50

scalarq1=@quantile(Y1,.2,1)

scalarq2=@quantile(Y1,.5,1)

二、研究统计量的分布

命令行

程序窗口

程序

从均匀分布/指数分布/任意分布中每次抽取100个样本数据，计算其均值。然后重复1000次，

查看其均值的直方图。

Workfilemyu11000

for!k=1to1000

seriesx{!k}

smpl1100

x{!k}=@rexp(3)

scalarxm=@mean(x{!k})

smpl11000

vector(1000)v

v(!k)=@mean(x{!k})

next

deletex??x???

mtos(v,y)

y.hist

三、估计圆周率

workfilemyworku1000

seriesx

seriesy

x=@runif(0,1)

y=@runif(0,1)

genrz=4*((x-0.5)^2+(y-0.5)^20.5^2)

scatxy

scalarzm=@mean(z)

四、一元回归的模拟：抽样3次，估计3次

workfile