异方差实验报告

合集下载

异方差性实验报告doc

异方差性实验报告doc

异方差性实验报告篇一:计量经济学上机实验报告(异方差性)提示:打包保存时自己的文件夹以“学号姓名”为文件夹名,打包时文件夹内容包括:本实验报告、EViews工作文件。

篇二:Eviews异方差性实验报告实验一异方差性【实验目的】掌握异方差性问题出现的来源、后果、检验及修正的原理,以及相关的Eviews操作方法。

【实验内容】以《计量经济学学习指南与练习》补充习题4-16为数据,练习检查和克服模型的异方差的操作方法。

【4-16】表4-1给出了美国18个行业1988年研究开发(R&D)费用支出Y与销售收入X的数据。

请用帕克(Park)检验、戈里瑟(Gleiser)检验、G-Q检验与怀特(White)检验来检验Y关于X的回归模型是否存在异方差性?若存在【实验步骤】一检查模型是否存在异方差性1、图形分析检验(1)散点相关图分析做出销售收入X与研究开发费用Y的散点相关图(SCATX Y)。

观察相关图可以看出,随着销售收入的增加,研究开发费用的平均水平不断提高,但离散程度也逐步扩大。

这说明变量之间可能存在递增的异方差性。

(2)残差图分析首先对数据按照解释变量X由小至大进行排序(SORT X),然后建立一元线性回归方程(LS Y C X)。

因此,模型估计式为: Y?187.507?0.032*X ----------(*) ?(0.17)(2.88) R2=建立残差关于X的散点图,可以发现随着X增加,残差呈现明显的扩大趋势,表明存在递增的异方差。

2、Park检验建立回归模型(LS Y C X),结果如(*)式。

生成新变量序列: GENR LNE2 = LOG(RESID^2)GENR LNX = LOG(X)生成新残差序列对解释变量的回归模型(LS LNE2 C LNX)。

从下图所示的回归结果中可以看出,LNX的系数估计值不为0且能通过显著性检验,即随机误差项的方差与解释变量存在较强的相关关系,即认为存在异方差性。

实验四异方差性的检验与处理

实验四异方差性的检验与处理

实验四异方差性的检验与处理集团标准化工作小组 [Q8QX9QT-X8QQB8Q8-NQ8QJ8-M8QMN]实验四 异方差性的检验及处理(2学时)一、实验目的(1)、掌握异方差检验的基本方法; (2)、掌握异方差的处理方法。

二、实验学时:2学时 三、实验要求(1)掌握用MATLAB 软件实现异方差的检验和处理; (2)掌握异方差的检验和处理的基本步骤。

四、实验原理1、异方差检验的常用方法(1) 用X-Y 的散点图进行判断(2). 22ˆ(,)(,)e x e y 或的图形 ,),x )i i y i i ((e 或(e 的图形)(3) 等级相关系数法(又称Spearman 检验)是一种应用较广的方法,既可以用于大样本,也可与小样本。

检验的三个步骤 ① ˆt t y y=-i e②|i x i i 将e 取绝对值,并把|e 和按递增或递减次序排序,计算Spearman 系数rs ,其中:21ni i d =∑s 26r =1-n(n -1)③ 做等级相关系数的显着性检验。

n>8时,/2(2),t t n α>-反之,若||i i e x 说明与之间存在系统关系,异方差问题存在。

(4) 帕克(Park)检验帕克检验常用的函数形式:若在统计上是显着的,表明存在异方差性。

2、异方差性的处理方法: 加权最小二乘法 如果在检验过程中已经知道:222()()()i i i ji u Var u E u f x σσ===则将原模型变形为:121(i i p pi iy x x uf xβββ=+⋅++⋅+在该模型中:即满足同方差性。

于是可以用OLS估计其参数,得到关于参数12,,,pβββ的无偏、有效估计量。

五、实验举例例101i i iy x u=++若用线性模型,研究不同收入家庭的消费情况,试问原数据有无异方差性如果存在异方差性,应如何处理解:(一)编写程序如下:(1)等级相关系数法(详见文件)%%%%%%%%%%%%%%% 用等级相关系数法来检验异方差性 %%%%%%%%[data,head]=xlsread('');x=data(:,1); %提取第一列数据,即可支配收入xy=data(:,2); %提取第二列数据,即居民消费支出yplot(x,y,'k.'); % 画x和y的散点图xlabel('可支配收入x(千元)') % 对x轴加标签ylabel('居民消费支出y(千元)') % 对y轴加标签%%%%%%%% 调用regres函数进行一元线性回归 %%%%%%%%%%%%xdata=[ones(size(x,1),1),x]; %在x矩阵最左边加一列1,为线性回归做准备[b,bint,r,rint,s]=regress(y,xdata);yhat=xdata*b; %计算估计值y% 定义元胞数组,以元胞数组形式显示系数的估计值和估计值的95%置信区间head1={'系数的估计值','估计值的95%置信下限','估计值的95%置信上限'};[head1;num2cell([b,bint])]% 定义元胞数组,以元胞数组形式显示y的真实值,y的估计值,残差和残差的95%置信区间head2={'y的真实值','y的估计值','残差','残差的95%置信下限','残差的95%置信上限'};[head2;num2cell([y,yhat,r,rint])]% 定义元胞数组,以元胞数组形式显示判定系数,F统计量的观测值,检验的P值和误差方差的估计值head3={'判定系数','F统计量的观测值','检验的P值','误差方差的估计值'};[head3;num2cell(s)]%%%%%%%%%%%%% 残差分析 %%%%%%%%%%%%%%%%%%figure;rcoplot(r,rint) % 按顺序画出各组观测值对应的残差和残差的置信区间%%% 画估计值yhat与残差r的散点图figure;plot(yhat,r,'k.') % 画散点图xlabel('估计值yhat') % 对x轴加标签ylabel('残差r') % 对y轴加标签%%%%%%%%%%%% 调用corr函数计算皮尔曼等级相关系数res=abs(r); % 对残差r取绝对值[rs,p]=corr(x,res,'type','spearman')disp('其中rs为皮尔曼等级相关系数,p为p值');(2)帕克(park)检验法(详见文件)%%%%%%%%%%%%%%% 用帕克(park)检验法来检验异方差性 %%%%%%%[data,head]=xlsread(''); %导入数据x=data(:,1);y=data(:,2);%%%%%% 调用regstats函数进行一元线性回归,linear表带有常数项的线性模型,r表残差ST=regstats(y,x,'linear',{'yhat','r','standres'});scatter(x,.^2) % 画x与残差平方的散点图xlabel('可支配收入(x)') % 对x轴加标签ylabel('残差的平方') %对y轴加标签%%%%%%% 对原数据x和残差平方r^2取对数,并对log(x)和log(r^2)进行一元线性回归ST1=regstats(log(.^2),log(x),'linear',{'r','beta','tstat','fstat'})% 输出参数的估计值% 输出回归系数t检验的P值% 输出回归模型显着性检验的P值(3)加权最小二乘法(详见文件)%%%%%%%%%%% 调用robustfit函数作稳健回归 %%%%%%%%%%%%[data,head]=xlsread(''); % 导入数据x=data(:,1);y=data(:,2);% 调用robustfit函数作稳健回归,返回系数的估计值b和相关统计量stats[b,stats]=robustfit(x,y) %调用函数作稳健回归% 输出模型检验的P值%%% 绘制残差和权重的散点图 %%%%%%%plot,,'o') %绘制残差和权重的散点图xlabel('残差')ylabel('权重'(二)实验结果与分析:第一步::用OLS方法估计参数,并保留残差(1)散点图图可支配收入(x)居民消费支出(y)散点图因每个可支配收入x的值,都有5个居民消费收入y与之对应,所以上述散点图呈现此形状。

计量经济学实验二

计量经济学实验二

实验二〔一〕异方差性【实验目的】掌握异方差性的检验及处理方法【实验内容】建立并检验我国制造业利润函数模型【实验步骤】【例1】表1列出了1998年我国主要制造工业销售收入与销售利润的统计资料,请利用统计软件Eviews建立我国制造业利润函数模型。

一、检验异方差性⒈图形分析检验⑴观察销售利润〔Y〕与销售收入〔X〕的相关图(图1):SCAT X Y图1 我国制造工业销售利润与销售收入相关图从图中可以看出,随着销售收入的增加,销售利润的平均水平不断提高,但离散程度也逐步扩大。

这说明变量之间可能存在递增的异方差性。

⑵残差分析首先将数据排序〔命令格式为:SORT 解释变量〕,然后建立回归方程。

在方程窗口中点击Resids按钮就可以得到模型的残差分布图〔或建立方程后在Eviews工作文件窗口中点击resid对象来观察〕。

图2 我国制造业销售利润回归模型残差分布图2显示回归方程的残差分布有明显的扩大趋势,即说明存在异方差性。

⒉Goldfeld-Quant检验⑴将样本安解释变量排序〔SORT X〕并分成两部分〔分别有1到10共11个样本合19到28共10个样本〕⑵利用样本1建立回归模型1〔回归结果如图3〕,其残差平方和为。

SMPL 1 10LS Y C X图3 样本1回归结果⑶利用样本2建立回归模型2〔回归结果如图4〕,其残差平方和为。

SMPL 19 28 LS Y C X图4 样本2回归结果⑷计算F 统计量:12/RSS RSS F ==,21RSS RSS 和分别是模型1和模型2的残差平方和。

取05.0=α时,查F 分布表得44.3)1110,1110(05.0=----F ,而44.372.2405.0=>=F F ,所以存在异方差性⒊White 检验⑴建立回归模型:LS Y C X ,回归结果如图5。

图5 我国制造业销售利润回归模型⑵在方程窗口上点击View\Residual\Test\White Heteroskedastcity,检验结果如图6。

异方差性实验报告doc

异方差性实验报告doc

异方差性实验报告doc异方差性实验报告篇一:计量经济学上机实验报告(异方差性)提示:打包保存时自己的文件夹以“学号姓名”为文件夹名,打包时文件夹内容包括:本实验报告、EViews工作文件。

篇二:Eviews异方差性实验报告实验一异方差性【实验目的】掌握异方差性问题出现的来源、后果、检验及修正的原理,以及相关的Eviews操作方法。

【实验内容】以《计量经济学学习指南与练习》补充习题4-16为数据,练习检查和克服模型的异方差的操作方法。

【4-16】表4-1给出了美国18个行业1988年研究开发(R&D)费用支出Y与销售收入X的数据。

请用帕克(Park)检验、戈里瑟(Gleiser)检验、G-Q检验与怀特(White)检验来检验Y关于X的回归模型是否存在异方差性?若存在【实验步骤】一检查模型是否存在异方差性1、图形分析检验(1)散点相关图分析做出销售收入X与研究开发费用Y的散点相关图(SCATX Y)。

观察相关图可以看出,随着销售收入的增加,研究开发费用的平均水平不断提高,但离散程度也逐步扩大。

这说明变量之间可能存在递增的异方差性。

(2)残差图分析首先对数据按照解释变量X由小至大进行排序(SORT X),然后建立一元线性回归方程(LS Y C X)。

因此,模型估计式为: Y?187.507?0.032*X ----------(*) ?(0.17)(2.88) R2=建立残差关于X的散点图,可以发现随着X增加,残差呈现明显的扩大趋势,表明存在递增的异方差。

2、Park检验建立回归模型(LS Y C X),结果如(*)式。

生成新变量序列: GENR LNE2 = LOG(RESID^2)GENR LNX = LOG(X)生成新残差序列对解释变量的回归模型(LS LNE2 C LNX)。

从下图所示的回归结果中可以看出,LNX的系数估计值不为0且能通过显著性检验,即随机误差项的方差与解释变量存在较强的相关关系,即认为存在异方差性。

异方差的检验与修正

异方差的检验与修正

西安财经学院本科实验报告学院(部)统计学院实验室 313 课程名称计量经济学学生姓名学号 1204100213 专业统计学教务处制2014年12 月 15 日《异方差》实验报告开课实验室:313 2014年12月22第六部分异方差与自相关4. 在本例中,参数估计的结果为:2709030.01402097.01402.728X X Y ++=Λ(2.218) (2.438) (16.999)922173.02=R D.W.=1.4289 F=165.8853 SE=395.2538三.检查模型是否存在异方差 1.图形分析检验 (1)散点相关图分析分别做出X1和Y 、X2和Y 的散点相关图,观察相关图可以看出,随着X1、X2的增加,Y 也增加,但离散程度逐步扩大,尤其表现在X1和Y.这说明变量之间可能存在递增的异方差性。

在Graph/scatter 输入log(x2) e^2,结果如下:(2)残差相关图分析建立残差关于X1、X2的散点图,可以发现随着X 的增加,残差呈现明显的扩大趋势,表明模型很可能存在递增的异方差性。

但是否确实存在异方差还应通过更进一步的检验。

2.GQ 检验首先在主窗口Procs菜单里选Sort current page命令,输入排序变量x2,以递增型排序对解释变量X2进行排序,然后构造子样本区间,分别为1-12和20-31,再分别建立回归模型。

(1)在Sample菜单里,将区间定义为1—12,然后用OLS方法求得如下结果(2)在Sample菜单里,将区间定义为20—31,然后用OLS方法求得如下结果则F的统计量值为:6699.834542929948192122===∑∑iieeF在05.0=α下,式中分子、分母的自由度均为9,查F分布表得临界值为:18.3)9,9(05.0=F,因为F=8.6699>18.3)9,9(05.0=F,所以拒绝原假设,表明模型确实存在异方差。

异方差实验报告

异方差实验报告

异方差实验报告引言异方差(heteroscedasticity)是指随着自变量的变化,因变量的方差也随之变化的现象。

在统计分析中,假设方差是恒定的是很常见的,但在实际应用中,许多变量的方差是不恒定的,需要进行异方差处理。

本实验旨在通过模拟数据和实际数据来探究异方差的影响并了解异方差检验方法。

实验设计本实验分为两个部分。

第一部分使用模拟数据,提供了不同阶段下的异方差数据集。

第二部分使用实际数据,通过观察数据的模式来判断是否存在异方差。

实验方法模拟数据在模拟数据部分,我们生成了四个数据集,每个数据集都包含一个自变量和一个因变量。

为了模拟异方差,我们设定了不同的标准差,并与自变量呈一定的关系。

具体参数如下:•数据集1:使用正态分布生成自变量和因变量,因变量的标准差为自变量的两倍。

•数据集2:自变量为正态分布,因变量为自变量的2次方,并加入了一个随机误差项,使得方差在自变量变大时也会变大。

•数据集3:自变量为均匀分布,因变量为自变量的指数函数,并加入了一个随机误差项,使得方差在自变量变大时也会变大。

•数据集4:自变量为正态分布,因变量为自变量的对数,并加入了一个随机误差项,使得方差在自变量变大时也会变大。

实际数据在实际数据部分,我们使用了一份销售数据。

该数据包含了不同日期下的产品销售量和价格。

我们首先观察数据的散点图,并通过直观感受来猜测是否存在异方差。

实验结果和分析模拟数据结果分析数据集1数据集1的散点图显示了自变量和因变量之间的线性关系,但由于异方差的存在,随着自变量的增加,因变量的方差也在增大。

这说明了异方差对回归结果的影响。

数据集2数据集2的散点图显示了自变量和因变量之间的非线性关系。

由于自变量的增大,因变量的方差也在增大。

这与模型中设定的异方差关系一致。

数据集3数据集3的散点图显示了自变量和因变量之间的指数关系。

随着自变量的增大,因变量的方差也在增大,符合预期的异方差模式。

数据集4数据集4的散点图显示了自变量和因变量之间的对数关系。

计量经济学实验5 异方差

计量经济学实验5 异方差

具体步骤是:
1 .选择普通最小二乘法估计原模型,得到随机误差 项的近似估计量 û t; 2.建立 1/| û t | 的数据序列; 3.选择加权最小二乘法,以 1/| û t |序列作为权,进
行估计得到参数估计量。实际上是以 1/| û t |乘原模型的两
边,得到一个新模型,采用普通最小二乘法估计新模型。
以不必把它们全包括在内。无交叉项选项仅使用解释变
量平方进行检验回归。
例:人均家庭交通及通讯支出(CUM)和可支配收入(IN ) 的回归方程的 White 异方差检验的结果:
该结果F 统计量和 Obs*R2 统计量的P值均很小,表明 拒绝原假设,即残差存在异方差性。
利用加权最小二乘法消除异方差
1.方差已知的情形 假设有已知形式的异方差性,并且有序列w,其值与误差标 准差的倒数成比例。这时可以采用权数序列为w 的加权最小二乘 估计来修正异方差性。对加权自变量和因变量最小化残差平方和 得到估计结果 :
四、实验原理与操作
异方差性检验
1. 图示检验法 (1) 用X-Y的散点图进行判断 观察是否存在明显的散点扩大、缩小或复杂型趋势(即 不在一个固定的带型域中)
2的散点图进行判断 (2)X - û i 首先采用OLS方法估计模型,以求得随机误差项的估计量 (注意,该估计量是不严格的),我们称之为“近似估计量”,用 2 表示。于是有 û i
5.随机误差项服从0均值、同方差的正态分布。即
ui ~
N (0, )
2
i=1,2,…,N
当随机误差项满足假定1 ~ 4时,将回归模型”称为 “标准回归模型”,当随机误差项满足假定1 ~ 5时,将回 归模型称为“标准正态回归模型”。如果实际模型满足不 了这些假定,普通最小二乘法就不再适用,而要发展其他 方法来估计模型。

异方差的检验与修正

异方差的检验与修正

西安财经学院本科实验报告学院(部)统计学院实验室313课程名称计量经济学学生姓名学号1204100213专业统计学教务处制2014年12 月15 日《异方差》实验报告五、实验过程原始记录(数据、图表、计算等) 一.选择数据1.建立工作文件并录入数据File\New\workfile, 弹出Workfile create 对话框中选择数据类型.Object\new object\group,按向上的方向键,出现两个obs 后输入数据.中国内地2006年各地区农村居民家庭人均纯收入与消费支出 单位:元城市 y x1 x2 城市 y x1 x2 北京 5724。

5 958.3 7317。

2 湖北 2732。

5 1934。

6 1484。

8 天津 3341。

1 1738.9 4489 湖南 3013。

3 1342.6 2047 河北 2495。

3 1607。

1 2194。

7 广东 3886 1313。

9 3765.9 山西 2253.3 1188。

2 1992.7 广西 2413。

9 1596。

9 1173。

6 内蒙古 2772 2560.8 781.1 海南 2232。

2 2213。

2 1042.3 辽宁 3066。

9 2026。

1 2064。

3 重庆 2205。

2 1234.1 1639。

7 吉林 2700.7 2623。

2 1017。

9 四川 2395 1405 1597.4 黑龙江 2618。

2 2622.9 929.5 贵州 1627。

1 961。

4 1023。

2 上海 8006 532 8606.7 云南 2195.6 1570。

3 680。

2 江苏 4135.2 1497。

9 4315.3 西藏 2002。

2 1399.1 1035.9 浙江 6057。

2 1403.1 5931。

7 陕西 2181 1070。

4 1189。

8 安徽 2420。

9 1472。

8 1496。

3 甘肃 1855.5 1167。

  1. 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
  2. 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
  3. 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。

《计量经济学》实训报告
实训项目名称异方差的检验及修正
实训时间 2011年12月13日
实训地点
班级
学号
姓名
实训(实践) 报告
实训名称异方差的检验及修正
一、实训目的
深刻理解异方差性的实质、异方差出现的原因、异方差的出现对模型的不良影响(即异方差的后果),掌握估计和检验异方差性的基本思想和修正异方差的若干方法;能够运用所学的知识处理模型中的出现的异方差问题,并要求初步掌握用EViews处理异方差的基本操作方法。

二、实训要求
使用教材第五章的数据做异方差的图形法检验、Goldfeld-Quanadt检验与White检验,使用WLS法对异方差进行修正。

三、实训内容
1、用图示法、戈德菲尔德、white验证法,验证该模型是否存在异方差。

2、用加权最小二乘法消除异方差。

四、实训步骤
练习题5.8数据1998年我国重要制造业销售收入和销售利润的数据
Y—销售利润,x—销售收入
1. 用OLS方法估计参数,建立回归模型:ls y c x
回归结果如下:
Y=12.036+0.1044x;
S = (19.5178) (0.00844)
T= (0.6167) (12.3667)
R^2=0.8547 S.E.=56.9037
2.检验是否存在异方差
(1) 图形检验:残差图形scat x e2
结果表明:
残差平方e2对解释变量的x的散点图主要分布在图形的下方,大致看出残差平方随X 的变动呈增大的趋势,因此,模型很可能出现异方差。

(2)戈德菲尔德-夸特检验
首先,对变量进行排序,在这个题目中,我选择递增型排序,这是y与x将以x按递增型排序。

然后构造子样本区间,建立回归模型。

在本题目中,n=28,删除中间的1/4,的观测值,即大约8个观测值,剩余部分平分得两个样本区间:1—10和19-28,他们的样本个数均为10。

用OLS方法得到前10个数的样本结果(ls y c x):
用OLS方法得到后10个数的样本结果(ls y c x):
接着,根据戈德菲尔德检验得到F统计量:(两个残差平方和相除,大的除以小的)F=63769.67/2577.969=24.736。

最后,进行判断。

在a=0.05,分子分母的自由度均为10-1-1=8,查F分布表,得到临界值F(8,8)=3.44,因为F=24.736 > 3.44,所以拒绝原假设。

结果表明:模型存在异方差。

(3)White检验
结论表明:查看P的概率小于0.05,存在明显的异方差。

3.消除异方差
使用加权最小二乘法(WLS)
选用权数Wl=1/x、W2=1/x2、W3=1/sqr(x),分别选用各个权数进行估计检验权数为w1时:
权数为w2时:
权数为w3时:
结果表明:w2权数是最好的。

五、实训分析、总结
1.可以看出运用加权最小二乘法消除了异方差之后,参数的t检验均显著,可绝系数大幅度提高,F检验也显著
2.销售收入每相差1元钱,销售利润就相差0.106892。

六、实训报告评价与成绩。

相关文档
最新文档