モーメント生成関数を特性関数より優先する場合は?


9

ましょう確率空間であり、およびletランダムベクトルです。LETの分布である、上のボレル測度。(Ω、F、P)バツ:Ω→RんPバツ=バツ∗PバツRん

  • 特性関数の関数で に対して定義された(確率変数は、すべてのに対してで制限されます)。これはのフーリエ変換です。バツ
    φバツ(t)=E[e私t⋅バツ]=∫Ωe私t⋅バツdP、
    t∈Rんe私t⋅バツL1(P)tPバツ
  • のモーメント生成関数(mgf)は、関数 すべての、上記積分が存在します。これはのラプラス変換です。バツ
    Mバツ(t)=E[et⋅バツ]=∫Ωet⋅バツdP、
    t∈Rん Pバツ

すでに、特性関数がRんどこにでも定義されていることがわかりますが、mgfにはバツに依存するドメインがあり、このドメインは{0}可能性があります(これは、例えば、コーシー分布の確率変数の場合に起こります) )。

これにも関わらず、特徴的な関数とmgfは多くの特性を共有しています。例えば:

  1. もしバツ1、…、バツん独立しており、次いで、
    φバツ1+⋯+バツん(t)=φバツ1(t)⋯φバツん(t)
    すべてについてt、および
    Mバツ1+⋯+バツん(t)=Mバツ1(t)⋯Mバツん(t)
    mgfが存在する すべてのt n(t)。
  2. 2つのランダムベクトルバツとYと同じ分布を有する場合にのみφバツ(t)=φY(t)すべてのためt。この結果のmgfアナログは、0のある近傍のすべてのtについてMバツ(t)=MY(t)場合、XとYは同じ分布を持つということです。t 0バツY
  3. 一般的な分布の特性関数とmgfの形式はよく似ています。例えば、もしバツ〜Nん(μ、Σ)(ん平均値と次元正常μと共分散行列Σ)は、
    φバツ(t)=exp⁡(私μ⋅t−12t⋅(Σt))
    と
    Mバツ(t)=exp⁡(μ⋅t−12t⋅(Σt))。
  4. いくつかの穏やかな仮定が成り立つ場合、特性関数とmgfの両方を微分してモーメントを計算できます。
  5. レヴィの連続性定理は、対応する特性関数の収束を使用して、確率変数のシーケンスが別の確率変数に分布で収束する時期を決定する基準を与えます。mgfの対応する定理があります(Curtiss 1942、定理3)。

特性関数とmgfは同じ目的で使用されることが多く、特性関数は常に存在するがmgfは常に存在しないという事実を考えると、mgfよりも特性関数を使用することを好むはずです。

質問。

  1. 特性関数よりもmgfの方が便利な例は何ですか?
  2. 特性関数ではできないmgfで何ができるのでしょうか?

1
この質問の鍵は、終わりに近い「紹介」という言葉ではありませんか?基本的な微積分への最小限の露出のみを想定し(そして快適ではない)、多くの場合それさえ想定しないコースに複素数の分析を含む何かを導入することは教育的に意味がありますか?
— whuber

1
@whuberそれも私が考えたものですが、教育についての質問をしたくないので、おそらく最後の段落を削除する必要があります
— Artem Mavrin

回答:


3

それは良い質問ですが、幅広い質問なので、言っておくべきことについてすべて言うことは約束できません。簡単に言えば、ライバルのテクニックは、何ができるかではなく、どれほどきちんとできるかという点で異なります。

特性関数は複素数の役割があるため、特に注意が必要です。学生が複素数について知っている必要があるということさえありません。関係する微積分には微妙な落とし穴があるということです。たとえば、変数シフトの置換で二乗を完了するだけで正規分布のMGFを取得できますが、多くのソースは、特性関数を使用するアプローチが不用意にふりをするのと同じくらい簡単です。ガウス積分の有名な正常化は、上の統合については何も述べていないので、それは、ない私c+Rとc∈R∖{0}。ああ、コンターに注意すれば、積分を評価することができます。実際には、もっと簡単な方法があります。これは、N(0、1)分布の特性関数φ(t)を満たしますφ˙=−tφ。しかし、MGFアプローチはさらに単純であり、学生が早い段階で必要とする分布のほとんどは、線分(例:ラプラス)または半線(例:ガンマ、幾何学的、負の二項)、またはR全体(例:ベータ、二項、ポアソン、ノーマル)。いずれにせよ、それは瞬間を研究するには十分です。

MGFでしかできないことはないと思いますが、目の前のタスクで最も簡単な方法を使用します。ポアソン分布のモーメントを計算する最も簡単な方法は何ですか?私はそれが再び異なる技術を使用することだと主張したい、確率生成関数G(t)=Etバツ=exp⁡λ(t−1)。次に、下降ポッホハンマー記号(バツ)kはE(バツ)k=G(k)(1)=λk。一般に、離散分布の場合はPGFを使用し、PDFの末尾に境界があるか超指数関数的減衰を伴う連続分布の場合はMGFを使用し、実際に必要な場合は特性関数を使用する価値があります。

そして、あなたが求めている質問によっては、代わりに、キュムラント生成関数を使用することが賢明であると考えるかもしれません。それは、MGFまたはCFの対数として定義されている場合です。例えば、私は最高のためのキュムラントのログMGFの定義という運動としてそれを残しておきますん Exp⁡(1)のIID与えκメートル=(メートル−1)!Σk=1んk−メートル、(それぞれ、平均と分散の非常に簡単な計算を提供κ1及びκ2)あなたは瞬間の点でそれらを書いた場合に比べ。


2
cfはRの複素数値関数の積分として定義されているため、「統合」についてのあなたの意見は理解できません。 輪郭積分と見なす必要はありません。複素数で不快な人にとっては、とにかく実数積分のペアと見なすことができます。mgfがどの点でも「単純」であるかは不明です。実際、cfは、収束について心配する必要がないという意味でより単純です。ic+R,R.
— whuber

1
私が何を意味するか@whuberです。∫R12πexp⁡(−x22+itx)dx=∫−it+R12πexp⁡(−y22−t22)dt
— JG、

私も同じくらい疑った。しかし、それはcf自体の固有の機能ではなく、積分を評価することを選択する方法の単なるアーティファクトではありませんか?
— whuber

@whuber問題は、多くのソースが、MGFの場合と同じくらい簡単に置換が機能するふりをすることですが、そうではありません。
— JG、

1
なぜそれができないのかについて少し詳しく教えていただけませんか?私はこの特定のケースで何も問題がないと思います。そして一般的に、上の元の積分は収束するので、この種の置換で問題が発生することはありません。R
— whuber

3

確率変数にすべてのモーメントがある場合、MGFが存在し、一般に少なくとも証明の特性関数と同じくらい役立ちます。

あなたの質問に答えるために、MGFが存在する場合、それはX関連する多くの極値計算の基礎を提供します。(用である最も単純なt≥0)、

P(X>r)=P(etX>etr)≤MX(t)/etr.

ここで、rhsはtにわたって最小化できます。奇妙なことに、この限界は、まれなイベントの見積もりを取得するために私たちが知っている数少ない簡単な方法の1つです。これの一般的な領域は、大きな偏差理論です。そこでは、より良い(よりタイトな)境界を取得するために大量の作業を行う必要があります。これの一般的な例は、Sn=X1+⋯+Xnを見ているため、X1のMGFが存在する場合、P (| S n − E [ X ] | > n r )を示すことができます。P(|Sn−E[X]|>nr)n指数関数的に減衰します。これは、より一般的にはCramerの定理として知られています。

これについての簡単なメモをいくつか示します。


1
最初の段落のすべては、私が間違っていると思う最後の文を除いて、質問ですでに言及されています。たとえば、対数正規分布のすべてのモーメントが存在しますが、そのmgfは正の実数に対して未定義です。回答の2番目の部分は、特徴的な機能のアナログがないと思われるmgfのアプリケーションを強調しているため、非常に役立ちます
— Artem Mavrin
弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.