我想在Julia中使用来自R中lme4库的lmer函数获得相同的结果.请在下面找到使用R的内置mtcars数据集的示例
library(lme4)
data<-mtcars
data$vs<-as.factor(data$vs)
data$am<-as.factor(data$am)
data$gear<-as.factor(data$gear)
str(data)
model <- lmer(mpg ~ cyl:gear + hp:am + (1|gear:am), data = data)
我已经找到了lmm()来自Julia的MixedModels包的函数,它应该能够运行相同的结果,但是我不知道如何从lmer()函数的第一个参数重写公式lmm().特别是interaction(:)运算符.
我将以简短的例子回答你的问题.