如何在R中为动物心率数据框绘制单张目标散点图？

问题分析 #

你遇到的问题根源在于数据格式和绘图函数的默认行为：

• 直接调用plot(data)时，因为你的数据是宽格式（每行对应一个动物，每列对应一个时间点的心率），R默认的plot函数会生成散点图矩阵（所有列两两组合绘制散点图），因此出现大量子图。
• 之前使用melt时未指定id.vars，导致数据转换为长格式时结构错误，无法生成预期的散点图。

解决方案步骤 #

1. 正确读取并整理数据 #

先读取数据，给列赋予有意义的名称，方便后续处理：

# 读取CSV数据（header=FALSE表示第一行不是列名）
data <- read.csv('1.csv', header = FALSE)
# 重命名列：V1为动物ID，V2到V27对应26个时间点的心率
colnames(data) <- c("animal_id", paste0("time_", 1:26))

2. 将宽格式转换为长格式 #

使用reshape2的melt函数，指定animal_id为标识变量，把每个时间点的列转换为行记录：

library(reshape2)
# 转换为长格式
data_long <- melt(
  data, 
  id.vars = "animal_id",          # 保留动物ID作为分组变量
  variable.name = "time_point",   # 新列：记录原列的时间点标识
  value.name = "heart_rate"       # 新列：记录对应心率值
)
# 将time_point转换为数值型（方便作为x轴）
data_long$time_num <- as.integer(gsub("time_", "", data_long$time_point))

3. 绘制目标散点图 #

根据需求选择以下两种绘图方式：

方式1：基础绘图包（无需额外安装）

绘制所有动物心率随时间变化的散点图，可区分不同动物：

# 绘制所有动物的心率散点
plot(heart_rate ~ time_num, 
     data = data_long,
     col = as.factor(data_long$animal_id),  # 不同动物用不同颜色区分
     xlab = "时间点（从刺激前到刺激后）",
     ylab = "心率",
     main = "动物受刺激前后心率变化")
# 添加图例（区分不同动物）
legend("bottomright",
       legend = levels(as.factor(data_long$animal_id)),
       col = 1:length(levels(as.factor(data_long$animal_id))),
       pch = 1,
       cex = 0.6)  # 缩小图例字号，避免遮挡图表

方式2：ggplot2（更美观、灵活）

如果需要更美观的可视化效果，使用ggplot2包：

install.packages("ggplot2")  # 首次使用需安装
library(ggplot2)

ggplot(data_long, aes(x = time_num, y = heart_rate, color = animal_id)) +
  geom_point(alpha = 0.7) +  # 设置点的透明度，避免重叠遮挡
  labs(
    x = "时间点（从刺激前到刺激后）",
    y = "心率",
    title = "动物受刺激前后心率变化散点图"
  ) +
  theme_minimal()

关键说明 #

• 宽格式转长格式是核心：只有将数据转换为动物ID-时间点-心率的长格式，才能绘制时间点x心率的目标散点图。
• 若你的V2到V27列实际对应具体的时间数值（而非序号），可以直接用原列名对应的数值替换time_num，比如提取列名中的数字或使用实际时间值。

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