mc2d (version 0.1-17)

plot.mc: 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

"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, ...)
"plot"(x, ...)
"plot"(x, ...)
"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 enveloppe (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 enveloppes 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.

Value

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

Details

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

For "VU" and "U" mcnodes, 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 enveloppe.

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 enveloppe corresponding to quantiles 0.025 and 0.975
plot(xVUM3, lim=c(0.025, 0.975)) 
## only one enveloppe 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)