如何使用matlab对时间序列进行ADF检验?|adftest函数(获取不同显著性下的统计结果)

news2024/12/25 9:02:07

ADF检验

迪基富勒检验(ADF检验)是一种常见的统计检验,用于检验给定时间序列是否平稳。在分析序列的平稳性时,它是最常用的统计检验之一。matlab中提供了函数adftest可以完成该检验,本文重点介绍该函数的用法。

Matlab|adftest

[h,pValue,stat,cValue,reg] = adftest(y)
[h,pValue,stat,cValue,reg] = adftest(y,param1,val1,param2,val2,...)

输入参数介绍
y - 时间序列数据的向量。最后一个元素是最近的观察结果。表示缺失值的nan会被删除。
lags- 可以是非负整数的标量或者向量,延时,默认为0。
model- 模型:‘AR’,‘ARD’,或’TS’,默认AR模型。
test- 检验方法,t1表示标准的t检验,t2表示非标准的t检验,'F’表示F检验。默认是t1检验。
alpha- 显著性水平,取值为0.001到0.999。默认为0.05。

输出参数介绍
h:
如果结果h=0,表示拒绝原假设,数据不平稳
如果结果h=1,表示不拒绝原假设,数据平稳
pValue:
检验统计量的p值
stat:
统计值
reg: 回归统计的系数

%   reg - Structure of regression statistics from the OLS estimation of 
%       coefficients in the alternative model. The number of records is
%       equal to the number of tests. Each record has the following fields:
%
%       num         Length of the input series y, with NaNs removed
%       size        Effective sample size, adjusted for lags, difference*
%       names       Regression coefficient names			
%       coeff       Estimated coefficient values
%       se          Estimated coefficient standard errors
%       Cov         Estimated coefficient covariance matrix
%       tStats      t statistics of coefficients and p-values
%       FStat       F statistic and p-value
%       yMu         Mean of y, adjusted for lags, difference*
%       ySigma      Standard deviation of y, adjusted for lags, difference*
%       yHat        Fitted values of y, adjusted for lags, difference*
%       res         Regression residuals
%       DWStat      Durbin-Watson statistic
%       SSR         Regression sum of squares
%       SSE         Error sum of squares
%       SST         Total sum of squares
%       MSE         Mean squared error
%       RMSE        Standard error of the regression
%       RSq         R^2 statistic
%       aRSq        Adjusted R^2 statistic
%       LL          Loglikelihood of data under Gaussian innovations
%       AIC         Akaike information criterion
%       BIC         Bayesian (Schwarz) information criterion
%       HQC         Hannan-Quinn information criterion

实战

通常使用
1.使用默认参数的结果

h = adftest(Y);  
disp(h);

2.如果想在不同显著性下获得统计结果,则需要对显著性参数进行调整。下图表示显著性为1%,%5,10%情况下的统计结果

%% ADF检验
% 如果结果h=1,表示不拒绝原假设,数据平稳
disp('统计结果**********************')
[h,pValue,stat,cValue,reg] = adftest(data,'alpha',[0.01,0.05,0.1]);
fprintf('显著性|%6.6f\n',h);
fprintf('检验统计量|%6.6f\n',stat);
fprintf('p-value|%6.6f\n',pValue);
fprintf(' 1%%|%6.6f\n',cValue(1));
fprintf(' 5%%|%6.6f\n',cValue(2));
fprintf('10%%|%6.6f\n',cValue(3));

在这里插入图片描述

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

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

相关文章

创建表和管理表

文章目录基础知识一条数据存储的过程标识符命名规则MySQL中的数据类型创建和管理数据库创建数据库使用数据库修改数据库删除数据库创建表创建方式1创建方式2查看数据表结构修改表追加一个列修改一个列重命名一个列删除一个列重命名表删除表清空表内容拓展拓展1:阿里…

setState的使用+React更新机制+events+受控和非受控组件

setState是异步更新 总结: 1.setState设计为异步,可以显著的提升性能 如果每次调用 setState都进行一次更新,那么意味着render函数会被频繁调用,界面重新染,这样效率是很低的;最好的办法应该是获取到多个更…

单片机开发---ESP32S3移植lvgl+触摸屏

书接上文 《单片机开发—ESP32-S3模块上手》 本章内容 熟悉一下ESP32S3的开发,修改范例程序的lvgl,使之能够匹配现在的显示屏。 具体工作大概为通过SPI接口连接一块SPI串口屏幕,并且适配lvgl,最后加上触摸屏作为输入。 屏幕 …

【计算机网络】第一章 计算机网络结构

文章目录第一章 体系结构1.1 计算机网络概述1.1.1 计算机网络的概念1.1.2 计算机网络的组成1.1.3 计算机网络的功能1.1.4 计算机网络的分类*1.1.5 计算机网络的标准化工作1.1.6 计算机网络的性能指标1.2 计算机网络体系结构与参考模型1.2.1 计算机网络分层结构1.2.2 计算机网络…

#8链表的中间结点#

链表的中间结点 1题目链接 链接 2思路 思路1:遍历一遍 计数 然后/2 再遍历一遍 思路2:slow fast指针 slow指针一次走1步 fast指针一次走2步 当fast为空的时候 slow的位置就是中间结点 奇数个: 1 2 3 4 5 fast走完第三次为空 slow走完第三次就是3 偶数个: 1 2 3 4 5 6 fast走完…

智慧型物业管理系统功能解析

随着当前社会经济的发展与科技发达,物业管理系统化已经成为常态了。尤其是随着智慧物业管理系统功能越来越多,人们对智慧物业管理系统的依赖就更明显了。毕竟系统真的可以给生活带来很多的便利之处: 业主可通过该系统查询自己住房的详细信息…

“揾”钱,最紧要系稳

我是腾讯安全的樊自磊。我们团队在腾讯主要负责金融风控产品,解决相关产品交付和服务维护工作,像国内知名大型国有银行、城商行、互联网金融公司等,都是我们的服务对象。今年春节,我和我的的同事们都在深圳为金融行业的网络安全进…

连续多输入多输出对象最优控制

连续多输入多输出对象最优控制 控制对象:平面二自由度机械臂 动力学模型: M ( q ) q + C ( q , q ) + G ( q ) =

【HDRP】自动生成的光照探针——Probe Volume

HDRP中,增加了Probe Volume,可代替旧版的光照探针Light Probe Group。 使用此功能的物体,不再需要光照贴图。 一、优缺点比较 详细说明可查看官方说明。 1.Probe Volume按像素而不是按对象发光,这意味着 HDRP 可以更准确地照亮…

软件著作权申请材料

(一)按要求填写的软件申请表; (二)软件的鉴别材料; 1、软件的操作手册:图文并茂的详细介绍软件的各功能,文档应不超过60页,超过60页应当删除中间内容,保留前后30页; 2、软件的源代码:每页不少于50行&…

Go第 15 章 :单元测试

Go第 15 章 :单元测试 15.1 先看一个需求 在我们工作中,我们会遇到这样的情况,就是去确认一个函数,或者一个模块的结果是否正确, 如: 15.2 传统的方法 15.2.1 传统的方式来进行测试 在 main 函数中&a…

k8s集群调度、亲和性、污点和容忍、pod状态、排障步骤

目录 一、调度约束 二、创建资源工作流程 三、Scheduler调度过程 1.Scheduler调度中考虑的问题 2.调度过程的步骤 3.预算策略(predicate)的常见算法 4.优选策略(priorities)的优先级 四、Pod 调度到指定Node节点 1.nodeNa…

微服务的服务拆分与远程调用

​ 哈喽大家好呀!好久不见甚是想念,给大家拜个年啦~应该不晚吧(ಥ_ಥ) 放假在家确实是容易躺平,有心而无力呀哈哈哈哈。但是闲着也是闲着,最近学了学微服务相关知识,马上也快毕业了就更到抓紧了 今天我来说说关于微服…

医疗影像容积重建-物体坐标系和渲染图像坐标系

1 物体坐标系(mm为单位,并三方向都是各向同性的情况)将所有数据的正中间(这里所有的数据,是指各体素按照实际位置在空间排列),作为物体坐标系的原点,以图像的image orientation X方向&#xff0…

Android开机动画

1.动画的位置 android开机动画位置在/system/media目录下,/system/media/bootanimation.zip 解压是如下的文件内容 part0,part1,part2存放的都是要播放开机动画的帧图片,desc.txt写的是播放方式,desc文件内容类似如下…

天云数据Hubble数据库荣获2022年度IT168技术卓越奖

在2022年即将结束之际,IT168再次启动“技术卓越奖”评选。由行业CIO/CTO大咖、技术专家及IT媒体三方联合评选,评判标准代表了用户和媒体声音。经过多方评审,天云数据Hubble数据库荣获2022年度IT168技术卓越奖。 IT168主编任朝阳说&#xff1a…

SAP FICO 内部订单详细解析

内部订单 内部订单是用来对企业内部某项工作或者任务编制计划、归集成本、结算的载体。比如市场推广活动、内部团队活动、研发项目、投资项目、在建工程项目等。 内部订单需要区别于销售订单、采购订单和生产订单。 销售订单和采购订单是企业与外部单位以合同或者契约为纽带&…

Socket缓冲区,可读可写条件

一 socket缓冲区 每个 socket 被创建后,都会分配两个缓冲区,输入缓冲区和输出缓冲区。 write()/send() 并不立即向网络中传输数据,而是先将数据写入缓冲区中,再由TCP协议将数据从缓冲区发送到目标机器。一旦将数据写入到缓冲区&…

从0到0.1学习 maven(一:概述及简单入门)

文章目录概述从没用的有趣小知识开始Maven是什么?为什么用Maven?使用与配置入门目录结构编写pom使用Archetype生成骨架小彩蛋:直接运行包含main的jar概述 从没用的有趣小知识开始 Maven可以翻译成“知识的累积”,“专家/行家”。…

数据质量管理深入浅出

质量是生活中最常关注的话题,我们都期望享用高质量的商品与服务,且企业也不断加大质量管理的投入,为了更好的用户体验。 在企业数字化转型浪潮下,传统手段已无法应对数字化转型中的数据质量管理需求,我们需要探索出一条…