R: Plots Results of a Monte Carlo Simulation

plot.mc {mc2d} | R 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 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.

### 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.

### 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 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)

[Package mc2d version 0.1-18 Index]