added estimated probability

modified to use estimated probability using estimated probability = 1/[1 + exp(-a - BX)].