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)
A "Topconfects" class object, as returned from limma_confects, edger_confects, etc.
Number if items to show.
c(lower, upper) limits on x-axis.
A ggplot2 object. Working non-interactively, you must print() this for it to be displayed.
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.
# 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)