保存 FittedModel 有问题

Question

在 Mathematica 7 和 8 中，我们发现 Mathematica 无法正确保存并获取拟合例程的结果。我希望这里有人知道：此错误的最佳实践解决方法是什么？

我们的任务是使用 NonlinearModelFit 拟合大量数据集，它为每个数据集返回一个 FittedModel。这些拟合需要一段时间，我们尝试使用 Save 和 DumpSave 将结果存储到文件中。后续分析使用 Get 将 FittedModel 拉回。Mathematica

执行 Save & 操作。没有警告，但返回的 FittedModel 有问题。我们调用 Get 并从磁盘加载 FittedModel，将结果命名为“foo”。尝试 foo[10] 使用输入 10 评估拟合模型。我们看到的错误是 foo["BestFitParameters"] 应该返回一个规则列表，而是将字符串插入到函数中，因为它是像 10 这样的数字。 [{"BestFitParameters","BestFitParameters"}] 应该在列表中给出两次规则，但对于有问题的 foo 保持未评估状态。

除此之外，还有一个奇怪且误导性的行为，具体取决于我们是否删除/清除 foo 或退出内核：

• 将 FittedModel 计算为 foo 将 FittedModel
• 保存或转储保存 foo 到磁盘
• 删除和/或清除 foo
• 从磁盘获取 foo
• foo 工作正常

但是如果我们重新启动内核，它将失败

• 将 FittedModel 计算为 foo
• 保存或转储 将 foo 保存到磁盘
• 终止并重新启动内核
• 从磁盘获取 foo
• foo 有问题

有人以前见过这种行为吗？

这个错误背后有一个很好的解释吗？

有没有好的解决方法？

Answer 1

这是一个已知的错误，幸运的是有一个简单的解决方法。在评估 Get 以读取保存的 FittedModel 之前，评估虚拟拟合模型代码。

In[1]:= NonlinearModelFit[Range[5], a x, x, a];

In[2]:= Get[FileNameJoin[{$HomeDirectory, "Desktop", "bignlm.mx"}]];

In[3]:= AbsoluteTiming[nlm["BestFitParameters"]]

Out[3]= {17.6010000, {a -> 1.45015, b -> 2.33999}}

如果跳过此初始评估，则由于无法自动加载实现该功能的所有依赖代码，拟合模型似乎无法正常工作。 In[1] 的评估基本上强制加载。