问题标签 [poly]

您不想绘制每个特定的多边形，而是希望概括该过程，以便您只需编写一次 sub 。编写一个将写入和定位多边形的子程序。在这种情况下，您可以将星形设置为任何您想要的大小，但您应该根据中心定义坐标。子程序应该将多边形的中心作为参数，因此星的每个顶点都需要基于中心点。

我正在使用lmR 上的函数调整响应面。在输出中，使用fit[["fitted.values"]]，报告的值与实际值一致。但是，当我使用为每个效果生成的系数 ( Intercept, x1, x1*x2) 并进行手动计算时，呈现的值非常不一致。

该fitted.values函数如何进行计算？

下面是部分代码和输出：

的输出head(data)：

的部分输出fit[["fitted.values"]]：

总结系数：

我所做的是如下应用这些系数，只考虑重要的系数：

例如，使用 head(data) 输出的第一个值：

这导致 7897.79 与估计的fitted.values67.02不同

我究竟做错了什么？

我正在使用lm()包含多项式的训练数据集。当我提前使用子集时，与在函数调用中使用参数[ ]相比，我得到不同的系数。为什么？subsetlm()

^{由reprex 包于 2021-12-26 创建(v2.0.1)}

我已经拟合了一个具有方差结构的二次模型，该结构允许每个因子水平的不同方差水平，并且我无法预测仅包含 2 个条目的新数据集。这是一个可重现的示例：

然而，预测在一个lm框架上运行良好：

为什么会这样？的包维护者之一emmeans似乎认为这可能与缺少信息有关attr(, “predvars”) （请参阅我们在此处的讨论https://github.com/rvlenth/emmeans/issues/133）

我已将此情况报告给 Bates 博士（nlme联系人），但我想我也会接触到更广泛的社区。提前致谢

poly_model 和 lin_model 产生不同的系数。

通过检查该图，看来 lin_model 系数是正确的。为什么 poly_model 似乎会产生不同的系数？这些图重叠，两者的 R2 相同。

所以我必须用SML写一个小程序->>

一个名为 'p0.sml' 的文件，其中包含一个名为 的函数，该函数接受值的epoly参数 alist到，以及一个真实值。该列表包含形式为 的多项式的系数，其中使用的实数是传递给函数的参数。您的实现必须接受系数列表作为第一个参数，并将值作为第二个参数。您的函数必须返回由传递给它的参数指定的多项式的值。reala0anxa0 + a1x + a2x 2 + … + anx nxxx

这是我到目前为止所拥有的，但由于 as 的语法错误而无法编译。“错误：在 AS 发现语法错误”。如果您有任何指示将不胜感激。

我正在尝试使用该effects包在线性混合模型中创建效果图。我指定型号

对于这个模型，我可以按照我的预期生成分析结果，summary或者anova当我尝试查看特定效果时：

返回一个错误：

（对于plot函数，错误是不同的，但可能是由错误引起的allEffects）

我理解，根据对这个问题和这个问题的回答，“数值溢出”可能是poly术语的问题。但是，我不清楚这意味着什么或如何克服这个问题。

我还在这篇文章和另一篇关于 lme4 的文章中看到我再也找不到了，我可能需要更新软件包，所以我更新了“效果”和“lme4”以试图解决这个问题，但无济于事。

因此，如果由于“数字溢出”而发生此错误，我该如何解决这个问题？或者如果不是数字溢出发生了什么，我该如何解决这个问题？

我使用的数据的一个子集dput是：

谢谢你的帮助！

编辑！！我terms(fit1)按照评论中的建议运行，结果如下：