問題タブ [model-fitting]

現在、次のプロットが生成されています。

データ：

コード：

各コロニーの成長率 (指数関数的成長を想定) を導出しようとしています。これを行う1つの方法は次のとおりだと思います。

1. 指数曲線を各コロニーのデータポイントに合わせてから
2. 曲線の方程式から成長率を抽出します。

これを達成するために nls を使用するか、指数曲線のパラメーターを 3 つのデータポイントから非常に簡単に導き出すことができるようにすることをお勧めします。

残念ながら、私はまだこのアドバイスを受け入れる準備ができていません。どちらの方法もわかりません。

ありがとう！

私はRでいくつかのデータ分析を行っており、データを3パラメーターのワイブル分布に適合させる方法を見つけようとしています。2パラメーターのワイブルでそれを行う方法を見つけましたが、3パラメーターでそれを行う方法を見つけるのに不足しています。

パッケージのfitdistr関数を使用してデータを適合させる方法は次のとおりです。MASS

x[[6]]は私のデータのサブセットであり、yはフィッティングの結果を保存している場所です。

関数f（x）を表す2つのベクトルと、別のベクトルf（a x + b）、つまりf（x）のスケーリングおよびシフトされたバージョンがあります。最適なスケールとシフトファクターを見つけたいと思います。

*最良-最小二乗誤差、最尤法などを使用します。

何か案は？

例えば：

* f（x）は元に戻せない場合があることに注意してください...

ありがとう、

オハド

「基本的なフィッティング」ツールを使用する場合、「フィッティング」が行われると、特定のポイントで値を評価/推定する機会があります。作図部分までしか再現できませんでした。「評価」関数をプログラムで再現して、特定のポイントの値を見積もり、コードで使用できるようにする方法がわかりません。今のところこれを実現できる唯一の方法は、GUIを使用することです。つまり、Figureウィンドウのメインメニューから「ツール>>基本的なフィッティング」を実行します。

自分自身を十分に明確にしているのかどうかはわかりませんが、さらに情報が必要な場合は遠慮なく質問してください。

たとえば、次の2つの変数の散布図があります。

そして、これら 2 つの変数の関係によりよく適合する関数を見つけたいと思います。

linear正確には、 、、 の 3 つのモデルのフィッティングを比較したいと思いexponentialますlogarithmic。

各関数を自分の値に当てはめ、それぞれのケースで尤度を計算し、AIC 値を比較することを考えていました。

しかし、どこから、どのように始めればよいのか、よくわかりません。これに関する可能な限りの助けをいただければ幸いです。

事前にどうもありがとうございました。

ティナ。

このようなデータセットで

このデータを以下のように定義されたモデルに当てはめたい

このページ ( http://www.inside-r.org/packages/cran/FME/docs/modCost ) のリファレンスを参考にして、次のコードを作成しました。

ただし、次の警告メッセージが表示されました。

どこに問題があるのか​​誰か教えてください。それともどうすれば手っ取り早く？ありがとう。

以前の投稿のフォローアップとして、 (Amelia で作成された) 多重代入データを使用して、回帰のために zelig で統計モデルの相対的な品質の通常の尺度を取得する方法を学習することに興味があります。

質問

(1) AIC、F統計量、およびMIデータの自由度の値を取得する方法を知っている人はいますか?

(2)個々のデータセットの要約によって生成された測定値の単純な平均を取ることができると今井幸助が答えた同様の質問を見つけました。「おそらく」は、私を少し疑わしくさせます。これについて何か考えはありますか？

どうもありがとう！！

コスタキスの紙の解を再現しようとしています。この論文では、de Heligman-Pollard モデルを使用して、要約された死亡表を完全な生命表に拡張します。モデルには、適合する必要がある 8 つのパラメーターがあります。著者は修正された Gauss-Newton アルゴリズムを使用しました。このアルゴリズム (E04FDF) は、コンピュータ プログラムの NAG ライブラリの一部です。Levenberg Marquardt は、同じ一連のパラメーターを生成するべきではありませんか? LM アルゴリズムのコードまたはアプリケーションの何が問題になっていますか?