pylimma.goana_trend

pylimma.goana_trend(index_de, covariate, n_prior=10, plot=False)[source]

Estimate per-gene DE probability from a covariate.

Port of R limma’s goanaTrend (goanaTrend.R). Used internally by goana() and kegga() for length / abundance bias correction. Stand-alone use is supported and validated against R.

Return type:

ndarray