LevyDistribution[μ,σ]
位置母数 μ,分散パラメータ σ のレヴィ(Lévy)分布を表す.
LevyDistribution
LevyDistribution[μ,σ]
位置母数 μ,分散パラメータ σ のレヴィ(Lévy)分布を表す.
詳細
- レヴィ分布における値
の確率密度は
に比例する. » - レヴィ分布LevyDistribution[0,σ]は逆ガンマ関数の特別ケース(
,
)である. » - LevyDistributionでは,μ は任意の実数で σ は任意の正の実数である.
- LevyDistributionでは,μ と σ は単位次元が等しい任意の数量でよい. »
- LevyDistributionはMean,CDF,RandomVariate等の関数とともに使うことができる.
予備知識
- LevyDistribution[μ,σ]は,区間
でサポートされ,実数 μ(「位置母数」と呼ばれる)および正の実数 σ(「分散母数」と呼ばれる)でパラメータ化された統計分布を表す.レヴィ分布の確率密度関数(PDF)は,大体において単一の「峰」(最大値)を持つ単峰性であるが,その全体的な形状(高さ,
の線近くでの集中度,最大値の水平位置)は μ と σ の値で決定される.これに加え,PDFの裾部は,PDFが
の大きい値について指数的というよりむしろ代数的に現象するという意味で「太い」(この現象は分布のSurvivalFunctionを分析することで数量的に正確にできる).レヴィ分布は,レヴィ分布に従う独立した確率変数の任意の線形結合が再びレヴィ分布に従うという意味でいわゆる安定分布(StableDistribution)であり,安定分布のある種の部分集合を説明するために使われる,より一般的な用語の「パレート・レヴィ分布」および「レヴィ(α)安定分布」と混同してはならない. - レヴィ分布はフランス人の数学者Paul Lévyに因んで名付けられた.レヴィ分布の歴史は(1920年代にレヴィによってはじめて詳しく説明された)他の安定分布のより一般的なクラスと切り離すことはできないが,簡約されたレヴィ分布のPDFについては,1919年にオランダ人の天文学者Holtsmarkによって重力場のランダム振動についての叙述の一部としてなはじめて書かれている.安定分布族は,当初は現実世界への応用が欠けているとして数学者や確率論の研究者に無視されたが,ブラウン運動(いわゆるレヴィフライト)を含むさまざまな概念の分析に用いられ,現在では金融工学や数理ファイナンスの分野の重要なツールとなっている.さらに,レヴィ分布は,地球磁気学,素粒子物理学,暗号学,信号解析,生物学のような分野における数多くの重要な現象の説明に使われている.
- RandomVariateを使って,レヴィ分布から,1つあるいは複数の機械精度あるいは任意精度(後者はWorkingPrecisionオプションを介す)の擬似乱数変量を得ることができる.Distributed[x,LevyDistribution[μ,σ]](より簡略な表記では xLevyDistribution[μ,σ])を使って,確率変数 x がレヴィ分布に従って分布していると宣言することができる.このような宣言は,Probability,NProbability,Expectation,NExpectation等の関数で使うことができる.
- レヴィ分布の確率密度関数および累積分布関数は,PDF[LevyDistribution[μ,σ],x]およびCDF[LevyDistribution[μ,σ],x]を使って得られる.平均,中央値,分散,原点の周りのモーメント,中心モーメントは,それぞれMean,Median,Variance,Moment,CentralMomentを使って計算することができる.しかし,LevyDistributionは裾部が長いため,
番目の原点の周りのモーメントと中心モーメント(平均と分散を含む)はすべての
について
であり,そのモーメント母関数(MomentGeneratingFunction)はIndeterminateである. - DistributionFitTestを使って,与えられたデータ集合がレヴィ分布と一致するかどうかを検定することが,EstimatedDistributionを使って与えられたデータからパラメトリックレヴィ分布を推定することが,FindDistributionParametersを使ってデータをレヴィ分布にフィットすることができる.ProbabilityPlotを使って記号レヴィ分布のCDFに対する与えられたデータのCDFのプロットを生成することが,QuantilePlotを使って記号レヴィ分布の変位値に対する与えられたデータの変位値のプロットを生成することができる.
- TransformedDistributionを使って変換されたレヴィ分布を表すことが,CensoredDistributionを使って上限値と下限値の間で切り取られた値の分布を表すことが,TruncatedDistributionを使って上限値と下限値の間で切断された値の分布を表すことができる.CopulaDistributionを使ってレヴィ分布を含む高次元分布を構築することが,ProductDistributionを使ってレヴィ分布を含む独立成分分布の結合分布を計算することができる.
- LevyDistributionは他の多くの分布と密接な関係がある.例えば,LevyDistributionはStableDistributionの一つの例である(LevyDistribution[μ,σ]は厳密にStableDistribution[0,1/2,1,μ+σ,σ]である).したがって,CauchyDistribution,LandauDistribution,NormalDistribution等の他の安定分布に関連している.LevyDistributionは,InverseGammaDistributionとPearsonDistribution両方の特殊ケースであり,(TransformedDistributionを介して)NormalDistributionの変換として得ることができる.LevyDistributionは,HalfNormalDistribution,ChiSquareDistribution,GammaDistribution,InverseGaussianDistributionとも関連している.
例題
すべて開く すべて閉じる例 (4)
Plot[Table[PDF[LevyDistribution[1, σ], x], {σ, {1 / 2, 1, 2}}]//Evaluate, {x, 1.1, 3}, Filling -> Axis, PlotRange -> All]Plot[Table[PDF[LevyDistribution[μ, 1.5], x], {μ, {-1 / 2, 1, 2}}]//Evaluate, {x, -1, 5}, Filling -> Axis, PlotRange -> All]PDF[LevyDistribution[μ, σ], x]Plot[Table[CDF[LevyDistribution[1, σ], x], {σ, {1 / 2, 1, 2}}]//Evaluate, {x, 1.1, 3}, Filling -> Axis, PlotRange -> All]Plot[Table[CDF[LevyDistribution[μ, 1.5], x], {μ, {-1 / 2, 1, 2}}]//Evaluate, {x, -1, 5}, Filling -> Axis, PlotRange -> All]CDF[LevyDistribution[μ, σ], x]Mean[LevyDistribution[μ, σ]]Variance[LevyDistribution[μ, σ]]Median[LevyDistribution[μ, σ]]スコープ (6)
data = RandomVariate[LevyDistribution[1, 2], 10 ^ 4];Show[
Histogram[data, {1, 10, 0.5}, "PDF"],
Plot[PDF[TruncatedDistribution[{1, 10}, LevyDistribution[1, 2]], x], {x, 1.1, 10}, PlotStyle -> Thick, PlotRange -> All]]sample = RandomVariate[LevyDistribution[1, 3], 10 ^ 3];edist = EstimatedDistribution[sample, LevyDistribution[μ, σ]]サンプルの密度ヒストグラムを推定分布の確率密度関数と比較する:
Show[Histogram[sample, {1, 20, 1}, "PDF"], Plot[PDF[TruncatedDistribution[{1, 20}, edist], x], {x, 1, 20}, PlotStyle -> Thick]]Moment[LevyDistribution[μ, σ], r]Expectation[x ^ r, xLevyDistribution[0, 1], Assumptions -> r < 1 / 2]Plot[Table[HazardFunction[LevyDistribution[μ, 2], x], {μ, {-1, 2, 3}}]//Evaluate, {x, -2, 10}, Filling -> Axis]Plot[Table[HazardFunction[LevyDistribution[2, σ], x], {σ, {1, 2, 3}}]//Evaluate, {x, 2.1, 6}, Filling -> Axis]HazardFunction[LevyDistribution[μ, σ], x]Plot[Table[Quantile[LevyDistribution[1, σ], q], {σ, {0.5, 1, 2}}]//Evaluate, {q, 0, 1}, Filling -> Axis]Quantile[LevyDistribution[μ, σ], q]母数でQuantityを一貫して使うとQuantityDistributionが与えられる:
geomMagneticReversalTime𝒟 = LevyDistribution[["0.01 millenia"], ["200 millenia"]]Quartiles[geomMagneticReversalTime𝒟]アプリケーション (1)
ファンデルワールス(van der Waals)スペクトル断片の半値幅を求める:
profile[ω_] = PDF[LevyDistribution[0, 1], ω]Plot[profile[ω], {ω, 0.1, 5}, Filling -> Axis]ω0 = ArgMax[profile[ω], ω]ω /. NSolve[profile[ω] == profile[ω0] / 2, ω, Reals]Differences[%]NProbability[ω > ω0, ωLevyDistribution[0, 1]]特性と関係 (7)
レヴィ分布は平行移動と正の因子によるスケーリングの下では閉じている:
TransformedDistribution[a * u + b, uLevyDistribution[μ, σ]]TransformedDistribution[x + y, {xLevyDistribution[Subscript[μ, 1], Subscript[σ, 1]], yLevyDistribution[Subscript[μ, 2], Subscript[σ, 2]]}]
LevyDistribution[0,σ]はInverseGammaDistributionの特殊ケースである:
PDF[LevyDistribution[0, σ], x]PDF[InverseGammaDistribution[1 / 2, σ / 2], x]Simplify[% - %%]レヴィ分布はタイプ5のPearsonDistributionの特殊ケースである:
PDF[PearsonDistribution[5, 1, -μ - (σ/3), (2/3), -(4/3)μ, (2/3)μ^2], x]//FullSimplifyPDF[LevyDistribution[μ, σ], x]//PowerExpandSimplify[% - %%]レヴィ分布はNormalDistributionを変換したものである:
TransformedDistribution[u ^ -2, uNormalDistribution[]]TransformedDistribution[(u - μ) ^ -2, uNormalDistribution[μ, σ]]レヴィ分布はStableDistributionである:
PDF[LevyDistribution[μ, σ], x]PDF[StableDistribution[0, 1 / 2, 1, μ + σ, σ], x]% - %%//FullSimplify考えられる問題 (2)
LevyDistributionは μ が実数ではない場合は定義されない:
Mean[LevyDistribution[I, 5]]LevyDistributionは σ が正の実数でなければ定義されない:
Mean[LevyDistribution[0, -5]]記号出力に有効ではないパラメータを代入すると無意味な結果になる:
CDF[LevyDistribution[μ, σ], 1] /. {μ -> 0, σ -> 1 + I}N[%]おもしろい例題 (1)
累積分布関数の等高線を持つ μ のさまざまな値についての確率密度関数:
dist = LevyDistribution[μ, 1.5];cdf = Function[{x, μ}, Evaluate[CDF[dist, 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[dist, x], {x, -3.5, 20}, {μ, -3, 3}, PlotTheme -> "Marketing", MeshFunctions -> {cdf}, Mesh -> {ql}, MeshStyle -> GrayLevel[0.8], MeshShading -> cl, AxesLabel -> Automatic, BaseStyle -> Opacity[0.9], ImageSize -> 400, PlotPoints -> 50, PlotRange -> {0, 0.4}], BarLegend["Rainbow", ql, LegendLabel -> "prob"]]関連するガイド
-
▪
- 裾の重い分布 ▪
- パラメトリック統計分布
テキスト
Wolfram Research (2008), LevyDistribution, Wolfram言語関数, https://reference.wolfram.com/language/ref/LevyDistribution.html (2016年に更新).
CMS
Wolfram Language. 2008. "LevyDistribution." Wolfram Language & System Documentation Center. Wolfram Research. Last Modified 2016. https://reference.wolfram.com/language/ref/LevyDistribution.html.
APA
Wolfram Language. (2008). LevyDistribution. Wolfram Language & System Documentation Center. Retrieved from https://reference.wolfram.com/language/ref/LevyDistribution.html
BibTeX
@misc{reference.wolfram_2026_levydistribution, author="Wolfram Research", title="{LevyDistribution}", year="2016", howpublished="\url{https://reference.wolfram.com/language/ref/LevyDistribution.html}", note=[Accessed: 02-November-2026]}
BibLaTeX
@online{reference.wolfram_2026_levydistribution, organization={Wolfram Research}, title={LevyDistribution}, year={2016}, url={https://reference.wolfram.com/language/ref/LevyDistribution.html}, note=[Accessed: 02-November-2026]}