BElikelihood (version 1.1)

plot.proLikelihood: Plot method for proLikelihood object

Description

This function generates a plot of a standardized profile likelihood after running the proLikelihood() function.

Usage

# S3 method for proLikelihood
plot(x, textx, texty = 0.9, textsize = 3, ...)

Value

ggplot2 object, a plot of the standardized profile likelihood with the maximum likelihood estimate and 1/8 and 1/32 likelihood intervals printed inside the plot.

Arguments

x

proLikelihood object

textx

numeric value, position (x-axis) of label for the maximum likelihood estimate and the 1/8 and 1/32 likelihood intervals.

texty

numeric value, position (y-axis) of label the maximum likelihood estimate and the 1/8 and 1/32 likelihood intervals.

textsize

numeric value text size of the label.

...

unused

Details

The function generates a plot of the standardized profile likelihood (the profile likelihood relative to the maximum) with the maximum likelihood estimate and 1/8 and 1/32 likelihood intervals for the parameter of interest (mean difference, total standard deviation ratio or within-subject standard deviation ratio depending on the ‘method’) printed inside the plot.

Examples

# \donttest{
data(dat)
cols <- list(subject = 'subject', formula = 'formula', y = 'AUC')
p4a <- averageBE(dat, colSpec = cols, xlength = 50)
p4t <- totalVarianceBE(dat, colSpec = cols, xlength = 50)
p4w <- withinVarianceBE(dat, colSpec = cols, xlength = 50)
plot(p4a)
plot(p4t)
plot(p4w)
# three period case
dd3 <- dat[dat$period < 4,]
p3a <- averageBE(dd3, colSpec = cols, xlength = 50)
plot(p3a)
# two period case
dd2 <- dat[dat$period < 3,]
p2a <- averageBE(dd2, colSpec = cols, xlength = 50)
plot(p2a)
# }