Skip to content

Conversation

@KrisThielemans
Copy link
Collaborator

WIP to fix #1472. Currently only handling Hessian.

accumulate_sub_Hessian_times_input_without_penalty used divide_and_truncate,
but this is incorrect as its threshold strategies are not appropriate
for (measured * forward_input)/forward_current^2. Most obvious example
is if forward_input contains negatives. We now use 0/x=0 (first thresholding
measured data to be non-negative).

UCL#1461
@KrisThielemans KrisThielemans self-assigned this Jul 25, 2024
@KrisThielemans KrisThielemans changed the title Poisson ll truncation change Poisson log-likelihood truncation/thresholding strategy Jul 25, 2024
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Projects

None yet

Development

Successfully merging this pull request may close these issues.

inconsistency between Poisson loglikelihood value, gradient (when using multiplicative factors) and Hessian due to thresholding

1 participant