powered by
Internal function, logdiff
logdiff(x, y)
numeric
max(x) + log(exp(x - max(x,y)) - exp(y-max(x,y)))
This function calculates the log of the difference of the exponentiated logs taking out the max, i.e. insuring that the difference is not negative