Mean response and SE for each replicate data set in an emaxsimB object
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
## 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)
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.