plot.LP {MCPMod} | R Documentation |
Plot method for LP objects
Description
Graphically displays an LP object.
Usage
## S3 method for class 'LP'
plot(x, line = TRUE, type = NULL, spldf = 5, ...)
Arguments
x |
LP object as obtained from a call to the |
line |
Logical indicating whether the power values should be smoothed. |
type |
One of "LP1", "LP2" or "both", availability depending on whether
the corresponding values have been calculated in the call to
the |
spldf |
Numeric determining the degrees of freedom for the smoothing spline
which is plotted if |
... |
Additional arguments. |
Details
The function produces a trellis display of the loss in power for different values
of the standardized model parameter. A smoothing spline (with spldf
degrees of freedom) is fit to these
points to give a smooth impression of the loss in power curve.
For models with two prior parameters
a trellis display is shown with the number of panels equal to len[2]
.
The number of points on which the power is evaluated is equal to len[1]
in each panel, where len
is an argument of the LP
function.
See Also
Examples
## Not run:
doses <- c(0,10,25,50,100,150)
models <- list(linear=NULL, emax=c(25),
logistic=c(50,10.88111), exponential=c(85),
betaMod=matrix(c(0.33,2.31,1.39,1.39),byrow=TRUE,nrow=2))
# Examples from JBS paper, p.654
LPobj <- LP(models, model = "emax", type = "both", paramRange = c(10,70),
doses = doses, base = 0, maxEff = 0.4, sigma = 1, n = 60,
alpha = 0.05, len = 15, scal = 200)
plot(LPobj)
plot(LPobj, line = FALSE, type = "LP1")
plot(LPobj, type = "LP1", spldf = 9)
## End(Not run)