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.

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

Force non-zero step size for general optimization after boundary fit in fast=true #201

Closed
wants to merge 3 commits into from

Conversation

palday
Copy link
Member

@palday palday commented Sep 27, 2019

Resolves #94, Closes #95.

For boundary fits in GLMM, the initial step for general optimization in fast=false includes a zero in (at least) one dimension. This forces that to be nonzero by replacing zeros with the smallest nonzero value.

Note: this will not address models where everything is zero (including fixed effects), but that's another can of worms.

@codecov-io
Copy link

codecov-io commented Sep 27, 2019

Codecov Report

Merging #201 into master will increase coverage by 1.97%.
The diff coverage is 100%.

Impacted file tree graph

@@            Coverage Diff             @@
##           master     #201      +/-   ##
==========================================
+ Coverage   92.96%   94.93%   +1.97%     
==========================================
  Files          18       18              
  Lines        1294     1166     -128     
==========================================
- Hits         1203     1107      -96     
+ Misses         91       59      -32
Impacted Files Coverage Δ
src/MixedModels.jl 100% <ø> (ø) ⬆️
src/generalizedlinearmixedmodel.jl 93.78% <100%> (+2.59%) ⬆️
src/utilities.jl 73.52% <0%> (-10.48%) ⬇️
src/linearmixedmodel.jl 92.15% <0%> (-4.2%) ⬇️
src/remat.jl 95.32% <0%> (-0.74%) ⬇️
src/varcorr.jl 100% <0%> (ø) ⬆️
src/linalg.jl 97.89% <0%> (+2.01%) ⬆️
src/linalg/rankUpdate.jl 97.82% <0%> (+7.82%) ⬆️
src/randomeffectsterm.jl 96.15% <0%> (+23.23%) ⬆️
... and 1 more

Continue to review full report at Codecov.

Legend - Click here to learn more
Δ = absolute <relative> (impact), ø = not affected, ? = missing data
Powered by Codecov. Last update e003fd8...4b277e5. Read the comment docs.

@palday
Copy link
Member Author

palday commented Sep 27, 2019

Stung by floating point comparison. I'll fix that.

@palday palday force-pushed the glmm_boundary_genopt branch 2 times, most recently from 923e233 to e5aae50 Compare October 1, 2019 08:03
@palday
Copy link
Member Author

palday commented Oct 1, 2019

Optimization gives different results on different platforms, but neither is the 'correct' answer, see #206.

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.

Gamma models convergence checks and run time Watch for boundary cases from fast=true in glmm fits
2 participants