【数据挖掘】时间序列教程【一】

第一章 说明

        对于时间序列的研究,可以追溯到19世纪末和20世纪初。当时,许多学者开始对时间相关的经济和社会现象进行研究,尝试发现其规律和趋势。其中最早的时间序列研究可以追溯到法国经济学家易贝尔(Maurice Allais)和英国经济学家詹姆斯·克拉克(James Clark)的研究。随着时间序列分析方法的不断发展和应用,时间序列研究逐渐成为了统计学、经济学、金融学、工程学等领域重要的研究方向。

第二章、时态数据的结构

时间序列数据为

  • 根据时间编制索引的观测值或测量值

  • 而不是Xi,我们用Xt表示

为什么这会让事情变得不同?

  • 时间索引具有特殊的排序。

  • 随时间测量的数据不可交换,这是我们在索引数据时通常假设的我我.

  • 时间也可以有其特殊的含义,代表其他未观察到的变量。

        需要明确的是,时间序列数据的一个关键属性是,它与通常分析的其他类型的数据区分开来,我们不认为我们可以随机排列数据的索引并以相同的分布对数据进行建模。数据有排序。此外,数据独立的较强假设通常不适用。

        时间序列数据的一个有趣且可能令人不安的特征是,原始形式的数据提供的真实信息很少。从某种意义上说,原始数据是最无用的数据形式。因此,绘制或汇总原始数据通常无法提供对正在发生的事情或为什么发生的很多见解。但是,由于时间索引具有如此特殊的意义,我们可以使用时间索引将时间序列数据分解为不同时间尺度的变化。时间尺度分析的正式方法有时称为傅里叶分析或光谱分析,但也有一些非正式方法也很有用。

        考虑时间序列数据的另一种方法是,时间序列实际上表示在不同时间尺度上变化的时间序列的混合。分析时间序列数据的部分工作是

  1. 挑选时间尺度的混合并描述它们之间的区别

  2. 根据经验属性或手头的科学问题确定感兴趣的时间尺度

2.1 示例:空气污染与健康

        例如,我们可能有兴趣研究长期暴露于环境空气污染如何影响您的预期寿命。例如,一些研究表明,与生活在更清洁的城市相比,一生生活在污染更严重的城市会使您的预期寿命减少多达 6 个月。在考虑如何解决这个问题以及如何分析数据时,我们主要感兴趣的是比较城市之间的长期平均污染水平,也许是几十年。我们不太可能关心某一天甚至一个月的污染水平有多高。

        另一方面,许多研究表明,空气污染的短期峰值会增加一个城市的心血管和呼吸系统疾病的死亡人数和住院人数。在这种情况下,我们可能有兴趣将空气污染的日常变化与住院或死亡率的日常变化进行比较。总体长期平均污染水平没有什么意义。

        考虑以下 10–1987 年期间密歇根州底特律的颗粒物 (PM1999) 数据的时间序列图。

        有人可能会问一个看似简单的问题:底特律的空气污染在1987年至1999年期间有所改善吗?事实上,在这段时间内,污染水平总体上略有下降,但持续不断。

实际上,当我们查看拟合的简单线性回归模型结果时,我们看到斜率的系数为负。

# A tibble: 2 x 5term        estimate std.error statistic   p.value<chr>          <dbl>     <dbl>     <dbl>     <dbl>
1 (Intercept) 48.4      1.67         28.9  9.59e-170
2 date        -0.00157  0.000184     -8.54 1.77e- 17

但是,在查看上面的图时,很难不注意到定期发生的极端峰值。粗略阅读该图显示了PM10水平达到100 \(mu\)g/m\(^3\)的天数。因此,底特律的PM10随着时间的推移而下降,但我们在某些日子里仍然会经历高水平。情况是否有所改善?

当然,答案与我们考虑数据的时间尺度有关。在长期的时间尺度上,事情似乎正在减少,因此趋势平稳。然而,在短期时间尺度上,我们仍然会看到大的峰值。没有一个答案;答案取决于时间尺度。

从政策角度来看,我们可以采用不同的策略来影响长期和短期时间尺度的空气污染。为了长期改变污染水平,我们可能会尝试将当地经济从基于化石燃料的能源转变为更可再生、污染更少的能源。这样的计划可能会产生重大影响,但需要大量时间才能实施。为了应对污染的短期波动,我们可能会实施交通禁令或有针对性的基于源的干预措施等政策,以缓解短期高峰。

现在假设我们想看看底特律的死亡率和空气污染之间是否存在任何关联。我们可以制作一个简单的散点图来查看是否存在简单的关联。

现在,这个散点图是我们在没有时间序列数据时可能会制作的。但是,由于我们确实有时间序列数据,因此我们应该立即开始根据不同的时间尺度变化来思考问题。我们关心的是污染和死亡率之间的长期关联,还是关心短期关联?

上图中显示的整体关联可以用简单的线性回归模型进行量化。

# A tibble: 2 x 5term        estimate std.error statistic    p.value<chr>          <dbl>     <dbl>     <dbl>      <dbl>
1 (Intercept)  46.0      0.226      204.   0         
2 pm10          0.0275   0.00564      4.88 0.00000108

两者之间似乎存在正相关,表明空气污染水平的增加与死亡率的增加有关。但是,我们能做更多的事情来获得更多的洞察力吗?

让我们计算PM10的年平均值,并对年度死亡总数求和,并制作这些年度汇总统计数据的散点图。

从这个图中我们可以看到,这种关联似乎相当强(当然,只有 13 个数据点)。当我们为这些数据拟合线性模型时,我们得到以下内容。

# A tibble: 2 x 5term        estimate std.error statistic       p.value<chr>          <dbl>     <dbl>     <dbl>         <dbl>
1 (Intercept)   10388.     987.      10.5  0.00000000725
2 pm10            190.      29.4      6.47 0.00000579   

从一年到下一年,年平均PM10变化一个单位,与190.4例死亡的变化有关。

现在,我们可以将PM10的每日偏差与其年平均值进行比较,并查看该偏差与每日死亡率之间的关联。

# A tibble: 2 x 5term        estimate std.error statistic p.value<chr>          <dbl>     <dbl>     <dbl>   <dbl>
1 (Intercept)  46.9      0.116      404.    0     
2 pm10dev       0.0142   0.00574      2.47  0.0136

这里的关联要小得多,但当然我们只关注PM10的每日变化,而不是年度变化。我们预计PM10从一天到第二天的一个单位变化会带来大量死亡。

人们可能会想知道:哪个估计是正确的?是平均PM10与死亡率之间的关联,还是平均PM10与死亡率之间的关联?答案是两者都是“正确的”,但每个都回答了不同的问题。日平均值着眼于短期变化,可以解释为代表污染的“急性”影响,而年平均值可能会重新反映空气污染水平的“慢性”影响。

在不同时间尺度上观察关联时要考虑的另一个问题是,在这个时间尺度上存在的混杂因素是什么?在观察PM10的逐年变化时,可能存在许多混杂因素,这些混杂因素也因PM10和死亡率而异。在观察PM10的每日变化时,每年平稳变化的相同混杂因素可能不会引起关注。但是,可能还有其他混杂因素每天都需要考虑。

2.2 固定变化与随机变化

大多数时间序列书籍倾向于将时间序列想象为仅由随机现象组成,而不是固定和随机现象的混合。因此,建模通常侧重于时序模型的随机方面。然而,世界上许多实时序列是由我们可能认为的固定和随机变化组成的。

  • 温度数据具有非“随机”的昼夜和季节性成分

  • 空气污染数据可能会根据交通或通勤模式产生星期的影响

虽然有时很容易将一切都视为随机的,但当我们缺乏对真正潜在现象的观察时,这通常是一根拐杖。此外,当某些东西是固定的时,将它视为随机将导致违反我们通常做出的平稳性假设(见下文)。

根据应用程序的性质,对固定或随机的相同现象进行建模可能是有意义的。换句话说,这取决于。

  • 在生物医学和公共卫生应用中,我们通常处理完全观察到的数据集,并试图解释“发生了什么?

  • 我们正在描述过去,也许是对未来做出推断

  • 在金融或控制系统应用中,我们可能会根据过去对未来事件进行预测。从过去数据中看似固定的事情将来可能会发生变化,因此我们可能希望允许模型“适应”未知的未来模式。

考虑以下 1990–1992 年马里兰州巴尔的摩的日平均温度图。正如人们从温度数据中预期的那样,有一个强烈的季节性模式,在夏季达到顶峰,在冬季出现低谷。

现在,这种所谓的季节性模式是固定的还是随机的?历史告诉我们,季节性模式是相当可预测的。我们通常不相信夏天可以冻僵,冬天也可以达到90度(F)。

讨论此问题的更正式方式可能是使用以下模型。设 \(y_1, y_2, \dots\) 是巴尔的摩每天 \(t\) 的温度值,并考虑以下模型,

\[ y_t = \mu + \varepsilon_t, \]

其中 \(\varepsilon\) 是期望值 \(\mu\) 和观测值 \(y_t\) 之间的随机偏差。在没有任何计算机的帮助下,我们可能会查看上面的图并估计\(\mu\)大约是50-55度。但是现在,假设您的工作是预测 \(t\) 的任何值的 \(\varepsilon_t\) 的值。很明显,如果 \(t\) 在年中下降,则很可能是 \(\varepsilon_t > 0\),如果 \(t\) 在年初或年底附近下降,那么很可能是 \(\varepsilon_t < 0\)。因此,我们只需知道 \(t\) 的值即可获得有关偏差 \(\varepsilon\) 的重要信息。换句话说,序列\(\varepsilon_t\)中嵌入了一个固定的季节性效应,我们可能很难将其视为“随机”。

但现在考虑以下模型。

\[ y_t = y_{t-1} + \varepsilon_t。 \]

此模型将 \(y_t\) 的值预测为与时间 \(t-1\) 处的值的偏差。所以今天的值等于昨天的值加上一个小的偏差。现在,假设您的工作是预测 \(\varepsilon_t\) 的值。这有点难,对吧?如果我知道昨天是70度,我确定今天会比70度温暖吗?还是更冷?如果我知道昨天是20度,我确定今天会更暖和还是更冷?在此模型中,偏差\(\varepsilon_t\)可能看起来更“随机”或更难预测。没有一个固定的规则说今天的温度总是比昨天的温度更暖(或更冷)。

考虑下面的不同时间序列,它显示了股票代码为 SPY 的交易所交易基金的加权中间交易价格。该基金追踪美国股票的标准普尔500指数。请注意,x 轴上的时间刻度以微秒为单位。

与温度时间序列相比,该图看起来不那么规则,并且没有可识别的模式。此外,在微秒级,我们可能不太熟悉这种股票价格可能存在的固定模式。有经验的交易者可能知道,在一天中的给定时间,在几十万微秒的窗口内,这种模式总是会出现。

然而,对于金融,有一种被称为有效市场假说的理论认为这种固定模式不应该存在。如果存在这样一个固定的模式,它将代表一个套利机会,或者一个没有风险的赚钱机会。例如,在上面的图中,我们可以在 2 万微秒的时间买入股票,然后在 5 万微秒左右卖出,以获得轻松的利润。如果这种模式每天都存在,我们可以告诉我们的经纪人每天执行这笔交易以获得微薄的利润。然而,随着这种模式的消息泄露到市场上,越来越多的人会开始和我同时买入,和我同时卖出。这将在买入时提高价格,在卖出时降低价格,最终获利机会将消失。

有效市场假说表明,这种固定模式的存在极不可能。因此,将此类数据建模为随机数据可能更有意义,而不是固定数据。这表明可以采用不同的建模策略和不同类型的模型。我们不会在这里详细讨论这些类型的模型。

2.3 时间序列分析的目标

人们希望从时间序列分析中得到什么?正在回答哪些问题?

2.3.1 预测

鉴于过去和现在,未来会是什么样子(及其不确定性)?

  • 鉴于过去10年的季度每股收益,苹果公司下个季度的每股收益是多少?

  • 鉴于过去200年的全球平均气温,未来100年全球平均气温是多少?

2.3.2 过滤

鉴于过去和现在的观察,我应该如何更新我对自然真实状态的估计?

  • 鉴于我目前对航天器位置和速度的估计,我应该如何根据新的陀螺仪和雷达测量更新我对位置和速度的估计?

  • 鉴于美国每月失业数据的历史以及我对当前失业水平的估计,我应该如何根据劳工统计局发布的最新数据修改我的估计?

  • 考虑到捐赠回报的历史,当年的回报,以及每年花费捐赠价值的目标百分比的需要,大学应该在下一个财政年度从捐赠基金中花费多少?

2.3.3 时间尺度分析

给定一组观察到的数据,哪些时间尺度的变化主导或解释了数据中的大部分时间变化。

  • 马里兰州巴尔的摩的温度观测是否有强烈的季节性周期?

  • 环境空气污染与死亡率之间的关联主要是由污染水平的年度大幅变化还是短期峰值驱动的?

2.3.4 回归建模

给定两种现象的时间序列,它们之间有什么关联?

  • 每日空气污染水平与心脏住院每日数值之间有什么关联?

  • 一个国家失业率的变化与国内生产总值的变化之间的滞后(以月为单位)是多少?

  • 大飓风发生后两周内发生的累计超额死亡人数是多少?

2.3.5 平滑

给定一个完整的(嘈杂的)数据集,我可以推断出过去自然的真实状态吗?

  • 给定一个噪声测量信号,我能从数据中重建真实信号吗?

  • 现在我的宇宙飞船已经绕月飞行了,它离月球最近的距离是多少?

  • (待续)

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

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

相关文章

设计模式学习之工厂方法模式

设计模式系列往期文章 设计模式学习之策略模式设计模式学习之策略模式在前端的应用设计模式学习之简单工厂模式 在上一篇文章中我们学习了简单工厂模式——这是工厂模式中最简单的一种模式&#xff0c;通过工厂类提供的方法创建类&#xff08;可以类比为产品&#xff09;&…

界面开发框架Qt新手入门教程:如何使用Calendar组件创建日历(二)

Qt 是目前最先进、最完整的跨平台C开发工具。它不仅完全实现了一次编写&#xff0c;所有平台无差别运行&#xff0c;更提供了几乎所有开发过程中需要用到的工具。如今&#xff0c;Qt已被运用于超过70个行业、数千家企业&#xff0c;支持数百万设备及应用。 本文中的CalendarWi…

主动式和被动式电容笔哪个好用?苹果平板平替笔排行

被动式电容笔与主动式电容笔最大的不同在于主动式电容笔具有更广泛的应用领域&#xff0c;可以与不同种类的电容屏幕进行匹配。随着人们对其了解的不断深入&#xff0c;其应用也日益广泛。除此之外&#xff0c;平替电容笔的技术&#xff0c;也在不断的改进和提高&#xff0c;逐…

ES批量上传数据 - Python操作ES

代码 # -*- coding: utf-8 -*- # Author : markadcimport randomfrom elasticsearch import Elasticsearch, helpersurl http://localhost:9200/ client Elasticsearch(url) index testsome [] for i in range(1, 10001):docu {uid: i, age: random.randint(1, 100)} …

【网络】协议的定制与Json序列化和反序列化

文章目录 应用层初识TCP协议通讯流程定制协议再谈协议网络版本计算器Protocal.hppCalServerCalClient Json的安装 应用层 我们程序员写的一个个解决我们实际问题, 满足我们日常需求的网络程序, 都是在应用层 初识TCP协议通讯流程 建立链接和断开链接 基于TCP协议&#xff0c…

[Unity][疑难杂症]Mac电脑打开工程unity attempt to write a readonly database

报错的根本原因&#xff1a; 你Unity项目所在的硬盘&#xff0c;格式不是MacOS或者APFS格式&#xff0c;我的是ExFAT所以报这个错。 可能出现的问题&#xff1a; 1、删除unity项目中的c#之类的文件就会报错&#xff0c;并且强制退出 2、新项目压根打不开&#xff0c;在生成资源…

图片:前端展示图像(img 、picture、svg、canvas )及常用图片格式(PNG、JPG、JPEG、WebP、GIF、SVG、AVIF等)

一、浏览器网页展示图片方法 1.1、HTML <img> 标签 <!DOCTYPE html> <html><head><title>图片展示</title></head><body><h1>图片展示</h1><img src"example.jpg" alt"Example Image" w…

IPO观察丨健身器材行业加“数”前行,康力源如何重构竞争壁垒?

近年来&#xff0c;健身器材市场的火热&#xff0c;引来了资本的广泛关注。数据显示&#xff0c;2021年中国健身器材市场规模为546.5亿元&#xff0c;同比增长15%&#xff0c;预计2024年或可达到799.6亿元。在这股浪潮之下&#xff0c;多家健身器材企业递出了IPO申请。比如&…

合同数智化如何助力地产企业实现变革“突围”?

从稳步发展到求新求变&#xff0c; 数智化成破局关键 近年来&#xff0c;随着宏观经济政策调整&#xff0c;在中央房住不炒的大基调下&#xff0c;房地产逐步回归居住属性。在这样的大背景下&#xff0c;针对不同类型的房地产企业&#xff0c;国家出台了不同的数字化转型指导文…

Linux——1初识linux

目录 1.1 硬件和软件 1.2 初识Linux 1.2.1 Linux的诞生 1.2.2 LInux内核 1.2.3 Linux发行版 1.3 虚拟机介绍 1.4 VMware WorkStation 安装 1.5 在VMware上安装Linux 1.6 远程连接Linux系统 1.6.1 图形化、命令行 1.6.2 FinalShell 1.1 硬件和软件 我们所熟知的计算…

基于安卓的初中英语语法APP设计与实现(源码+文档+PPT)

智能手机的普及使移动学习成为一种全新的学习方式,越来越受到研究者与学习者的关注。调查显示,移动学习对英语学习特别是应用英语的学习有良好的效果。基于Android设计开发的英语语法APP包含“练习题册”、“语法分析”、“搜索语法”三大功能模块,能够帮助学生利用碎片时间学习…

发起投票平台投票吧网络投票平台网络投票平台

小程序投票活动如何做&#xff1f;很多企业在运营当中&#xff0c;都会通过投票活动来进行推广&#xff0c;从而达到吸粉、增加用户粘度等效果。而此类投票活动&#xff0c;通过小程序就可以实现&#xff0c;操作简单。 我们现在要以“促进诚信经营”为主题进行一次投票活动&am…