第11章

11.1干预分析

library(TSA)
win.graph(width = 4.875,height = 2.5,pointsize = 8)
data(airmiles)
plot(log(airmiles),ylab = 'log(airmiles',xlab = 'year')

#美国航空的每月客运里程：1996年1月~2005年5月

as.vector(diff(diff(window(log(airmiles),end = c(2001,8)),12))),lag.max = 48)

#干预期（1-B）（1-B12）log（航空客运里程）的样本ACF　　

air.ma1 <- arimax(log(airmiles),order = c(0,1,1),seasonal = list(order = c(0,1,1),period = 12),xtransf = data.frame(I911 = 1*(seq(airmiles)==69),I911 = 1*(seq(airmiles)==69)),transfer = list(c(0,0),c(1,0)),xreg = data.frame(Dec96 = 1*(seq(airmiles)==12),Jan97 = 1*(seq(airmiles)==13),Dec02 = 1*(seq(airmiles)==84)),method = 'ML')
air.ma1

　　#对数化航空客运里程的干预模型的估计

plot(log(airmiles),ylab = 'Log(airmiles)')
points(fitted(air.ma1))

　　#对数化的航空客运里程与拟合值

Nine11p <- 1*(seq(airmiles)==69)
plot(ts(Nine11p*(-0.0949)+filter(Nine11p,filter = 0.8139,method = 'recursive',side = 1)*(-0.2715),frequency = 12,start = 1996),ylab = '9/11',type = 'h')
abline(h = 0)

　　#9.11事件对航空客运量序列所造成影响的估计

11.2异常值

　　指的是一些不规则的观测值，其出现可能源自测量误差与复制误差其中之一，或者两者都有可能，也可能源于基础过程发生了短期性变化。

　　对于时间序列来说可识别的异常值有两种，可加异常值与新息异常值，简记为AO/IO.

转载于:https://www.cnblogs.com/xxupup/p/10588650.html

《时间序列分析及应用.R语言》第十一章阅读笔记相关推荐

时间序列分析及应用r语言pdf_R语言时间序列分析（十一）：指数平滑法
作者:黄天元,复旦大学博士在读,热爱数据科学与开源工具(R),致力于利用数据科学迅速积累行业经验优势和科学知识发现,涉猎内容包括但不限于信息计量.机器学习.数据可视化.应用统计建模.知识图谱等,著有& ...
时间序列分析及应用r语言pdf_R语言：时间序列经典分析法(二)
题记:本文是个人的读书笔记,仅用于学习交流使用.本文将深入研究时间序列技术. 01 解决什么问题? 前面一章,介绍了时间序列中涉及到的基本概念,本章将在此基础上介绍如何对时间序列的资料进行分析,怎么选 ...
R语言使用timeROC包计算无竞争情况下的生存资料多时间AUC值、R语言使用timeROC包的plotAUCcurve函数可视化多时间生存资料的不同标记物情况下对应的AUC曲线、并进行对比
R语言使用timeROC包计算无竞争情况下的生存资料多时间AUC值.R语言使用timeROC包的plotAUCcurve函数可视化多时间生存资料的不同标记物情况下对应的AUC曲线.并进行对比目录
R语言实战前四章学习小结-支付宝账单的简单作图
学习了R语言的前四章后,决定做个小项目练练手,翻了下知乎,发现有人用自己的支付宝账单做过相关的简单分析,所以我也打算用这个练练手首先从支付宝首页下载了2017年1月到6月的支付宝账单,对其中的付款时 ...
R语言与抽样技术学习笔记（Jackknife）
R语言与抽样技术学习笔记(Randomize,Jackknife,bootstrap) Jackknife算法 Jackknife的想法在我很早的一篇博客<R语言与点估计学习笔记(刀切法与最小二 ...
R语言与函数估计学习笔记（函数模型的参数估计）
R语言与函数估计学习笔记毫无疑问,函数估计是一个比参数估计要复杂得多的问题,当然也是一个有趣的多的问题.这个问题在模型未知的实验设计的建模中十分的常见,也是我正在学习的内容的一部分. 关于函数估计我 ...
R语言基础入门（学习笔记通俗易懂版）
文章目录 R语言预备知识获取工作目录设置工作目录注释变量名的命名赋值变量的显示查看与清除变量函数帮助文档查询函数安装R包文件的读取文件的输出软件的退出与保存 R语言语法向量 ...
【R语言】第二章上机实验记录
书籍:<R语言与数据挖掘> 1. 实验目的了解R语言中数据类型的判别及转换函数.应用方法了解R语言中对数据结构操作的函数.应用方法了解R语言中读写数据文件的方法 2.实验内容掌握读 ...
c语言第七章函数笔记,我的C语学习笔记-C语言教程(七).doc
我的C语学习笔记- C语言教程(七) C语言教程---第一章: C语言概论 C语言教程---第二章: 数据类型.运算符.表达式 C语言教程---第三章: C语言程序设计初步 C语言教程---第四章: ...

《时间序列分析及应用.R语言》第十一章阅读笔记

第11章

11.1干预分析

11.2异常值

《时间序列分析及应用.R语言》第十一章阅读笔记相关推荐

最新文章

热门文章