plot.sim: Plot method for simulation 'sim' objects

In kkholst/lava: Latent Variable Models

Plot method for simulation 'sim' objects

Description

Density and scatter plots

Usage

## S3 method for class '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)

kkholst/lava documentation built on Feb. 22, 2024, 4:07 p.m.