《药物评价实验设计与统计学基础》P31 所介绍的算法过于简化,仅仅是大概的理解,具体算法是不准确的。而且不能处理概率是0或100%的情况。
现代probit回归,不是这么简单计算的,只能通过非线性最小二乘法或最大似然法逐步迭代逼近,具体见Finney的《Probit Analysis》。
所以实际应用时,可以采用楼上建议的SAS里面的 PROCEDURE Probit(看SAS手册,这个算法引用了Finney的《Probit Analysis》里面的从Bliss发展而来的算法),SPSS的命令行手册也有Probit过程,同样引用了Finney的经典文献。
R的drc包也是同样引用了同样的Finney的经典文献。