e0<-10
dose<-1
led50<-log(0.5)
lambda<-2
target<- -1.5
emax<-solveEmax(target,dose,led50,lambda,e0)
emax
dose1<-solveDose(target,led50,lambda,emax,e0)
dose1
emaxfun(dose=dose1,parm=c(led50,lambda,emax,e0)) - e0
Run the code above in your browser using DataLab