生成两个具有预定义相关性的负二项分布随机变量

假设我有一个负二项分布变量 X1 和 NB(mu=MU1,size=s1) 和一个负二项分布变量 X2 和 NB(mu=MU2,size=s2)。我拟合了一个负二项式回归，以从我的数据中估计 Mu 和大小

我可以使用 R 中的 rnbinom() 函数从此分布生成随机抽取。

X1model<-rnbinom(n=1000,mu=MU1fitted,size=s1fitted)
X2model<-rnbinom(n=1000,mu=MU2fitted,size=s2fitted)

那些抽签现在是独立的。但是，我如何从这些分布中提取，以便它们表现出预定义的相关性 r，这是我观察到的原始数据 X1、X2 之间的相关性。

这样：

cor(X1,X2,method="spearman") = r = cor(X1model,X2model,method="spearman")

-甚至更好地从具有任意预设相关性的那些中抽取 r