Question

我是Julia编程语言的新手，但是，我正在使用Linear Mixed Effects Model，我发现很难将fixed和random效果估算保存在{{ 1}}文件。

可以找到示例代码：

.csv

有关如何获取固定（using MixedModels @time modelOutput = fit(lmm(Y~ A + B + (0 + A | group), data))）和随机（fixef(modelOutput)）效果的参考资料，但是使用我遇到错误的DataFrame。

感谢任何建议。

Answer 1

好的，我实际上是花时间为你做这件事。 CoefTable是statmodels here中定义的类型。根据这些信息，我们可以从CoefTable实例中提取相关信息，如下所示：

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

这将给出nvar-by-4 DataFrame，然后您可以使用writetable("output.csv",df)

按前面所述写入csv

Answer 2

我遇到许多问题，无法获得公认的工作答案；自那时以来，朱莉娅已经发展了很多。我主要根据jglmm R package的代码重写了它，并从其他来源进行了一些改编/整合……

"""
    outfun(m, outfn="output.csv")

output the coefficient table of a fitted model to a file
"""
outfun = function(m, outfn="output.csv")
    ct = coeftable(m)
    coef_df = DataFrame(ct.cols);
    rename!(coef_df, ct.colnms, makeunique = true)
    coef_df[!, :term] = ct.rownms;
    CSV.write(outfn, coef_df);
end

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

2 个答案: