plot {cosa}

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 class '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 class '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 class '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

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