无法使用 rmvbin 从 R 中的二项式分布生成相关随机数

Question

我正在尝试从 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?

我找不到任何理由。

Answer 1

这是一个数学/统计“问题”，而不是 R 问题（从某种意义上说，这不是问题，而是模型的结果）

简短版：对于二元二进制数据，边际概率与可观察到的相关性之间存在联系。如果您对边际概率 $p_A$ 和 $p_B$ 以及同时概率 $p_{AB}$ 进行一些无聊的处理，您就会看到它。换句话说：边际概率限制了允许的相关范围（反之亦然），而您在通话中违反了这一点。

对于二元高斯随机变量，边际和相关性是分开的，可以相互独立地指定。

这个问题可能应该转移到 stats exchange。