如何从联合区间均匀生成随机数？R语言是否有内置函数？

通用思路 #

要从多个联合区间中均匀生成随机数，核心是保证每个区间被选中的概率和其长度成正比，具体步骤如下：

• 计算每个区间的长度，再求和得到所有区间的总长度；
• 以「单区间长度/总长度」为概率，随机选中一个目标区间；
• 在选中的区间内，生成标准均匀分布的随机数。

R语言实现（针对区间[-1.5, -0.5] ∪ [0.5, 1.5]） #

R语言没有直接处理多联合区间均匀随机数的内置函数，但可以通过简单逻辑实现，以下提供两种方案：

方案一：通用型实现（适配任意数量/长度的区间） #

这个方法不管区间数量多少、长度是否一致都能使用：

# 定义目标联合区间
intervals <- list(c(-1.5, -0.5), c(0.5, 1.5))
# 计算每个区间的长度
interval_lengths <- sapply(intervals, function(x) x[2] - x[1])
# 计算各区间的抽样概率
probabilities <- interval_lengths / sum(interval_lengths)

# 封装成生成n个随机数的函数
sample_uniform_union <- function(n) {
  # 随机选择要生成数的区间
  selected_intervals <- sample(1:length(intervals), n, replace = TRUE, prob = probabilities)
  # 在选中的区间内生成均匀随机数
  sapply(selected_intervals, function(i) {
    runif(1, min = intervals[[i]][1], max = intervals[[i]][2])
  })
}

# 示例：生成10个符合要求的随机数
sample_uniform_union(10)

方案二：针对本例的简化实现 #

由于本例中两个区间长度相等（均为1），总长度为2，可通过坐标变换简化操作：

# 生成n个目标随机数的函数
sample_symmetric_union <- function(n) {
  # 先生成[0,2)区间的均匀随机数
  base_num <- runif(n, 0, 2)
  # 将[0,1)映射到[-1.5,-0.5]，[1,2)映射到[0.5,1.5]
  ifelse(base_num < 1, base_num - 1.5, base_num - 0.5)
}

# 示例：生成10个符合要求的随机数
sample_symmetric_union(10)

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