fields (version 16.3)

envelopePlot: Add a shaded the region between two functions to an existing plot

Description

This function shades the region vertically between two functions, specified as pairs of x and y vectors, and draws the functions in a darker shade. More formally, it shades all points (x,y) such that f1(x) < y < f2(x) or f2(x) < y < f1(x). When both functions have the same group of x values, the x values only need to be set once but y2 needs to be passed in by name. If the two functions intersect, the vertical space between the functions will be shaded on both sides, as implied in the definition above.

Usage

envelopePlot(x1, y1, x2 = x1, y2,
                         col ="thistle1" , lineCol = "thistle3", ...)

Arguments

x1

The x coordinates for the first function (or possibly both functions).

y1

The y coordinates for the first function.

x2

The x coordinates for the second function.

y2

The y coordinates for the second function.

col

The color to make the filling between the functions.

lineCol

The color to make the lines representing the functions.

...

Additional arguments to the base R function polygon

Author

Matt Iverson

Examples

x <- seq(0, 2*pi,, 100)
y1 <- cos(x)
y2 <- sin(x)
plot(x, y1, type="l")
envelopePlot(x, y1, y2=y2)

x1 <- c(0, 0.5, 1)
y1 <- c(0, 2, 1)
x2 <- c(0, 1)
y2 <- c(-1, 0)
plot(x1, y1, type="l", ylim=c(-1, 2))
envelopePlot(x1, y1, x2, y2)