問題タブ [nls]

次のコードを再現しようとしています ( nls 関数はうまく機能しません) が、for ループ、sprintf、および as.formula() を使用して追加の実装を行い、特定のスペクトルのピーク数に応じて変数を追加します。ピーク間で一貫性を持たせるために、各ピークの変数名をベクトル化したので、ピーク番号 1 には関連する 'alfa[1]'、'peak[1]'、'height[1]' があります。

これまでのところ、期待される式が得られました。

parそれにもかかわらず、同じシステムを回線に複製しようとすると、いくつかの問題が発生します。次のように表示されます。

しかし、代わりに、すべての文字列を折りたたんでas.formula後でコマンドを使用すると、次のようになりました。

折りたたんだ文字列を出力すると文字列は想定通りなので、どうにかして as.formula コマンドに繋げられるのではないかと考えています(つまり、流用されたコマンドではないかもしれません)

ggplot2 を使用して、いくつかのデータをプロットし、非線形曲線を当てはめようとしています。すでに定義した nls オブジェクトで stat_function を使用したいのですが、結果として切り捨てられた曲線が生成されます。次のページを見ましたが、これまでのところ解決策が見つかりませんでした。

R - 線形領域の ggplot2 外挿回帰線

stat_function() を使用して R の ggplot で多数のカスタム関数をプロットする

http://docs.ggplot2.org/0.9.3/stat_function.html

ggplot の curve() と同等

私は ggplot2 の経験があまりないので、単純なものが欠けている場合は申し訳ありません。助けていただければ幸いです。ここに私のサンプルデータ:

ggplot2 でプロットを作成するためにこれまでに持っているものは次のとおりです。

これらはすべてうまく機能しますが、nls オブジェクトを作成して stat_function で使用しようとすると、切り捨てられた曲線が得られます。曲線を x 軸まで延長したいと思います。また、エラーが生成されます。

私が欲しいのは、基本的な r 機能を使用した次のようなものですが、ggplot2 によって提供されるすべてのベルとホイッスル、および最終的には各処理 (Tx) に対して 1 つの nls 曲線を備えています。

これは、nls と LPPL の両方が私にとってかなり新しいものであるという点で、これまで R で行った中で最も困難なことです。

以下は、私が取り組んできたスクリプトの一部です。df は、S&P 500 の終値である日付と Y の 2 つの列で構成されるデータ フレームです。関連性があるかどうかはわかりませんが、日付は 2003 年 1 月 1 日から 2007 年 12 月 31 日までです。

実行すると、次のメッセージが表示されます。

LPPL の公式は、この PDF ファイルの 5 番目の画面にあります。

私がどこで間違っているか知っていますか？これは別のモデルでは正しく機能していたので、新しい方程式のコードを変更しました。この投稿のこのコードについては、jlhoward の功績によるもので、R で nls を使用して研究を再作成します。

ご協力ありがとうございました。

jlhoward のコメントによると、df.rda はここからダウンロードできます: https://drive.google.com/file/d/0B4xAKSwsHiEBb2lvQWR6T3NzUjA/edit?usp=sharing

y=a^(b^x)以下のデータに合わせて使用​​したいのですが、

非線形最小二乗法を使用すると、

エラーが発生します:初期パラメーター推定値での特異勾配行列。結果はおおよそ a = 1.1466、b = 0.6415 であり、a=1、b=0.5 と定義しているため、最初のパラメーター推定値に問題はないはずです。

他のトピックで、曲線を変更すると便利であると読みました。みたいなことを考えていたのですlog y=log a *(b^x)が、関数指定の扱いがわかりません。何か案が？

3 つのモデルのカーブ フィットを r 二乗値で比較したいと考えています。nlsおよびdrcパッケージを使用してモデルを実行しました。ただし、どちらのパッケージも r-2 乗値を計算していないようです。ただし、「残差標準誤差」と「残差二乗和」が得られます。

これら 2 つは、モデルの適合度を比較するために使用できますか?

nls() 関数で情報を合わせようとしています。式は次のとおりです。

しかし、Rは怒っています:)

Rを使用して次の式に当てはめたいデータがあります。

Z(u,w)=z0*F(w)*[1-exp((-b*u)/F(w))]

ここで、z0 と b は定数で、F(w), w=0,...,9 は、F(0)=1 および u=1,...,50 の w に依存する減少ステップ関数です。

Z(u,w) は、50x10 行列 (行の横に u=50,...,1、列に沿って w=0,...,9) の形式で観測されたデータのセットです。たとえば、それほど素晴らしいことは説明していませんが、Z(42,3) は 9 行下、4 列目の要素になります。

F(0)=1 を使用すると、次のコードで最初の列 (つまり w=0) だけを使用して b と z0 の推定値を取得できました。

次に、各列を調べて見つけた b と z0 の値を使用して、w=1,...,9 の F(w) を見つけました。

ただし、b と z0 は最初の列だけでなく、すべてのデータに適合させる必要があるため、12 個のパラメーターすべて (b、z0、および F(w) の 10 個の値) を一度に推定する方法を見つけたいと考えていました。

これを行う方法を知っている人はいますか？すべての助けをいただければ幸いです。

ありがとうジェームズ

Rでグラフ内の複数の線の「一般的な」方程式を見つける方法はありますか? 次のグラフは、私が達成したかったことの例です。t の主要な線のみを示しましたが、それらの主要な線の間に隠れた線もあります。私が望むのは、線の「一般的な」方程式を取得することです。どうすればいいですか？

このRの問題を解決する方法を知っている人はいますか? 以下のデータの x=5 のように、関係の変化点を見つけることです。

SAS は見事にうまく機能します。

R は失敗します。

信号が明確な単純なケースなので、R が機能するはずです。Rで何が間違っているのだろうか（https://stats.stackexchange.com/questions/7527/change-point-analysis-using-rs-nlsのコードを使用しました）。誰でも提案/解決策はありますか?

ありがとう！

ウィレム