doselev<-c(0,5,25,50,100,350)
n<-c(78,81,81,81,77,80)
### population parameters for simulation
e0<-2.465375
ed50<-67.481113
dtarget<-100
diftarget<-9.032497
emax<-solveEmax(diftarget,dtarget,log(ed50),1,e0)
sdy<-8.0
pop<-c(log(ed50),emax,e0)
dose<-rep(doselev,n)
meanlev<-emaxfun(dose,pop)
y<-rnorm(sum(n),meanlev,sdy)
testout<-fitEmax(y,dose,modType=4)
sigma(testout)
Run the code above in your browser using DataLab