dose.resp.ss<-function(L,alpha,beta,q=0.05){ ea<-exp(alpha) w0<-function(a,b){return(ea/(1+ea)^2)} w1<-function(a,b){return(b*(ea-ea^2)/((1+ea)^3))} w2<-function(a,b){return((b^2)*ea*(1-4*ea+(ea^2))/((1+ea)^4))} N<-((qnorm(1-q/2)^2)*(w0(alpha,beta)+1.5*w2(alpha,beta))) / (L^2*(w0(alpha,beta)+1.5*w2(alpha,beta))* (w0(alpha,beta)+0.5*w2(alpha,beta))-w1(alpha,beta)^2) return(N) }