Explication intuitive de la densité de la variable transformée?

37

Supposons que soit une variable aléatoire avec pdf f X ( x ) . Alors la variable aléatoire Y = X 2 a le pdfXfX(x)Y=X2

fY(y)={12y(fX(y)+fX(y))y00y<0

Je comprends le calcul derrière cela. Mais j'essaie de trouver un moyen de l'expliquer à quelqu'un qui ne connaît pas le calcul. En particulier, j'essaie d'expliquer pourquoi le facteur 1y apparaît devant. Je vais tenter le coup:

Supposons que X ait une distribution gaussienne. Presque tout le poids de son pdf est entre les valeurs, par exemple, 3 et 3. Mais que les cartes à 0 à 9 pour Y . Ainsi, le poids lourd dans le pdf pour X a été étendue à travers une gamme plus large de valeurs dans la transformation à Y . Ainsi, pour que fY(y) soit un véritable pdf, le poids très lourd doit être réduit du facteur multiplicatif 1y

Comment ça sonne?

Si quelqu'un pouvait fournir une meilleure explication ou un lien vers celui d'un document ou d'un manuel, je l'apprécierais beaucoup. Je trouve cet exemple de transformation variable dans plusieurs livres de statistiques mathématiques d'introduction et de probabilité. Mais je ne trouve jamais d'explication intuitive avec cela :(

lowndrul
la source
Je pense que votre explication est correcte.
highBandWidth
2
L'explication est juste, mais elle est purement qualitative: la forme précise du facteur multiplicatif reste un mystère. Le pouvoir -1/2 apparaît simplement comme par magie. Ainsi, à un certain niveau, vous devez faire la même chose que Calcul: trouver le taux de changement de la fonction racine carrée.
whuber

Réponses:

37

Les PDF sont des hauteurs mais ils sont utilisés pour représenter une probabilité par surface. Il est donc utile d’exprimer un fichier PDF de manière à nous rappeler que la surface est égale à la hauteur multipliée par la base.

Initialement, la hauteur pour toute valeur x est donnée par le PDF fX(x) . La base est le segment infinitésimal dx , d'où la distribution (c'est-à-dire la mesure de probabilité opposée à la fonction de distribution ) est en réalité la forme différentielle, ou "élément de probabilité".

PEX(x)=fX(x)dx.

Ceci, plutôt que le PDF, est l’objet avec lequel vous voulez travailler, tant sur le plan conceptuel que pratique, car il inclut explicitement tous les éléments. les éléments nécessaires pour exprimer une probabilité.

Lorsque nous ré-exprimons x en termes de y=x2 , les segments de base dx sont étirés (ou pincés): en quadrillant les deux extrémités de l'intervalle de x à x+dx nous voyons que la base de la zone y doit être un intervalle de longueur

dy=(x+dx)2x2=2xdx+(dx)2.

Puisque le produit de deux infinitésimaux est négligeable par rapport aux infinitésimaux eux-mêmes, nous concluons

dy=2xdx, whence dx=dy2x=dy2y.

Ceci étant établi, le calcul est trivial car il suffit de brancher la nouvelle hauteur et la nouvelle largeur:

PEX(x)=fX(x)dx=fX(y)dy2y=PEY(y).

Parce que la base, en termes de y , est dy , quel que soit le multiplie doit être la hauteur, que l' on peut lire directement sur le moyen terme comme

12yfX(y)=fY(y).

Cette équation PEX(x)=PEY(y) est effectivement une loi de conservation de la surface (= probabilité).

Two pdfs

Ce graphique montre avec précision les parties étroites (presque infinitésimales) de deux PDF liés par y=x2 . Les probabilités sont représentées par les zones ombrées. En raison de la compression de l'intervalle [0.32,0.45] via la quadrature, la hauteur de la région rouge ( y , à gauche) doit être proportionnellement étendue pour correspondre à la surface de la région bleue ( x , à droite).

whuber
la source
2
J'aime les infinitésimaux. C'est une merveilleuse explication. Penser en termes de , ce qui peut être vu clairement émerger de la dérivée de la transformation, est beaucoup plus intuitif que de penser en 2x . Je pense que c’était là où était mon problème. y
lowndrul
@whuber, je crois que vous la première ligne devrait être ? Est-ce ce que vous entendez par pdf X ( x ) ? PS: aussi curieux de savoir ce que vous pensez de ma réponse (ci-dessous). P(X(x,x+dx))=fx(x)dxpdfX(x)
Carlos Cinelli
@Carlos It's a little more rigorous to express the idea in the way I did at the outset: the PDF is what you multiply the Lebesgue measure dx by in order to get the given probability measure.
whuber
@whuber but if the pdf is what you multiply then it is the term fX(x), not the product fx(x)dx as you wrote, right? It is not clear why you call the product fX(x)dx a pdf.
Carlos Cinelli
1
@Carlos: thank you; now I see your point. I made some edits to address it.
whuber
11

How about, if I manufacture objects that are always square and I know the distribution of the side lengths of the squares; what can I say about the distribution of the areas of the squares?

In particular, if I know the distribution of a random variable X, what can I say about Y=X2? One thing that you can say is

FY(c)=P(Yc)=P(X2c)=P(cXc)=FX(c)FX(c).

So a relationship is established between the CDF of Y and CDF of X; what is the relationship between their PDFs? We need calculus for that. Taking the derivatives of both sides gives you the results you wanted.

schenectady
la source
2
(+1) Although this is not a full answer, it presents a good way to go about finding fY and clearly shows why it is a sum of two pieces, one for each square root.
whuber
1
I don't get why pdf(x) = f(x)dx. What about pdf(x) dx = f(x), density = prob mass/interval...what i'm getting wrong?
Fernando
2

Imagine we have a population and Y is a summary of that population. Then P(Y(y,y+Δy)) is counting the proportion of individuals that have variable Y in the range (y,y+Δy). You can consider this as a "bin" of size Δy and we are counting how many individuals are inside that bin.

Now let us re-express those individuals in terms of another variable, X. Given that we know that Y and X are related as Y=X2, the event Y(y,y+Δy) is the same as the event X2(x2,(x+Δx)2) which is the same as the event X(|x|,|x|+Δx) or X(|x|Δx,|x|). Thus, the individuals that are in the bin (y,y+Δy) must also be in the bins (|x|,|x|+Δx) and (|x|Δx,|x|). In other words, those bins must have the same proportion of individuals,

P(Y(y,y+Δy))=P(X(|x|,|x|+Δx))+P(X(|x|Δx,|x|))

Ok, now let's get to the density. First, we need to define what a probability density is. As the name suggests, it is the proportion of individuals per area. That is, we count the share of individuals on that bin and divide by the size of the bin. Since we have established that the proportions of people are the same here, but the size of the bins have changed, we conclude the density will be different. But different by how much?

As we said, the probability density is the proportion of people in the bin divided by the size of the bin, thus the density of Y is given by fY(y):=P(Y(y,y+Δy))Δy. Analogously, the probability density of X is given by fX(x):=P(X(x,x+Δx))Δx.

From our previous result that the population in each bin is the same we then have that,

fY(y):=P(Y(y,y+Δy))Δy=P(X(|x|,|x|+Δx))+P(X(|x|Δx,|x|))Δy=fX(|x|)Δx+fX(|x|)ΔxΔy=ΔxΔy(fX(|x|)+fX(|x|))=ΔxΔy(fX(y)+fX(y))

That is, the density fX(y)+fX(y) changes by the factor ΔxΔy, which is the relative size of stretching or squeezing the bin size. In our case, since y=x2 we have that y+Δy=(x+Δx)2=x2+2xΔx+Δx2. If Δx is tiny enough we can ignore Δx2, which implies Δy=2xΔx and ΔxΔy=12x=12y, and that is why the factor 12y shows up in the transformation.

Carlos Cinelli
la source