Hi,

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:

typeof(dataframe[:X_categorical])

CategoricalArrays.CategoricalArray{Int64,1,UInt32,Int64,CategoricalArrays.CategoricalValue{Int64,UInt32},Union{}}

If anybody might have any insight I would be so appreciative.

-Sarah