def getNegLL(self, t, mu, sqrtVar, controls, cases):
z = (mu-t) / sqrtVar
logProbControls = stats.norm(0,1).logcdf(-z)
logProbCases = stats.norm(0,1).logcdf(z)
ll = logProbControls[controls].sum() + logProbCases[cases].sum()
return -ll
评论列表
文章目录