陈涛涛的博客分享 http://blog.sciencenet.cn/u/ttchen86

博文

利用R进行时间序列Mann-Kendall突变检测

已有 5491 次阅读 2020-4-24 21:48 |个人分类:R学习|系统分类:科研笔记| 突变检测

library(tidyverse)

library(pheno)

library(rkt)

library(cowplot)


x <- runif(100,min = -3,max = 3)

seqMK.series <- seqMK(x) %>% as.data.frame() %>%

  mutate(date=seq_along(prog))


G <- ggplot(data = seqMK.series, aes(x = date, y = prog)) +

  geom_line(color = "red") +

  geom_line(aes(y = retr), color = "blue") +

  geom_hline(yintercept = c(-1.96, 1.96), linetype = "dashed") +

  labs(x="Date or Year",y= "Kendall's normalized tau's") +

  cowplot::theme_cowplot(font_size = 8,line_size = 0.5) 


G.list <- plot_grid(G,G,G,G,nrow = 2,labels = letters)


Taotao::preview.ggplot(G.list,width = 17,height = 12)


突变点检测.jpg



https://m.sciencenet.cn/blog-3427939-1230000.html

上一篇:利用R绘制指定区域的DEM图
下一篇:利用R进行通径分析

0

该博文允许注册用户评论 请点击登录 评论 (0 个评论)

数据加载中...

Archiver|手机版|科学网 ( 京ICP备07017567号-12 )

GMT+8, 2024-6-3 05:39

Powered by ScienceNet.cn

Copyright © 2007- 中国科学报社

返回顶部