We recommend using dirichlet_multinomial_anova_mc instead.

dirichlet_multinomial_glm(x, y, concShape = 1.0001, concRate = 1e-04)

Arguments

x

[samples] x [covariates] matrix

y

[samples] x [introns] matrix of intron usage counts

concShape

Gamma shape parameter for concentration parameter

concRate

Gamma rate parameter for concentration parameter