Dirichlet multinomial GLM likelihood ratio test for a single cluster (single overdispersion parameter)

Source: R/dm_glm.R

We recommend using dirichlet_multinomial_glm_anova_mc instead.

dirichlet_multinomial_anova(
  xFull,
  xNull,
  y,
  concShape = 1.0001,
  concRate = 1e-04
)

Arguments

xFull

[samples] x [covariates] matrix for the alternative model
xNull

[samples] x [covariates] matrix for the null model
y

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

Gamma shape parameter for concentration parameter
concRate

Gamma rate parameter for concentration parameter