質問をすることでしか得られない、回答やアドバイスがある。

15分調べてもわからないことは、質問しよう!

ただいまの
回答率

89.53%

fitdistでNumericリストをわたしているのにエラーになる

解決済

回答 1

投稿

  • 評価
  • クリップ 0
  • VIEW 107

canadanac

score 5

ある統計量をガンマ関数でフィッチングさせてみようと思い、fitdistrplusというパッケージで以下のようにfitdistを使いましたが、エラーになってしまいます。フィットさせる数字になにか制限があるのでしょうか。

$summary(vx)
Min. 1st Qu.  Median    Mean 3rd Qu.    Max. 
806    2029    3063    3262    4286    7483 
$vx
[1] 2431 2497 6571 4327 1403 6479 2519 4404 2563 3397 1302 2432 4229 1823 2444
[16]  806 2685 4136 3379 2318 3283 2536 3966 4437 3093 3814 4450 1395 6912 6188
[31] 6322 2444 2294 4321 4271  891 2327 1604 1474 1187 1549 2605 1025  939 5819
[46] 7483 7063 3542 3576 1623 3995 3672 1348 1408 6837 4902 1394 3123 1178 3294
[61] 5083 3677 5189 4710 4286 3633 5361 4844 1748 1096 6005 5797 3552 2641 3537
[76] 2479 2818 3703 2256 2007 1629 3063 6511 3561 2025 1521 1936 4771 2151 4943
[91] 3424 2892 2855 2607 2083 3732 2590 2486 4056 2531 4294 1865 1458 2198 2012
[106] 1183 2223 1877 1595 1961 3260 3089 4408 1429 2029 3630 2255 3260 2693 7153
[121] 4093 4311 3934 3177 4110 4216 7067 2891 4552 2802 3082 4718 1226 2794 2133
[136] 4530 5564 4532 2391 1237 1193
>$ fit=fitdist(vx,"gamma")
<simpleError in optim(par = vstart, fn = fnobj, fix.arg = fix.arg, obs = data,     gr = gradient, ddistnam = ddistname, hessian = TRUE, method = meth,     lower = lower, upper = upper, ...): non-finite finite-difference value [2]>
Error in fitdist(vx, "gamma") : 
the function mle failed to estimate the parameters, 
with the error code 100

  • 気になる質問をクリップする

    クリップした質問は、後からいつでもマイページで確認できます。

    またクリップした質問に回答があった際、通知やメールを受け取ることができます。

    クリップを取り消します

  • 良い質問の評価を上げる

    以下のような質問は評価を上げましょう

    • 質問内容が明確
    • 自分も答えを知りたい
    • 質問者以外のユーザにも役立つ

    評価が高い質問は、TOPページの「注目」タブのフィードに表示されやすくなります。

    質問の評価を上げたことを取り消します

  • 評価を下げられる数の上限に達しました

    評価を下げることができません

    • 1日5回まで評価を下げられます
    • 1日に1ユーザに対して2回まで評価を下げられます

    質問の評価を下げる

    teratailでは下記のような質問を「具体的に困っていることがない質問」、「サイトポリシーに違反する質問」と定義し、推奨していません。

    • プログラミングに関係のない質問
    • やってほしいことだけを記載した丸投げの質問
    • 問題・課題が含まれていない質問
    • 意図的に内容が抹消された質問
    • 過去に投稿した質問と同じ内容の質問
    • 広告と受け取られるような投稿

    評価が下がると、TOPページの「アクティブ」「注目」タブのフィードに表示されにくくなります。

    質問の評価を下げたことを取り消します

    この機能は開放されていません

    評価を下げる条件を満たしてません

    評価を下げる理由を選択してください

    詳細な説明はこちら

    上記に当てはまらず、質問内容が明確になっていない質問には「情報の追加・修正依頼」機能からコメントをしてください。

    質問の評価を下げる機能の利用条件

    この機能を利用するためには、以下の事項を行う必要があります。

回答 1

checkベストアンサー

0

the function mle failed to estimate the parameters,

vxの値がデカすぎるのではないでしょうか。
fitdistに渡すのは確率変数などではないですか?
fitdistは使ったことないので間違ってたらすみません。

> fitdist(vx/10000, "gamma")

Fitting of the distribution ' gamma ' by maximum likelihood 
Parameters:
       estimate Std. Error
shape  4.179747  0.4792592
rate  12.813815  1.5611958

投稿

  • 回答の評価を上げる

    以下のような回答は評価を上げましょう

    • 正しい回答
    • わかりやすい回答
    • ためになる回答

    評価が高い回答ほどページの上位に表示されます。

  • 回答の評価を下げる

    下記のような回答は推奨されていません。

    • 間違っている回答
    • 質問の回答になっていない投稿
    • スパムや攻撃的な表現を用いた投稿

    評価を下げる際はその理由を明確に伝え、適切な回答に修正してもらいましょう。

  • 2020/02/05 18:01

    アドバイスありがとうございます!
    fitdist(vx/1000, "gamma")でも大丈夫だったので確率変数でなくてもよいみたいですが、値が大きすぎたようです。
    ありがとうございました!

    キャンセル

15分調べてもわからないことは、teratailで質問しよう!

  • ただいまの回答率 89.53%
  • 質問をまとめることで、思考を整理して素早く解決
  • テンプレート機能で、簡単に質問をまとめられる