Map scRNA-seq counts to a known covariance structure

CountMap(sigma, ngene, n, a1 = 3, b1 = 1, a20 = 2, b20 = 3, a30 = 1, b30 = 10)

Arguments

sigma

a covairance matrix

ngene

number of genes

n

number of samples

a20

the hyperparameter for generating beta, where beta~Gamma(a20,b20)

b20

the hyperparameter for generating beta, where beta~Gamma(a20,b20)

a10

the hyperparameter for generating alpha, where alpha~Gamma(a10,b10)

b10

the hyperparameter for generating alpha, where alpha~Gamma(a10,b10)

a3

the hyperparameter for generating pij, where pij~beta(a3, b3)

b3

the hyperparameter for generating pij, where pij~beta(a3, b3)

Value

a list of i) raw count data matrix; ii) adjacency matrix; iii) the true expression mean level theta matrix; iv) zij: the dropout matrix; v) sigma: the covariance matrix; vi) precision: the precision matrix.