Simulating single cell and bulk RNA-seq data accounting for gene-gene correlations

Question

[ali_karimnezhad]

I am thinking of simulating some data to evaluate a single cell proportion estimation method. To do so, I am looking to simulating two datasets: 1) a single cell gene expression count matrix for some given cell type proportions 2) a simulated bulk RNA-seq matrix for the same genes. In this simulation, I want to be able to control gene-gene correlations so that, for example, all genes in both simulated datasets (single cell and bulk RNA-seq) have strong/weak correlation. I appreciate for ideas.