Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Jump to bottom

dwelltime: add standard errors based on hessian approximation #690

Open
wants to merge 2 commits into
base: main
Choose a base branch
from
Open

dwelltime: add standard errors based on hessian approximation #690

wants to merge 2 commits into from

Commits on Sep 5, 2024

  1. dwelltime: add asymptotic standard errors 

    When calculating the asymptotic uncertainty interval, we should take into account that we actually impose a linear sum constraint (otherwise, the amplitudes will have indeterminate confidence intervals).

One could add the constraint explicitly to the Hessian by simply adding a large penalty term to the relevant derivatives. Considering that the sum constraint is of the form (1 - sum(a_i)) ** 2, this would result in adding a constant term d^2f/daidaj = -c with c large to all amplitude terms.

What is ugly is that we would need to choose this constant as large as possible without incurring numerical issues. This is why it is preferable to project onto the null space and then calculate the result back instead.
    @JoepVanlier
    JoepVanlier committed Sep 5, 2024
    Configuration menu
    Copy the full SHA
    516e86e View commit details
    Browse the repository at this point in the history

  2. profiles: provide the option to plot w/stderr

    @JoepVanlier
    JoepVanlier committed Sep 5, 2024
    Configuration menu
    Copy the full SHA
    fcc014a View commit details
    Browse the repository at this point in the history