数据如下👇

非生长季相关性矩阵数据.xlsx

其中挺多细节值得学习

# 载入包
pacman::p_load(tidyverse,ggthemes,showtext,rio)

# 导入数据
data_result <- import("H:\\\\数据处理\\\\2022-05-04_new\\\\非生长季\\\\非生长季相关性矩阵数据.xlsx")

#散点描边与散点形状有关，shape=16无描边属性，shape=21有描边属性和填充属性
font_add("ST_TNR", regular = "G:\\\\字体\\\\宋体和time new roman\\\\ST_TNR.ttf")
showtext_auto()

ggplot() +
  geom_point(data = data_result, aes(x_c, y_c, size = abs(cor), fill = p_d), color = "#999999", shape = 21 ) +
  geom_point(data = data_result[which(data_result$Neg_Pos == "Pos"), ], aes(x_c, y_c, size = abs(cor), color = p_d), shape = 16) +
  scale_fill_manual(values = (c( "#01334E", "#0254BA", "#3877C6", "#9BC2E3" ))) +
  scale_color_manual(values = (c( "#691B11", "#C62E20", "#FC3850", "#FEAEB9" ))) +
  theme_bw() +
  theme(panel.grid.major.x = element_blank(), panel.grid.minor.x = element_blank(), axis.text.x = element_text(angle = 45, hjust = 1, family = "ST_TNR" ), axis.text.y = element_text(family = "ST_TNR"), legend.margin = margin(0, unit = "pt"), text = element_text(family = "ST_TNR")) +
  labs(x = "景观格局指数", y = "大气污染物") +
  guides(size = guide_legend(title = "相关系数", order = 1, override.aes = list(shape = 21, color = "black", fill = "white")), fill = guide_legend(title = "负相关", order = 3, override.aes = list(shape = 21, size = 5, color = "white")), color = guide_legend(title = "正相关", order = 2, override.aes = list(size = 5))) +
  scale_y_discrete(label = c("CO", expression(NO[2]), expression(O[3]), expression(PM[10]), expression(PM[2.5]), expression(SO[2])))

ggsave("F:/YR/本科论文修改/地理科学/编修1/图件_修改/PDF/图12 非生长季景观格局指数与大气污染物浓度相关性.pdf",width = 8,height = 4,units = "in")