如何计算生成相关二元变量的边际概率

Question

假设我想创建50个长度为100的二进制变量,每个变量彼此相关.

因为我用指定的rho创建了一个相关矩阵:

cor.mat <- matrix(0.2,nrow=50, ncol=50)
diag(cor.mat) <- 1

接下来我使用rmvbin:

library(bindata)
rmvbin<-rmvbin(100, margprob=x, bincorr=cor.mat)

但是,我不确定如何计算margprob参数.有人可以帮忙吗？

它应该是每行和每列概率之和的向量吗？

Answer 1

margprob应该只是任何单个二元变量为 1 的概率的重复向量，与其余变量无关；称此值为p。假设变量分布相同（给定您的相关矩阵似乎就是这种情况）， margprob=rep(p,50)。

它不应该是每行和每列的概率之和的向量，因为相关矩阵不能用于确定边际概率。如果您无法确定随机变量的边际概率是多少，则必须为问题提供更多背景信息，这将是一个更适合 math.stackexchange.com 的问题。