問題タブ [non-linear-regression]

以下のような非線形回帰があります。

y_t= (シータ) sum(((1-シータ)^k) (x_(tk))) + アルファ * (ガンマ-z_t)

x_t と z_t は独立変数であり、theta と alpha と gamma はパラメーターです。パラメータを推定するために、次の再帰関数を書きました。

シータ値を知って回帰モデルの最初の文を計算します。このモデルのシータやその他のパラメータを推定するにはどうすればよいですか?

モデルの非線形性と標準誤差帯域による制限された 3 次スプラインを含む回帰の予測線をプロットしたいと思います。予測されたポイントを取得できますが、線とエラー バンドをプロットするかどうかはわかりません。ggplot が推奨されますが、ベース グラフィックスでも問題ありません。ありがとう。

ドキュメントの例を次に示します。

モデルの予測値のプロット:

関数を使用してnlsLM、非線形回帰を当てはめています。nlsLMモデル オブジェクトから帽子の値とクックの距離を抽出するにはどうすればよいですか?

nlsまたは関数を使用して作成されたオブジェクトを使用してnlreg、観測値のハット値とクック距離を抽出する方法は知っていますが、 を使用してそれらを取得する方法がわかりませんnslLM。

誰でもこれについて私を助けることができますか? ありがとう！

(2 次元) 線分のセットがあります。線分に平行な 2 次曲線を当てはめたいと思います。f(x,y) = ax^2 + b xy + cy^2 + dx + ey + f = 0 のような暗黙の関数を使用してこれを行いたいと思います。また、曲線がで開始する必要があるため、f を決定することができます。

これまでに試したことは次のとおりです。線分に垂直な線を計算しました。次に、線のグラデーションに垂直な曲線の交点にグラデーションを設定したいと思います（したがって、線分に等しくなります）。残念ながら、2 つの問題があります。1) 2 つの解 (交点) があるかもしれません。2) f の係数がわからない限り、交点はわかりません。最終的に、解き方がわからない非線形方程式系が得られます。

ここまでは、線形方程式系を解くために特異値分解を使用しました。

最初の X が 36 個の機能である一連の機能 X から変数 y を予測しようとしています。これに関して 2 つの質問があります。

1. 多項式機能の作成中にブール属性 (0,1) を処理する方法は? たとえば、それらを二乗しても意味がありません。

私がこれまでに持っているコード:

2. 多項式回帰の特徴選択を行う方法は? 36 個の変数に対して次数 2 の多項式機能を作成すると、X のサイズが大幅に増加するためです。たとえば、MSE に基づいて最適なモデルを返す選択を実行する方法はありますか?

予測/適合値を取得するために次を実行しています:

ただし、これは私に単一の値しか与えません。すべての予測値を取得するにはどうすればよいですか?

1. DRC パッケージの曲線がプロットに正しくフィットしない
2. Broken = TRUE グラフィカル パラメータがプロットの外に表示される
3. type = "confidence" を適用するとエラーが発生する

こんにちは、みんな、

R パッケージの DRC で drm をプロットするのに問題があります。

私のデータは 3 つの変数で構成されています。濃度 (mg/L でのアッセイで使用される化合物の濃度の数値ベクトル)、除草剤 (化合物の種類の文字ベクトル)、阻害 (試験処理ごとの阻害の数値ベクトル)。

用量反応曲線は正常に機能し、個々のプロットとモデルは機能します。

問題は、化合物のうちの 2 つの曲線が、その除草剤の最大阻害を超え続けていることです。

y 軸を再スケーリングしましたが、問題は解決しません。

データフレームを再編成して、代わりに別の除草剤が一番上になるように (つまり、A または B の代わりに D または E)、継続する不適切な曲線が最初の 2 つの除草剤に影響を与えるようにしました。

前もって感謝します！

私のコードを以下に示します。エラー (Q 3) は最後の行です。

編集、マイデータ: