如何用gganimate制作样本量影响样本均值与标准误的动画？

问题说明 #

原animation包代码出现致命错误，核心原因是函数参数传递不匹配（plo函数仅定义了samplesize和iter两个参数，但调用时传入了三个参数），且基础绘图的循环逻辑和gganimate的数据驱动模式不兼容。以下是用gganimate实现相同效果的替代方案：

实现代码 #

# 首次运行先安装所需包
install.packages(c("tidyverse", "gganimate", "gifski"))

# 加载依赖包
library(tidyverse)
library(gganimate)

# 设定分布与动画参数
mu <- 10.0
sd <- 1.0
sample_sizes <- c(10, 50, 100, 500, 1000, 10000)
iterations <- 100  # 每个样本量的迭代次数

# 生成所有动画帧的数据集（设定随机种子保证结果可复现）
set.seed(123)
animation_data <- map_dfr(sample_sizes, function(n) {
  map_dfr(1:iterations, function(i) {
    # 生成正态分布样本
    sample_data <- rnorm(n, mean = mu, sd = sd)
    # 计算当前样本的均值和标准差
    sample_mean <- mean(sample_data)
    sample_sd <- sd(sample_data)
    # 整理为规范数据框
    tibble(
      x = sample_data,
      sample_size = n,
      iteration = i,
      mean_val = sample_mean,
      sd_val = round(sample_sd, 2)
    )
  })
})

# 绘制基础图形并添加动画逻辑
p <- ggplot(animation_data, aes(x = x)) +
  # 直方图：保持原代码的配色、分箱规则
  geom_histogram(
    breaks = seq(4, 16, 0.5),
    fill = rgb(0.1, 0.9, 0.1, 0.2),
    color = "grey",
    alpha = 0.8
  ) +
  # 添加样本均值竖线
  geom_vline(aes(xintercept = mean_val), color = "red", linewidth = 2) +
  # 固定坐标轴范围，保证动画切换时的视觉连贯性
  xlim(c(5, 15)) +
  ylim(0, max(floor(animation_data$sample_size / 3))) +
  # 动态标题：显示当前样本量、均值、标准差
  labs(
    x = "x (高斯样本)",
    title = "N = {closest_state$sample_size}, xbar = {round(closest_state$mean_val, 2)}, s = {closest_state$sd_val}"
  ) +
  # 控制动画帧切换逻辑：每个样本量下迭代100次
  transition_states(
    state = interaction(sample_size, iteration),
    transition_length = 0,
    state_length = 0.12  # 对应原代码的帧间隔参数
  ) +
  ease_aes("linear")

# 批量导出不同样本量对应的GIF
walk(sample_sizes, function(n) {
  subset_data <- filter(animation_data, sample_size == n)
  p_subset <- p %+% subset_data
  anim_save(
    filename = paste0("N=", n, ", mu=", mu, ", sd=", sd, ".gif"),
    animation = p_subset,
    width = 480,
    height = 320,
    renderer = gifski_renderer(loop = FALSE)
  )
})

核心要点 #

• 采用数据驱动模式生成所有帧的数据集，完全适配gganimate的设计逻辑，避免了基础绘图循环的错误
• 用transition_states控制动画帧切换，state_length参数对应原代码的interval，精确控制每帧停留时间
• 固定坐标轴范围，保证不同样本量动画切换时的视觉一致性
• 用walk函数批量导出GIF，和原代码的循环输出逻辑保持一致

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