plot.emaxsimobj: Plot dose response from a data set generated by emaxsim

In clinDR: Simulation and Analysis Tools for Clinical Dose Response Modeling

Plot dose response from a data set generated by emaxsim

Description

Plot of population dose response curve, dose group means with CIs, predictive intervals, and the model-based estimated dose response curve.

Usage

## S3 method for class 'emaxsimobj'
plot(
    x, xlim, xat=NULL, ylim, xlab, ylab,
    plotDif=FALSE,
    plotResid=FALSE,
    clev = 0.9, 
    plotPop=c('m','3','4'),     
    negC = FALSE, 
    logScale=FALSE,
    predict=TRUE, 
    plot=TRUE, ...)

Arguments

x	Extracted data object from emaxsim
xlim	x-axis limits
xat	The points at which tick-marks are to be drawn. Errors occur if the points are outside the range of xlim. By default (when NULL) tickmark locations are computed.
ylim	y-axis limits
xlab	x-axis label
ylab	y-axis label
plotDif	When TRUE, the difference with placebo is plotted.
plotResid	When TRUE, residuals (dose group means) are plotted.
clev	Level for confidence intervals
plotPop	Plot population dose response curve when plotPop='m' using linear interpolation between population means, when PlotPop='3' or '4', using the population Emax parameters that must be supplied with the data generator (see FixedMean or RandEmax). If the Emax parameters are not available, linear interpolation is used.
negC	If the ED50<lower ED50 limit, TRUE causes the Emax model to be plotted in addition to the alternative model selected.
logScale	If TRUE, log scale is used for dose.
predict	When TRUE, predictive intervals are plotted with grey errorbars in addition to the confidence intervals.
plot	Return plotting output without plotting.
...	Other plot parameters (not used).

Value

ggplot object is returned

Author(s)

Neal Thomas

See Also

emaxsim, print.emaxsimobj, summary.emaxsimobj, update.emaxsimobj

Examples

## Not run: 
## emaxsim changes the random number seed

nsim<-50
idmax<-5
doselev<-c(0,5,25,50,100)
n<-c(78,81,81,81,77)

### population parameters for simulation
e0<-2.465375 
ed50<-67.481113 

dtarget<-100
diftarget<-9.032497
emax<-solveEmax(diftarget,dtarget,log(ed50),1,e0)

sdy<-7.967897
pop<-c(log(ed50),emax,e0)    
meanlev<-emaxfun(doselev,pop)  

###FixedMean is specialized constructor function for emaxsim
gen.parm<-FixedMean(n,doselev,meanlev,sdy)  
D1 <- emaxsim(nsim,gen.parm)
e49<-D1[49]

plot(e49,clev=0.8)  

## End(Not run)

clinDR documentation built on Aug. 9, 2023, 9:08 a.m.