predict.emaxsimB | R Documentation |
Return warning and explanation that only predicted values at doses included in the study are available. The code needed to obtain predicted values at other doses is indicated.
## S3 method for class 'emaxsimB'
predict(object,
dose, dref=0, ...)
object |
Output of |
dose |
Vector (can be a single value) of doses where dose response curve is to be evaluated. |
dref |
A reference dose (0 by default) for contrasts, but other values can be specified. If specified, a single reference value must be given. |
... |
Optional arguments are not used. |
No output.
Neal Thomas
emaxsimB
, summary.emaxsimB
,
plot.emaxsimB
## Not run:
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,seed=12357,mcmc=mcmc,check=FALSE)
predict(D1,dose=20)
## End(Not run)
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.