I’m trying to find a way to run a Linear Mixed Model with certain fixed effects using the :
μ = @formula(T ~ X_cts + X_categorical)
mixed_effect_model = MixedEffectModel(μ, Σ, IdentityLink(), NormalResponse(0.000001))
But for whatever reason Julia doesn’t seem to like my X_categorical variable…
When I do the @formula without this categorical variable and just the continuous variable, μ = @formula(T ~ X_cts), I have no issues.
I have the following:
If anybody might have any insight I would be so appreciative.