Create a ggplot2 object showing the confect, effect, and average expression level of top features in a Topconfects object.

confects_plot(confects, n = 50, limits = NULL)

Arguments

confects

A "Topconfects" class object, as returned from limma_confects, edger_confects, etc.

n

Number if items to show.

limits

c(lower, upper) limits on x-axis.

Value

A ggplot2 object. Working non-interactively, you must print() this for it to be displayed.

Details

For each gene, the estimated effect is shown as a dot. The confidence bound is shown as a line to positive or negative infinity, showing the set of non-rejected effect sizes for the feature.

Examples

# Generate some random effect sizes with random accuracies n <- 100 effect <- rnorm(n, sd=2) se <- rchisq(n, df=3)^-0.5 # Find top confident effect sizes confects <- normal_confects(effect, se) # Plot top confident effect sizes confects_plot(confects, n=30)