Plot the output of emaxsimB
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.
## S3 method for class 'emaxsimB' plot(x, id = x$idmax, plotDif= TRUE, ...)
x |
Output of |
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 |
ggplot object is returned
Neal Thomas
## 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)
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.