3
在Julia中可能存在LMM中随机效应的相互作用吗?Julia中混合效应模型中随机斜率的相互作用
这给出了一个错误
model = fit!(lmm(@formula(response ~ 1 + A*B + (1+A*B|sub)), data)
ERROR: MethodError: no method matching getindex(::DataFrames.DataFrame, ::Expr)
拆开方面也没有帮助。
model = fit!(lmm(@formula(response ~ 1 + A*B + (1+A+B+A&B|sub)), data)
这工作
mode2 = fit!(lmm(@formula(response ~ 1 + A*B + (1+A+B|sub)), data)
注意,没有任何问题,当你有固定效应的相互作用。
谢谢!这似乎工作。我正在使用Julia v0.6.0,并且能够获得模型的完整输出。它不想用'fast = true'运行,但给出了和R中lmer大致相同的结果(一般来说,我看到一些变化,所以不能确定这是由于重新定义了计算还是正常行为) 。不确定Cholesky分解,但与R相比,随机效应的相关性是关闭的。 – mike