r - 从二项分布生成相关随机数

Question

我试图找到一种方法来从几个二项分布中生成相关的随机数。

我知道如何使用正态分布(使用 MASS::mvrnorm )，但我没有找到适用于二项式响应的函数。

Answer 1

您可以使用 copula 生成相关制服。包，然后使用 qbinom函数将它们转换为二项式变量。这是一个快速示例:

library(copula)

tmp <- normalCopula( 0.75, dim=2 )
x <- rcopula(tmp, 1000)
x2 <- cbind( qbinom(x[,1], 10, 0.5), qbinom(x[,2], 15, 0.7) )

现在 x2是一个矩阵，其中 2 列代表 2 个相关的二项式变量。