小议非二元响应变量的EC50

其实解决这个问题的R包还是很多的，MASS里就有一个dose.p来计算LD50，我们来仔细看看其中的代码：

dose.p <- function(obj, cf = 1:2, p = …