I am trying to fit a fairly complicated mixed model using Julia and the MixedModels package, on a dataset that has ~3.5 million rows (see screen shot for model specification and errors).

The model is nested, with the random effect grouping factors V*, S* and M* representing 3 different crossed effects, where * represents the level of the hierarchy.

I simplified the model to troubleshoot, keeping only the V1, S1, and M1 terms, but got the same error. Interestingly, the model will fit without errors if I keep V1 and S1, V1 and M1, or S1 and M1 - just not all three. The model will run with V1, S1 and M1 using lme4 in R though. Being new to Julia, I don’t know how else to go about troubleshooting this error message, so I would appreciate any help anybody can offer.

I am running Julia v0.5.0 on Windows 7, using MixedModels v0.7.5.