問題タブ [piecewise]

次の手順に従って、segmentedパッケージで行った 1 つのブレークポイントを使用して区分線形回帰の結果をプロットしました。

そして、ブレークポイントのある2行を表示しないこのプロットを取得します!!! 正しいプロットを取得するには、コードをどのように変更すればよいですか?

これは私のデータセットです: (ChH には 11 個の欠損データがあります)

と：

区分的に連続した実験データがあり、各部分は線形に適合する必要があります。ただし、手動でノットを特定するのは簡単ではないため、ノットが正確にどこにあるのか (つまり、勾配が変化しているポイント) を知らずにフィットさせたいと考えています。

これまでのところ、このタスクにpy-earth を使用するというアドバイスを見つけましたが、それを実装する方法を理解できませんでした。つまり、X と Y の変数のリストだけがあれば、そのような区分的線形回帰をどのように実行できるかということです。誰かがそれを行う方法についてアドバイスをもらえますか?

UPD :結局のところ、問題は配列の形式が異なるためでした。私の配列の "X=numpy.array([X]).T" で解決され、py_earth が機能するようになりました。ただし、これは「大雑把」で、数ノットに対して 1 本の線を示しています。区分線形回帰の他の解決策を提案できる人はいますか?

それぞれが基本関数のシフトされたコピーである一連の区分関数が必要ですが、シンピーでそれを行う方法がわかりません。次に例を示します。

ただし、最後の行でエラーが発生します。

Piecewise を呼び出し可能にして、 fsがxとs の関数であることを指定するにはどうすればよいですか?

こんにちは、これらの値を区分線形関数に適合させる方法を見つけようとしています。この質問を読みましたが、先に進むことができません ( Python で区分線形フィットを適用する方法は? )。この例では、2 セグメントの場合に区分関数を実装する方法を示します。ただし、図のように 3 セグメントの場合に行う必要があります。

私はこのコードを書きました：

しかし、これは出力です

なにか提案を？問題は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])特に 2 番目に あると思いlambdaます。

編集1：

AL が提供するソリューションとは異なるデータを使用しようとすると、良い結果が得られません。

私はこの結果を得ます：

と

と

次のコードを実行します。

私はプロットを取得します：

しかし、トンの値をプロットしようとすると、ゼロとは異なります。たとえば、次のようになります。

プロットは、すべての xon 値を整数に丸めているようです:

私のコードの何がそのような動作を引き起こしますか?

このコードは正しく実行されます:

しかし、関数が区分的になると、例えば:

次のエラーが表示されます。

ファイル "//anaconda/lib/python2.7/site-packages/sympy/core/relational.py"、行 > 103、ゼロ以外で TypeError("\n%s の真値を判別できません" % self)

TypeError: ton <= t の真理値を特定できません

私が理解している限り、エラーは と の両方が正と負になる可能性があるという事実によるものtonですt。それが正しいか？正の積分限界を設定してtも、エラーは消えません。指定された区分関数の積分を計算するにはどうすればよいですか?

更新:機能する関数の更新バージョン: