R 语言实现 Superposed Epoch Analysis 时序叠加分析

一、安装并引用 “burnr” 包

Superposed Epoch Analysis 可以使用 “burnr” 包 来实现，首先安装 “burnr” 包：

命令行输入：

install.packages('burnr')

 其次将 burnr 包添加到引用：

library(burnr)

二、准备输入数据

Superposed Epoch Analysis 的输入数据格式如下：

1、 准备时间序列数据放入EXCEL,格式如下：

 之后，读入R，并修改行名称为年份：

data.runoff=read.csv("data_dir")
row.names(data.runoff)=c(1557:2005)

结果如下：

 2、准备矢量格式的事件年数据

data.event=c(1580,1600,1641,1660,1815,1883,1902,1912,1991)

三、执行SEA 并将结果可视化

 运行代码：

data.res=sea(pure,year,4,4)

执行结果格式如下： 

 数据可视化：

plot(res[["actual"]]$lag,res[["actual"]]$mean,type = 'l')

参考文献：

https://www.rdocumentation.org/packages/burnr/versions/0.4.0/topics/run_sea