1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
|
\name{twoord.plot}
\alias{twoord.plot}
\title{ Plot with two ordinates }
\description{
Two sets of values are displayed on the same plot with different ordinate
scales on the left and right.
}
\usage{
twoord.plot(lx,ly,rx,ry,data=NULL,main="",xlim=NULL,lylim=NULL,rylim=NULL,
mar=c(5,4,4,4),lcol=1,rcol=2,xlab="",lytickpos=NA,ylab="",ylab.at=NA,
rytickpos=NA,rylab="",rylab.at=NA,lpch=1,rpch=2,
type="b",xtickpos=NULL,xticklab=NULL,halfwidth=0.4,axislab.cex=1,
do.first=NULL,xaxt="s",...)
}
\arguments{
\item{lx,ly,rx,ry}{y and optional x values for the plot}
\item{data}{an optional data frame from which to obtain the above values}
\item{main}{Title for the plot}
\item{xlim}{optional x limits as in \samp{plot}}
\item{lylim,rylim}{optional y limits for the left and right axes
respectively}
\item{mar}{optional margin adjustment, defaults to \samp{c(5,4,4,4)}}
\item{lcol,rcol}{colors to distinguish the two sets of values}
\item{xlab}{X axis label as in \samp{plot}}
\item{lytickpos}{Optional positions for the left axis tick labels.}
\item{ylab}{Left Y axis label as in \samp{plot}}
\item{ylab.at}{Optional position for the left Y axis label}
\item{rytickpos}{Optional positions for the right axis tick labels.}
\item{rylab}{Right Y axis label}
\item{rylab.at}{Optional position for the right Y axis label}
\item{lpch,rpch}{plot symbols to distinguish the two sets of values}
\item{type}{as in \samp{plot}}
\item{xtickpos}{Optional positions for x-axis tick labels.}
\item{xticklab}{Optional labels for x-axis. Useful for things like dates.}
\item{halfwidth}{Half the width of the bars in user units. The bars are
centered on successive integers if no \samp{x} values are supplied.}
\item{axislab.cex}{Character expansion for the axis labels and tick labels.}
\item{do.first}{Optional command(s) that will be executed immediately after
the blank plot is displayed.}
\item{xaxt}{Whether to display the x-axis - "n" = no.}
\item{...}{additional arguments passed to \samp{plot} and \samp{points}.}
}
\value{nil}
\details{
\samp{twoord.plot} automates the process of displaying two sets of
values that have different ranges on the same plot. It is principally
useful in illustrating some relationship between the values across the
observations. It is assumed that the \samp{lx} and \samp{rx} values
are at least adjacent, and probably overlapping.
It is best to pass all the arguments \samp{lx, ly, rx, ry}, but the
function will attempt to substitute sensible x values if one or two
are missing.
If at least one of the \samp{type} arguments is "bar", bars will be plotted
instead of points or lines. It is best to plot the bars first (i.e. relative
to the left axis) if the other type is points or lines, as the bars will
usually obscure at least some of the points or lines. Using NA for the color
of the bars will partially correct this. If both types are to be bars,
remember to pass somewhat different x values or the bars will be overplotted.
Note that more values can be added to the plot using \samp{points} or
\samp{lines}, but remember that these will be plotted relative to the left
ordinate.
The \samp{do.first} argument is useful for adding a background color or grid
to the plot as shown in the first two examples.
}
\note{
There are many objections to the use of plots with two different
ordinate scales, and some of them are even sensible and supported by
controlled observation. Many of the objections rest on assertions that the
spatial arrangement of the values plotted will override all other
evidence. Here are two:
The viewer will assume that the vertical position of the data points
indicates a quantitative relationship.
To some extent. It is probably not a good idea to have the spatial
relationship of the points opposed to their numerical relationship. That
is to say, if one set of values is in the range of 0-10 and the other
20-100, it is best to arrange the plot so that the latter values are
not plotted below the former. See the second example, which illustrates
a method for separating the two series and offsetting the axes.
The viewer will assume that an intersection of lines indicates an
intersection of values.
If the visual elements representing values can be arranged to avoid
intersections, so much the better. Many people have no trouble
distinguishing which visual elements are linked to which axis as long as
they are both coded similarly, usually with colors and/or symbols. In the
special case where there is an underlying relationship between the two
such as the probability of that value occurring under some conditions, it
may help to mark the point(s) where this occurs.
It may be useful to consider \samp{gap.plot} or separate plots as an
alternative.
}
\author{Jim Lemon (thanks to Christophe Dutang for the idea of using bars
and lines in the same plot, Clair Crossupton for pointing out that
dates on the x-axis weren't very good, Jacob Kasper for the axis
character expansion and Ye Lin for finally motivating me to add the
do.first argument.)}
\seealso{\link{plot}}
\examples{
xval1 <- seq.Date(as.Date("2017-01-02"),
as.Date("2017-01-10"), by="day")
xval2 <- seq.Date(as.Date("2017-01-01"),
as.Date("2017-01-15"), by="day")
going_up<-seq(3,7,by=0.5)+rnorm(9)
going_down<-rev(60:74)+rnorm(15)
twoord.plot(2:10,going_up,1:15,going_down,xlab="Sequence",
ylab="Ascending values",rylab="Descending values",lcol=4,
main="Plot with two ordinates - points and lines",
do.first="plot_bg();grid(col=\"white\",lty=1)")
axis.Date(1,xval2)
# now separate the lines
twoord.plot(2:10,going_up,1:15,going_down,xlab="Sequence",
lylim=range(going_up)+c(-1,10),rylim=range(going_down)+c(-10,2),
ylab="Ascending values",ylab.at=5,rylab="Descending values",
rylab.at=65,lcol=4,main="Plot with two ordinates - separated lines",
lytickpos=3:7,rytickpos=seq(55,75,by=5),
do.first="plot_bg();grid(col=\"white\",lty=1)")
twoord.plot(2:10,going_up,1:15,going_down,xlab="Sequence",
lylim=range(going_up)+c(-1,10),rylim=range(going_down)+c(-10,2),
type=c("bar","l"),ylab="Ascending values",ylab.at=5,
rylab="Descending values",rylab.at=65,
main="Bars on left axis, lines on right axis",
lytickpos=3:7,rytickpos=seq(55,75,by=5),
lcol=3,rcol=4,do.first="plot_bg()")
twoord.plot(2:10,going_up,1:15,going_down,xlab="Sequence",
lylim=c(-3,8),rylim=c(50,100),type=c("l","bar"),
ylab="Ascending values",rylab="Descending values",
lytickpos=3:7,rytickpos=seq(55,75,by=5),ylab.at=5,rylab.at=65,
main="Lines on left axis, bars on right axis",
lcol=3,rcol=4,do.first="plot_bg(\"yellow\")")
# histogram with density curve superimposed
xhist<-hist(rnorm(100),plot=FALSE)
xdens<-dnorm(seq(-3,3,by=0.05))
twoord.plot(xhist$mids,xhist$counts,seq(-3,3,by=0.05),
xdens,type=c("bar","l"),lcol=4,rcol=2,ylab="Counts",
rylab="Density",main="Histogram and density curve",
halfwidth=0.2,lylim=c(0,max(xhist$counts)+1),rylim=c(0,0.45),lwd=2)
}
\keyword{misc}
|
