问题标签 [piecewise]

我按照以下步骤绘制了分段线性回归的结果，其中包含我通过segmented包完成的一个断点：

我得到了这个图，它没有显示带有断点的两条线！！！我应该如何更改我的代码以获得正确的情节？

这是我的数据集：（ChH 有 11 个缺失数据）

和：

我有一个实验数据，它是分段连续的，每个部分都应该线性拟合。但是，我想在不知道结到底在哪里（所以斜率变化的点）的情况下拟合它，因为手动确定它们并不容易。

到目前为止，我找到了使用 py-earth完成此任务的建议，但不明白如何实现它；也就是说，只有 X 和 Y 的变量列表，如何能够进行这种分段线性回归。有人可以给我建议如何做到这一点吗？

UPD：原来，问题是由于数组格式不同。我的数组的“X=numpy.array([X]).T”解决了它，现在 py_earth 正在工作。然而，它是“粗糙”的，显示一个线几个结。任何人都可以为分段线性回归提出其他解决方案吗？

目标是实现一个分段表达式，当 n 为偶数时给出 0，当 n 为奇数时给出 1。一种方法是使用下面的 floor 函数：

但是，这会返回错误的输出：

正确的输出应该是：

另一种实现相同目的的方法是使用：

这将返回正确的输出。有没有办法使用原始代码中的 floor 函数来实现这一点？

我想要一系列分段函数，每个函数都是基本函数的移位副本，但我不知道如何在 sympy 中做到这一点。这是一个例子：

但是，最后一行给了我一个错误：

如何使分段可调用并指定fs是x和s的函数？

假设是具有 n 维v的向量，我想进行以下分段操作：aEigen::VectorXd

• vby的分段乘法a，即向量(a[1]*v[1], ..., a[n]*v[n])和
• 的分段平方v，即向量(v[1]*v[1], ..., v[n]*v[n])。

是否Eigen提供了上述操作的方法，还是我需要手动实现它们？当然也有非常简单的，但我希望它们尽可能快地运行。

我试图编写一个代码，它是一个更大程序的一部分，它将z在s. 但是，当我运行代码时，我只得到z=0，或者如果 lastelse被忽略，则代码返回零向量。

有人知道我在哪里犯错了吗？我已使用此来源的方法 1 。任何帮助将不胜感激，我几个月来一直在努力完成这项工作。

为了进一步参考，这解决了我的问题。谢谢@丹！

我想做一个具有固定效应（有和没有交互）的分段多元回归。SAS 是否有任何程序来处理此任务。我读过 PROC NLIN 可以进行分段回归，但是我只能找到关于单个预测变量的分段回归的教程。如果您能指出有关此问题的任何示例，我将不胜感激。教程的链接也很好。

我目前正在使用 SAS EG 6.1。

非常感谢，

麦

嗨，我想弄清楚如何用分段线性函数拟合这些值。我已阅读此问题，但无法继续前进（如何在 Python 中应用分段线性拟合？）。在此示例中，展示了如何为 2 段案例实现分段函数。但我需要在如图所示的三段式案例中进行操作。

我写了这段代码：

但这是输出

有什么建议吗？我相信问题return np.piecewise(x , [x <= x0, (x>= x1)] , [lambda x:k0*x + y0-k0*x0, lambda x:k1*(x-(x1+x0))-y1, lambda x:k2*x + y1-k2*x1])尤其在第二个 lambda.

编辑1：

如果我尝试使用 AL 提供的解决方案不同的数据，我不会得到好的结果。

我得到这个结果：

和

和

运行以下代码：

我得到了情节：

但是当我尝试绘制 ton 的值时，wchich 不同于零，例如：

该图似乎将所有 xon 值四舍五入为整数：

我的代码中的什么导致了这种行为？

此代码运行正确：

但是当函数变得分段时，例如：

我收到以下错误：

文件“//anaconda/lib/python2.7/site-packages/sympy/core/relational.py”，行> 103，非零 引发TypeError（“无法确定\ n％s的真值”％自我）

TypeError：无法确定 ton <= t 的真值

据我了解，错误是由于两者都ton可以t是积极的和消极的。这是正确的吗？如果我为错误设置正积分限制t不会消失。如何计算给定分段函数的积分？

更新：该功能的更新版本，有效：