ols残差_涨知识丨OLS原理的矩阵方法很难?Just So So

53678685d1622fc92e8f7f88d1cb9716.png

对计量经济学初学者而言,OLS原理的矩阵表示通常令人“发怵”。其原因主要在于,至少在财经类课程体系中,关于矩阵微分的先行课程是缺失的。鉴于计量经济学的进阶课程大多采用矩阵语言,笔者认为有必要专文论述如何“搞掂”关于OLS原理的矩阵方法,以降低后续学习的门槛。

一、从OLS的基本原理谈起

对于多元回归模型(1):19758e5c5566e0eaca1dfbf4f2b627a1.pngOLS原理就是,选择参数估计值以使得残差平方和最小,即:d38c8b75f270396b6e4a7ad02f804db4.png若定义目标函数为Q,则由上述最优化问题的一阶条件可形成一个包括k+1个正规方程的方程组。 3a7cc9c25eb54e1c626f3c5f3af03e93.png求解上述正规方程组(3),即获得各个参数的OLS估计量。现在若我们引入向量与矩阵定义:3525c07ac38d79bbe0ceb284eb8d5f7a.png则多元回归模型(1)可表示为:8a73407bdbfd10b6f77b200d051bb521.png最优化问题(2)可表示为:25f95b107c8b3e1e653e3c0ca3bad8fa.png正规方程组(3)可表示为:e3a448b18bfcabf8a4970eb94da24ac1.png

二、矩阵微分规则的引出与应用

我们考察式(6)。在这里,a4889f61ba12e202c61730188c4acb26.png0是k+1维列向量。用式(6)来描述正规方程组(3),看似十分平凡,但其实隐含了一个关于矩阵微分的一般规则:一个标量对一个m维列向量求导,等价于该标量对这个m维列向量中的每一个元素求导,其求导结果是一个m维列向量。这是一个简单而重要的规则,接下来我们将反复利用此规则。最优化问题(5)的目标函数Q可进一步展开成:bccbe70d31f5b0eb94cdfa787cffb2a9.png由于标量Q只可能被分解成标量,式(7)中最后一个等号右边的四项均为标量,并且有:7251e8fa5f7055af15f77b96110064ce.png根据式(8),我们需依次解决四个问题:

(一)ad179ecc2a6414d622ef9d10b9c2dfe5.png

标量8b6e86fd287ffd9e8d85d976343ff4fb.png,其不是dcafa0c168ea78c5847c6f1a458d1ea8.png中任何元素的函数。因此,有:66a3c2aede430ee5fdde521f366d2aeb.png从形式上看,式(9)与我们在《微积分》课程中所熟悉的微分规则c29623f150970e1ac7866927885d1af4.png是一致的,其中8d977c05840bf65ab1568f9b6d0e88ab.png为常数,151172689d7ed101a51dc614a2b0cbe6.png为变量。这里的0是标量,而式(9)中的0是k+1维列向量。

(二)081774ba0834179ee478f3481b8eea9d.png

由于c7cb1d7b4ec60f5a6b2925f5d128b46c.png为标量,而dcafa0c168ea78c5847c6f1a458d1ea8.png为k+1维列向量,我们可迅速判断fa383c9f911abd0cfeab70a56b9cfa4c.png为k+1维行向量。若定义:b5e7f996bc26b2013c18ad43c892ce76.png,则0591e1a0f68225cfc07e6cbceddcea9a.png。显然有:86dfcbfddecf8af1895e603ffcb1602f.png从形式上看,式(10)与我们在《微积分》课程中所熟悉的微分规则c773bab052beea011b8665424557975a.png是一致的。关键的差别在于,在式(10)中,fa383c9f911abd0cfeab70a56b9cfa4c.png不能同8d977c05840bf65ab1568f9b6d0e88ab.png那样,被直接置于等号右边——为了满足矩阵微分规则,我们还需对其进行转置处理,以使其变为一个列向量。

(三)3482e98b9cc433f153ed8bb18d37b848.png

敏锐的读者会发现,由于标量的转置等于标量本身,即有:1190934128f65da592e8ce144bc72945.png故有:bab2cfdce6c30558846a7f9251da4d63.png当然,我们还可将标量f9c78f75f118798dbdfa479055558f41.png中的因式55d429a40f162fb02c64ee1ba63145a3.png定义为列向量229558986e59b841c8f14cf5a013066b.png,从而有:f3d9dfe81ca2c2f16bb0423edca21d82.png。因此,a7b123a06f90ffa745b3544d708a0694.png。从形式上看,式(12)与我们在《微积分》课程中所熟悉的微分规则c773bab052beea011b8665424557975a.png是一致的。重要的是,在这里,55d429a40f162fb02c64ee1ba63145a3.png能同8d977c05840bf65ab1568f9b6d0e88ab.png那样,被直接置于等号右边。如此处理满足矩阵微分规则,原因在于55d429a40f162fb02c64ee1ba63145a3.png是一个列向量。

(四)d44302cda15adb2bd8b0f0c8d67d1f60.png

现在我们碰到了棘手的问题——在42e74201b8e2a051f1af8b1e586bfc14.png中,dcafa0c168ea78c5847c6f1a458d1ea8.png出现了两次。但非常幸运的是,对问题(一)、(二)与(三)的讨论已暗示,这个求导结果也应该与传统的微分规则具有一定的一致性。关于函数积的微分规则表明:ded237f1681cb6b8cd07d8a372515e8a.png参照式(13)中第二个等号右边的表达式,我们可以猜测:f1d146b682ab7138180de1c3bb628ced.png在这里,我们很容易注意到:2827b206e2228a4ed16dc9a875b699db.png为列向量;14a07337070fdbf83be069d897cc3320.png为行向量。为了满足矩阵微分规则,我们需要对行向量14a07337070fdbf83be069d897cc3320.png取转置,以将其转化为列向量。在回答问题(一)、(二)与(三)时,我们对相应的矩阵微分规则进行了具体的验证。当然我们也可以验证式(14)是成立的,但由于比较复杂,在此略去。三、OLS估计量的“三步”记忆法及启示将式(9)、(10)、(12)、(14)带入式(8),并结合式(6),有:e5d630d4deefdc21ea388702eb342bce.png进而有:91b2237a160af20de096d7dbc904f7eb.png假定34cdb16a499fd3b4bc7dd0f8b3f6871a.png888918bd2a5874b501635ec2367dde49.png逆存在,则有:d72ad230e0b3954a21ea351e4e266342.png

(一)“三步”记忆法

我们可通过如下“三步”来记忆式(17):Step1:918f9b113c800bcee32f80a5a9aa1265.png。注意,等式两边左乘ec47400532cbf3788145d5ed0a39a7d6.png而不是559909f25a2127021e9fc737300cf0d0.png。考虑矩阵的维数,08f5b9a3b6a42eef76ed96dbfa90b33d.png显然是无意义的。Step2:eeabd5a49beb3b457d12e62d9f59f3b6.png。前提是91b4bc8bc2f4e9eb5c79198525835970.png存在。Step3:省略c96c50eba9239b72ed6f44d6dfdb9acb.png,则有2a1c8b66079e23870967044b8eecc9d6.png。既然两者近似相等,那么就以0fe06ac8e1fb620f940dfbdeeab82ecc.png作为9e0b42389c83d97fb293cc6b830c672a.png的估计量。

(二)启示

如果c96c50eba9239b72ed6f44d6dfdb9acb.png能省略,那么意味着1b9ed5865dfac0134c81085b1c1a1afa.png。但这有何依据呢?为了回答此问题,我们来考察列向量f59528b710b24788c1bf515f913723b7.png25abca5b22c9beb04a9f624cf0e979c3.png由式(18)可知,若式(19)成立,则省略c96c50eba9239b72ed6f44d6dfdb9acb.png就是比较合理的。545640bb14b8dcce3984147915b55d27.png那么,式(19)意味着什么呢?很容易发现,其意味着:第一,误差项的样本均值近似为零;第二,误差项与任何一个解释变量均近似地样本不相关。现在的问题是,上述两个结论成立吗?答案是,若“误差项期望值等于零”与“误差项与任意解释变量均不相关”这两大假定成立,则上述两个结论至少在大样本下是成立的。原因在于,这两大假定是两个总体矩条件,上述两个结论其实是相应的样本矩条件,而根据矩估计原理,样本矩是对总体矩的一致估计。根据上述讨论,我们可以获得两大启示:第一,给定上述两个假定成立,OLS估计本质上是矩估计的特例;第二,如果上述两个假定不成立,那么OLS估计量0fe06ac8e1fb620f940dfbdeeab82ecc.png就不会是对真实参数向量9e0b42389c83d97fb293cc6b830c672a.png比较“靠谱”的估计。反过来这意味着,上述两大假定成立,对于保证OLS估计量具有良好性质至关重要。四、如何保证2b64e03882e60d57155809e047da47f8.png存在?2b64e03882e60d57155809e047da47f8.png存在,表明34cdb16a499fd3b4bc7dd0f8b3f6871a.png是一个的满秩方阵,亦即e858b413fa1be01ca68a420824aba4cf.png×e858b413fa1be01ca68a420824aba4cf.png245b423f7b8adf5854b7c9765a74fb60.png。按照矩阵理论有:4405819c099c57a8abbcd164e44a98d0.png,故这进一步意味着,559909f25a2127021e9fc737300cf0d0.png作为一个的矩阵,秩等于k+1,亦即559909f25a2127021e9fc737300cf0d0.png必须列满秩。559909f25a2127021e9fc737300cf0d0.png满足列满秩假定,意味着构成559909f25a2127021e9fc737300cf0d0.png的k+1个列向量线性无关——这k+1个列向量中的任何一个向量,均不能是其余列向量的线性组合。若此假定被违背,则出现完全共线情况,此时91b4bc8bc2f4e9eb5c79198525835970.png不存在,OLS法失效。在此我们列举一个559909f25a2127021e9fc737300cf0d0.png不满足列满秩假定的例子。对于模型(20):727e1530f4ff7d2bf76f68b07d430ead.png假设a4c4ee8b20510960eeb235be57390479.png,则矩阵559909f25a2127021e9fc737300cf0d0.png中的第一个列向量是后两个列向量的线性组合,故三个列向量完全共线,559909f25a2127021e9fc737300cf0d0.png不具有列满秩性质。此时一个与原模型等价的新模型是:db1c5d851039f0ae1842e26155a3d386.png在这里,f3587e4421697a5b1457901b70f9f021.png为任意常数。1bdfe120c3f3c601cb6705315b20c6cd.png。现在我们不妨问这样一个问题,如果真能够将fcf0bb2e0b6849aca0bb4ebe7828c488.png710d7695cfcd1ae4874390abee24492e.png2b5d763ae63ed3425082c1f37da9ffbd.png进行回归,那么回归结果所估计的到底是式(20)还是式(21)呢?显然,我们无法确定。用计量经济学术语来讲,就是当559909f25a2127021e9fc737300cf0d0.png不满足列满秩假定时,模型(20)或者(21)是无法被识别的。值得指出的是,559909f25a2127021e9fc737300cf0d0.png不满足列满秩假定的一个特殊例子是,样本容量小于待估计参数的数量。例如,对于模型(20),其有三个参数需要估计。假定我们仅有两个观测值,那么559909f25a2127021e9fc737300cf0d0.png将是一个2×3的矩阵,其秩最大为2,故559909f25a2127021e9fc737300cf0d0.png不满足列满秩假定。其实,从直觉上很容易理解,模型(20)的样本回归方程代表一个平面,而要确定一个平面,至少需要3个点(观测值)。五、回到一元线性回归模型对于一元回归模型(22):529c1116d7c335bcecd8a22195bdd0f1.png此时,矩阵559909f25a2127021e9fc737300cf0d0.png由列向量3112a1630bfff6d9677905654487d6b3.pnga168f55319988059f2baf709a82a4c35.png构成,矩559909f25a2127021e9fc737300cf0d0.png列满秩假定成立表明:ee213068cb33cf0813dac77c8eead5c3.png,其中14bd50822383321a8a4c0b4bfb5cfc2c.png为任意常数。亦即,变量26e22925d2d9a98bd68104b2bd5aa49c.png的N次观测值不能为一个常数。对于一元线性回归模型,斜率估计量的公式为:e8186394d17808bfc6564268cd295dc3.png显然,若变量26e22925d2d9a98bd68104b2bd5aa49c.png的N次观测值为一个常数,则77a76b1ede1bd6747ed03138d10f3b53.png,而这是一个不定型。我们从直觉上很容易理解,当变量26e22925d2d9a98bd68104b2bd5aa49c.png的N次观测值为一个常数时,由于缺乏对照,变量26e22925d2d9a98bd68104b2bd5aa49c.pngc2ef52097c49e4a1af2670169937ffd9.png的影响是根本无法被识别的。

·END·

ed1a15806f4c70b960e2167400ecafec.png

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/561608.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

asr语音转写_搜狗智能录音笔C1正式上市 语音转文字准确率达95%

中新网3月19日电 3月18日,搜狗旗下全新智能硬件——搜狗智能录音笔C1在搜狗商城、京东等平台正式发售。在当天的产品沟通会上,搜狗公司CTO杨洪涛表示,搜狗智能录音笔C1集成了搜狗的核心AI技术,能够有效提高用户的工作和生活效率。…

730阵列卡支持多大硬盘_华为1000多手机哪款好?推荐只此一款!华为品牌性价比最高千元机...

自荣耀品牌独立运营之后,华为品牌就主攻高端机型,其低端机型均交给荣耀品牌来运营。其实,华为品牌中端机型中,还是有几款性价比表现极为不俗的,比如中端性价比最高机型——华为nova 5Z,华为nova 5Z性价比直…

jsp 中div居中_让div在屏幕中居中(水平居中+垂直居中)

html代码<div class"main"><h1>测试</h1></div>方法一&#xff1a;div使用绝对布局&#xff0c;设置margin:auto,并设置top、left、right、bottom的值相等即可&#xff0c;不一定要都是0。方法二&#xff1a;绝对布局&#xff0c;用margin来…

使用eolink优雅地进行API接口管理

为什么使用eolink? 我们都知道在一个项目团队中是由很多角色组成的&#xff0c;例如&#xff1a;业务>产品>设计>前端>后端>测试等。每个角色各司其职&#xff0c;一起合作完成项目的生命周期。而前端与后端的沟通则是主要通过接口来实现&#xff0c;通过接口…

没有icon_ICON设计干货来啦~

​ICON设计是UI设计中非常重要的一个环节&#xff0c;它能够帮助我们美化画面&#xff0c;传递情绪&#xff0c;快速导航&#xff0c;而且不受语言限制&#xff0c;那今天就来给大家分享下ICON的设计干货。首先图标可以分为3大类&#xff1a;功能性图标 , 应用图标 , 创意图标。…

基于JavaSpringBoot+Vue+uniapp微信小程序实现校园宿舍管理系统

博主介绍&#xff1a;✌全网粉丝20W,csdn特邀作者、博客专家、CSDN新星计划导师、java领域优质创作者,博客之星、掘金/华为云/阿里云/InfoQ等平台优质作者、专注于Java技术领域和毕业项目实战✌ &#x1f345;文末获取联系&#x1f345;精彩专栏推荐订阅&#x1f447;&#x1f…

实战 es6_腾讯云 Elasticsearch 实战篇(二十二) ES6.8权限使用配置

前言|在前面的章节中我们讲了开源架构ELK、腾讯云Elasticsearch产品系列。我们也知道了&#xff0c;在构建腾讯云ES集群的时候&#xff0c;我们选择的6.8.2白金版具备充分的安全的机制来保证数据和访问的安全。那么&#xff0c;它到底是如何实现的呢&#xff1f;我们今天就来简…

基于Java+Springboot+Vue+elememt社区疫情返乡管控系统设计实现

博主介绍&#xff1a;✌全网粉丝20W,csdn特邀作者、博客专家、CSDN新星计划导师、java领域优质创作者,博客之星、掘金/华为云/阿里云/InfoQ等平台优质作者、专注于Java技术领域和毕业项目实战✌ &#x1f345;文末获取联系&#x1f345;精彩专栏推荐订阅收藏&#x1f447;&…

策略 python_Python版商品期货多品种均线策略

完全移植自「CTP商品期货多品种均线策略」&#xff0c;由于Python版本商品期货策略还没有一个多品种的策略&#xff0c;所以就移植了JavaScript版本的「CTP商品期货多品种均线策略」。提供一些Python商品期货多品种策略的设计思路、例子。不论JavaScript版本还是Python版本&…

mc用云服务器搭建_最全的云服务器架设我的世界私服教程,不看后悔哦!

我们将通过一个一步步的讲述通过新睿云的云服务器来向你展示如何搭建一个“我的世界(Minecraft)”服务器端。这是一个可以多人在线的游戏&#xff0c;能让您体会到极致的游戏乐趣。我们的建议是不要使用个人电脑&#xff0c;即使从技术角度来说你能做到&#xff0c;但个人电脑不…

基于Java+SpringBoot+vue等疫情期间在线网课管理系统详细设计实现

博主介绍&#xff1a;✌全网粉丝20W,csdn特邀作者、博客专家、CSDN新星计划导师、java领域优质创作者,博客之星、掘金/华为云/阿里云/InfoQ等平台优质作者、专注于Java技术领域和毕业项目实战✌ &#x1f345;文末获取联系&#x1f345;精彩专栏推荐订阅收藏&#x1f447;&…

智能家居中语音识别算法研究_语音识别研究获进展

中国科学院自动化研究所智能交互团队在环境鲁棒性、轻量级建模、自适应能力以及端到端处理等几个方面进行持续攻关&#xff0c;在语音识别方面获新进展&#xff0c;相关成果将在全球语音学术会议INTERSPEECH2019发表。现有端到端语音识别系统难以有效利用外部文本语料中的语言学…

自动驾驶技术越来越火,浅谈一些对百度Apollo开放平台8.0的看法和认知

最近几年来&#xff0c;包括我在内的很多开发者们&#xff0c;都投身进入了智能汽车领域&#xff0c;其实不论是初学者还是大佬&#xff0c;接触智能汽车开发行业后都需要学习很多知识和既能。 我在最近一段时间接触到了百度的Apollo开放平台&#xff0c;上手的感觉很好用&…

命令行输入mysql不行_MySQL命令行无法插入中文数据

下面的方法是针对window系统的注意&#xff1a;1 下面的my.ini文件在有些版本的Mysql中的安装目录找不到&#xff0c;比如我用的版本5.5.60-log就找不到&#xff0c;想要找到这个文件&#xff0c;需要设置显示window系统隐藏文件(不懂的自行百度&#xff0c;因系统的版本不同而…

mysql悲观锁关键字_MySQL悲观锁 select for update实现秒杀案例(jfinal框架)

MySQL悲观锁 select for update实现秒杀案例(jfinal框架)发布时间&#xff1a;2018-08-17作者&#xff1a;laosun阅读(4287)为了方便测试&#xff0c;博主使用最新的jfinal框架&#xff0c;里边的东西就修改了一下a_little_config.txt&#xff0c;配置数据库链接的&#xff0c;…

pl sql 连接mysql_PL/SQL 连接mysql步骤

下面就将PL/SQL的配置说明一下。一、安装Oracle客户端&#xff0c;让后配置 安装目录下面的C:\ORACLE\instantclient_11_2\NETWORK\ADMIN 的 tnsnames.ora文件&#xff0c;如下&#xff1a;修改成自己的配置名 (DESCRIPTION (ADDRESS (PROTOCOL TCP)(HOST 服务器IP)(PORT…

mysql 创建表时提示错误代码_MySQL创建表时遇到的错误

学习MySQL第一个程序就遇到了bug。创建表时需注意的事项加以总结&#xff1a; 1、先创建数据库&#xff1a; create database 数据库名&#xff1b;然后连接要用的数据库 use 数据库名&#xff1b;下面就可以执行建表语句 了。2、错误提示&#xff1a;ERROR 1064 (4…

根据工序画出aoe网_这些金刚网纱窗竟然含“毒”!选错就得病

窗户进行装修我们通常会安装纱窗&#xff0c;这样可以阻止一些蚊虫或者灰尘进入家里。金刚网纱窗是业主们常常挑选的一种纱窗。不过最近听别人说金刚网纱窗是含毒的&#xff0c;这不禁让安装这种纱窗的业主感到吃惊&#xff0c;而这消息是真是假呢我们看看专家的说法&#xff0…

防火墙阻止网页连接网络连接_win7电脑打不开网页怎么办 win7电脑打不开网页解决步骤【图文步骤】...

我们在使用电脑的过程当中&#xff0c;经常都会浏览一些网页&#xff0c;不过最近有位win7系统用户在使用电脑浏览网页的时候却遇到了网页打不开网页的问题&#xff0c;这导致了很多事情都无法操作&#xff0c;其实网页打不开的问题&#xff0c;一般都跟网络有关&#xff0c;那…

mysql mysqldumpslow_慢日志分析工具—mysqldumpslow 和 mysqlsla

前提&#xff1a;分析mysql性能的时候会查看数据库的哪些sql语句有问题&#xff0c;效率低。这就用到了数据库的慢查询&#xff0c;作用就是&#xff1a; 它能记录下所有执行超过long_query_time时间的SQL语句&#xff0c;帮你找到执行慢的SQL&#xff0c;方便我们对这些SQL进行…