Генерация коррелированных случайных чисел из биномиальных распределений в R

Я пытаюсь найти способ генерировать коррелированные случайные числа из нескольких биномиальных распределений.

Я знаю, как это сделать с нормальными распределениями (используя mvrnorm), но не нашел функции, применимой к биномиальным.