# 加载ggplot2包
library(ggplot2)
# 自定义文字大小
axis_title_size <- 18
axis_text_size <- 12
label_text_size <- 3.8
title_size <- 18
facet_label_size <- 14
legend_text_size <- 14 # 调整图例文本大小
legend_key_size <- unit(1, "cm") # 调整图例键大小
# 自定义框的高度
tile_height <- 0.6
# 自定义分面间距
facet_spacing <- unit(0.3, "cm")
# 创建数据框
sites <- c("Site 1", "Site 2", "Site 3")
weeks <- c("Week 4", "Week 6", "Week 8", "Week 10")
materials <- c("B73", "Mo17")
df <- expand.grid(site = sites, time = weeks, material = materials, replicate = c("Rep 1", "Rep 2", "Rep 3"))
# 创建新标签映射
label_map <- function(site, week, material) {
site_number <- gsub("Site ", "", site)
week_number <- gsub("Week ", "", week)
material_letter <- ifelse(material == "B73", "B", "M")
paste0(material_letter, site_number, "W", week_number)
}
# 将映射应用到数据框，生成新的标签列
df$label <- mapply(label_map, df$site, df$time, df$material)
# 绘制图形
p <- ggplot(df, aes(x = interaction(time, material, sep = "-", lex.order = TRUE), y = replicate)) +
geom_tile(aes(fill = material), color = "white", width = 0.9, height = tile_height) +
geom_text(aes(label = label), color = "black", size = label_text_size, vjust = 0.5, hjust = 0.5) +
scale_fill_manual(values = c("B73" = "#00BFFF", "Mo17" = "#FF4500")) +
facet_wrap(~ site, ncol = 1, scales = "free_y") +
labs(title = "Field situation diagram", x = "Growth Stage", fill = "Variety") +
theme_minimal() +
theme(axis.text.x = element_text(angle = 45, hjust = 1, size = axis_text_size),
axis.text.y = element_text(size = axis_text_size),
axis.title.x = element_text(size = axis_title_size),
axis.title.y = element_text(size = axis_title_size),
plot.title = element_text(size = title_size, hjust = 0.5),
panel.grid.major = element_blank(),
panel.grid.minor = element_blank(),
strip.text = element_text(size = facet_label_size),
panel.spacing = facet_spacing,
legend.text = element_text(size = legend_text_size),
legend.key.size = legend_key_size)
# 打印图形
print(p)