if (FALSE) {
## emaxsim changes the random number seed
nsim<-50
idmax<-5
doselev<-c(0,5,25,50,100)
n<-c(78,81,81,81,77)
### population parameters for simulation
e0<-2.465375
ed50<-67.481113
dtarget<-100
diftarget<-9.032497
emax<-solveEmax(diftarget,dtarget,log(ed50),1,e0)
sdy<-7.967897
pop.parm<-c(log(ed50),emax,e0)
meanlev<-emaxfun(doselev,pop.parm)
###FixedMean is specialized constructor function for emaxsim
gen.parm<-FixedMean(n,doselev,meanlev,sdy)
D1 <- emaxsim(nsim,gen.parm)
print(D1,c(31,50),digits=2,id=4)
print(D1,c(1,20))
D1 ### implicitly calls print with default parameter settings
}
# \dontshow{
## emaxsim changes the random number seed
nsim<-3
doselev<-c(0,5,25,50,100)
n<-c(78,81,81,81,77)
### population parameters for simulation
e0<-2.465375
ed50<-67.481113
dtarget<-100
diftarget<-9.032497
emax<-solveEmax(diftarget,dtarget,log(ed50),1,e0)
sdy<-7.967897
pop.parm<-c(log(ed50),emax,e0)
meanlev<-emaxfun(doselev,pop.parm)
###FixedMean is specialized constructor function for emaxsim
gen.parm<-FixedMean(n,doselev,meanlev,sdy)
D1 <- emaxsim(nsim,gen.parm,nproc=1)
print(D1,c(2,3),digits=2,id=4)
# }
Run the code above in your browser using DataLab