当前位置: 代码迷 >> 综合 >> R 语言实现 Superposed Epoch Analysis 时序叠加分析
  详细解决方案

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

热度:85   发布时间:2024-01-09 04:44:15.0

一、安装并引用 “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 

  相关解决方案