Quand les approximations de la série Taylor aux attentes des fonctions (entières) convergent-elles?


10

Prenez une espérance de la forme pour une variable aléatoire univariée et une fonction entière (c.-à-d., L'intervalle de convergence est la ligne réelle entière)E(f(X))Xf()

J'ai une fonction de génération de moment pour et je peux donc facilement calculer des moments entiers. Utilisez une série de Taylor autour de puis appliquez l'attente en termes d'une série de moments centraux, = f (\ mu) + \ sum_ {n = 2 } ^ {\ infty} \ frac {f ^ {(n)} (\ mu)} {n!} E \ left [(x - \ mu) ^ n \ right] Tronque cette série, E_N (f (x) ) = f (\ mu) + \ sum_ {n = 2} ^ {N} \ frac {f ^ {(n)} (\ mu)} {n!} E \ left [(x - \ mu) ^ n \droite] XμE(x)

E(f(x))=E(f(μ)+f(μ)(xμ)+f(μ)(xμ)22!+)
=f(μ)+n=2f(n)(μ)n!E[(xμ)n]
EN(f(x))=f(μ)+n=2Nf(n)(μ)n!E[(xμ)n]

Ma question est la suivante: dans quelles conditions sur la variable aléatoire (et tout élément supplémentaire sur f() également) l'approximation de l'attente converge-t-elle lorsque j'ajoute des termes (c'est-à-dire limNEN(f(x))=E(f(x)) ).

Puisqu'il ne semble pas converger pour mon cas (une variable aléatoire de poisson et f(x)=xα ), existe-t-il d'autres astuces pour trouver des attentes approximatives avec des moments entiers lorsque ces conditions échouent?



@Jonathan Merci. Voir mes modifications maintenant que c'est devenu plus clair. Très utile, même si je ne pouvais pas vraiment le casser. De cela, il apparaît qu'une condition suffisante pour que cela fonctionne est que ma variable aléatoire soit fortement concentrée? Bien que j'aie du mal à comprendre exactement comment utiliser l'inégalité de Hoeffding, etc. pour comparer à ces notes.
jlperla

Que voulez-vous dire "une variable aléatoire de poisson et "? Est-ce un cas ou deux, et quel est le pdf? f(x)=xα
Carl

@Carl Il y a quelques années, mais si je me souviens bien, la variable était pour certains avec PDF de en.wikipedia.org/wiki/Poisson_distribution . Que était la fonction sur laquelle je répondais aux attentes. iexPoisson(λ)λf(x)E(f(x))
jlperla

Je ne sais pas ce que vous demandez. Que diriez-vous que les moments supérieurs de la distribution de Poisson autour de l'origine sont des polynômes de Touchard dans : où les {accolades} désignent des nombres de Stirling du second type? mkλ
mk=i=0kλi{ki},
Carl

Réponses:


1

Par votre hypothèse que est analytique réel, Converge presque sûrement (en fait sûrement) vers .f

yn=f(μ)+f(μ)(xμ)+f(μ)(xμ)22!++f(n)(μ)(xμ)nn!
f(x)

Une condition standard dans laquelle, comme la convergence implique la convergence des attentes, c'est-à-dire est que comme pour certains tels que . (Théorème de convergence dominé.)

E[f(x)]=E[limnyn]=limnE[yn],
|yn|yyE[y]<

Cette condition serait remplie si la série de puissances converge absolument comme, c'est-à-dire et

y=n0|f(n)(μ)||xμ|nn!<a.s.
E[y]<.

Votre exemple d'une variable aléatoire de Poisson et , , suggérerait que l'intégrabilité ci-dessus du critère de limite absolue est la plus faible possible, en général.f(x)=xααZ+


-1

L'approximation convergera si la fonction f (x) admet une expansion en série de puissance, c'est-à-dire que toutes les dérivées existent. Il sera également pleinement atteint si les dérivés d'un seuil spécifique et supérieur sont égaux à zéro. Vous pouvez vous référer à Populis [3-4] et à Stark et Woods [4].


"Il sera également pleinement atteint si les dérivés d'un seuil spécifique et supérieur sont égaux à zéro." Si les dérivées existent et sont égales à zéro, n'est-ce pas une autre façon de dire polynôme?
Accumulation

Ce n'est pas vrai. Lorsque "tous les dérivés existent" au moment de l'expansion de la série de puissance, la série de puissance n'a pas besoin de converger n'importe où. (L'exemple standard est la série Maclaurin de ) Un autre est que même lorsque la série converge à un moment donné, elle n'a pas besoin de converger partout. Un exemple simple est la série Maclaurin deLorsque cela se produit, la convergence dépend des détails de la variable aléatoire. Par exemple, supposons que ait une distribution de Student t et considéronsFinalement, n'existe même pas! e1/x2.1/(1x).X
1/(1X)=1+X+X2++Xn+.
E(Xn)
whuber
En utilisant notre site, vous reconnaissez avoir lu et compris notre politique liée aux cookies et notre politique de confidentialité.
Licensed under cc by-sa 3.0 with attribution required.