plot.emaxsimB {clinDR}

Plot the output of emaxsimB

Description

A Q-Q plot of the posterior mean of the mean dose response at a specified dose minus the population value divided by the posterior SD of the mean difference.

Usage

## S3 method for class 'emaxsimB'
plot(x, id = x$idmax, plotDif= TRUE, ...)

Arguments

x

Output of emaxsimB

id

Index of the dose to be assessed (placebo index=1).

plotDif

If true (default), the estimates and population values are differences with placebo. IF false, absolute dose response values are used.

...

Optional parameters passed to the plotting function

Value

ggplot object is returned

Author(s)

Neal Thomas

See Also

emaxsimB, print.emaxsimB, summary.emaxsimB

Examples

## Not run: 
## emaxsimB changes the random number seeds
nsim<-50
idmax<-5
doselev<-c(0,5,25,50,100)
n<-c(78,81,81,81,77)
Ndose<-length(doselev)

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

dtarget<-100
diftarget<-2.464592
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<-FixedMean(n,doselev,meanlev,sdy)  

prior<-emaxPrior.control(epmu=0,epsca=30,difTargetmu=0,
		difTargetsca=30,dTarget=100,p50=50,sigmalow=0.1,
		sigmaup=30,parmDF=5)
mcmc<-mcmc.control(chains=1,warmup=500,iter=5000,seed=53453,propInit=0.15,adapt_delta = 0.95)

D1 <- emaxsimB(nsim,gen, prior, modType=3,mcmc=mcmc,check=FALSE)

plot(D1,id=3)

## End(Not run)