問題タブ [gamma-distribution]

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

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

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

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

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

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

R で optim() を使用して積分を含む可能性を解き、2 つのパラメーターのヘッセ行列を取得するのに問題があります。アルゴリズムは収束しましたが、optim() で hessian=TRUE オプションを使用するとエラーが発生します。エラーは次のとおりです。

integrate(integrand1, lower = s1[i] - 1, upper = s1[i]) のエラー: 非有限関数値

NAの警告メッセージもありました

これが私のコードです：

ご協力いただきありがとうございます！

ガンマ累積分布関数のポイントごとのプロファイルの可能性に基づく信頼区間の構築を手伝ってくれる人はいますか?

推定されたガンマ パラメータの信頼区間に基づくプロファイル尤度の構築については認識していますが、推定された cdf については認識していません。理論とRプログラミングを使ったコーディングを手伝ってください!!

前もって感謝します。

Gamma(alpha=scale=29, beta=size=3) をプロットしようとしていますが、エラー TypeError: _parse_args() got an unexpected keyword argument 'size' for y1. ガンマ分布のscipy ドキュメントに従っていますが、ガンマ分布にはサイズ パラメータがあるため、このエラーが発生する理由がわかりません。

エラーメッセージ