MKpower (version 1.0)

volcano: Volcano Plots

Description

Produce volcano plot(s) for simulations of power and type-I-error of tests.

Usage

# S3 method for sim.power.ttest
volcano(x, alpha = 1, shape = 19, 
                                    hex = FALSE, bins = 50, ...)
# S3 method for sim.power.wtest
volcano(x, alpha = 1, shape = 19, 
                                    hex = FALSE, bins = 50, ...)

Value

Object of class gg and ggplot.

Arguments

x

object of class sim.power.ttest.

alpha

bleding factor (default: no blending.

shape

point shape used.

hex

logical, should hexagonal binning be used.

bins

number of bins used for hexagonal binning.

...

further arguments that may be passed through).

Author

Matthias Kohl

Details

The plot generates a ggplot2 object that is shown.

Missing values are handled by the ggplot2 functions.

References

Wikipedia contributors, Volcano plot (statistics), Wikipedia, The Free Encyclopedia, https://en.wikipedia.org/w/index.php?title=Volcano_plot_(statistics)&oldid=900217316 (accessed December 25, 2019).

For more sophisticated and flexible volcano plots see for instance: Blighe K, Rana S, Lewis M (2019). EnhancedVolcano: Publication-ready volcano plots with enhanced colouring and labeling. R/Bioconductor package. https://github.com/kevinblighe/EnhancedVolcano.

See Also

volcano

Examples

res1 <- sim.power.t.test(nx = 5, rx = rnorm, rx.H0 = rnorm, 
                        ny = 10, ry = function(x) rnorm(x, mean = 3, sd = 3), 
                        ry.H0 = function(x) rnorm(x, sd = 3))
volcano(res1)
# \donttest{
## low number of iterations to reduce computation time
res2 <- sim.power.wilcox.test(nx = 6, rx = rnorm, rx.H0 = rnorm,
                      ny = 6, ry = function(x) rnorm(x, mean = 2), 
                      ry.H0 = rnorm, iter = 100, conf.int = TRUE)
volcano(res2)
# }