如何在ggplot分面图中为对应子图添加专属回归直线

核心问题原因 #

你的代码中beta_mean数据框没有包含Continent列，导致ggplot无法将回归直线与分面匹配，因此每个子图都会绘制所有6条直线。

解决步骤 #

1. 为beta_mean添加大陆匹配列 #

首先需要将merged_data中的大陆因子水平与beta_mean的系数行一一对应，添加Continent列：

# 按merged_data中大陆的因子水平顺序，为beta_mean添加Continent列
beta_mean$Continent <- levels(merged_data$Continent)

这一步确保每一行系数对应一个大陆，与你Stan模型中beta_c的输出顺序保持一致。

2. 修改ggplot代码 #

无需额外调整geom_abline的映射，因为当beta_mean包含Continent列后，ggplot的分面机制会自动在每个子图中筛选出对应大陆的系数，只绘制单条回归直线：

ggplot(merged_data, 
       aes(x = log_wealth_per_capita,  
           y = log_total_deaths_per_million)) +
  geom_text(aes(label = iso3), size = 3) + 
  geom_abline(data = beta_mean, aes(intercept = Intercept, slope = Slope), 
              color = "red") +  # 现在会自动匹配对应分面的直线
  facet_wrap(~ Continent) +  
  labs(x = "Log(Wealth per Capita)", 
       y = "Log(Total Deaths per Million)") +  
  ggtitle("Scatterplots of Wealth per Capita vs. Total Deaths per Million, by Continent") +  
  theme_minimal()  

验证说明 #

执行上述步骤后，每个分面子图只会显示对应大陆的那条回归直线，因为beta_mean中的Continent列与分面变量完全匹配，ggplot会在每个分面中仅保留当前大陆的系数行用于绘制直线。

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