Matlab を使用して、poly1 曲線を使用して 2 つの座標 (x、y) にデータを当てはめています。

問題は、フィッティング ラインを長くする方法が見つからないことです。

(180, 930) から (191, 944) までが必要ですが、代わりに、Matlab はデータの近くにフィッティング ラインを描画します。これは、これら 2 つの座標の間にあります。

私を助けることができるfitコマンド（またはcftoolのいくつかの設定）への引数はありますか？

さらに、cftool で「軸の制限を調整する」オプションを試しましたが、まったく役に立ちませんでした。

既出の質問を検索しましたが、これに関連するものは見つかりませんでした。私はこのプログラムに慣れていないので、これがばかげた質問だったらごめんなさい、よろしくお願いします、ジョバンニ

編集:最初の画像のコードは次のとおりです。

2 番目のコードは次のとおりです。

ここで 2 つの適合を見つけることができます。最初の適合はトリミングされていないように見えますが、2 番目の適合は次のhold onとおり です。 docs.google.com/file/d/0B749BCu7mnZHeDVTOGRuSkktUmc/edit?usp=sharing

私が持っているデータ セットに MATLAB の曲線を当てはめようとしていて、最小二乗最小化を取得しようとしています。曲線はガウス関数と指数関数の畳み込みであるため、for ループを使用しました。曲線には 5 つの変数があります。

counts は、私が入力した唯一のデータです。

この曲線を最小化して変数 (a、b、sigma、t0、lambda) が可能な限り最適に適合する関数を見つけるのに苦労しています。可能な限りカウントします。どんな助けでも大歓迎です！

私はこれに空白を描いています。

私は2つのデータセットを持っています:

d2にどれだけうまく適合するかを示す、ある種の統計値、おそらくr値のようなものを取得したいと思いますd1。

ここでまた問題が発生します。私はいくつかのサンプル データを持っています。それらに j という名前を付けましょう。データのサンプルごとに、データに適合できる理論モデルがあります。簡単にするために、フィッティング関数に f(sig1, µ1, sig2, m2) という名前を付けます。もちろん、各サンプルのデータはわずかに異なるため、各サンプルから異なる sig1、µ1、sig2、µ2 を取得します。ここで問題が発生します。測定ごとに µ1 が同じである必要がありますが、事前にはわかりません。もちろん、すべてのサンプル j を適合させることができます。たとえば、μ1 の平均を取り、それらすべてを再度適合させますが、それは最良の近似値ではありません。したがって、すべての測定に対して同時に µ1 グローバルを適合させ、単一の測定ごとに他のすべてのパラメータを適合させたいと考えています。私の質問が今明確になったことを願っています。

これは、問題のもっともらしい方法について尋ねた相互検証済みの質問への続きです。この質問はよりプログラミング指向であるため、SO に投稿します。

バックグラウンド

既知の日付が 1 年以上にわたる曲線があります。この曲線の y 値は、毎日の気温と塩分の記録から計算されたd18O 値の予測です。また、炭酸カルシウムからなるシェルの d18O 値も測定しました。これらの値は、距離軸に沿って測定されます。最初と最後の測定は、曲線の開始と終了とほぼ同時に (正確ではありませんが) 行われます。

d18O 値は、未知の確率誤差内で曲線の予測値と一致することが知られています。測定値の x 軸を変更することによって (または少なくともインデックスを曲線のインデックスと一致させることによって)、測定値を曲線に最適に適合させたいと考えています。このようにして、測定値の日付を推定し、年間の殻の成長率をさらに推定することができます。成長率は変動すると予想され、成長の停滞（成長の停止）が生じる可能性があります。ただし、測定値間の増加は > 0 でなければなりません(制約)。

サンプルデータ

データセットの例 (curveおよびmeasured) は次のとおりです。

...そして、これがどのように見えるかです：

実際の問題

の x 軸を変更することにより、measRが適合する数学的/統計的方法を見つけたいと考えています。さらに、適合した「x軸」を相互に比較するために、ある種の適合度統計を取得したいと考えています（異なる制約で複数のモデルを実行する場合）。私は「x 軸モデル」を成長モデルと呼んでいます。値間の距離が > 0 でなければならないことを指定して、フィッティングを制限したい。つまり、 の値は の値の後に発生する必要がある。また、成長率 (つまり、隣接する 2 つのインデックス ポイント間の最大距離) を制限できるようにしたいと考えています。これを実証するために、手動で行います。curvemeasmeasMeasindex == 2index == 1

まず、meas(赤い破線)のインデックスのいくつかをcurve(黒い線) のインデックスに固定する必要があります。最初と最後のポイントに加えて、最も高い値を持つポイントを固定することを選択します。

次に、これらの固定点の間で直線的な成長を仮定します。成長はcurve指数に沿って進みます。Curveには 1 日 1 つの値があります。したがって、成長は毎日の規模になります。

この後、インデックスを日付に置き換えて結果をプロットします。

上のプロットは、3 つのアンカー ポイントに基づく結果の適合を示しています。以下のプロットは、1 日間隔で時間に沿ったモデル化された成長を示しています。3 月初めの成長曲線の曲がりは、私には理解できない面白い数学的アーティファクトです (パッケージのna.approx関数による)。zoo

私は何を試しましたか

以前の質問から、ダイナミック タイム ワーピングが解決策になる可能性があることがわかりました。また、dtw 関数を含むR パッケージも見つけました。良い。実際、動的タイムワーピングは、その質問のサンプルデータセットで機能しましたが (制約の設定を除く)、 (と呼ばれる)curveよりもはるかに多くのデータポイントがあるこのデータセットでは機能しません。measpoints前の質問で）。スペースを節約するために、ここにコード/図をコピーしません。その質問に対する私の答えで、私が試したことを見ることができます。問題は、最も単純なパターンを除いて、どのステップ パターンもこれらのタイプのデータを処理できないことです。最も単純なステップ パターンは、測定値を曲線に数回一致させます。これは避けたいことです。測定ポイントごとに定義された日付が必要だからです。また、測定ポイント間で成長率が 0 を超えなければならないという制約を設定するのは難しいようです。

質問

私の質問は 2 つあります。まず、動的タイム ワーピングよりも問題を解決するためのより良い方法はありますか? 第二に、Rで実際にこれを行うにはどうすればよいですか? .

EDITS 9. 2013年12月私は質問をより明確にしようとしました。