plot.LP {MCPMod} | R Documentation |
Graphically displays an LP object.
## S3 method for class 'LP': plot(x, line = TRUE, type = NULL, spldf = 5, ...)
x |
LP object as obtained from a call to the LP function |
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 LP function. |
spldf |
Numeric determining the degrees of freedom for the smoothing spline
which is plotted if line==TRUE . Note that spldf should
be larger than 1 and smaller than len (default: 5). |
... |
Additional arguments. |
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.
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)