Skip to content

[WIP] Add support for uncertainty in ratios#61

Open
glouppe wants to merge 23 commits into
diana-hep:masterfrom
glouppe:histogram-uncertainty
Open

[WIP] Add support for uncertainty in ratios#61
glouppe wants to merge 23 commits into
diana-hep:masterfrom
glouppe:histogram-uncertainty

Conversation

@glouppe

@glouppe glouppe commented Apr 25, 2017

Copy link
Copy Markdown
Contributor

No description provided.

Comment thread carl/ratios/classifier.py
(std_p / p[:, 1]) ** 2 -
2 * (-std_p ** 2) / (p[:, 0] *
p[:, 1]))) ** 0.5
# nb: cov(p, 1-p) = -var(p) = -std^2(p)

Copy link
Copy Markdown
Contributor Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Instead, derive std_r as dr/dp * std_p

Copy link
Copy Markdown
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I just tried it both ways and got the same answer. Maybe good to double check.

Copy link
Copy Markdown
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

not sure if there are numerical issues.
pretty sure this simplifies to
std_r = std_p / ( p[:,1]**2)

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment

Labels

None yet

Projects

None yet

Development

Successfully merging this pull request may close these issues.

2 participants