时间序列的预处理之纯随机性检验
目录
1.纯随机序列的定义
2.性质
3.纯随机性检验
1.纯随机序列的定义
- 纯随机序列也称为白噪声序列,满足如下性质:
2.性质
- 纯随机性(无记忆性)
- 方差齐性
举例,随机生成1000个白噪声序列
用正态分布序列 rnorm(数量,均值,方差),如下为,1000个标准正态的分布图形
a<-rnorm(1000)
x<-ts(a)
plot(x)
时序图如下:
自相关图:
acf(x)
返回:
如图,可以看出自相关系数基本分布在二倍标准差之间,但由于数据是随机的,所以有一定的误差
3.纯随机性检验
Bartlett定理:如果一个时间序列是纯随机的,得到一个观察期数为n的观察序列,那么该序列的延迟非零期的样本自相关系数将近似服从均值为0,方差为序列观察期数倒数的正态分布
原假设:
备择假设:
Q统计量(Box和Pierce):服从自由度为m的卡方分布,其对大样本检验效果较好
LB统计量(Box和Ljung):Q统计量的修正,现在使用较普遍
拒绝域:
R语言白噪声检验
Box.test(x,type=,lag=6)
其中:
type='Box-Pierce' Q统计量 ,默认
type='Ljung-Box' LB统计量
举例1:对上面随机生成的白噪声序列进行6阶和12阶的LB统计量
for(i in 1:2)print( Box.test(x,type='Ljung-Box',lag=6*i))
返回:
举例2:对1900年到1998年全球7级以上地震法伤次数序列进行平稳性和纯随机性检验
读取数据
a<-read.table('D:/桌面/E2_5.csv',sep=',',header=T)
a
返回:
选择变量序列
x<-ts(a$number,start=1900)
x
返回:
绘制时序图:
plot(x)
绘制自相关图:
acf(x)
6阶LB统计量:
Box.test(x,type='Ljung-Box',lag=6)
返回:
时间序列的预处理之纯随机性检验相关推荐
- 时间序列分析学习笔记:时间序列的预处理(平稳性检验、纯随机性检验)
1 时间序列预处理 2 平稳性检验 2.1 特征统计量(概率分布的意义) 2.2 时间序列的概率分布 2.3 概率分布族应用的局限性 2.4 特征统计量(均值.方差) 2.5 平稳时间序列的定义(严平 ...
- 平稳性检验(描述性)与纯随机性检验
这篇博客主要记录人大出版<应用时间序列分析>第二章的笔记.本章主要介绍进行时序分析前的预处理,即平稳性检验与纯随机性检验. 平稳性检验(描述性) 平稳性检验的方法分为描述性方法与计量性方法 ...
- 序列的平稳性与纯随机性检验,模型的有效性,参数的显著性,最优模型准则AIC,SBC
问题:某城市过去63年中每年的降雪量数据(单位:mm). 判断该序列的平稳性与纯随机性 如果序列平稳且非白噪声,选泽适当模型拟合序列发展 利用拟合模型,预测该城市未来5年的降雪量 data <- ...
- 时间序列的预处理——平稳性检验的R语言实现(二)
这一部分是时间序列预处理R语言的实现. 目标是将课本和上课知识点整合. 老师是用一节课讲完的,本篇文章只做了平稳性检验~~~ 下一篇再写纯随机性检验 全部代码 #input data yield &l ...
- 时间序列的平稳性检验与随机性检验
1. 时间序列的定义 1.1 什么是时间序列 在统计研究中,常用按时间顺序排列的一组随机变量X1,X2,⋯ ,Xt,⋯X _ { 1 } , X _ { 2 } , \cdots , X _ ...
- 时间序列分析——基于R | 第2章 时间序列的预处理习题代码
时间序列分析--基于R | 第2章 时间序列的预处理习题 1.考虑序列{1,2,3,4,5,-,20} 1.1判断该序列是否平稳 x <- seq(1,20);x ## [1] 1 2 3 4 ...
- 二、时间序列的预处理
一般情况下,拿到一个观察值序列之后,首先要对它的平稳性和纯随机性进行检验,这两个重要的检验称为序列的预处理.根据检验的结果可以将序列分为不同的类型,对不同类型的序列我们会采用不同的分析方法. 一.平稳 ...
- ch02_时间序列的预处理
第二章 时间序列的预处理 文章目录 2.1 平稳性检测 2.1.1 特征统计量 2.1.1.1 概率分布 2.1.1.2 特征统计量 2.1.2 平稳时间序列的定义 2.1.3 平稳时间序列的统计性质 ...
- 时间序列平稳性分析和白噪声检验
文章目录 一.时间序列平稳性 1.ADF检验 1.1.ADF检验原理 1.2.ADF的python实现 2.ACF和PACF 二.白噪声检验 原始的负荷时间序列曲线 一.时间序列平稳性 时间序列分析之 ...
最新文章
- http://blog.sina.com.cn/s/blog_458f3c010100n4st.html
- b+tree数据结构可视化_数据结构: B+Tree及其应用
- 【Java报错】记录一次调用递归方法导致的 StackOverFlowError 及如何重构递归代码避免栈溢出
- 在VS2010 中使用subversion 进行代码的分支与合并
- Spring中的Spring JSR-250 注释之@Resource
- 视觉SLAM笔记(24) 图像基础操作
- 南林计算机课程表,南京林业大学课程表.doc
- latex中pdflatex与xelatex的区别
- 2022-2028全球与中国移动卫星通信市场现状及未来发展趋势
- 计算机 取得高级权限,win10获取system权限,win7获取最高权限
- face_recognition 人脸识别
- Windows批处理文件bat学习(一)
- TreeUtil(树形结构转换)
- 苏州大学计算机学院导师洪宇介绍,苏州大学计算机科学与技术学院硕导介绍:马小虎...
- Python 如何实时绘制数据
- socket通信函数
- 元境技术助力元宇宙营销 联合发起商广协元宇宙营销研究院
- STM32F407使用USB作虚拟串口出现黄感叹号
- 【php毕业设计】基于php+mysql+apache的课程网站管理系统设计与实现(毕业论文+程序源码)——课程网站管理系统
- 转载百度百科 python
热门文章
- 360浏览器怎么关闭全屏看大图提示条
- acer电脑设置u盘启动方法
- JEECG 缓存用法
- 一般java,【面试题】java一般
- java实现把数据写入到Excel并下载
- 操作痕迹包括那些_高级消防设施操作员专题之:走近气体灭火系统
- 电子科大16秋计算机应用基础在线作业一,电子科大16秋《计算机应用基础(本科)》在线作业3...
- python 整数输出 d f_如何将数字(10,11,12,13,14,15)分配给Python 3中的字母(A,B,C,D,E,F)?...
- 消息订阅与发布(pubsub)
- php 物联网应用,蜂窝物联网的概念以及应用