在R语言survey包的svyplot()中为气泡填充颜色的方法

嘿，我刚好在使用survey包处理分层样本可视化时踩过类似的坑，让我来帮你解决气泡填充颜色的问题！

首先，你遇到的pch参数不生效、无法设置填充色的核心原因是：你可能用了不支持填充的符号类型。R里的pch符号分两类：一类是空心符号（比如1-10），它们没有填充属性；另一类是可填充符号（21-25），这类符号可以分别用bg设置填充色，col设置轮廓色。

结合survey包的分层样本数据，我给你一套完整的解决方案：

1. 准备数据与调查对象 #

先模拟一组分层抽样数据，创建survey设计对象：

library(survey)

# 模拟分层抽样数据
set.seed(123)
sample_data <- data.frame(
  stratum = rep(c("城市", "乡村", "郊区"), each = 50),
  x = rnorm(150),
  y = rnorm(150),
  sampling_weight = runif(150, 1, 10)  # 抽样权重，范围1-10
)

# 创建分层抽样调查设计
svy_design <- svydesign(
  ids = ~1,
  strata = ~stratum,
  weights = ~sampling_weight,
  data = sample_data
)

2. 绘制带填充色的权重气泡图 #

使用svyplot函数，指定支持填充的pch=21（填充圆形），然后用bg设置填充色，col设置轮廓色，同时用sqrt(权重)来映射气泡大小（保证面积和权重成正比）：

# 基础填充色版本
svyplot(y ~ x, design = svy_design,
        pch = 21,  # 关键：选可填充的圆形符号
        bg = "#4287f5",  # 填充颜色（这里用蓝色）
        col = "black",  # 气泡轮廓颜色
        cex = sqrt(sampling_weight)/2,  # 面积正比于权重，所以用sqrt转换
        main = "分层样本权重气泡图",
        xlab = "变量X", ylab = "变量Y"
)

3. 进阶：按分层设置不同填充色 #

如果想给不同分层的气泡设置不同颜色，只需要把bg参数替换成分层对应的颜色向量即可：

# 定义分层对应的颜色
stratum_colors <- c("城市" = "#4287f5", "乡村" = "#f54242", "郊区" = "#42f56f")

# 按分层着色的气泡图
svyplot(y ~ x, design = svy_design,
        pch = 21,
        bg = stratum_colors[sample_data$stratum],  # 按分层分配填充色
        col = "black",
        cex = sqrt(sampling_weight)/2,
        main = "按分层着色的权重气泡图",
        xlab = "变量X", ylab = "变量Y"
)

关键点总结 #

• 必须使用pch=21-25这类可填充符号，其中21是最常用的填充圆形；如果用1-10的空心符号，bg参数会完全失效
• 填充色用bg参数设置，col参数控制的是气泡的轮廓颜色
• 气泡面积要和权重成正比，所以cex需要用sqrt(权重)计算（因为面积是半径的平方，sqrt转换后，cex的平方刚好和权重成正比）

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