不能使用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} $做一些无聊的玩杂耍，你可以看到它。换句话说：边际概率对允许的相关范围进行了限制（反之亦然），并且您在通话中违反了这一点。

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

问题应该转移到统计交换。