plot.mc: Plots Results of a Monte Carlo Simulation

View source: R/plot.mc.R

plot.mcR Documentation

Plots Results of a Monte Carlo Simulation

Description

Plots the empirical cumulative distribution function of a ‘⁠mcnode⁠’ or a ‘⁠mc⁠’ object ("0" and "V" nodes) or the empirical cumulative distribution function of the estimate of a ‘⁠mcnode⁠’ or ‘⁠mc⁠’ object ("U" and "VU" nodes).

Usage

## S3 method for class 'mc'
plot(x, prec=0.001, stat=c("median", "mean"), lim=c(0.025, 0.25, 0.75,
	  0.975), na.rm=TRUE, griddim=NULL, xlab=NULL, ylab="Fn(x)",
	  main="", draw=TRUE, paint=TRUE, xlim=NULL, ylim=NULL, ...)
## S3 method for class 'mcnode'
plot(x, ...)
## S3 method for class 'plotmc'
plot(x, ...)
## S3 method for class 'mccut'
plot(x, stat=c("median", "mean"), lim=c(0.025, 0.25, 0.75, 0.975),
	  griddim=NULL, xlab=names(x), ylab="Fn(x)", main="",
	  draw=TRUE, ...)

Arguments

x

a ‘⁠mcnode⁠’ or a ‘⁠mc⁠’ objects

prec

the precision of the plot. 0.001 will provide an ecdf from the 0.000, 0.001, .002, ..., 1.000 quantiles.

stat

the function used for estimates (2D ‘⁠mc⁠’ or ‘⁠mcnode⁠’). By default the median.

lim

a vector of numbers (between 0 and 1) indicating the envelope (2D ‘⁠mc⁠’ or ‘⁠mcnode⁠’) . Maybe ‘⁠NULL⁠’ or empty.

na.rm

Should NA values be discarded

griddim

a vector of two integers, indicating the size of the grid of the graph. If ‘⁠NULL⁠’, the grid is calculated to produce a "nice" graph.

xlab

vector of labels for the x-axis. If ‘⁠NULL⁠’, use the name of the node.

ylab

vector of labels for the y-axis.

main

vector of main titles of the graph.

draw

Should the plot be drawn?

paint

Should the envelopes be filled?

xlim

x coordinate range. ‘⁠xlim⁠’ is either a vector of length 2, used for each graph, or a list of vectors of length 2, whose ith element is used for the ith graph. By default, the data range is used as ‘⁠xlim⁠’.

ylim

y coordinate range. ‘⁠ylim⁠’ is either a vector of length 2, used for each graph, or a list of vectors of length 2, whose ith element is used for the ith graph. By default, the data range is 0-1.

...

further arguments to be passed to ‘⁠plot.stepfun⁠’.

Details

⁠plot.mcnode⁠’ is a user-friendly function that send the ‘⁠mcnode⁠’ to ‘⁠plot.mc⁠’.

For ‘⁠"VU"⁠’ and ‘⁠"U"⁠’ ‘⁠mcnode⁠’s, quantiles are calculated using quantile.mc within each of the ‘⁠nsu⁠’ simulations (i.e. by columns of each ‘⁠mcnode⁠’). The medians (but may be the means using ‘⁠stat="mean"⁠’) calculated from the ‘⁠nsu⁠’ values are plotted. The 0.025 and 0.975 quantiles, and the 0.25 and 0.75 quantiles (default values of ‘⁠lim⁠’) of these quantiles are used as the envelope.

Value

A ‘⁠plot.mc⁠’ object, list of the quantiles used to plot the draw.

References

Cullen AC and Frey HC (1999) Probabilistic techniques in exposure assessment. Plenum Press, USA, pp. 81-155.

See Also

ecdf, plot, quantile.mc

Examples

data(total)

plot(xVUM3)
## only one envelope corresponding to quantiles 0.025 and 0.975
plot(xVUM3, lim=c(0.025, 0.975)) 
## only one envelope not painted
plot(xVUM3, lim=c(0.025, 0.975), paint=FALSE) 

def.par <- par(no.readonly = TRUE)
par(mar=c(4, 4, 1, 1))
plot(total)
par(def.par)




mc2d documentation built on June 22, 2024, 10:54 a.m.