r - 如何在R中格式化SPEI包的硬编码绘图功能的x轴?
问题描述
我正在使用 SPEI 包及其 32 年的每月样本数据。我想修改 x 轴标签以反映年份而不是数字。但是,硬编码绘图功能不允许我这样做。我厌倦了提取SPEI$fitted
数据并尝试使用 ggplot 复制相同的图,但没有成功。这是示例代码
install.packages("SPEI")
library(SPEI)
data("wichita")
wichita$PET=hargreaves(Tmin=wichita$TMIN, Tmax = wichita$TMAX, lat = 37.64)
wichita$BAL=wichita$PRCP - wichita$PET
SPEI_12=spei(wichita[,"BAL"],12)
plot.spei(SPEI_12, main = 12-Month SPEI)
解决方案
我不太了解这个plot.spei
功能,所以我使用了ggplot2
.
基本上,我用拟合值的 ts 构建了一个数据框,并为正 ( ) 或负 ( ) 值创建了颜色/填充条件。pos
neg
library(zoo)
library(tidyverse)
DF <- zoo::fortify.zoo(SPEI_12$fitted)
DF <- DF %>%
dplyr::select(-Index) %>%
dplyr::mutate(Period = zoo::as.yearmon(paste(wichita$YEAR, wichita$MONTH), "%Y %m")) %>%
na.omit() %>%
dplyr::mutate(sign = ifelse(ET0_har >= 0, "pos", "neg"))
ggplot2::ggplot(DF) +
geom_bar(aes(x = Period, y = ET0_har, col = sign, fill = sign),
show.legend = F, stat = "identity") +
scale_color_manual(values = c("pos" = "darkblue", "neg" = "red")) +
scale_fill_manual(values = c("pos" = "darkblue", "neg" = "red")) +
scale_y_continuous(limits = c(-3, 3),
breaks = -3:3) +
ylab("SPEI") + ggtitle("12-Month SPEI") +
theme_bw() + theme(plot.title = element_text(hjust = 0.5))
编辑:另一个想法。
DF2 <- DF %>%
tidyr::spread(sign, ET0_har) %>%
replace(is.na(.), 0)
ggplot2::ggplot(DF2) +
geom_area(aes(x = Period, y = pos), fill = "blue", col = "black") +
geom_area(aes(x = Period, y = neg), fill = "red", col = "black") +
scale_y_continuous(limits = c(-3, 3),
breaks = -3:3) +
ylab("SPEI") + ggtitle("12-Month SPEI") +
theme_bw() + theme(plot.title = element_text(hjust = 0.5))
推荐阅读
- java - 写入输出到控制台的字符串并进行编辑
- c# - 创建索引失败
- python - 使用 sklearn.feature_extraction.text 时出现 KeyError
- swift - 更改选项卡时清除视图数据 SwiftUI
- python-2.7 - 具有多层的 TensorFlow 训练函数
- c++ - 在 Windows 10 C++ 中使用 include sys/resource.h 时出错
- typo3 - 复制页面时,slug 字段未正确更新
- opencart - 如何修复 Opencart 3 页面加载错误?
- php - 记录未添加到数据库codeigniter 4
- jbehave - 无法在没有 delemeters 的情况下转换 Jbehave 命名参数