使用Julia语言的线性混合效果模型中随机效果组的标签

时间:2018-08-31 16:29:49

标签: julia mixed-models

我是Julia编程语言的新手,但是,我正在拟合线性混合效果模型,并且发现使用ranef()函数很难获得随机效果“组名”的标签。

可以找到示例代码:

$using MixedModels
$@time modelOutput = fit(LinearMixedModel(Y~ A + B + (0 + A | group), data))

我可以使用(ranef(modelOutput))效果看到每个组的随机效果估计,但是无法获得组中每个级别的标签。我很难将随机效应估计值映射到特定的组级别,因为我大约有60个级别

任何建议都值得赞赏。

0 个答案:

没有答案