• ベストアンサー

ワイブル分布について

エクセルでワイブル分布を求めました。 求めたワイブル分布を元にランダムに その分布に入る値を求めたいのですが どのようにしたらよいのでしょうか よろしくお願いします。

  • conc
  • お礼率91% (11/12)

質問者が選んだベストアンサー

  • ベストアンサー
  • shkwta
  • ベストアンサー率52% (966/1825)
回答No.1

Excelの=WEIBULL(X,α,β,TRUE)では、ワイブル累積分布関数を F(X)=1 - EXP(-((X/β)^α))) で計算するようになっています。 そこで、この逆関数を考えます。 X=β*((-ln(1-F(X))^(1/α)) このF(X)を、0以上1未満の一様乱数 RAND() で置き換えて、 X=β*((-ln(1-RAND())^(1/α)) とすると、得られたXはワイブル分布に従います。

conc
質問者

お礼

ありがとうございました うまくいきました 逆関数は発想になかったので 大変参考になりました

関連するQ&A

  • ワイブル分布

    物の故障などに使われる、ワイブル分布 について教えて下さい。よろしくしくお願いします。

  • ワイブル分布の確率密度関数と累積分布と関係

    初歩的な質問で申し訳ないのですが、どうしてもわからないので 質問させていただきました。 ワイブル分布で、故障率をプロットしたいのですが、 このときエクセルのワイブル関数で確率密度関数と累積分布をプロットすると、 以下のような数字になります x 累積分布関数       x 確率密度関数 0 0.0%       0 0.0% 0.3 0.1%      0.3 0.7% 0.6 0.8%      0.6 5.4% 0.9 4.0%      0.9 17.5% 1.2 12.2%      1.2 37.9% 1.5 27.1%      1.5 61.5% 1.8 48.1%      1.8 75.7% 2.1 70.3%      2.1 68.7% 2.4 87.4%      2.4 43.5% 2.7 96.4%      2.7 17.8% 3 99.4%       3 4.3% 3.3 99.9%      3.3 0.5% 3.6 100.0%      3.6 0.0% 確率密度関数の値を累積したものが累積分布になると思っていたのですが、 累積分布の値はそのような数字になりません。 確率密度関数の値を累積したものが累積分布にならないのはなぜでしょうか。 それぞれの使い方が違うのでしょうか。 そうであれば故障率としてはどちらを使えばいいのでしょうか。 本当に初歩的な質問で申し訳ございませんが、ご教授いただきたくお願い申し上げます。

  • ワイブル分布より裾の厚い分布は何がありますか?数値

    ワイブル分布より裾の厚い分布は何がありますか?数値は0<x<無限の範囲です。

  • ワイブル分布の標準偏差

    統計初心者です。 ワイブル確率紙のプロットからmの値と平均寿命、標準偏差が求められるようですが、標準偏差はどのように使うのでしょうか?正規分布していないものの標準偏差とは?何でしょうか? 平均寿命±3シグマで99.7%がその範囲に入る??(正規分布ではないので違うと思いますが?) よろしくご教授願います。

  • ワイブル分布を使った推定

    ある製品が販売後t年経っても廃棄されずに残っている率(残存率Y(t)と呼ぶことにします)のデータがt=1からt=10まであります。これを元に、t=11からt=50までの残存率を推定したいと思っています。 ワイブル分布を説明した http://ja.wikipedia.org/wiki/%E3%83%AF%E3%82%A4%E3%83%96%E3%83%AB%E5%88%86%E5%B8%83 を見て、ワイブル・プロットという方法を使うと次のようにできると考えたのですが、妥当なやり方でしょうか。 ・t=1~10の時のln(t)を説明変数、ln(ln(1/Y(t)))を被説明変数とした単回帰を行い、回帰式を求める ・回帰式にt=11~50を代入して、ln(ln(1/Y(t)))の推定値を求める(これをZ(t)と書くことにする) ・1/(exp(exp(Z(t)))を計算することにより、Y(t)を求める。 Excelで作ってみたところ、それっぽい曲線はできたのですが。 よろしくお願いいたします。

  • 指数分布族としてのワイブル分布

    指数分布族は対数尤度関数がLL(θ,φ)={xθ-A(θ)}/φ + c(x,φ)という共通のかたちをもつそうです(θはcanonical parameter,φはdispersion parameterと呼ばれるものです)。ワイブル分布の確率密度関数が f = k λ^(-k) x^(k-1) exp[-(x/λ)^k]と書けるときに、対数尤度関数は LL = k ln(x/λ) + ln(k/x) - (x/λ)^k になりますが、このときに、θ、φ、A(θ)、c(x,φ)はそれぞれどのように書けるのでしょうか?(c(x,φ)はexplicitに書けない??)ご存知のかたがいらっしゃいましたらお教えいただけると助かります。 念のため上の式をTeX形式で書くと f=k\lambda^{-k}x^{k-1}\exp\left[-\left(x/\lambda\right)^k\right] LL=k\ln\left(x/\lambda\right)+\ln\left(k/x\right)-\left(x/\lambda\right)^{k} となります。

  • 指数分布とは

    こんにちは。 ワイブル分布で、m=1のときは指数分布に従うとあります。 指数分布とは、時間と共に不良率は減少していくと思いますがなぜそうなるのですか? おしえてください

  • ワイブル関数累積分布関数でのt=ηの点はなんと説明しますか。

    ワイブル関数累積分布関数F(t)=1-exp(-(t/η)^m)でm=1、t=ηの点は63%ですがこの点をなんと説明すればよいですか。tを機械の寿命とすればηは何ですか。わかっているようで、わからない。なさけない。教えてください。

  • ワイブル分布・確率密度関数について

    ワイブル分布に関する確率密度関数の導出が理解できずに困っています。 まず、2母数ワイブル分布の確率密度関数を f (N |α,β) =(α/β)*(N/β)^(α-1 )*exp{ -(N/β)^α} ・・・ (1) N : 寿命 α : 形状母数(既知とする) β : 尺度母数(未知とする) とします。ここで、試験数 n 、N={N1,N2,・・・Nn}の結果が 得られた場合、βの最尤推定値Βは、 Β = { (ΣNi^α)/n }^(1/α) ・・・ (2) [ i = 1→n ] となります。Βは確率変数で、Βの確率密度関数を次のように表せるそうです。 f (Β) = (n^n / Γ(n) )*(α/β)*(Β/β)^(α*n-1)*exp{ -n*(Β/β)^α} ・・・ (3) Γ(n) : ガンマ関数 (2)式の導出は何とかわかったのですが、(3)がなぜこのような式に なるのかわかりません。(3)の導出方法をどなたか教えてください。 よろしくお願いします。 ※式の誤解がないように不要な括弧を多用しており、見にくいかもしれません。 ソースは下記のPDFファイルの3ページ目(表紙・目次除く)にありますので、 そちらの方が見やすいと思います。 また、下記の資料によると、「変数2n(Β/β)^αが自由度2nのカイ二乗分布に 従う」と書かれていますが、これも良くわかりません。できれば、こちらの方も よろしくお願い致します。 http://repository.tksc.jaxa.jp/dr/prc/japan/contents/AA0065137000/65137000.pdf 「宇宙航空研究開発機構研究開発資料―疲労寿命安全率に関するノート」

  • ワイブル確率紙

    ワイブル確率紙にはグラフ真ん中上部に○印があります。これはどういうものでしょうか?単なるグラフの書き始めの点なのでしょうか?それともやはりここにあるべき理由があるのでしょうか?教えてください。宜しくお願いします。

専門家に質問してみよう