-
Notifications
You must be signed in to change notification settings - Fork 47
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
the end of zerocorr!
and checking consistent ReMat eltypes
#406
Conversation
Codecov Report
@@ Coverage Diff @@
## master #406 +/- ##
=========================================
Coverage ? 94.30%
=========================================
Files ? 23
Lines ? 1632
Branches ? 0
=========================================
Hits ? 1539
Misses ? 93
Partials ? 0
Continue to review full report at Codecov.
|
Co-authored-by: Alex Arslan <ararslan@comcast.net>
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.
This discussion makes me think that removing zerocorr!
altogether might be the best course.
I've been thinking that too. I wonder if we should remove it from the package but include it (or something like it) in the docs as an example of how advanced users can set restrictive variance structures. |
Sounds like a plan. Do you have time to do that or would you prefer that I do it? |
I haven't added the docs yet, but I've marked |
I probably won't get around to a vignette on |
~ the end of
zerocorr!` and checking consistent ReMat eltypes
~ the end of
zerocorr!` and checking consistent ReMat eltypeszerocorr!
and checking consistent ReMat eltypes
@dmbates Since I changed something in the |
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.
Looks good. Thanks for doing this.
@@ -78,6 +78,18 @@ function LinearMixedModel( | |||
for (i, x) in enumerate(Xs) | |||
if isa(x, AbstractReMat{T}) | |||
push!(reterms, x) | |||
elseif isa(x, ReMat) # this can occur in weird situation where x is a ReMat{U} | |||
# avoid keeping a second copy if unweighted | |||
z = convert(Matrix{T}, x.z) |
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.
I tend to write this construction as T.(x.z)
but this is quite reasonable too.
Closes #396.