cosa (version 2.1.0)

plot: Power and MDES Curves

Description

Plots statistical power or minimum detectable effect size curves with (1-\(\alpha\))x100 % confidence interval for the design of interest.

Usage

# S3 method for power
plot(x, score = NULL, ypar = "mdes",  xpar = NULL,
          xlim = NULL, ylim = NULL,
          xlab = NULL, ylab = NULL,
          main = NULL, sub = NULL,
          locate = FALSE, benchmark = NULL, ...)
# S3 method for mdes
plot(x, score = NULL, ypar = "mdes",  xpar = NULL,
          xlim = NULL, ylim = NULL,
          xlab = NULL, ylab = NULL,
          main = NULL, sub = NULL,
          locate = FALSE, benchmark = NULL, ...)
# S3 method for cosa
plot(x, score = NULL, ypar = "mdes",  xpar = NULL,
          xlim = NULL, ylim = NULL,
          xlab = NULL, ylab = NULL,
          main = NULL, sub = NULL,
          locate = FALSE, benchmark = NULL, ...)

Arguments

x

an object returned from functions in cosa package.

score

vector or list; an empirical score variable or an object with class 'score' returned from the inspect.score() function.

ypar

character; "mdes" or "power" on y axis.

xpar

character; one of the sample sizes on x axis.

xlim

limits for xpar.

ylim

limits for ypar.

xlab

x axis label.

ylab

y axis label.

main

title for the plot.

sub

subtitle for the plot.

locate

logical; TRUE locates parameter values for design x on the plot.

benchmark

benchmark line.

...

other graphical parameters to pass to plot.new().

Examples

# NOT RUN {
d1 <- mdes.bcrd3r2(rho2 = .10, rho3 = .20, omega3 = .30,
                   n1 = 20, n2 = 44, n3 = 50)
plot(d1, xpar = "n3", xlim = c(30, 100))
# }