問題タブ [non-linear-regression]

次の非線形モデルのパラメーターを暗黙の形式で推定しようとしています。その式は y = f(x;theta) + u ではなく、u = g(x,y;theta) で与えられ、u はランダム項です。 x と y は、説明変数と被説明変数です。g から f を取得することはできません (g は y で反転できない、および/または g は u で加法的ではないため)。以下の NLS と GMM はどちらも実行されません。NLS および GMM のドキュメントでこの問題について何も見つけることができなかったので、どうすればよいのだろうか。このトピックに関するヘルプをいただければ幸いです。

より具体的には、g の次の仕様を検討しました (ただし、NLS と GMM は実行されません)。

助けてください。以下を使用して、かなり適切な非線形回帰曲線をデータに適合させることができました。

次に、次を使用してプロットを生成します。

95% 信頼区間の上限と下限に点線を引きたいと思います。線形回帰分析のコードを書くことはできましたが、現在の非線形回帰に苦労しています。次のコード:

線形回帰の場合、これにより、fit、lwr、upr の 3 つの列を持つデータフレームが得られます。私の非線形回帰では、何もできない数値のリストが表示されます。助言がありますか？非線形回帰の上限線と信頼線を取得する簡単な方法はありますか?

R で次のデータの等温線モデルを当てはめたいと思います。最も単純な等温線モデルは、ここで指定されたラングミュア モデルです。モデルはページの下部で指定されています。エラーをスローする私のMWEを以下に示します。Isotherm モデル用の R パッケージがあるのだろうか。

編集済み

一部の非線形モデルは、線形モデルに変換できます。私の理解では、非線形モデルの推定値とその線形モデル形式の間には 1 対 1 の関係があるかもしれませんが、対応する標準誤差は互いに関連していません。この主張は本当ですか？線形性に変換して非線形モデルをフィッティングする際に落とし穴はありますか?

多数の非線形モデルをデータセットに当てはめようとしていますが、モデル パラメータ値をすべて正の値に制限する必要があります。

lsqcurvefit 関数を試してみましたが、うまくいきました。ただし、モデルの比較には AICc スコアも必要です (明らかに、調整された R-2 乗スコアは適切ではありません)。lsqcurvefit はそのようなスコアを計算しないようです (?)。

そこで、AICc を返す fitnlm を試してみました。ただし、今回はパラメーター値をバインドできません (その結果、意味のない負のパラメーター値がいくつか得られました)。

ユーザーが同時に AICc スコアを計算し、同時にパラメーター値をバインドできるようにする MATLAB 関数を知っている人はいますか? lsqcurvefit 出力から AICc を計算する方法などの代替ソリューションも非常に高く評価されます。

非線形モデルのパラメータ推定に GA パッケージの Genetic Algorithm 関数を使用しています。すべての変数に多重共線性の問題があるシミュレーションデータを使用します。GA 関数を使用すると、GA 関数で推定したパラメーターに偏りがあることがわかりました。なぜそれが起こったのか説明してもらえますか？多重共線性データを使用して GA 関数を使用してパラメーターを推定するソリューションは何ですか? ありがとう。

私は長期的に定期的に収集された気象データを分析しています (ほとんどのデータで 15 ～ 60 分ごと)。気温や日射量の影響には日周期があります。私は、放射線が雲によって遮られていない場合、1 年のうちの任意の日の太陽放射線への平均露出を説明しようとしています。私は数年分のデータにアクセスでき、R に入力したデータをその年の日付に従って平均化できます。雲のない日の平均日射量を説明するために、平均を行う前にいくつかのデータを捨てる必要があります。

どうやらグラフを載せるほどの評判はないようですが、雲のない日の放射パターンのグラフは放物線のような形をしています。曇りの日は、複数のピークを持つ曲線によって識別できます。2 次回帰の R^2 値を使用して、2 種類の曜日を区別できます。

(編集 -- すべての放射線データと日付/時刻は、1 つのテキスト ファイルの 2 つの列で報告されます。以下のデータを日付ごとに分けて、私が分析しようとしているパターンを読者が簡単に確認できるようにしています。データを共有してパターンを表示するためのより洗練された方法を知りません。)

R^2 の違いは、曇りの日と晴れの日を区別する方法として使用できます。より大きなデータ セット内で毎日このプロセスを繰り返す方法を見つけることができれば.

日付と時刻、またはすべての日の放射線測定値が単一の列内で報告される単一のデータ フレームから複数の二次回帰を行う方法はありますか。

理想的には、2 つの列を持つテーブルで終了したいと考えています。1 つの列には年間通算日が含まれ、2 番目の列には 2 次回帰分析からの R^2 値が含まれます。複数の R^2 または調整済みの R^2 のいずれかが機能すると思います (ただし、R^2 の 2 つのバージョンの違いを十分に理解していないため、一方を優先して使用するよう説得することはできませんでした)もう一方。)

二次回帰分析から R^2 値のみを報告する方法や、分析しているデータの日数だけ二次回帰を繰り返す方法がわかりません。私は潜在的に 10 年間のデータを見ているので、単一の表で分析結果を分析して報告できることは、どの日のデータを使用できるかを分類する素晴らしい方法です。