无法使用rmvbin从R中的二项分布生成相关随机数
Can't generate correlated random numbers from binomial distributions in R using rmvbin
我正在尝试从 R 中的二项分布中获取相关随机数的样本。我尝试使用 rmvbin，它在某些概率下运行良好：
> rmvbin(100, margprob = c(0.1,0.1), bincorr=0.5*diag(2)+0.5)
而下一个非常相似的调用会引发错误：
> rmvbin(100, margprob = c(0.01,0.01), bincorr=0.5*diag(2)+0.5)
Error in commonprob2sigma(commonprob, simulvals) :
Extrapolation occurred ... margprob and commonprob not compatible?
我找不到任何理由。
这是一个 math/stats "problem" 而不是 R 问题（从某种意义上说，这不是问题而是模型的结果）
简短版本：对于双变量二进制数据，可以观察到边际概率和相关性之间存在 link。如果您对边际概率 $p_A$ 和 $p_B$ 以及同时概率 $p_{AB}$ 进行一些无聊的处理，您就会看到它。换句话说：边际概率限制了允许相关性的范围（反之亦然），而您在调用中违反了这一点。
对于双变量高斯随机变量，边际和相关性是分开的，可以相互独立指定。
这个问题可能应该转移到 stats exchange。
我正在尝试从 R 中的二项分布中获取相关随机数的样本。我尝试使用 rmvbin，它在某些概率下运行良好：
> rmvbin(100, margprob = c(0.1,0.1), bincorr=0.5*diag(2)+0.5)
而下一个非常相似的调用会引发错误：
> rmvbin(100, margprob = c(0.01,0.01), bincorr=0.5*diag(2)+0.5)
Error in commonprob2sigma(commonprob, simulvals) :
Extrapolation occurred ... margprob and commonprob not compatible?
我找不到任何理由。
这是一个 math/stats "problem" 而不是 R 问题（从某种意义上说，这不是问题而是模型的结果）
简短版本：对于双变量二进制数据，可以观察到边际概率和相关性之间存在 link。如果您对边际概率 $p_A$ 和 $p_B$ 以及同时概率 $p_{AB}$ 进行一些无聊的处理，您就会看到它。换句话说：边际概率限制了允许相关性的范围（反之亦然），而您在调用中违反了这一点。
对于双变量高斯随机变量，边际和相关性是分开的，可以相互独立指定。
这个问题可能应该转移到 stats exchange。