Description Usage Arguments Details Value Methods Note Author(s) References See Also Examples

Generic function for producing kernel density plots.

1 2 3 4 | ```
simDensityplot(x, ...)
## S4 method for signature 'SimResults'
simDensityplot(x, true = NULL, epsilon, NArate, select, ...)
``` |

`x` |
the object to be plotted. For plotting simulation results, this
must be an object of class |

`true` |
a numeric vector giving the true values. If supplied, reference lines are drawn in the corresponding panels. |

`epsilon` |
a numeric vector specifying contamination levels. If supplied, the values corresponding to these contamination levels are extracted from the simulation results and plotted. |

`NArate` |
a numeric vector specifying missing value rates. If supplied, the values corresponding to these missing value rates are extracted from the simulation results and plotted. |

`select` |
a character vector specifying the columns to be plotted. It
must be a subset of the |

`...` |
additional arguments to be passed down to methods and eventually
to |

For simulation results with multiple contamination levels or missing value rates, conditional kernel density plots are produced.

An object of class `"trellis"`

. The
`update`

method can be used to update
components of the object and the `print`

method (usually called by default) will plot it on an appropriate plotting
device.

`x = "SimResults"`

produce kernel density plots of simulation results.

Functionality for producing conditional kernel density plots was added in version 0.2. Prior to that, the function gave an error message if simulation results with multiple contamination levels or missing value rates were supplied.

Andreas Alfons

Alfons, A., Templ, M. and Filzmoser, P. (2010) An Object-Oriented Framework for
Statistical Simulation: The **R** Package simFrame. *Journal of
Statistical Software*, **37**(3), 1–36. URL
http://www.jstatsoft.org/v37/i03/.

`simBwplot`

, `simXyplot`

,
`densityplot`

,
`"SimResults"`

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 | ```
#### design-based simulation
set.seed(12345) # for reproducibility
data(eusilcP) # load data
## control objects for sampling and contamination
sc <- SampleControl(size = 500, k = 50)
cc <- DARContControl(target = "eqIncome", epsilon = 0.02,
fun = function(x) x * 25)
## function for simulation runs
sim <- function(x) {
c(mean = mean(x$eqIncome), trimmed = mean(x$eqIncome, 0.02))
}
## run simulation
results <- runSimulation(eusilcP,
sc, contControl = cc, fun = sim)
## plot results
tv <- mean(eusilcP$eqIncome) # true population mean
simDensityplot(results, true = tv)
#### model-based simulation
set.seed(12345) # for reproducibility
## function for generating data
rgnorm <- function(n, means) {
group <- sample(1:2, n, replace=TRUE)
data.frame(group=group, value=rnorm(n) + means[group])
}
## control objects for data generation and contamination
means <- c(0, 0.25)
dc <- DataControl(size = 500, distribution = rgnorm,
dots = list(means = means))
cc <- DCARContControl(target = "value",
epsilon = 0.02, dots = list(mean = 15))
## function for simulation runs
sim <- function(x) {
c(mean = mean(x$value),
trimmed = mean(x$value, trim = 0.02),
median = median(x$value))
}
## run simulation
results <- runSimulation(dc, nrep = 50,
contControl = cc, design = "group", fun = sim)
## plot results
simDensityplot(results, true = means)
``` |

simFrame documentation built on May 30, 2017, 12:25 a.m.

Embedding an R snippet on your website

Add the following code to your website.

For more information on customizing the embed code, read Embedding Snippets.