lmomco (version 2.3.1)

add.log.axis: Add a Polished log Axis to a Plot

Description

This function provides special support for adding superior looking base-10 logarithmic axes relative to R's defaults, which are an embarassment. The Examples section shows an overly elaborate version made by repeated calls to this function with a drawback that each call redraws the line of the axis so deletion in editing software might be required. This function is indexed under the “lmomco functions” because of its relation to add.lmomco.axis and is not named add.lmomcolog.axis because such a name is too long.

Usage

add.log.axis(make.labs=FALSE, logs=c(2, 3, 4, 5, 6, 7, 8, 9), side=1,
             two.sided=FALSE, label=NULL, ...)

Arguments

make.labs

A logical controlling whether the axis is labled according to the values in logs.

logs

A numeric vector of log-cycles for which ticking and (or) labeling is made. These are normalized to the first log-cycle, so a value of \(3\) would spawn values such as \(\cdots, 0.03, 0.3, 3, 30, \cdots\) through a range exceeding the axis limits. The default anticipates that a second call to the function will be used to make longer ticks at the even log-cycles; hence, the value 1 is not in the default vector. The Examples section provides a thorough demonstration.

side

An integer specifying which side of the plot the axis is to be drawn on, and argument corresponds to that in the Axis function. The axis is placed as follows: 1=below, 2=left, 3=above, and 4=right.

two.sided

A logical controlling whether the side oppose of side also is to be drawn.

label

The label (title) of the axis, which is placed by a call to function mtext, and thus either the xlab or ylab arguments for plot should be set to the empty string "".

...

Additional arguments to pass to Axis.

Value

No value is returned. This function is used for its side effects.

See Also

add.lmomco.axis

Examples

# NOT RUN {
par(mgp=c(3,0.5,0)) # going to tick to the inside, change some parameters
X <- 10^sort(rnorm(65)); pp <- pp(X) # generate synthetic data
plot(qnorm(pp), X, xaxt="n", yaxt="n", xlab="", ylab="", log="y",
     xlim=c(-2,3), ylim=c(0.01, 600), pch=6, yaxs="i", col=4)
add.lmomco.axis(las=2, tcl=0.5, side.type="RI", otherside.type="NPP")
# Logarithmic axis: the base ticks to show logarithms
add.log.axis(side=2,      tcl=0.8*abs(par()$tcl), two.sided=TRUE)
#                   the long even-cycle tick, set to inside and outside
add.log.axis(logs=c(1),   tcl=-0.5*abs(par()$tcl), side=2, two.sided=TRUE)
add.log.axis(logs=c(1),   tcl=+1.3*abs(par()$tcl), side=2, two.sided=TRUE)
#                   now a micro tick at the 1.5 logs but only on the right
add.log.axis(logs=c(1.5), tcl=+0.5*abs(par()$tcl), side=4)
#                   and only label the micro tick at 1.5 on the right
add.log.axis(logs=c(1.5), side=4, make.labs=TRUE, las=3)
#                   add the bulk tick labeling and axis label.
add.log.axis(logs=c(1, 2, 4, 6), side=2, make.labs=TRUE, las=1, label="QUANTILE")
par(mgp=c(3,1,0)) # restore defaults
# }