plot-Samples-ModelEff-method: Plot the fitted dose-effcacy curve using a model from...

In crmPack: Object-Oriented Implementation of CRM Designs

Plot the fitted dose-effcacy curve using a model from ModelEff class with samples

Description

Plot the fitted dose-effcacy curve using a model from ModelEff class with samples

Usage

## S4 method for signature 'Samples,ModelEff'
plot(
  x,
  y,
  data,
  ...,
  xlab = "Dose level",
  ylab = "Expected Efficacy",
  showLegend = TRUE
)

Arguments

x	the Samples object
y	the ModelEff model class object
data	the Data object
...	not used
xlab	the x axis label
ylab	the y axis label
showLegend	should the legend be shown? (default)

Value

This returns the ggplot object for the dose-efficacy model fit

Examples

## we need a data object with doses >= 1:
data <-DataDual(x=c(25,50,25,50,75,300,250,150),
                y=c(0,0,0,0,0,1,1,0),
                w=c(0.31,0.42,0.59,0.45,0.6,0.7,0.6,0.52),
                doseGrid=seq(25,300,25),
                placebo=FALSE)
##plot the dose-efficacy curve with samples using the model from 'ModelEff' 
##class e.g. 'Effloglog' class model
##define the model (see Effloglog example)
Effmodel<-Effloglog(Eff=c(1.223,2.513),Effdose=c(25,300),nu=c(a=1,b=0.025),data=data,c=0)
## define the samples obtained using the 'Effloglog' model (see details in 'Samples' example)
##options for MCMC
options<-McmcOptions(burnin=100,step=2,samples=200)
## samples must be of 'Samples' class
samples <- mcmc(data=data,model=Effmodel,options=options)
## plot the fitted dose-efficacy curve including the 95% credibility interval of the samples
## 'x' should be of 'Samples' class and 'y' of 'ModelEff' class
plot(x=samples,y=Effmodel,data=data)

crmPack documentation built on Sept. 3, 2022, 1:05 a.m.