我是Julia编程语言的新手,但是,我正在拟合线性混合效果模型,并且发现使用ranef()
函数很难获得随机效果“组名”的标签。
可以找到示例代码:
$using MixedModels
$@time modelOutput = fit(LinearMixedModel(Y~ A + B + (0 + A | group), data))
我可以使用(ranef(modelOutput))
效果看到每个组的随机效果估计,但是无法获得组中每个级别的标签。我很难将随机效应估计值映射到特定的组级别,因为我大约有60个级别
任何建议都值得赞赏。