我已经在Julia 1.0.4中使用MixedModels和StatsModels包创建了线性混合模型,因为R中的lme4似乎没有在我的数据集中运行。但是,现在有了模型,我还无法找到将其应用于测试数据的方法,因此无法评估模型的准确性。是否有朱莉娅(Julia)替代允许新数据的R预测()方法?我已经尝试过使用Julia的predict()方法,但是它只能预测与我创建模型时使用的数据相同的数据。
我就像下面的示例一样创建了模型:
mm1 = fit(LinearMixedModel, @formula(Y ~ 1 + (1 | G)), dyestuff)
但是似乎没有一种方法可以允许
predict(mm1, test_data)
在茱莉亚工作。我也尝试过使用RCall将mm1变量发送到R并以这种方式使用predict()方法,但是@rput似乎无法使用这种变量。谢谢!