seas.var.plot: Plot seasonal normals of a given variable

In seas: Seasonal Analysis and Graphics, Especially for Climatology

Plot seasonal normals of a given variable

Description

Plot seasonal normals of a variable using boxplots.

Usage

seas.var.plot(x, var, width = 11, start = 1, rep = 0, start.day = 1,
              col, ylim, add.alt, alt.ylab, main, ylab, ylog, ...)

Arguments

x	a data.frame with Date and var columns of data; x can also have id or name attributes to help give a title for the plot
var	a variable; a column name in x; this can also have attributes of units and long.name to help give a title for the y-axis
width	size of bin; see mkseas
start	starting bin number; e.g., if width="mon" and start=5, the plot will start on "May" at the left-hand side
rep	repetition of the bins in the boxplot
start.day	when width is numeric, this is the starting day of the year for the fist bin, or it can be a Date to specify a month and day (year is ignored)
col	colour for the boxplots; the default is "lightgrey"
ylim	c(min, max) range for y-axis
add.alt	this adds an alternative axis, and is specified by c(slope,inter); for example, if the primary measure is in °C, a secondary scale in K would be c(1,273.15), or in °F would be c(5/9,32); if ylog=TRUE, then this can also be TRUE to display the log10 transformed values of var on the alternative axis
alt.ylab	label for the alternate y-axis (the primary y-axis label is set through attributes for var in x)
main	title for plot; if it is missing, then it will automatically be generated
ylab	y-axis label; if it is missing, then it will automatically be generated
ylog	used to log10 transform values of var for the boxplots; this has a similar but different affect than specifying par(ylog=TRUE) before this function
...	ignored

Details

Shows normals of a seasonal variable using boxplots.

Value

Returns values from boxplot statistics on the variable.

Note

This function was formerly named plot.seas.var, but required renaming as it is not an S3 method.

Author(s)

Mike Toews

See Also

seas.var.plot, plot.seas.norm, year.plot.

Use mksub to make a subset of x.

Examples

opar <- par(no.readonly=FALSE)
on.exit(par(opar))
data(mscdata)
dat <- mksub(mscdata, id=1108447)

seas.var.plot(dat, var="t_max", col="tomato",
  add.alt=c(5/9, 32), alt.ylab="F")
abline(h=0)

par(opar) # reset graphics parameters

seas.var.plot(dat, var="t_min",
  start=18, rep=16)

pdat <- dat[dat$precip > 0,]
attr(pdat$precip, "long.name") <- "precipitation intensity"
attr(pdat$precip, "units") <- "mm/day"

par(ylog=TRUE)
seas.var.plot(pdat, var="precip", col="azure")
title(sub="These boxplots are simply plotted on a log-y scale")

par(opar)

seas.var.plot(pdat, var="precip", col="azure", ylog=TRUE)
title(sub="These boxplots are based on log-transformed values")

seas.var.plot(pdat, var="precip", col="azure", ylog=TRUE, add.alt=TRUE)
title(sub="The actual axis for graph is on the right-side")

seas documentation built on May 2, 2022, 5:08 p.m.