susie_plot_changepoint: Plot changepoint data and susie fit using ggplot2

In susieR: Sum of Single Effects Linear Regression

Plot changepoint data and susie fit using ggplot2

Description

Plots original data, y, overlaid with line showing susie fitted value and shaded rectangles showing credible sets for changepoint locations.

Usage

susie_plot_changepoint(
  s,
  y,
  line_col = "blue",
  line_size = 1.5,
  cs_col = "red"
)

Arguments

s	A susie fit generated by susie_trendfilter(y,order = 0).
y	An n-vector of observations that are ordered in time or space (assumed equally-spaced).
line_col	Color for the line showing fitted values.
line_size	Size of the lines showing fitted values
cs_col	Color of the shaded rectangles showing credible sets.

Value

A ggplot2 plot object.

Examples

set.seed(1)
mu = c(rep(0,50),rep(1,50),rep(3,50),rep(-2,50),rep(0,300))
y = mu + rnorm(500)
# Here we use a less sensitive tolerance so that the example takes
# less time; in practice you will likely want to use a more stringent
# setting such as tol = 0.001.
s = susie_trendfilter(y,tol = 0.1)

# Produces ggplot with credible sets for changepoints.
susie_plot_changepoint(s,y) 

susieR documentation built on March 7, 2023, 6:11 p.m.