問題タブ [fitdistrplus]

正規分布に合わせて次のコードを使用しています。「b」のデータセットのリンク (大きすぎて直接投稿できません) は次のとおりです。

bのリンク

ヒストグラムをプロットした後、正規分布がうまく適合するはずです。

xy.coords(x, y, xlabel, ylabel, log) のエラー:

'x' はリストですが、コンポーネント 'x' と 'y' がありません

当てはめた分布をプロットしようとすると上記のエラーが発生し、要約統計量でさえ f1n がオフになっています。

助けていただければ幸いです。

異なるデータセットを使用して、qq プロットを 1 つずつ取得することしかできません。

2 つの qq プロットを同じ座標系に配置するにはどうすればよいですか?

fitdistrplus を使用してワイブル分布をさまざまなデータ セットに適合させ、生データを CDF にプロットするコードがあります。

この計算から適合の残差にアクセスして、+/- 5 疑似シグマ スクリーンを適用し、外れ値の数を特定してカウントしたいと考えています。これらは外因性の欠陥母集団と見なされるためです。

残差のベクトル、リスト、または行列を取得するにはどうすればよいですか。fitdistr は resid() 関数を使用します。fitdistrplus はこの関数を継承する必要がありますが、resid(zFIT) を呼び出しても機能しないようです (ここで、zFIT は私の fitdist オブジェクトです.....)

何か案は？どうもありがとう

いくつかのデータがあり、それらを分布に当てはめたとします。 x がサンプル外のデータ ポイントである場合、gammaの間隔確率を見つける方法は?Pr(1 < x <= 1.5)

ガンマ分布を 338 要素で構成されたデータセットに、固定の低しきい値(私は R を使用しています)で適合させたいと考えています。下限を表すために、位置を追加して 3 つのパラメーターを持つガンマを使用することを考えました。これが私のコードです：

コードを実行するたびに、同じエラーが発生しました。

そこで、ガンマを 2 つのパラメーターでフィッティングして得られた値を使用して、開始値を変更しようとしました。これにより、次の結果が得られました。

しかし、それでもうまくいきません..2 つのパラメーターを持つガンマでさえうまくいかない場合は理解できますが、そうではなく、自分で説明することはできません。さらに、QQ プロットと 2 つのパラメーターを持つガンマの ecdf はあまり良くありません...しかし、低しきい値に関してスケーリングされた元のデータセットの分布に適合すると、適合は完璧に見えます:

しかし、そのようにすることが正しいかどうかはわかりません...パラメータは非常に異なり、データに関連する返品期間を計算するために必要です! そのため、位置パラメーターを使用してガンマについて考えました。

psデータが多すぎるため添付しませんでしたが、それらの概要を報告できます。

以下のコードを実行しています。つまり、負の二項分布の AIC 値をベクトルに保存します。

次のエラー メッセージが表示されます。

fitdist(a, distr = "nbinom") のエラー: 関数 mle がパラメーターの推定に失敗し、エラー コード 100 が追加されました: 警告メッセージ: 1: In sqrt(diag(varcovar)): NaN が生成されました 2: In sqrt (1/diag(V)) : 生成された NaN 3: cov2cor(varcovar) : diag(.) には 0 または NA エントリがありました。非有限の結果は疑わしい

また、出力には 62 個の値しか含まれていませんが、336 セットの分布があるため、336 個の値を期待しています。