如何优化24小时双处理组每两小时孢子计数环形柱状图？

问题根源分析 #

1. 时间变量转换错误：原代码dc$Time=as.numeric(levels(dc$Time))[dc$Time]仅适用于有序因子且levels与数值完全匹配的场景，若因子顺序混乱或为无序因子，会导致数值映射错误。
2. 极坐标起始角度未调整：默认coord_polar(start=0)会让0刻度在右侧，而非顶部。
3. 刻度与柱子位置未对齐：未设置合适的刻度间隔与柱子偏移，导致时段范围不清晰。

分步解决方案 #

1. 正确转换时间变量

将Time列转为数值型，优先用字符转数值的方式，避免因子索引错误：

# 若Time是因子类型
dc$Time <- as.numeric(as.character(dc$Time))

# 读取数据时直接指定类型更稳妥
dc <- read_excel('Spore Data for R.xlsm', "Test", col_types = c("numeric", "numeric", "character"))
# 注：col_types按实际列顺序调整，确保Time列设为numeric

2. 优化绘图代码

调整极坐标起始角度、刻度设置与柱子布局，让0/24刻度置顶，柱子对应时段范围：

library(ggplot2)
library(readxl)

# 读取并预处理数据
dc <- read_excel('Spore Data for R.xlsm', "Test")
dc$Time <- as.numeric(as.character(dc$Time))

ggplot(dc, aes(x = Time, y = Mean.Treatment, fill = Treatment)) +
  # 设置柱子宽度与分组间距，确保两组不重叠
  geom_col(position = position_dodge(width = 1), width = 0.8) +
  # 极坐标起始角度设为-pi/2（即顶部），可选顺时针排列刻度
  coord_polar(start = -pi/2, direction = 1) +
  # 设置x轴每2小时一个刻度，扩展范围避免柱子贴边
  scale_x_continuous(
    breaks = seq(0, 24, by = 2),
    limits = c(-1, 25),
    expand = c(0, 0)
  ) +
  # 优化主题提升可读性
  theme_minimal() +
  theme(
    axis.text.x = element_text(size = 10, face = "bold"),
    legend.position = "bottom",
    panel.grid = element_blank()
  ) +
  labs(x = "采样时段（小时）", y = "孢子计数均值", fill = "处理组")

关键参数说明 #

• coord_polar(start = -pi/2)：将极坐标的0度起始点从右侧移至顶部（对应-90度弧度）。
• position_dodge(width = 1)：让黑暗/光照组的柱子在同一时间刻度旁左右错开，避免重叠。
• scale_x_continuous(limits = c(-1,25))：扩展x轴范围，防止首尾时段的柱子被截断。
• direction = 1：可选参数，让时间刻度按顺时针方向排列，符合日常时钟阅读习惯。

内容的提问来源于stack exchange，提问作者Oakley Robertson