不同起止时间与频率的时序数据绘图方法咨询

先解决plot函数的报错 #

你遇到的Error in xy.coords(x, y) : 'x' and 'y' lengths differ，核心原因是没有按name分组绘图，或者某组数据里存在NA/重复值导致timestamp和value的有效长度不匹配。

快速排查与修复 #

1. 先检查每组数据的完整性：

library(dplyr)
# 按name分组，查看每组的时间戳、数值数量及是否有NA
data %>%
  group_by(name) %>%
  summarise(
    时间戳数量 = n_distinct(timestamp),
    数值数量 = n(),
    存在缺失值 = any(is.na(value))
  )

如果某组的时间戳数量≠数值数量，说明有重复时间戳或缺失值，可通过distinct(timestamp, .keep_all = TRUE)去重，或tidyr::fill(value)补全缺失值。

2. 用基础plot函数正确分组绘图：

# 获取所有唯一分组名称
unique_names <- unique(data$name)
# 设置多图布局（根据分组数量调整行列）
par(mfrow = c(2, 2))
for (name in unique_names) {
  # 提取当前分组的数据
  subset_data <- data[data$name == name, ]
  # 确认长度一致后绘图
  if (nrow(subset_data) == length(subset_data$value)) {
    plot(subset_data$timestamp, subset_data$value, 
         main = name, xlab = "时间", ylab = "数值", type = "l")
  } else {
    warning(paste("分组", name, "数据长度不匹配，跳过绘图"))
  }
}
# 恢复默认布局
par(mfrow = c(1,1))

推荐更适合时序绘图的R库 #

1. ggplot2（最常用，适合静态图） #

分组绘图和子图展示非常便捷，自动适配不同起止时间和采样频率：

library(ggplot2)
# 所有分组在同一张图（按颜色区分）
ggplot(data, aes(x = timestamp, y = value, color = name)) +
  geom_line() +
  labs(x = "时间", y = "数值", title = "各分组时序变化") +
  theme_minimal()

# 每个分组单独画子图（x轴自动适配各自时间范围）
ggplot(data, aes(x = timestamp, y = value)) +
  geom_line(color = "steelblue") +
  facet_wrap(~name, scales = "free_x") +
  labs(x = "时间", y = "数值", title = "各分组独立时序图") +
  theme_minimal()

2. plotly（交互式绘图，适合Shiny App） #

支持缩放、悬停查看数值，完美适配Shiny的交互需求：

library(plotly)
# 基于ggplot2转交互式图
p <- ggplot(data, aes(x = timestamp, y = value, color = name)) +
  geom_line() +
  labs(x = "时间", y = "数值", title = "交互式时序图")
ggplotly(p)

# 原生plotly语法
plot_ly(data, x = ~timestamp, y = ~value, color = ~name, type = "scatter", mode = "lines")

3. xts/zoo（专门处理时序数据） #

适合对时序数据做清洗、对齐操作后绘图：

library(xts)
# 按分组转换为xts对象
xts_list <- lapply(unique_names, function(n) {
  sub <- data[data$name == n, ]
  xts(sub$value, order.by = sub$timestamp)
})
names(xts_list) <- unique_names

# 合并后绘图
plot.xts(do.call(merge, xts_list), 
         main = "各分组时序图", 
         col = 1:length(unique_names))
legend("topright", legend = unique_names, col = 1:length(unique_names), lty = 1)

Shiny App适配建议 #

可以做一个简单的交互式选择组件，让用户选择查看指定分组的时序图：

library(shiny)
library(ggplot2)

# 预处理数据（提前导入并转换时间格式）
data$timestamp <- as.POSIXct(data$timestamp)

ui <- fluidPage(
  selectInput("name_select", "选择分组", choices = unique(data$name)),
  plotOutput("time_plot")
)

server <- function(input, output) {
  output$time_plot <- renderPlot({
    subset_data <- data[data$name == input$name_select, ]
    ggplot(subset_data, aes(x = timestamp, y = value)) +
      geom_line(color = "steelblue") +
      labs(x = "时间", y = "数值", title = paste("分组", input$name_select, "的时序变化")) +
      theme_minimal()
  })
}

shinyApp(ui, server)

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