問題タブ [nls]

条件付き nls を R 2.15.1 に適合させようとしています。R 2.13 では同じコードが正常に動作していましたが、R 2.15.1 ではエラーがスローされます。

エラーは

編集済み

私は次のモデルに適合しています：

Rの非線形回帰を解決しようとしていますが、構文エラーメッセージが表示されます。私はデバッグを試みましたが、同僚を募集しましたが無駄になりました。私はRを初めて使用するため、手順が正しく設定されているかどうかわかりません。

どんな洞察も素晴らしいでしょう-

ありがとう

toy_data_set：

対数関係にあると思われるデータポイント (x と y) がいくつかあります。

3ここで、グラフに適合し、他のデータポイント (または)を推測できる基になる関数を見つけたいと思います82。私はについて読みましたがlm、nls実際にはどこにも行きません。

最初に、プロットに最も似ていると思われる関数を作成しました。

その後、以下を使用してフィッティング モデルを生成しようとしましたnls。

ここから何をすべきかについて、誰かが私を正しい方向に向けることができますか?

ファローアップ

あなたのコメントを読み、もう少しグーグルで調べた後、 の開始パラメータを調整するとa、突然モデルが収束しました。bc

非線形最小二乗法で推定するモデルがあります。モデルは、合計 4 つの係数に依存する独自の式を持つ 5 つの係数を持つ主な式を持つように指定されています。最初の式自体は線形であり、モデルを非線形にする係数の個々の式を含めることです。

R の nls 関数は、係数のない式を要求します。その場合、それらを定義する方法がわかりません。

物事を明確にするために、以下に具体的な式を掲載します。代数が含まれる場合は、インデックスに「_」を使用し、そのようなインデックスのコンテンツに「{}」ブラケットを使用します。構文は非常に TEX に似ており、変数名は 1 文字にすぎないため、乗算が発生するすべての場所にわざわざ "*" を配置しませんでした。

データ セット全体に対して非線形最小二乗回帰を実行し、そのデータ セットのいくつかのサブセットに対して回帰を繰り返す必要があります。単一のサブセットに対してこれを行うことができます。例 (ここで、y は一般的なロジスティック方程式であり、x は 1 から 20 までのベクトルです):

しかし、これを 3:20、4:20、5:20 などで実行したいので、for ループを試しました。

しかし、「xj [i]のエラー：無効な添字タイプ 'リスト'」が表示されます。nls() を 20 回コピーして貼り付ける必要はありません。どんな助けでも大歓迎です。

ヘルプ (predict.nls) で述べたように、se.fit=TRUE の場合、予測の標準誤差を計算する必要があります。ただし、次のコードはそれを表示せず、予測のみを表示します。

コードの何が問題になっていますか? ありがとうございました！

Rでモデルを推定し、係数の1つを他の係数よりも小さくする制約を追加したいと思います。どうやってやるの？

加重モデルと非加重モデルの偏差を AIC 値と比較すると、いくつかの不一致が発生します。一般的な例 (「nls」から):

これは重み付けされていない適合です。「nls」のコードでは、「nls」がベクトルを生成することがわかりますwts <- rep(1, n)。

加重フィットの場合:

ここでは、2 つの反復で 8 つの濃度のそれぞれに増加する重みを割り当てます。

今deviance私は得る：

加重適合の逸脱度が大幅に高い（適合が悪い）ことを教えてくれます。

今AIC（またはBIC）で私は得る

これは、2番目の適合が桁違いに優れていることを示しています（AICが低い）。なんでそうなの？

AIC教科書にある残差二乗和に基づいて定義すると

私は得る

stats:::AIC.logLikこれは、対数尤度に基づく上記の AIC ( ) と同じ値です。

しかし

これはより高く、2 番目のモデルのより高い逸脱にも完全に適合します。

誰でも私を啓発できますか？加重フィッティングの場合、「nls」モデルの標準 AIC 実装は適用されませんか?

乾杯、アンドレイ

nlsループにフィットするために残差のブートストラップを実行したいと思います。私はを使用nlsBootし、計算時間を短縮するために、それを並行して実行したいと思います（現時点ではWindows 7システムで）。これが私の問題を再現するいくつかのコードです：

これは環境の問題だと思います。問題のコードを確認したところ、呼び出しnlsBootで無名関数を使用したことが原因のようです。lapply

nlsBoot並列ループで使用する方法はありますか？または、関数を変更する必要がありますか？（for代わりにループを使用することもできますlapply。）

Oracle 11g r2 を使用しており、言語ソートを使用してソート順を設定するように DB を構成しようとしています。やった

Oracleを再起動した後、これらのパラメーターを確認しました：

それは私に正しい値を示しています。しかし、並べ替えると

結果が正しくない順序で表示されます。つまり、最初に数字、次に文字です。

しかし、もし私がするなら

それは私に正しい順序を示しています。

システムの変更が正しい結果を表示しない理由を知っている人はいますか?