我想将多项式系数附加到 data.frame，如下所示：

我想知道是否有这样的简洁方法

谢谢

我正在做一个带有 spml 函数（splm 包）的空间 durbin 模型，我已经将一个多项式（2 级）作为回归量包括在内，显然我也需要包括它的空间滞后。不幸的是，函数 slag 和 lag.listw 不适用于 poly 函数，这给了我错误：

没有适用于“炉渣”的方法应用于“c（'poly'，'matrix'）”类的对象

和

lag.listw(dist.listw, poly(GDP.PCAP, 2)) 中的错误：对象长度不同

手动插入两个变量被排除在外，因为该模型不适用于奇异性问题。

如果有人可以帮助我，我将不胜感激。

谢谢西尔维娅

我想将包含多项式的列添加到数据框（DF）。

背景：我需要在 glmnet 设置中使用多项式。我不能在 glmnet() 估计命令中直接调用 poly()。我收到一个错误，可能是因为我的“Xtrain”数据包含因子。我的解决方法是将我的 Xtrain DF 分成两部分，一部分包含所有因素（不需要转换），另一部分包含其余部分，即。数字列。

现在我想将带有多项式的列添加到我的数字 DF 中。这是我的问题的一个最小示例。

但是，我无法弄清楚如何“cbind”结果。我最终想要的是一个包含 x, x^2, y, y^2 的 DF。顺序无所谓。但是，理想情况下，我也会有列标签（以识别多边形）。例如像这样：

谢谢...干杯！

我需要绘制 x 和 y 之间的关系，其中 x 的多项式预测 y。这是使用 poly() 函数完成的，以确保多项式是正交的。

考虑到线性、二次和三次项，我如何绘制这种关系？问题是不同项的系数没有像 x 那样缩放。

我在下面提供了一些示例代码。我尝试将每个多项式的对比度值重新分配给 x。

该解决方案给出了不可能的预测值。

预先感谢您的帮助 ！

最良好的祝愿，埃里克

这是一个示例代码：

我有定义桥梁开口形状的数据。我开发的代码用水平线切割桥梁开口，就好像它是水位一样，包括桥墩和墙壁上的交点。我想确保形状始终是边界形状而不是对角线切割。当我使用排序并绘制数据时，我得到一个 Zig Zaw 形状......所以我从另一个进程获得的坐标序列是 A，我需要的序列是 B......？如何确保左侧和右侧的垂直点是避免 Zig Zaw 的序列？

我在 R 中使用正交回归来找到两个变量之间的关系。我想在最佳拟合线上的几个点中找到切线。

给了我一个非常意想不到的结果。我怀疑这是因为 poly 函数使用正交多项式。在上面的回归中我得到

返回我

这显然是错误的，“就像系数一样”（知道它们没有错）。

现在我想找到 x=2 和 x=7 的切线

如果我只看数字，我怀疑在 x=2 中，切线会是 6.5 之类的东西？(23-10) / (3-1)

因为它是二阶多元回归，所以输入我从回归中得到的汇总变量是没有意义的，因为它给出了毫无意义的结果

我试图更好地理解为什么 stat smooth 不会绘制我的多项式回归线，除非我的 x 变量（自变量）首先被分配为绘图之外的值（例如 x <- dataset$Salary）

数据集

返回错误的初始图

错误：'newdata' 有 80 行，但找到的变量有 10 行

有效的解决方案

据我了解

x <- dataset$Salary 与 dataset$Salary 没有什么不同，只是包含在一个值中。我唯一的想法是它与 poly() 如何查看 x（一个数字向量）与它如何将 dataset$Salary 视为提取的向量有关。一世

除此之外，我希望得到相同的结果，但事实并非如此。

我还尝试将 x 重命名为 t，它与第一张图所做的完全一样，所以我不明白如果 x 只是值的名称，为什么它如此重要。

我得到了以下回归模型，我正在尝试做一些预测，但我总是得到一个错误：

poly(Sepal.Width, 2, coefs = list(alpha = c(3.05733333333333, : object 'Sepal.Width' not found) 中的错误

我知道我可以简单地访问预测的拟合值。但是假设我的新数据如下所示：

如何在不出现此错误的情况下对这个新数据集进行预测？

我正在尝试使用step具有 1,400 个变量的 AIC（通过 ）运行逐步回归，但我的计算机只是冻结了。如果我包含 <300 个变量（运行 13 小时后），它会起作用。

在我运行逐步回归之前，有没有办法消除一些变量（如果 p 值 >.7）？

我的目标是在主题内跨时间使用多个 IV 运行二次函数。我遇到了一些代码，有点困惑。下面是我正在尝试运行的可重现示例。遵循代码将是我的问题。

我的问题如下：

1. 在非 poly() 中，lmer 语句time*IV2将给出时间和 IV2 相互作用的系数以及低阶系数时间和 IV2。我是否正确使用 poly() 语句不会将较低项放入模型中？

2. 我被教导说，如果你包括较高的术语，你也应该包括较低的术语。这对于 r 中的 poly() 函数仍然正确吗？

如果是这样使用任何一个都有意义

我会假设以上两个是等价的，但是，我错了，因为第二个给了我一个错误：

错误：删除列无法生成完整的列秩设计矩阵

3. 我删除了哪些列？

感谢您的帮助。我对 r 很陌生，所以我正在尽力了解这些功能的作用，而不是仅仅遵循一个食谱。