Displaying 1 result from an estimated 1 matches for "const2a".
Did you mean:
const2
2005 Jul 13
1
Fieller's Conf Limits and EC50's
...pha<-1-conf.level
zalpha.2 <- -qnorm(alpha/2)
gamma <- zalpha.2^2 * var.b1 / (b1^2)
eta = family(obj)$linkfun(p) #based on calcs in V&R's dose.p
EC50 <- (eta-b0)/b1
const1 <- (gamma/(1-gamma))*(EC50 + cov.b0.b1/var.b1)
const2a <- var.b0 + 2*cov.b0.b1*EC50 + var.b1*EC50^2 -
gamma*(var.b0 - cov.b0.b1^2/var.b1)
const2 <- zalpha.2/( (1-gamma)*abs(b1) )*sqrt(const2a)
LCL <- EC50 + const1 - const2
UCL <- EC50 + const1 + const2
conf.pts <- c(LCL,EC50,UCL)...