ErlangDistribution[k,λ]
形状母数 k,比率 λ のアーラン(Erlang)分布を表す.
ErlangDistribution
ErlangDistribution[k,λ]
形状母数 k,比率 λ のアーラン(Erlang)分布を表す.
詳細
- アーラン分布における値
の確率密度は,
では
に比例し
では0である. - ErlangDistribution[k,λ]は各 xiがExponentialDistribution[λ]に従う,TransformedDistribution[x1+…+xk,…]に等しい.
- ErlangDistributionでは,k は任意の正の整数でよく,λ は任意の正の実数でよい.
- ErlangDistributionでは,λ は任意の単位次元の数量でよく,k は無次元量でよい. »
- ErlangDistributionは,Mean,CDF,RandomVariate等の関数とともに用いることができる.
予備知識
- ErlangDistribution[k,λ]は,2つの値 k と λ によってパラメータ化された区間
上の連続統計分布を表す.この分布は,独立同分布に従う (それぞれがXiExponentialDistribution[λ]の分布に従う)k 個の確率変数
の和
の分布として定義される.ここで,k は正の整数であり,λ は関連する確率密度関数(PDF)の全体的な高さ/傾きを決定する実数値の「サービス率」である.アーラン分布のPDFは連続で単峰性か単調減少で,領域の下限境界に近付く潜在的な不連続性がある.これに加え,PDFの裾部は,PDFが
の大きい値について指数的にではなく代数的に減少するという意味で「厚い」(この動作は,分布のSurvivalFunctionを分析することで数量的に正確にできる). - アーラン分布は,デンマーク人の数学者かつ統計学者であるAgner Krarup Erlangによって,1900年代の初めに電話回数の頻度のモデル化のために提案された.続く数年でアーラン分布は,待ち行列理論,特に系の現在の状態を与えた場合の待ち行列の条件付き待ち時間のモデル化で定期的に利用されるツールとなった.例えば,すでに3人が並んでいる待ち行列に並んだ場合に,期待される待ち時間
は TErlangDistribution[3,λ]を満足する.ただし,λ は1人あたりの平均応対時間の逆数である. - RandomVariateを使ってアーラン分布から,1つあるいは複数の機械精度あるいは任意精度(後者はWorkingPrecisionオプションを介す)の擬似乱数変量を得ることができる.Distributed[x,ErlangDistribution[k,λ]](より簡略すると xErlangDistribution[k,λ])を使って,確率変数 x が,アーラン分布に従って分布していると宣言することができる.このような宣言は,Probability,NProbability,Expectation,NExpectation等の関数で使うことができる.
- 確率密度関数および累積分布関数は,PDF[ErlangDistribution[k,λ],x]およびCDF[ErlangDistribution[k,λ],x]を使って得られることがある.平均,中央値,分散,原点の周りのモーメント,中心モーメントは,それぞれMean,Median,Variance,Moment,CentralMomentを使って計算することができる.
- DistributionFitTestを使って,与えられたデータ集合がアーラン分布と一致するかどうかを検定することが,EstimatedDistributionを使って与えられたデータからパラメトリックアーラン分布を推定することが,FindDistributionParametersを使ってデータをアーラン分布にフィットすることができる.ProbabilityPlotを使って記号アーラン分布のCDFに対する与えられたデータのCDFのプロットを生成することが,QuantilePlotを使って記号アーラン分布の変位値に対する与えられたデータの変位値のプロットを生成することができる.
- TransformedDistributionを使って変換されたアーラン分布を表すことが,CensoredDistributionを使って上限値と下限値の間で切り取られた値の分布を表すことが,TruncatedDistributionを使って上限値と下限値の間で切断された値の分布を表すことができる.CopulaDistributionを使ってアーラン分布を含む高次元分布を構築することが,ProductDistributionを使ってアーラン分布を含む独立成分分布の結合分布を計算することができる.
- アーラン分布は,(その定義を通して)指数分布を含む他の多くの分布と関連している.ErlangDistributionは, GammaDistribution[k, 1/λ]とPearsonDistribution[3, 1,(-k + 1)/λ, 0, 1/λ, 0]のPDFが両方ともErlangDistribution[k,λ]のPDFと厳密に等しいという意味で,GammaDistributionとPearsonDistributionの両方の特殊ケースでもある.k がInfinityに近付くにつれてErlangDistribution[k,λ]はNormalDistribution[k/λ,
/λ]に収束し,ParetoDistributionはExponentialDistributionおよびErlangDistributionの商として得ることができる.さらに,ErlangDistributionは変換されたChiSquareDistributionに関連しており,ここから,NormalDistribution,LogNormalDistribution,RayleighDistribution,WeibullDistribution,MaxwellDistributionとの関係を継承している.
例題
すべて開く すべて閉じる例 (4)
Plot[Table[PDF[ErlangDistribution[k, .3], x], {k, {1, 3, 5}}]//Evaluate, {x, 0, 22}, Filling -> Axis]Plot[Table[PDF[ErlangDistribution[4, λ], x], {λ, {0.5, 1, 2}}]//Evaluate, {x, 0, 15}, Filling -> Axis, PlotRange -> All]PDF[ErlangDistribution[k, λ], x]Plot[Table[CDF[ErlangDistribution[k, .3], x], {k, {1, 3, 5}}]//Evaluate, {x, 0, 22}, Filling -> Axis]Plot[Table[CDF[ErlangDistribution[4, λ], x], {λ, {0.5, 1, 2}}]//Evaluate, {x, 0, 15}, Filling -> Axis, PlotRange -> All]CDF[ErlangDistribution[k, λ], x]Mean[ErlangDistribution[k, λ]]Variance[ErlangDistribution[k, λ]]Median[ErlangDistribution[k, λ]]スコープ (8)
data = RandomVariate[ErlangDistribution[4, 0.5], 10 ^ 4];Show[
Histogram[data, 15, "PDF"],
Plot[PDF[ErlangDistribution[4, 0.5], x], {x, 0, 20}, PlotStyle -> Thick]]sample = RandomVariate[ErlangDistribution[2, 5.5], 10 ^ 3];edist = EstimatedDistribution[sample, ErlangDistribution[k, λ]]サンプルの密度ヒストグラムを推定分布の確率密度関数と比較する:
Show[Histogram[sample, Automatic, "PDF"], Plot[PDF[edist, x], {x, 0, 1.5}, PlotStyle -> Thick]]DiscretePlot[Skewness[ErlangDistribution[k, λ]], {k, 1, 10}, PlotStyle -> PointSize[Medium]]Skewness[ErlangDistribution[k, λ]]Limit[Skewness[ErlangDistribution[k, λ]], k -> ∞]DiscretePlot[Kurtosis[ErlangDistribution[k, λ]], {k, 1, 8}, PlotStyle -> PointSize[Medium]]Kurtosis[ErlangDistribution[k, λ]]k が大きくなると,尖度はNormalDistributionの尖度に近付く:
Limit[Kurtosis[ErlangDistribution[k, λ]], k -> ∞]FormulaGrid[list_, type_] := Grid[...]FormulaGrid[Table[Moment[ErlangDistribution[k, λ], r], {r, 3}], M]Moment[ErlangDistribution[k, λ], r]FormulaGrid[Table[CentralMoment[ErlangDistribution[k, λ], r]//Factor, {r, 3}], CM]CentralMoment[ErlangDistribution[k, λ], r]FormulaGrid[Table[FactorialMoment[ErlangDistribution[k, λ], r]//FullSimplify, {r, 3}], FM]FormulaGrid[Table[Cumulant[ErlangDistribution[k, λ], r], {r, 3}], C]Cumulant[ErlangDistribution[k, λ], r]Plot[Table[HazardFunction[ErlangDistribution[10, λ], x], {λ, {1, 1.5, 2}}]//Evaluate, {x, 0, 40}, Filling -> Axis]Plot[Table[HazardFunction[ErlangDistribution[k, .3], x], {k, {2, 3, 5}}]//Evaluate, {x, 0, 40}, Filling -> Axis]HazardFunction[ErlangDistribution[k, λ], x]Plot[Table[Quantile[ErlangDistribution[k, 4], q], {k, {1, 3, 5}}]//Evaluate, {q, 0, 1}, Filling -> Axis]Plot[Table[Quantile[ErlangDistribution[10, λ], q], {λ, {.1, .3, .5}}]//Evaluate, {q, 0, 1}, Filling -> Axis]Quantile[ErlangDistribution[k, λ], q]母数でQuantityを一貫して使うとQuantityDistributionが与えられる:
time𝒟 = ErlangDistribution[3, Quantity[5, 1 / "Minutes"]]UnitConvert[time𝒟, "Seconds"]アプリケーション (5)
交通信号による遅れは平均0.5分で指数分布に従うと仮定する.ドライバーが同期されていない7つの信号がある道路を通るとし,すべての信号を通過する場合の遅れの分布を求める:
Mean[ExponentialDistribution[λ]] == Quantity[0.5, "Minutes"]trafficDelayDistribution = ErlangDistribution[7, Quantity[2, 1/"Minutes"]];Probability[d > Quantity[5, "Minutes"], dtrafficDelayDistribution]N[%]電話の通話時間は指数分布に従っていると仮定する.平均通話時間は3.7分である.連続する9回の通話時間が25分よりも長くなる確率を求める:
Mean[ExponentialDistribution[λ]] == Quantity[3.7, "Minutes"]NSolve[%, λ]phoneDurationDistribution = (ErlangDistribution[9, λ] /. First[%])NProbability[d > Quantity[25, "Minutes"], dphoneDurationDistribution]論理素子の遅延は指数分布に従い平均の遅延は
秒であると仮定する.組合せ論理回路内で最長の論理素子のシーケンスが6個であるとして,6個すべての素子を通過する際の遅延が
秒より長くなる確率を求める:
Mean[ExponentialDistribution[λ]] == Quantity[1/1000000000, "Seconds"]NSolve[%, λ]logicDelayDistribution = (ErlangDistribution[6, λ] /. First[%])Probability[d > Quantity[8.5 10^-9, "Seconds"], dlogicDelayDistribution]寿命にA,B,Cの3段階があるデバイスがある.各段階における時間は平均を10時間とする指数分布に従う.段階Cを過ぎるとデバイスは故障する.このデバイスの故障までの時間の分布を求める:
𝒟 = ErlangDistribution[3, Quantity[1/10, 1/"Hours"]];Plot[PDF[𝒟, Quantity[x, "Hours"]], {x, 0, 70}, Filling -> Axis, AxesLabel -> {"h"}]Mean[𝒟]このようなデバイスが少なくとも40時間使用に耐える確率を求める:
NProbability[x ≥ Quantity[40, "Hours"], x𝒟]独立した30個のデバイスの故障までの時間のシミュレーションを行う:
ListPlot[{RandomVariate[𝒟, 30], {{0, Mean[𝒟]}, {30, Mean[𝒟]}}}, Filling -> Axis, Joined -> {False, True}, AxesOrigin -> {0, 0}, AxesLabel -> Automatic]10個のデバイスで始まるシステムがある.デバイスのうち1個だけがアクティブで残りの9個は待機状態である.各デバイスの寿命は母数
でExponentialDistributionに従う.待機中のデバイスがある限り,1つのデバイスが故障するとすぐに別のデバイスに切り換えられる.この場合,このシステムの寿命は次の分布に従う:
𝒟 = ErlangDistribution[10, Quantity[1/300, 1/"Hours"]];Plot[SurvivalFunction[𝒟, Quantity[x, "Hours"]], {x, 0, 7000}, Filling -> Axis, AxesLabel -> {"h"}]Mean[𝒟]Probability[x ≥ Quantity[5000, "Hours"], x𝒟]N[%]このようなシステムを独立で30回実行した場合の寿命のシミュレーションを行う:
ListPlot[{RandomVariate[𝒟, 30], {{0, Mean[𝒟]}, {30, Mean[𝒟]}}}, Filling -> Axis, Joined -> {False, True}, AxesOrigin -> {0, 0}, AxesLabel -> Automatic]特性と関係 (9)
アーラン分布は正の因子によるスケーリングの下では閉じている:
TransformedDistribution[a * u, uErlangDistribution[k, λ]]TransformedDistribution[u + v, {uErlangDistribution[Subscript[k, 1], λ], vErlangDistribution[Subscript[k, 2], λ]}]CharacteristicFunction[ErlangDistribution[k, λ], t] ^ n//PowerExpandCharacteristicFunction[ErlangDistribution[ n k, λ], t]FullSimplify[% - %%]ErlangDistribution[k,λ]は k->∞で正規分布に収束する:
{μ, σ} = {Mean[ErlangDistribution[k, λ]], Sqrt[Variance[ErlangDistribution[k, λ]]]}Block[{λ = 0.01}, Table[Plot[{PDF[ErlangDistribution[k, λ], x], PDF[NormalDistribution[μ, σ], x]}, {x, μ - 2σ, μ + 2σ}, AxesOrigin -> {μ - 3σ, 0}, Ticks -> {{{μ - 2σ, HoldForm[μ - 2σ]}, {μ, HoldForm[μ]}, {μ + 2σ, HoldForm[μ + 2σ]}}, Automatic}, PlotLabel -> k], {k, {5, 10, 25, 100}}]]
ExponentialDistributionに従っている k 個の変数の総和はアーラン分布に従う:
CharacteristicFunction[ExponentialDistribution[λ], t] ^ kCharacteristicFunction[ErlangDistribution[k, λ], t]% - %% //FullSimplify[#, Element[k, Integers]]&TransformedDistribution[u + v + w, {u, v, w}ProductDistribution[{ExponentialDistribution[λ], 3}]]アーラン分布はタイプ3のピアソン分布(PearsonDistribution)の特殊ケースである:
PDF[PearsonDistribution[3, 1, (-k + 1/λ), 0, (1/λ), 0], x]//Refine[#, λ > 0]&PDF[ErlangDistribution[k, λ], x]//Refine[#, λ > 0]&% - %%//FullSimplify[#, λ > 0]&アーラン分布はGammaDistributionの特殊ケースである:
PDF[ErlangDistribution[k, λ], x]PDF[GammaDistribution[k, 1 / λ], x]//PowerExpand% - %%//SimplifyParetoDistributionはExponentialDistributionとErlangDistributionの商から求められる:
𝒟 = TransformedDistribution[u / v, {uExponentialDistribution[λ], vErlangDistribution[n, λ]}];PDF[𝒟, x]PDF[ParetoDistribution[1, n, 0], x]FullSimplify[% - %%, x ≠ 0]ParetoDistributionはErlangDistributionの商として求めることができる:
TransformedDistribution[u / v, {uErlangDistribution[1, a], vErlangDistribution[1, b]}]おもしろい例題 (1)
累積分布関数の等高線を持つ λ のさまざまな値についての確率密度関数:
cdf = Function[{x, λ}, Evaluate[CDF[ErlangDistribution[3, λ], x]]];
ql = {0.025, 0.10, 0.25, 0.5, 0.75, 0.90, 0.975};
cl = Table[ColorData["Rainbow"][q], {q, Join[{0.0}, ql]}];Legended[Plot3D[PDF[ErlangDistribution[3, λ], x], {x, 0, 4}, {λ, 0.2, 4}, PlotTheme -> "Marketing", MeshFunctions -> {cdf}, Mesh -> {ql}, MeshStyle -> GrayLevel[0.8], MeshShading -> cl, AxesLabel -> Automatic, PlotPoints -> 100, BaseStyle -> Opacity[0.9], ImageSize -> 400, PlotRange -> {0, 1.2}], BarLegend["Rainbow", ql, LegendLabel -> "prob"]]関連するガイド
-
▪
- 信頼性解析の分布 ▪
- 待ち行列過程 ▪
- 指数関連分布 ▪
- 通信システムの分布 ▪
- パラメトリック統計分布
テキスト
Wolfram Research (2010), ErlangDistribution, Wolfram言語関数, https://reference.wolfram.com/language/ref/ErlangDistribution.html (2016年に更新).
CMS
Wolfram Language. 2010. "ErlangDistribution." Wolfram Language & System Documentation Center. Wolfram Research. Last Modified 2016. https://reference.wolfram.com/language/ref/ErlangDistribution.html.
APA
Wolfram Language. (2010). ErlangDistribution. Wolfram Language & System Documentation Center. Retrieved from https://reference.wolfram.com/language/ref/ErlangDistribution.html
BibTeX
@misc{reference.wolfram_2026_erlangdistribution, author="Wolfram Research", title="{ErlangDistribution}", year="2016", howpublished="\url{https://reference.wolfram.com/language/ref/ErlangDistribution.html}", note=[Accessed: 11-August-2026]}
BibLaTeX
@online{reference.wolfram_2026_erlangdistribution, organization={Wolfram Research}, title={ErlangDistribution}, year={2016}, url={https://reference.wolfram.com/language/ref/ErlangDistribution.html}, note=[Accessed: 11-August-2026]}