gear (version 0.3.4)

plot.evgram: Plot evgram object

Description

Plots evgram object produced by evgram using the plot function.

Usage

# S3 method for evgram
plot(x, ..., split = FALSE, add_legend = TRUE, args_legend = list())

Arguments

x

An evgram object produced by the evgram function.

...

Additional arguments to pass the plot function to change aspects of the plot.

split

A logical value indicating whether, for a directional semivariogram, the directional semivariograms should be displayed in a single or split panels. Default is FALSE, for a single panel.

add_legend

A logical value indicating whether a legend should be included for a plot of directional semivariograms when split = FALSE.

args_legend

A list with arguments for legend.

Details

The arguments after ... are only considered if a directional semivariogram is provided, i.e., when x$ndir != 1.

split will split directional semivariograms into different plots automatically using autosize.

add_legend and args_legend are only used for a directional semivariogram when split = FALSE

See Also

xyplot, evgram

Examples

# NOT RUN {
data(co)
# omnidirectional example
v = evgram(Al ~ 1, co, ~ easting + northing)
plot(v)
plot(v, main = "semivariogram of Al")

# directional semivariograms overlaid
v2 = evgram(Al ~ 1, co, ~ easting + northing,
            angle = 22.5, ndir = 4)
plot(v2)
plot(v2, ylab = "semi-variance", pch = 2:5, type = "p")
plot(v2, lty = 2:5, type = "l", 
     args_legend = list(x = "bottomright",
                        legend = c("22.5", "67.5", "112.5", "157.5")))
# directional semivariograms split
plot(v2, split = TRUE)
plot(v2, split = TRUE, col = 2, pch = 3, type = "b",
     main = c("(a)", "(b)", "(c)", "(d)"))
# }