无季节效应的非平稳序列分析(一)

目录

Cramer分解定理(1961年提出)

差分

R语言函数 diff

例题:

过差分:

小结


Cramer分解定理(1961年提出)

任何一个时间序列 {x_t} 都可以分解为两部分的叠加:其中一部分是由多项式决定的确定性趋势成分,另一部分是平稳的零均值误差成分,即

Box和Jenkins用大量的案例分析证明了差分方法是一种非常简便、有效的确定性信息提取方法
Cramer分解定理在理论上保证了适当阶数的差分一定可以充分提取确定性信息
就有:

其中确定性部分看成常数c:

为什么确定性部分可以看成常数呢?看下面!!

举例如:

它的一阶差分为:

二阶差分为:

从二次差分以后的差分值都为2,所以确定性部分可以看成一个常数。

差分是什么哇?接下来我们就介绍一下差分

差分

差分运算的实质是使用自回归的方式提取确定性信息

 则可写成:

差分方式的选择:

  1. 序列蕴含显著的线性趋势,一阶差分就可以实现趋势平稳
  2. 序列蕴含着曲线趋势,通常低阶(二阶或三阶)差分就可以提取出曲线趋势的影响
  3. 对于蕴含着固定周期的序列进行步长为周期长度的差分运算,通常可以较好地提取周期信息

R语言函数 diff

1差分  diff(x)

2差分  diff(x,1,2)

k差分  diff(x,1,k)

d差分  diff(x,d,1) 或 diff(x,d)

一阶差分后再进行d步差分 diff(diff(x),d)

例题

例5.1 1964年-—1999年中国纱年产量序列蕴含着一个近似线性的递增趋势。对该序列进行一阶差分运算,考察差分运算对该序列线性趋势信息的提取作用
 

a<-read.table("D:/桌面/5_1.csv",sep=',',header=T)
a
x<-ts(a$output,start=1949) 
plot(x)  #plot(x,type='o') #绘制时序图
difx<-diff(x)  #1阶差分
plot(difx) #对1阶差分绘制时序图

返回:

如图,还具有趋势,需要二阶差分

difx2<-diff(x,1,2)
plot(difx2)

返回:

过差分

理论上,足够次的差分运算可以充分地提取原序列中的非平稳序列的确定性信息。
但应当注意的是,差分运算的阶数并不是越多越好。因为差分计算是一种对信息的提取、加工过程,每次差分都会有信息的损失
在实际应用中差分运算的阶数得适当,应当避免过度差分的使用

举例

假设序列:

接下来考察一阶差分和二阶差分后序列的平稳性与方差

一阶差分:平稳

二阶差分:平稳

虽然二阶差分也是平稳的,但其方差要比一阶差分的大,所以二阶差分就是过差分

小结

  • 差分提取确定性信息:

  • R语言中差分运算diff函数

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

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

相关文章

求一个任意实数c的算术平方根g的算法设计思想_算法复习第四篇——贪心法

公元2020年5月5日&#xff0c;距离算法考试仅剩4天。一、知识归纳1.设计思想只根据当前已有的信息就做出选择&#xff0c;而且一旦做出了选择&#xff0c;将来无论如何都不能更改不从整体最优考虑&#xff0c;所做的选择只是在某种意义上的局部最优这种选择并不总能获得整体最优…

安装百分之80卡住_关注丨男子翻越高铁站台丢命,家属向铁路部门索赔80万!法院这样判...

去年3月&#xff0c;一名男子翻越高铁站台被卡住致死引发广泛关注。事发后&#xff0c;其家属将铁路部门告上法庭&#xff0c;索赔80余万元。日前&#xff0c;法院宣判&#xff1a;死者杨某擅自闯入危险区域负全责&#xff0c;其父母要求铁路部门赔偿的诉请被驳回。事件还原201…

ARIMA模型的介绍

目录 ARIMA模型结构 随机游走模型(random walk) 举例&#xff1a; ARIMA模型的性质 小结 ARIMA模型结构 使用场合&#xff1a;差分平稳序列拟合 模型结构 则有&#xff1a; 随机游走模型(random walk) 模型结构&#xff1a; Karl Pearson(1905)在《自然》杂志上提问:假…

java斗地主游戏

使用java实现斗地主小游戏。 程序分析&#xff1a; 1、创建扑克牌的实体类Poker&#xff0c;设置四个参数&#xff1a;花色、数字、牌值&#xff08;判断大小&#xff09;、是否地主牌&#xff0c;实现getset方法和构造方法&#xff1b; 2、创建玩家的实体类Player&#xff0…

为什么总是封板又打开涨停_警惕!如果股票涨停板反复打开说明了什么?

股市禁忌“三进三出”&#xff1a;“三进”为&#xff1a;一是“高位进”&#xff0c;这就是我们常说的追涨。犯有这种错误的投资者&#xff0c;喜欢在股票已经有了很大的涨幅后买进&#xff0c;其买进的依据往往是股票已经放量突破了历史高点或突破了某个高点阻力位等&#xf…

python中表示空类型的是_python中怎么表示空值

首先了解python对象的概念 python中&#xff0c;万物皆对象&#xff0c;所有的操作都是针对对象的。 那什么是对象&#xff1f;5是一个int对象&#xff0c;‘oblong是一个str对象&#xff0c;异常也是一个对象&#xff0c;抽象一点是&#xff0c;人&#xff0c;猫&#xff0c;够…

动态背景登录

前端动态背景登录 效果演示 完整代码 首先我们要导入jquery.js包 html代码 <!doctype html> <html lang"zh"> <head> <meta charset"UTF-8"> <meta http-equiv"X-UA-Compatible" content"IEedge,chrome1&qu…

ARIMA模型的建模和预测

目录 基本过程&#xff1a; 1.Green 函数递推公式 2.ARIMA模型的预测 例题 小结 基本过程&#xff1a; 1.Green 函数递推公式 ARIMA模型 也可写成 确定 设 则可得Green函数递推公式 2.ARIMA模型的预测 原则&#xff1a;最小均方误差预测原理 传递形式&#xff1a; 预测…

python画画用哪库好_数据可视化哪款工具更好用?对比7款Python 数据图表工具的性能...

更多腾讯海量技术文章&#xff0c;请关注云加社区&#xff1a;https://cloud.tencent.com/developer/column 作者&#xff1a;CDA数据分析师 摘要&#xff1a;Python 的科学栈相当成熟&#xff0c;各种应用场景都有相关的模块&#xff0c;包括机器学习和数据分析。数据可视化是…

jsp论坛网站模版_网站被降权了?看看这些解决方法,或许有帮助哦

点击标题下「蓝色微信名」可快速关注在前面的文章中&#xff0c;笔者已经给大家介绍了网站被降权的因素。那么当遇到降权的问题该怎么办呢&#xff1f;接下来&#xff0c;笔者就结合本身经历过的网站降权来给大家介绍一下补救措施。1、造成网站首页降权的原因&#xff1a;服务器…

ARIMA模型之疏系数模型

目录 1.疏系数模型的定义 2.拟合ARIMA疏系数模型函数 例题&#xff1a; 小结 1.疏系数模型的定义 ARIMA(p,d,q)模型是指d阶差分后自相关最高阶数为p&#xff0c;移动平均最高阶数为q的模型&#xff0c;通常它包含pq个独立的未知系数: 如果该模型中部分自相关系数 ,1≤j<…

动画分享菜单

制作一个动画分享的菜单 效果演示 注意&#xff1a;在我们的程序开始之前我们一定要导入jquery.js架包并且找到如图所示的矢量图标。 完整代码 html代码 <!DOCTYPE html> <html> <head> <meta http-equiv"Content-Type" content"text/…

java 双向链表_23张图!万字详解「链表」,从小白到大佬

链表和数组是数据类型中两个重要又常用地基础数据类型&#xff0c;数组是连续存储在内存中的数据结构&#xff0c;因此它的优势是可以通过下标迅速的找到元素的位置&#xff0c;而它的缺点则是在插入和删除元素时会导致大量元素的被迫移动&#xff0c;为了解决和平衡此问题于是…

安庆机器人编程选_安庆市首届中小学机器人创客竞赛举行 460名选手比智慧拼创新 科技感满满!...

机器人熟练完成“奥林匹克火种”采集、点燃、传递等动作;“小狗皮皮的智能家园”、智能停车场、多功能智能小车等创意智能作品巧思新奇……11月28日&#xff0c;在2020年安庆市中小学机器人创客竞赛活动现场&#xff0c;可谓科技感满满&#xff0c;炫酷十足。本次比赛由市教育体…

有季节效应的非平稳序列分析之因素分解

目录 6.1 因素分解理论 6.2因素分解模型 6.2.1因素分解模型的选择 6.2.2趋势效应的提取 简单中心移动平均的良好属性 R语言中&#xff0c;使用filter函数可以做简单移动平均 6.2.3 季节效应的提取 6.2.4 X11季节调节模型 X11模型分析步骤&#xff1a; ​ 确定性…

python变量的使用_python – 如何在变量中使用冒号(:)

你想要一个 slice() object&#xff1a; index slice(0, 2) print(somelist[index]) slice()模拟您可以在[start&#xff1a;stop&#xff1a;stride]预订语法中指定的start,stop和stride值作为对象. 从文档&#xff1a; Return a 07001 object representing the set of indic…

for循环执行 mybatis_Mybatis中使用循环遍历

Mybatis中传参数的方法1.使用注解绑定&#xff0c;mapper.xml 对应方法 不需要指定 parameterType,(推荐使用注解绑定方式)方法接口&#xff1a;List selectByExample(Param("example") CalculateIdeacommissionsum example,Param("roleNameList") List ro…

指数平滑预测模型

目录 1.分类 2.简单指数平滑 简单指数平滑预测&#xff1a; 平滑系数的确定 R中实现&#xff1a; 3. Holt两参数指数平滑 4. Holt-Winters三参数指数平滑 5.ARIMA加法季节模型 R实现 例题&#xff1a; 1.分类 2.简单指数平滑 简单移动平均法向前预测1期&#xff1a; …

Java反射机制API

实现 Java 反射机制的类都位于 java.lang.reflect 包中&#xff0c;java.lang.Class 类是 Java 反射机制 API 中的核心类。 java.lang.Class 类 java.lang.Class 类是实现反射的关键所在&#xff0c;Class 类的一个实例表示 Java 的一种数据类型&#xff0c;包括类、接口、枚…

ARIMA乘法季节模型

目录 ARIMA乘法季节模型 例题1 例题 2 例题3 ARIMA乘法季节模型 序列的季节效应、长期趋势效应和随机波动之间有着复杂的相互关联性,简单的季节模型不能充分地提取其中的相关关系&#xff0c;这时常采用乘积季节模型。 例题1 我国1949-2008年年末人口总数&#xff08;单…