时间序列研(part10)--误差修正模型

article/2025/9/26 14:23:12

学习笔记,仅供参考,有错必纠


文章目录

    • 时间序列
      • 误差修正模型
        • F检验
        • 似然比(LR)检验
        • W检验
        • LM乘数检验
        • LR, W和LM检验
        • 自相关的LM检验



时间序列


误差修正模型


在用“一般到特殊”方法建立模型时的,首先应对初始模型(即对回归参数不加任何约束的动态分布滞后模型)的随机误差项进行异方差和自相关检验。对模型的其他检验都应建立在随机误差项是一个白噪声序列的基础之上。在检验约束条件是否成立的过程中逐步剔除不显著变量,化简模型,同时还要保持模型随机误差项的非自相关性和同方差性不被破坏。在这个过程中要用到许多统计量.

F检验

把样本数据取对数后建立回归模型,随机误差项一般不会存在异方差。对于随机误差项的一阶自相关检验可用DW统计量完成。对于ADL(1,1)模型,约束条件 α 1 = 0 , β 1 = 0 , β 0 = 0 \alpha_1 = 0, \beta_1 = 0, \beta_0 = 0 α1=0,β1=0,β0=0 α 1 + β 0 + β 1 − 1 = 0 \alpha_1 + \beta_0 + \beta_1 - 1 = 0 α1+β0+β11=0是否成立可用t检验完成。如果t统计量的绝对值大于临界值,则相应约束条件不成立,相应解释变量不能轻易地从模型中剔除掉。否则接受相应约束条件,从模型中剔除相应解释变量。

对于联合线性约束条件可用F检验完成。假定模型误差项服从正态分布,共有m个线性约束条件,则所用统计量是 :

在这里插入图片描述

在零假设“约束条件真实”条件下, F ∼ F ( m , T − k ) F \sim F(m, T-k) FF(m,Tk).

因为两个模型都是用OLS法估计的,所以可把被解释变量的总平方和(SST)分解为回归平方和 (SSR) 与误差平方和(SSE)两部分。

对于不加约束的模型有:

SST = SSRu + SSEu .

对于施加约束条件的模型有:

SST = SSRr + SSEr .

如果约束条件成立,那么在施加约束条件下求到的SSEr 不会比不加约束条件的SSEu大很多,用样本计算的F值不会很大。

若F值小于临界值,则约束条件是可接受的(真实的)。否则应该拒绝零假设。

注意,F检验的零假设是m个约束条件同时为零,备择假设是m个约束条件不同时为零。所以拒绝零假设并不排除有部分约束条件为零。应利用t检验进一步对每一个参数进行显著性判别。

比如对ADL(1,1)模型,检验联合约束条件 α 1 = β 1 = 0 \alpha_1 = \beta_1 = 0 α1=β1=0,则:

在这里插入图片描述


似然比(LR)检验

以上介绍的t检验和F检验只适用于对线性约束条件的检验。对于非线性约束条件, α 1 β 0 + β 1 = 0 \alpha_1 \beta_0 + \beta_1 = 0 α1β0+β1=0,则无法用t或F检验完成。下面介绍三种常用的检验方法,即:

似然比(LR)检验;

沃尔德(W)检验;

拉格朗日(lagrange)乘数(LM)检验。

这三种检验所用统计量都是利用极大似然估计法计算的。LR检验由内曼—皮尔逊(Neyman-Pearson 1928)提出,只适用于对线性约束的检验。W检验和LM检验既适用于对线性约束条件的检验,也适用于对非线性约束条件的检验。

LR检验的基本思路是如果约束条件成立则相应约束模型与非约束模型的极大似然函数值应该是近似相等的。

下式表示非约束模型的极大似然函数:

在这里插入图片描述

其中 β ^ \hat{\beta} β^ σ ^ 2 \hat{\sigma}^2 σ^2分别是对 β \beta β(参数集合), σ 2 \sigma^2 σ2的极大似然估计。

下式表示约束模型的极大似然函数:

在这里插入图片描述

其中 β ~ \tilde{\beta} β~ σ ~ 2 \tilde{\sigma}^2 σ~2分别是对 β \beta β(参数集合), σ 2 \sigma^2 σ2的极大似然估计。

定义似然比(LR)统计量为:

在这里插入图片描述

中括号内是两个似然函数之比(似然比检验由此而得名)。在零假设约束条件成立条件下:

在这里插入图片描述

其中m表示约束条件个数。用样本计算LR统计量。判别规则是:

L R < χ α 2 ( m ) LR < \chi^2_{\alpha}(m) LR<χα2(m) , 则接受零假设,约束条件成立。

L R > χ α 2 ( m ) LR > \chi^2_{\alpha}(m) LR>χα2(m) , 则拒绝零假设,约束条件不成立。

ADL(1,1)模型,检验联合约束条件 α 1 = β 1 = 0 \alpha_1 = \beta_1 = 0 α1=β1=0,则:

在这里插入图片描述

LR统计量只适用于对线性约束条件的检验。对非线性约束条件应该采用如下两种检验方法。

W检验

W检验的优点是只需估计无约束模型。当约束模型的估计很困难时,此方法尤其适用。W检验由沃尔德(Wald 1943)提出,适用于线性与非线性约束条件的检验。

W检验的原理是测量无约束估计量与约束估计量之间的距离。先举一个简单例子。比如对模型:

在这里插入图片描述

检验线性约束条件 β 2 = β 3 \beta_2 = \beta_3 β2=β3是否成立。W检验只需对无约束模型进行估计,因为对约束估计量 β ~ 2 \tilde{\beta}_2 β~2 β ~ 3 \tilde{\beta}_3 β~3来说,必然有 β ~ 2 − β ~ 3 = 0 \tilde{\beta}_2 - \tilde{\beta}_3 = 0 β~2β~3=0。如果约束条件成立,则无约束估计量 β ^ 2 − β ^ 3 \hat{\beta}_2 - \hat{\beta}_3 β^2β^3应该近似为零。如果约束条件不成立,则无约束估计量应该显著地不为零。关键是要找到一个准则,从而判断什么是显著地不为零。

在这里插入图片描述

其中:

在这里插入图片描述

∂ f ( β ^ ) ∂ β ^ \frac{\partial f(\hat{\beta})}{\partial \hat{\beta}} β^f(β^)表示用无约束估计量 代替后的偏导数矩阵,其中第i行第j列位置上的元素表示第i个约束条对第 j个无约束估计量的偏导数值. V a r ( β ^ ) Var(\hat{\beta}) Var(β^) β ^ ) \hat{\beta}) β^)的估计的方差协方差矩阵.

在约束条件成立条件下:

在这里插入图片描述

W统计量的具体表达式为:

在这里插入图片描述

在零假设 β 1 β 2 = β 3 \beta_1 \beta_2 = \beta_3 β1β2=β3成立条件下,W统计量近似服从 χ 2 ( 1 ) \chi^2(1) χ2(1)分布.

LM乘数检验

与W检验不同的是拉格朗日(Lagrange)乘数(LM)检验只需估计约束模型。所以当施加约束条件后模型形式变得简单时,更适用于这种检验。LM检验是由艾奇逊—西尔维(Aitchison-Silvey 1960)提出的。LM检验另一种表达式是由拉奥(Rao 1948)提出的,称为得分检验。

首先给出非约束模型的对数似然函数:

在这里插入图片描述

对于非约束极大似然估计量 β j ^ \hat{\beta_j} βj^必然有:

在这里插入图片描述

在这里插入图片描述

假定有两个约束条件f1(b)=0和f2(b)=0。为求这两个约束条件下的对数似然函数的极大似然估计量,应按拉格朗日乘数法则建立如下函数:

在这里插入图片描述

在这里插入图片描述

对于线性回归模型,通常是通过一个辅助回归式计算LM统计量的值。LM统计量与辅助回归式的可决系数 R 2 R^2 R2有直接联系,而辅助回归式的形式直接与被检验的约束条件有关。

LM检验的实际步骤如下:

在这里插入图片描述

在这里插入图片描述


  • 例子

下面介绍用LM辅助回归方法检验约束条件 β 2 + β 3 = 1 \beta_2 + \beta_3 = 1 β2+β3=1

在这里插入图片描述

LM检验的实际步骤:

在这里插入图片描述

  • 实例(湾湾制造业生产函数)

在这里插入图片描述

(1)用OLS法估计约束模型,计算残差序列 u t ^ \hat{u_t} ut^

在这里插入图片描述

(2)确定LM辅助回归式的解释变量。

(3)建立LM辅助回归式如下:

在这里插入图片描述

(4)用OLS法估计上式并计算可决系数 R 2 R^2 R2

在这里插入图片描述

(5)用第四步得到的 R 2 R^2 R2计算LM统计量的值:

在这里插入图片描述


LR, W和LM检验


对LR,W和LM检验方法的选择应以做实际计算时的难易程度而定。一般来说W和LM检验应优于LR检验,因为W和LM检验只需要估计一个模型即可,而LR检验需估计约束与非约束两个模型。对W 和LM检验方法的选择应以约束模型与非约束模型哪个更容易估计而定。应该注意,即使三种检验方法都可使用,它们的计算结果通常也是不相同的。因为三个统计量只是渐近相同,对于线性回归模型,在小样本条件下有如下关系成立。:
L M ≤ L R ≤ W LM \le LR \le W LMLRW
上式说明只有当 LM检验的结果为拒绝零假设(约束条件不成立)或者W检验的结果为接受零假设(约束条件成立)时,三种检验的结论才是一致的。

实际中,三种检验方法有可能得出相互不一致的结论。另外只有当用参数的样本估计值计算的约束条件完全成立时,即把参数估计值代入约束条件能准确成立时,式中的三个统计量才有完全相等的关系。

当对数似然函数中只含有一个参数 β \beta β时,LM, LR 和W三种检验的关系可用图表示:

在这里插入图片描述

自相关的LM检验

DW统计量只适用于一阶自相关检验,而对于高阶自相关检验并不适用。利用LM统计量可建立一个适用性更强的自相关检验方法,既可检验一阶自相关,也可检验高阶自相关。

考虑两种误差过程的模型:

在这里插入图片描述

在这里插入图片描述


http://chatgpt.dhexx.cn/article/Xm3fCdgH.shtml

相关文章

pyspider框架之Tripadvisor酒店数据爬取

需求 网站入口&#xff1a;www.tripadvisor.com 网页下端&#xff0c;遍历点开进入所有城市链接&#xff1a; 点击后进入该城市的所有hotel 代码 #!/usr/bin/env python # -*- encoding: utf-8 -*- # Created on 2018-09-06 11:16:59 # Project: trip_hotelfrom pys…

爬取携程和蚂蜂窝的景点评论数据\携程评论数据爬取\旅游网站数据爬取

本人长期出售超大量微博数据、旅游网站评论数据,并提供各种指定数据爬取服务,Message to YuboonaZhang@Yahoo.com。同时欢迎加入社交媒体数据交流群:99918768 前言 为了获取多源数据需要到各个网站获取一些景点的评论信息和图片,首先选到了携程和蚂蜂窝这两个网站,将一些…

Python爬虫实战--TripAdvisor爬虫

目标站点分析 目标URL&#xff1a;https://www.tripadvisor.cn/Attractions-g60763-Activities-New_York_City_New_York.html#FILTERED_LIST 明确内容&#xff1a; 在图中&#xff0c;我们明确爬去的内容为&#xff1a;title&#xff0c;comment, cate&#xff08;即绿色框中的…

如何用Python从TripAdvisor抓取数十万条酒店评论

我从TripAdvisor抓取一些酒店评论&#xff0c;然后发现了一种从它们那里刮掉数十万条酒店评论的好方法。 让我们假设&#xff0c;例如&#xff0c;我们要从大加那利岛刮掉酒店评论。如果转到TripAdvisor&#xff0c;我们将看到URL为&#xff1a; https://www.tripadvisor.com…

美通企业日报 | 猫途鹰联手携程打造中国顶级旅行平台;强生战略合作阿里旗下Lazada...

今日看点 TripAdvisor与携程集团联手打造面向中国出境旅行者的顶级旅行平台。全球领先的旅游平台TripAdvisor&#xff08;猫途鹰&#xff09;宣布其中国子公司已与携程集团达成战略合作&#xff0c;以打造中国顶级的旅行计划和预订网站&#xff0c;为渴望探索世界的中国旅行者提…

Python爬取TripAdvisor

直接上代码&#xff1a; #爬取tripadvisor纽约市酒店超值排名#引入requests 获取html文件&#xff0c;才能从html获取信息 import requests #利用BeautifulSoup解析文件&#xff0c;获取想要的到的数据 from bs4 import BeautifulSoup #这段代码只用在获取等待&#xff0c;避免…

猫途鹰公布2023年全球十大最佳旅行体验和十大顶级景点 | 美通社头条

美通社消息&#xff0c;旅游指南平台猫途鹰(Tripadvisor)公布2023年旅行者之选&#xff1a;最佳“必做之事”。 随着夏季旅游的全面展开&#xff0c;这些是猫途鹰评论家们最喜欢的来自世界各地的活动&#xff0c;为希望创造难忘时刻的旅行者提供一份明确的非凡体验清单&#xf…

猫途鹰公布2019年“旅行者之选”全球最佳海滩榜单

全球旅游规划和预订平台猫途鹰(TripAdvisor)公布2019年“旅行者之选”最佳海滩榜单。获奖海滩是基于过去12个月内全球上亿旅行者的评分和点评的数量及质量综合计算得出&#xff0c;巴西费尔南多迪诺罗尼亚群岛&#xff08;Fernando de Noronha&#xff09;的桑乔湾海滩&#xf…

【Python】代码:获取猫途鹰的London酒店信息:基于Scrapy框架和requests库

本文以代码分析的形式记录&#xff1a;利用Scrapy框架和requests库爬取tripadvisor(猫途鹰)多个城市的酒店信息&#xff0c;数据量300w条(1.09G)&#xff0c;运行时间约7h。多个城市与单个城市的操作类似&#xff0c;为避免代码过于冗长&#xff0c;本文仅以爬取London酒店的评…

利用 pyspider 框架抓取猫途鹰酒店信息

利用框架 pyspider 能实现快速抓取网页信息&#xff0c;而且代码简洁&#xff0c;抓取速度也不错。 环境&#xff1a;macOS&#xff1b;Python 版本&#xff1a;Python3。 1.首先&#xff0c;安装 pyspider 框架&#xff0c;使用pip3一键安装&#xff1a; pip3 pyspider 2.终端…

可怕的pyspider猫途鹰

1.启动pyspider 2.新建一个项目 3.代码 4. 注意事项&#xff1a;网址什么的都变了 5.存储到MongoDB&#xff0c; 注意这个地方我错了三次 6.在tableau可视化才发现错误的1,2 之后就能可视化了&#xff0c;本次实验是个半成品。后期会补充。 #!/usr/bin/env python # -*- enco…

爬虫-猫途鹰

from bs4 import BeautifulSoup import requests url https://www.tripadvisor.cn/ wb_data requests.get(url) soup BeautifulSoup(wb_data.text,lxml) for i in soup.select(li):if len(i.select(.ranking))>0:sorti.select(.ranking)[0].text #排名countryi.select(.c…

JS DOM 编程复习笔记--父元素、子元素和兄弟元素(三)

今天我们来复习DOM中的获取父元素、子元素和兄弟元素的API&#xff0c;它们主要有parentNode、firstChild、firstElementChild、lastChild、lastElementChild、childNodes、children、nextElementSibling、nextSibling、previousElementSibling、previousSibling等。 目录 获取…

jquery 在兄弟节点前、或兄弟节点后添加最新元素

使用 jquery 封装好的方法操作 dom&#xff0c;非常方便 1、在兄弟节点前添加最新元素 使用 before() 方法 演示代码如下 <!DOCTYPE html> <html><head><meta charset"UTF-8"><title></title><script src"js/jquery-…

CSS第一章:4.元素关系(兄弟、祖后代关系);关系选择器

总览 1.元素关系 2.关系选择器 一、元素关系 1.元素关系1&#xff1a;父子 父元素&#xff1a;直接包含子元素的元素叫做父元素 子元素&#xff1a;直接被父元素包含的元素叫做子元素 2.元素关系2&#xff1a;祖先后代 祖先元素&#xff1a;直接或间接包含后代元素的元素…

选择兄弟元素中的第几个元素

:nth-child(anb) 这个CSS伪类首先找到所有当前元素的兄弟元素&#xff0c;然后按照位置先后顺序从1开始排序&#xff0c;选择的结果为CSS伪类:nth-child括号中表达式&#xff08;anb&#xff09;匹配到的元素集合&#xff08;n0&#xff0c;1&#xff0c;2&#xff0c;3...&am…

兄弟元素选择器

兄弟元素选择器 语法1&#xff1a;前一个元素 后一个元素作用&#xff1a;选中一个元素后紧挨着的指定的兄弟元素。语法2&#xff1a;前一个元素 ~ 后边一类元素作用&#xff1a;选中后边的所有兄弟元素 举例1&#xff1a; <!DOCTYPE html> <html lang"en&qu…

Thinking -- CSS从根解决选择前一个兄弟元素

Thinking系列&#xff0c;旨在利用10分钟的时间传达一种可落地的编程思想。 开发中遇到这样一个诉求&#xff1a;特定class的元素单独占一行&#xff0c;现需要针对其前一个兄弟元素增加相应标识&#xff0c;以使其占据所在行的剩余所有空间。 换句话&#xff1a;就是如何选中…

JS 兄弟元素和兄弟节点

获取兄弟节点和兄弟元素 例子&#xff1a; //body内容<ul id"ul">wwwwww<!-- 我是注释 --><li>我是li标签1</li><li id"li2">我是li标签2</li><li>我是li标签3</li><li>我是li标签4</li><…

css-知识点(学习笔记)

一、鼠标悬停样式 hover伪类基础用法 在前端学习的初期&#xff0c;想必大家用的很多的属性之一就有hover吧,hover伪类元素使用的三种情况 1、 hover用于父子元素 父子元素直接在父元素的id或者class或者标签名后紧接着:hover 然后空格子元素的名字就可以了。 2、 hover用于…