plot.ebnm: Plot an ebnm object

In ebnm: Solve the Empirical Bayes Normal Means Problem

Plot an ebnm object

Description

Given one or more fitted ebnm object(s), produces a plot of posterior means vs. observations. If desired, a plot of cumulative distribution functions of fitted prior(s) can also be produced.

Usage

## S3 method for class 'ebnm'
plot(
  x,
  ...,
  incl_pm = TRUE,
  incl_cdf = FALSE,
  subset = NULL,
  remove_abline = FALSE
)

Arguments

x	The fitted ebnm object.
...	Additional ebnm objects to be included on the same plots.
incl_pm	Plot posterior means vs. observations?
incl_cdf	Plot the cumulative distribution functions?
subset	The subset of observations to include on the plot of posterior means vs. observations. Can be a numeric vector corresponding to indices of observations to plot, or a character vector if observations are named. If subset = NULL then all observations will be plotted.
remove_abline	To better illustrate shrinkage effects, the plot of posterior means vs. observations includes the line y = x by default. If remove_abline = TRUE, then this line will not be drawn.

Examples

theta <- c(rep(0, 100), rexp(100))
theta[1:50] <- 0
s <- 1
x <- theta + rnorm(200, 0, s)
pn.res <- ebnm_point_normal(x, s)
plot(pn.res)

pe.res <- ebnm_point_exponential(x, s)
plot(pn.res, pe.res)

# Customize plot:
library(ggplot2)
plot(pn.res, pe.res, remove_abline = TRUE) +
  theme_bw() +
  labs(x = "Simulated data")

ebnm documentation built on Oct. 13, 2023, 1:16 a.m.