不能用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}$做了一些无聊的杂耍，你就可以看到它。换句话说:边际概率限制了允许关联的范围(反之亦然)，而您在调用中违反了这一点。

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

这个问题可能应该转移到统计数据交换。