Density and scatter plots

# S3 method for sim
plot(
  x,
  estimate,
  se = NULL,
  true = NULL,
  names = NULL,
  auto.layout = TRUE,
  byrow = FALSE,
  type = "p",
  ask = grDevices::dev.interactive(),
  col = c("gray60", "orange", "darkblue", "seagreen", "darkred"),
  pch = 16,
  cex = 0.5,
  lty = 1,
  lwd = 0.3,
  legend,
  legendpos = "topleft",
  cex.legend = 0.8,
  plot.type = c("multiple", "single"),
  polygon = TRUE,
  density = 0,
  angle = -45,
  cex.axis = 0.8,
  alpha = 0.2,
  main,
  cex.main = 1,
  equal = FALSE,
  delta = 1.15,
  ylim = NULL,
  xlim = NULL,
  ylab = "",
  xlab = "",
  rug = FALSE,
  rug.alpha = 0.5,
  line.col = scatter.col,
  line.lwd = 1,
  line.lty = 1,
  line.alpha = 1,
  scatter.ylab = "Estimate",
  scatter.ylim = NULL,
  scatter.xlim = NULL,
  scatter.alpha = 0.5,
  scatter.col = col,
  border = col,
  true.lty = 2,
  true.col = "gray70",
  true.lwd = 1.2,
  density.plot = TRUE,
  scatter.plot = FALSE,
  running.mean = scatter.plot,
  ...
)

Arguments

x

sim object

estimate

columns with estimates

se

columns with standard error estimates

true

(optional) vector of true parameter values

names

(optional) names of estimates

auto.layout

Auto layout (default TRUE)

byrow

Add new plots to layout by row

type

plot type

ask

if TRUE user is asked for input, before a new figure is drawn

col

colour (for each estimate)

pch

plot symbol

cex

point size

lty

line type

lwd

line width

legend

legend

legendpos

legend position

cex.legend

size of legend text

plot.type

'single' or 'multiple' (default)

polygon

if TRUE fill the density estimates with colour

density

if non-zero add shading lines to polygon

angle

shading lines angle of polygon

cex.axis

Font size on axis

alpha

Semi-transparent level (1: non-transparent, 0: full)

main

Main title

cex.main

Size of title font

equal

Same x-axis and y-axis for all plots

delta

Controls the amount of space around axis limits

ylim

y-axis limits

xlim

x-axis limits

ylab

y axis label

xlab

x axis label

rug

if TRUE add rug representation of data to x-axis

rug.alpha

rug semi-transparency level

line.col

line colour (running mean, only for scatter plots)

line.lwd

line width (running mean, only for scatter plots)

line.lty

line type (running mean, only for scatter plots)

line.alpha

line transparency

scatter.ylab

y label for density plots

scatter.ylim

y-axis limits for density plots

scatter.xlim

x-axis limits for density plots

scatter.alpha

semi-transparency of scatter plot

scatter.col

scatter plot colour

border

border colour of density estimates

true.lty

true parameter estimate line type

true.col

true parameter colour

true.lwd

true parameter line width

density.plot

if TRUE add density plot

scatter.plot

if TRUE add scatter plot

running.mean

if TRUE add running average estimate to scatter plot

...

additional arguments to lower level functions

Examples

n <- 1000
val <- cbind(est1=rnorm(n,sd=1),est2=rnorm(n,sd=0.2),est3=rnorm(n,1,sd=0.5),
             sd1=runif(n,0.8,1.2),sd2=runif(n,0.1,0.3),sd3=runif(n,0.25,0.75))

plot.sim(val,estimate=c(1,2),true=c(0,0),se=c(4,5),equal=TRUE,scatter.plot=TRUE)

plot.sim(val,estimate=c(1,3),true=c(0,1),se=c(4,6),xlim=c(-3,3),
  scatter.ylim=c(-3,3),scatter.plot=TRUE)

plot.sim(val,estimate=c(1,2),true=c(0,0),se=c(4,5),equal=TRUE,
  plot.type="single",scatter.plot=TRUE)

plot.sim(val,estimate=c(1),se=c(4,5,6),plot.type="single",scatter.plot=TRUE)

plot.sim(val,estimate=c(1,2,3),equal=TRUE,scatter.plot=TRUE)

plot.sim(val,estimate=c(1,2,3),equal=TRUE,byrow=TRUE,scatter.plot=TRUE)

plot.sim(val,estimate=c(1,2,3),plot.type="single",scatter.plot=TRUE)

plot.sim(val,estimate=1,se=c(3,4,5),plot.type="single",scatter.plot=TRUE)


density.sim(val,estimate=c(1,2,3),density=c(0,10,10), lwd=2, angle=c(0,45,-45),cex.legend=1.3)