关于 一元线性回归、异方差性与Stata实现
第一次blog献给了CSDN与计量经济学,但是不得不说 markdown编辑器还行
(大抵约
莫这个经济学专业的少女年 所学习过的 语言的太少了呢,基本HTML一点都不会rua)
既计量经济学作业,遂顺手笔记耳
文章目录
- 1. 一元线性回归 Stata实现
- 2. 检验异方差
- 2.1检验方法:回归残差图
- 2.2检验方法:Breusch-Pagan / Cook-Weisberg test (for heteroskedasticity)
- 2.3 检验方法:White Test
- 2.4 Spearman’s correlation coefficient for ranked data 斯皮尔曼等级相关检验法
- 3. 建立存在异方差的估计模型
- 3.1稳健(Robust)方差的回归模型
- 3.2加权最小二乘法
- 第9章 论那个少年即金融工程学生的行为
1. 一元线性回归 Stata实现
1.1命令窗口:
. rename var1 x //把导入的数据的 变量名 改为 x (default : var1). rename var2 y. label values x 居民可支配收入 //把 变量x 贴一个标签 "居民可支配收入". label values y 居民消费支出 //把 变量y 贴一个标签 "居民消费支出". regress y x //进行OLS回归
1.2菜单栏:
“统计”—“线性模型及相关”—“线性回归”
2. 检验异方差
“所谓异方差, 是指我们的误差项u 对不同的个体是不同的,从统计上讲,就是违反了假设: 球形误差Cov(u,x)=0的基本假设”
( 即: 随机误差 和 自变量 存在一定关系 , “随机”误差 不是随机的 )
2.1检验方法:回归残差图
2.1.1命令窗口:
.rvpplot x //完成回归后,显示 “残差-自变量 图”
2.1.2菜单栏:
“图形”—“回归诊断图”—“残差对解释变量”
2.1.3结果:
∴ 显然 随机误差 随着x 的变化 而改变
※!!!也可以做 u^2-x (残差平方-自变量)的图像
(
btw.
)
2.2检验方法:Breusch-Pagan / Cook-Weisberg test (for heteroskedasticity)
2.2.1命令窗口:
. estat hettest //完成回归后,tests for heteroskedasticity// .estat :回归后的后验估计统计量
2.2.2菜单栏:
“统计”—“回归诊断”—“模型检定等”
2.2.2结果:
2.3 检验方法:White Test
2.3.1命令窗口:
. estat imtest, white // imtest : imformation matrix test
2.3.2菜单栏:
“统计”—“回归诊断”—“模型检定等”
2.3.3结果:
(貌似 White检验 被 C-T检验整合、修正了?)
P值=0.0659 > α=0.05 ∴接受H_0:存在异方差
2.4 Spearman’s correlation coefficient for ranked data 斯皮尔曼等级相关检验法
2.4.1命令窗口:
.regress y x //对原模型应用OLS.predict res , residuals //把回归的残差 保存到 变量 res 中.generate abs_res = abs(res) // 对残差求绝对值,然后保存在 变量abs_res中.sort abs_res // 对 残差的绝对值 排序(升序)// !!!注意:sort 命令会 扩展到整个数据集,所以重新导入了 已排序的自变量.generate sort_x .spearman sort_x abs_res //对已排序的自变量和残差绝对值 做Spearman Test//一定记得“绝对值”“排序”
2.4.2 结果
3. 建立存在异方差的估计模型
3.1稳健(Robust)方差的回归模型
3.1.1命令窗口:
.regress y x, vce(robust) // 或者 . regress y x ,robust
3.1.2菜单栏:
“统计”—“线性模型及相关”—“线性回归”"SE/Robust" —“稳健”
3.1.3结果:
3.2加权最小二乘法
3.2.1命令窗口:
.quietly regress y x // 进行OLS回归, 在regress,命令前加上quietly选项,不显示回归结果
.predict u ,residuals //用 predict命令生成残差 u
.generate sigma = abs(u) // 另(σ_i ) ̂=|u_i |,保存在 sigma 变量中.generate z = y/sigma
.generate x1 = 1/sigma
.generate x2 = x/sigma .regress z x1 x2 , noconstant
3.2.2结果
对回归结果进行 White检验
数据数据!
第9章 论那个少年即金融工程学生的行为
233333 真是hardcore呢,懂的人自然懂吧
出自 AN INQUIRY INTO THE NATURE AND CAUSES OF THE WEALTH OF NATURE
那个不被人理解的,谦逊而又质朴的大师
……
凡愿意细心探讨这个极重要科学的原来的人,都得对他十分留意
关于 一元线性回归、异方差性与Stata实现相关推荐
- STATA:陈强教授 计量经济学及stata应用第四章 一元线性回归
第四章 一元线性回归 一元回归的stata实例 (1)带常数项的回归 (2)不带常数项的回归 stata命令运行结果的储存与调用 4.9 总体回归函数与样本回归函数:蒙特卡洛模拟 附录(1) 随机数的 ...
- 一元线性回归决定系数_回归分析|笔记整理(1)——引入,一元线性回归(上)...
大家好! 新学期开始了,不知道大家又是否能够适应新的一学期呢?先祝所有大学生和中小学生开学快乐! 本学期我的专业课是概率论,回归分析,偏微分方程,数值代数,数值逼近,金融时间序列分析,应用金融计量学和 ...
- 统计学 一元线性回归
统计学 一元线性回归 回归(Regression):假定因变量与自变量之间有某种关系,并把这种关系用适当的数学模型表达出来,利用该模型根据给定的自变量来预测因变量 线性回归:因变量和自变量之间是线性关 ...
- TensorFlow基础3-机器学习基础知识(解析法实现一元线性回归、多元线性回归)
记录TensorFlow听课笔记 文章目录 记录TensorFlow听课笔记 一,机器学习基础 1.1一元线性回归 1.2解析法实现一元线性回归 1.3解析法实现多元线性回归 二,代码实现一元/多元回 ...
- 机器学习公开课笔记(1):机器学习简介及一元线性回归
初步介绍 监督式学习: 给定数据集并且知道其正确的输出应该是怎么样的,即有反馈(feedback),分为 回归 (Regressioin): map输入到连续的输出值. 分类 (Classificat ...
- UA MATH 571B 回归 QE练习题 一元线性回归理论
UA MATH 571B 回归 QE练习题 一元线性回归理论 2015/1/5 2015/5/5 2016/5/6 2017/1/5 2017/5/6 这是2015年1月第五题,2015年5月第五题, ...
- UA MATH571A 一元线性回归II 统计推断2
UA MATH571A 一元线性回归II 统计推断2 β0\beta_0β0的分布 拟合与预测 拟合值的区间估计 预测值的区间估计 数值例子:女性肌肉量与年龄的关系 β0\beta_0β0的分布 ...
- UA MATH571A 一元线性回归IV 模型诊断
UA MATH571A 一元线性回归IV 模型诊断 解释变量 解释变量的可视化 残差 残差的性质 Semistudentized Residual 残差的可视化 残差关于解释变量的图 残差关于拟合值的 ...
- UA MATH571A 一元线性回归III 方差分析与相关性分析
UA MATH571A 一元线性回归III 方差分析与相关性分析 ANOVA Table F检验 回归系数的F检验 F检验与t检验等价 广义线性检验方法 R2R^2R2 数值例子:女性肌肉量与年龄的关 ...
最新文章
- cutadapt 处理fastq文件
- ads无法启用状态服务器,NAC ADSSO 无法工作在Microsoft 2008服务器版本
- AlertDialog.Builder setCancelable用法
- 5.7 echo:显示一行文本
- PAT_B_1055_Java(25分)
- 实验11 编写子程序
- nw.js package一般设置
- 敏捷开发用户故事系列之八:验收标准
- sap产品图谱 - road to sap.pdf_蛇胆陈皮胶囊化学成分及指纹图谱研究
- 【Servlet笔记】Servlet入门
- javase二维数组笔记
- Qt 窗口嵌入桌面 窗口 放在桌面
- java项目启动就调用方法_Java项目启动时执行指定方法的几种方式
- 仿QQ项目(一):好友列表
- 【吐槽】Android 第一步,雀氏纸尿裤
- mongdb系列之最详细基础知识
- Live2D桌面动画模型,在桌面显示
- python通过jira的接口操作jira
- 洛谷P1149(NOIP2008) 火柴棒等式 (C语言 + 详细注释)
- 浙大计算机跨专业考研心路历程记录