-
Notifications
You must be signed in to change notification settings - Fork 2
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.
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
likelihood-ratio test #90
Conversation
Codecov ReportAll modified and coverable lines are covered by tests ✅
❗ Your organization needs to install the Codecov GitHub app to enable full functionality. Additional details and impacted files@@ Coverage Diff @@
## main #90 +/- ##
==========================================
+ Coverage 94.16% 94.23% +0.07%
==========================================
Files 10 10
Lines 2090 2117 +27
==========================================
+ Hits 1968 1995 +27
Misses 122 122 ☔ View full report in Codecov by Sentry. |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Nice one, this is exactly how I thought of that implementation
# Since N/p_c is an integer, alpha*p_c %% N == 0 iff alpha %% (N/p_c) == 0 | ||
# In other words, alpha must be a multiple of (N/p_1, N/p_2,...,N/p_C) | ||
# (alpha = k*(N/p_1) or k*(N/p_2) or ... or k*(N/p_C) for some integer k (including 0)) | ||
# However, alpha must be at most M, and a valid bound from above for integer k is max(p_1,...,p_C). |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Here M is not defined. I know it is defined in the paper but not really in the comments
possible_alphas <- unique(sort(alpha_matrix[alpha_matrix <= M])) | ||
|
||
# Recalculate the r_alpha vector using its definition directly. |
There was a problem hiding this comment.
Choose a reason for hiding this comment
The reason will be displayed to describe this comment to others. Learn more.
Just a nitpick but there is not recalculation
of the r_alfa vector, this is the first time the vector is calculated
No description provided.