使用changepoint包绘图返回NULL，无法用grid.arrange组合ggplot求助

问题核心在于changepoint包的plot()函数基于R基础绘图系统，仅即时渲染图像而不返回可操作的绘图对象（返回NULL），因此无法直接用gridExtra::grid.arrange与ggplot图表组合。以下提供两种可行解决方案：

方案1：用gridGraphics捕获基础绘图为grob对象 #

通过gridGraphics包将基础绘图转换为grid系统的grob对象，即可支持与ggplot图表组合：

# 安装并加载依赖包
install.packages("gridGraphics")
library(gridGraphics)
library(gridExtra)
library(changepoint)

# 原有分析代码
set.seed(1)
df <- data.frame(date_rec = seq(as.Date("2014-01-01"),as.Date("2015-12-31"), by=1),
                 number = sample(c(1:10),size = 730, replace=TRUE))

monthly_data <- aggregate(number ~ year(date_rec) + month(date_rec), data = df, FUN = sum)
df_ts <- ts(monthly_data$number, frequency = 12, start = c(2014, 1))
cpt_mean <- cpt.mean(df_ts, method = "AMOC", penalty = "AIC")

# 捕获变点图为可组合的grob对象
grid.echo({plot(cpt_mean)})
p_changepoint <- grid.grab()

# 与其他ggplot图表组合（示例）
# grid.arrange(p1, p2, p3, p_changepoint, ncol=2)

方案2：用ggplot手动重绘变点图（推荐） #

直接基于ggplot2重构变点分析图，生成原生ggplot对象，与其他图表组合更流畅：

library(ggplot2)
library(changepoint)
library(gridExtra)

# 原有分析代码
set.seed(1)
df <- data.frame(date_rec = seq(as.Date("2014-01-01"),as.Date("2015-12-31"), by=1),
                 number = sample(c(1:10),size = 730, replace=TRUE))

monthly_data <- aggregate(number ~ year(date_rec) + month(date_rec), data = df, FUN = sum)
# 转换为日期格式适配ggplot
monthly_data$date <- as.Date(paste(monthly_data$`year(date_rec)`, monthly_data$`month(date_rec)`, "01", sep = "-"))
df_ts <- ts(monthly_data$number, frequency = 12, start = c(2014, 1))
cpt_mean <- cpt.mean(df_ts, method = "AMOC", penalty = "AIC")

# 提取变点核心信息
cpt_pos <- cpts(cpt_mean)  # 变点对应的索引位置
mean_segments <- param.est(cpt_mean)$mean  # 各分段的均值

# 构建分段均值的辅助数据框
seg_df <- data.frame(
  start_date = c(monthly_data$date[1], monthly_data$date[cpt_pos]),
  end_date = c(monthly_data$date[cpt_pos - 1], monthly_data$date[nrow(monthly_data)]),
  mean_val = mean_segments
)

# 绘制ggplot版变点图
p_changepoint <- ggplot(monthly_data, aes(x = date, y = number)) +
  geom_line(color = "steelblue") +
  geom_point(size = 2, color = "darkblue") +
  # 添加分段均值线
  geom_segment(data = seg_df, aes(x = start_date, xend = end_date, y = mean_val, yend = mean_val),
               color = "red", linewidth = 1.2) +
  # 添加变点竖线
  geom_vline(xintercept = monthly_data$date[cpt_pos], color = "darkred", linetype = "dashed", linewidth = 1) +
  labs(x = "日期", y = "数值", title = "变点分析结果") +
  theme_minimal()

# 直接与其他ggplot对象组合（示例）
# grid.arrange(p1, p2, p3, p_changepoint, ncol=2)

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