Computing log liklihood given a probability function

[jbrea]

I think the easiest would be to implement your own log-likelihood function. Something like

function loglikelihood(n, m, T)
    prob = NTCP.(X_Data, n, m, T)
    sum(y == 1 ? log(p) : log(1 - p) for (y, p) in zip(Y_data, prob))
end