ecdf.Rd
makes a ecdf plot for enrichment test result returned by
enrichment_test
with test equals to kst.
ecdf(x, ...) # S3 method for EnrichmentKST ecdf( x, level, alternative = c("greater", "less"), colors = c("blue", "red"), ... )
x | EnrichmentKST. Returned by |
---|---|
... | other arguments not supported |
level | character. The level to plot the ecdf plot. Must be in the group. |
alternative | character. Either greater or less |
colors | character |
a ggplot2 object
enrichment_test
lpd = transform_by_sample(lipidome, function(x) log(x/sum(x))) design = model.matrix(~Treatment * Timepoint + Subject, data = lpd$pdata) fit = model_fit(lpd, design, "TreatmentMed:TimepointPre", "limma") en = enrichment_test(lpd, fit, "class", "kst") ecdf(en, level = "PC", alternative = "greater")