使用Julia语言在csv中导出线性混合效果模型输出

面具

我是Julia编程语言的新手,但是,我很适合,Linear Mixed Effects Model并且发现很难将fixedrandom效果估计值保存.csv文件中。

可以找到示例代码:

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

有关于如何获得固定(fixef(modelOutput))和随机(ranef(modelOutput))效果的可用参考,但是使用DataFrame却遇到了错误。

任何建议表示赞赏。

尼尔斯·古达(Nils Gudat)

好的,我实际上是花时间为您做的。ACoefTable是在statmodels 此处定义的类型有了这些信息,我们可以从CoefTable实例中提取相关信息,如下所示:

df = DataFrame(variable = ct.rownms,
               Estimate = ct.mat[:,1],
               StdError = ct.mat[:,2],
               z_val = ct.mat[:,3])

这将提供一个nvar-by-4 DataFrame,您可以按照前面所述使用以下方法将其写入csv:writetable("output.csv",df)

本文收集自互联网,转载请注明来源。

如有侵权,请联系 [email protected] 删除。

编辑于
0

我来说两句

0 条评论
登录 后参与评论

相关文章