Attente d'une fonction d'une variable aléatoire de CDF

9

Est-il possible de calculer l'espérance d'une fonction d'une variable aléatoire avec seulement le CDF du VR? Dis que j'ai une fonctiong(x) qui a la propriété g(x)dx et la seule information que j'ai sur la variable aléatoire est le CDF.

Par exemple, j'ai un scénario où il y a trois temporisateurs qui peuvent être modélisés comme des variables aléatoires exponentielles X1,X2,X3 avec paramètres de débit λ1,λ2,λ3respectivement. Pour chaque instant, je gagne une récompense en fonction d'une fonction de récompenseg(x). Autrement dit, ma récompense pour attendre le tempst peut être écrit comme 0tg(x)dx. cependant,g(x) connaît des rendements décroissants de sorte que la récompense marginale reçue d'attendre une seconde à t=0 est supérieur à une seconde par exemple t=27. Ce «jeu» se termine lorsque l'une des deux choses se produit. Soit les deux minuteriesX1 ou X2 doit sonner ou minuteries X1 ou X3doit sonner. J'essaie de trouver la récompense attendue de jouer à ce jeu.

Actuellement, je peux calculer le CDF de la variable aléatoire en modélisant le temps jusqu'à la fin du jeu, mais je ne sais pas comment utiliser ces informations lorsque ce dont j'ai vraiment besoin est une récompense associée à ce temps.

Jusqu'à présent, j'ai les variables aléatoires supplémentaires:

W12=max(X1,X2)W13=max(X1,X3)Z=min(W12,W13)
Soit également désigne le CDF de Le CDF de , peut s'écrire: Fi(x),i{1,2,3}XiZ
FZ(t)=F1(t)F2(t)+F1(t)F3(t)F1(t)F2(t)F3(t)

Je sais quand une variable aléatoire prend des valeurs non négatives, vous pouvez utiliser un raccourci pour calculer l'attente à l'aide du CDF. Autrement dit, . Y a-t-il quelque chose de similaire que je pourrais utiliser pour une fonction d'une variable aléatoire, ou est-il nécessaire de calculer d'abord le pdf de pour calculerE[X]=0F(Xx)dxZ0g(t)fz(t)dx

CoconutBandit
la source
2
Qu'entendez-vous par "uniquement" des informations? Le CDF vous dit tout sur le VR qui pourrait être lié aux attentes! Il semble que votre problème sous-jacent soit lié à la forme de calcul sous laquelle le CDF vous est remis. Veuillez expliquer votre situation. BTW, peut être indéfini ou infini même lorsque l'intégrale deest fini. E[g(X)]|g|
whuber
2
Je pense que vous recherchez une intégration par parties en.wikipedia.org/wiki/Integration_by_parts
seanv507

Réponses:

13

Lorsque est le CDF d'une variable aléatoire et est une fonction (mesurable), l'espérance de peut être trouvée comme une intégrale de Riemann-StieltjesFXgg(X)

E(g(X))=g(x)dF(x).

Cela exprime la loi du statisticien inconscient.

Si est également différentiable, écrivez et intégrez par parties pour donnergdF=d(1F)

E(g(X))=g(x)(1F(x))|+(1F(x))g(x)dx

à condition que les deux addends convergent. Cela signifie plusieurs choses, qui peuvent être simplement exprimées en brisant l'intégrale à une valeur finie définie telle que :0

  1. limxg(x)(1F(x)) et existent et sont finies. Si c'est le cas, le premier addend est la différence de ces deux.limxg(x)(1F(x))

  2. limtt0(1F(x))g(x)dx et existent et sont finis. Si c'est le cas, le deuxième ajout est la somme de ces deux.limt0t(1F(x))g(x)dx

Un bon endroit pour rompre l'intégrale est à n'importe quel zéro de , car - à condition que diminue finalement assez rapidement pour les grands--qui fait disparaître le premier addend, ne laissant que l'intégrale de contre la fonction de survie .gg|x|g1F

Exemple

L'espérance d'une variable non négative est obtenue en appliquant la formule à la fonction d'identité pour laquelle et en utilisant le fait que l'intégration peut commencer à zéro:Xg(x)=xg(x)=1

E(X)=x(1F(x))|0+0(1F(x))dx.

À condition que (c'est-à-dire que la fonction de survie n'a pas de queue trop lourde), la limite supérieure du premier terme disparaît. Sa limite inférieure disparaît évidemment. Il ne nous reste que l'intégrale, donnant l'expression dans la question.limxx(1F(x))=0

whuber
la source
Merci, cela ressemble exactement à ce que je voulais. J'ai juste besoin de lire sur mon intégration Riemann-Stieltjes maintenant.
CoconutBandit
Dans votre application, parce que est continuellement différenciable partout sauf à , vous pouvez diviser l'intégrale à en deux intégrales Riemann et ignorer complètement les complications. F00
whuber
Qu'entendez-vous par «complications»? De plus, dans votre deuxième point, devrait-il être ? Sinon, pourquoi le devenu ? t0(1F(x))g(x)dxt0(1F(x))g(x)dxg(x)g(x)
CoconutBandit
(1) Merci, ces nombres premiers devaient être là. (2) "Complications" fait référence à la nécessité de l'intégrale de Riemann-Stieltjes au lieu de l'intégrale de Riemann.
whuber
1
Celui-ci utilise trois règles de différenciation de base: la règle de somme, la règle de produit et le fait que les constantes n'ont aucune dérivée.
whuber