ezsim (version 0.5.5)

plot.summary.ezsim: Plot an summary.ezsim Object

Description

Plot the summary statistics for several estimators in the same plot. Summary statistics abd estimators are separated by colour and linetype. The longest scalars parameter will be the x-variable of the plot. The rest of the scalars parameters will be become the facets of the plot (see ggplot2).Banker parameters will not be shown in the graph.

Usage

"plot"(x, parameters_priority, ylab = "Summary Statistics", title, pdf_option, return_print, ...)

Arguments

x

An summary.ezsim Object

parameters_priority

Display priority of parameter. Any missed parameters will be sorted by length.

ylab

Label of y-axis

title

Title of the plot

return_print

If TRUE, return a list of ggplot2 object. If FALSE(default), all of the plot will be printed out.

pdf_option

A list of option pass to pdf. If it is not missing, the plot will export to a pdf file

...

unused

Value

Optional: a ggplot2 object

See Also

ezsim,summary.ezsim

Examples

## Not run: 
# ezsim_basic<-ezsim(
#     m             = 100,
#     run           = TRUE,
#     display_name  = c(mean_hat="hat(mu)",sd_mean_hat="hat(sigma[hat(mu)])"),
#     parameter_def = createParDef(list(n=seq(20,80,20),mu=c(0,2),sigma=c(1,3,5))),
#     dgp           = function() rnorm(n,mu,sigma),
#     estimator     = function(x) c(mean_hat = mean(x),
#                                  sd_mean_hat=sd(x)/sqrt(length(x)-1)),
#     true_value    = function() c(mu, sigma / sqrt(n-1))
# )
# ## Plot the summary ezsim
# plot(summary(ezsim_basic,c("q25","q75")))
# plot(summary(ezsim_basic,c("q25","q75"),subset=list(estimator='mean_hat')))
# plot(summary(ezsim_basic,c("median"),subset=list(estimator='sd_mean_hat')))
# ## End(Not run)