Closed ararslan closed 6 months ago
Looks like https://github.com/JuliaStats/MixedModels.jl/blob/cfd302310ba66ad5008e97738bd6710cf5a7509a/src/bootstrap.jl#L250 is problematic for rank deficient models since β_names has one name per fixed effect (which is ≤ the number of coefficients) but fixef!(βsc, m) is being called with the full coefficient vector as the first argument.
β_names
fixef!(βsc, m)
Looks like https://github.com/JuliaStats/MixedModels.jl/blob/cfd302310ba66ad5008e97738bd6710cf5a7509a/src/bootstrap.jl#L250 is problematic for rank deficient models since
β_names
has one name per fixed effect (which is ≤ the number of coefficients) butfixef!(βsc, m)
is being called with the full coefficient vector as the first argument.