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
|
\name{dispersion}
\alias{dispersion}
\alias{dispbars}
\title{ Display a measure of dispersion. }
\description{
Display lines or capped bars at specified points on a plot representing
measures of dispersion.
}
\usage{
dispersion(x,y,ulim,llim=ulim,intervals=TRUE,arrow.cap=0.01,arrow.gap=NA,
type="a",fill=NA,lty=NA,pch=NA,border=NA,col=par("fg"),display.na=TRUE,
...)
}
\arguments{
\item{x,y}{x and y position of the centers of the bars}
\item{ulim,llim}{The extent of the dispersion measures.}
\item{arrow.cap}{The width of the cap at the outer end of each bar
as a proportion of the width of the plot.}
\item{arrow.gap}{The gap to leave at the inner end of each bar.
Defaults to two thirds of the height of a capital "O".}
\item{intervals}{Whether the limits are intervals (TRUE) or absolute values
(FALSE).}
\item{type}{What type of display to use.}
\item{fill}{Color to fill between the lines if \samp{type} is not NULL and
\samp{fill} is not NA.}
\item{lty}{Line type for redrawing the lines if necessary.}
\item{pch}{Symbol for redrawing the points if necessary.}
\item{border}{Line type for drawing a border on the confidence band.}
\item{col}{Color for the lines or capped bars.}
\item{display.na}{Whether to display NA values as lines going off the plot.}
\item{...}{additional arguments passed to \samp{arrows} or \samp{lines}
depending upon \samp{type}.}
}
\details{
\samp{dispersion} displays a measure of dispersion on an existing plot.
Currently it will display either vertical lines with caps (error bars) or lines
that form a "confidence band" around a line of central tendency. If \samp{fill}
is not NA and \samp{type} is \samp{"l"}, a polygon will be drawn between the
confidence lines. Remember that any points or lines within the confidence band
will be obscured, so pass point and/or line types as in the second example.
The default behavior is to display an undefined dispersion (e.g. a variance with
only one observation) as a line going off the plot. If \samp{display.na} is FALSE,
NA values will not be displayed, allowing the user to show only upper or lower
dispersion limits.
The \samp{intervals} argument allows the user to pass the limits as either
intervals (the default) or absolute values.
If \samp{arrow.gap} is greater than or equal to the upper or lower
limit for a bar, \samp{segments} is used to draw the upper and
lower caps with no bars to avoid zero length arrows.
}
\value{nil}
\author{Jim Lemon}
\seealso{\link{arrows}, \link{segments},\link{lines}}
\examples{
disptest<-matrix(rnorm(200),nrow=20)
disptest.means<-rowMeans(disptest)
row.order<-order(disptest.means)
se.disptest<-unlist(apply(disptest,1,std.error))
plot(disptest.means[row.order],main="Dispersion as error bars",
ylim=c(min(disptest.means-se.disptest),max(disptest.means+se.disptest)),
xlab="Occasion",ylab="Value")
dispersion(1:20,disptest.means[row.order],se.disptest[row.order])
plot(disptest.means[row.order],main="Dispersion as confidence band",
ylim=c(min(disptest.means-se.disptest),max(disptest.means+se.disptest)),
xlab="Occasion",ylab="Value")
dispersion(1:20,disptest.means[row.order],se.disptest[row.order],type="l",
fill="#eeccee",lty=2,pch=1)
disptest2<-matrix(sample(c(TRUE,FALSE),200,TRUE),nrow=10)
disptest.prop<-rowMeans(disptest2)
disptest.ulim<-disptest.llim<-rep(NA,10)
for(i in 1:10) {
disptest.ulim[i]<-binciWu(disptest2[i,],20)
disptest.llim[i]<-binciWl(disptest2[i,],20)
}
plot(disptest.prop,main="Dispersion as binomial confidence intervals",
ylim=c(min(disptest.llim),max(disptest.ulim)),
xlab="Sample",ylab="Proportion")
dispersion(1:10,disptest.prop,disptest.ulim,disptest.llim,
interval=FALSE,lty=2,pch=1)
}
\keyword{misc}
|
