Biorąc pod uwagę n równomiernie rozłożonych rv, jaki jest PDF dla jednego rv podzielony przez sumę wszystkich nv.v?


10

Interesuje mnie następujący rodzaj przypadku: istnieją ciągłe zmienne losowe „n”, które muszą sumować się do 1. Co to byłby plik PDF dla każdej pojedynczej takiej zmiennej? Jeśli więc n=3 , to jestem zainteresowany rozkładem dla X1X1+X2+X3 , gdzieX1,X2iX3są rozmieszczone równomiernie. Średni oczywiście w tym przykładzie jest1/3, a średnia wynosi tylko1/n, a mimo to jest łatwe do dystrybucji symulacji w R, nie wiem, co rzeczywiste równanie PDF lub CDF jest.

Ta sytuacja jest związana z dystrybucją Irwin-Hall ( https://en.wikipedia.org/wiki/Irwin%E2%80%93Hall_distribution ). Tylko Irwin-Hall jest rozkładem sumy n jednolitych zmiennych losowych, podczas gdy chciałbym rozkład dla jednego z n jednolitych zmiennych rv podzielony przez sumę wszystkich n zmiennych. Dzięki.


1
Jeśli n ciągłych jednorodnych zmiennych losowych sumuje się do 1 , to przy n=3 , X1+X2+X3=1 a więc rozkład X1X1+X2+X3=X1jest taki sam jak rozkładX1, prawda?
— Dilip Sarwate,

1
Powinienem się poprawić: rozkłady N jednolitych nie sumują się do 1. Zakładam, że każdy z nich jest jednolity między 0 a 1, więc ich suma może wynosić od 0 do N. Myślę o wzięciu każdej zmiennej jednolitej i podzieleniu przez sumę wszystkich N zmiennych zmiennych, aby otrzymać zestaw N zmiennych losowych, które sumują się do 1 i mają wartość oczekiwaną 1 / N. Uwaga: usunąłem słowo „mundur” z pierwszego zdania. Rozkład, którego szukam, nie jest jednolity, ale pochodzi od podzielenia jednej z N zmiennych jednolitych przez sumę wszystkich N zmiennych zmiennych. Po prostu nie jestem pewien jak.
— user3593717,

W którym Xi są rozmieszczone w postępie geometrycznym, wektor zmiennych znormalizowanych ma rozkład dirichleta. Może to być interesujące samo w sobie, ale przeanalizowane może również zapewnić taktykę dla tego rodzaju sytuacji.
— przypuszcza

Odpowiedzi:


4

Punkty przerwania w domenie sprawiają, że jest nieco niechlujny. Proste, ale żmudne podejście polega na osiągnięciu ostatecznego wyniku. Dla n=3, niech Y=X2+X3, W=X2+X3X1,aT=1+W.NastępnieZ=1T=X1X1+X2+X3.

Wartości graniczne są 1 do 1 i 2 w W , 2 i 3, do T , i 1 / 3 i 1 / 2 do Z . Znalazłem pełny pdfY,W,T,1/31/2Z.

f(z)={     1(1−z)2 ,if 0≤z≤1/33z3−9z2+6z−13z3(1−z)2 ,if 1/3≤z≤1/2       1−z3z3 ,if 1/2≤z≤1

Plik cdf można następnie znaleźć jako

F(z)={           z(1−z) ,if 0≤z≤1/312+−18z3+24z2−9z+16z2(1−z) ,if 1/3≤z≤1/2        56+2z−16z2 ,if 1/2≤z≤1

+1 Nice. Also, your density agrees beautifully with simulation.
— Glen_b -Reinstate Monica

2

Let Y=∑i=2nXi. We can find the cdf of X1/∑i=1nXi by calculating

P(X1∑i=1nXi≤t)=P(X1≤t∑i=1nXi)=P((1−t)X1≤t∑i=2nXi)=P(X1≤t1−tY)=∫01P(x1≤t1−tY) dx1=∫01(1−FY(1−ttx1)) dx1=1−∫01FY(1−ttx1) dx1
We then differentiate and substitute the Irwin-Hall pdf to obtain the desired pdf:
f(t)=∫01fY(1−ttx1)⋅x1t2 dx1=1t2∫01∧(n−1)t1−t∑k=0⌊1−ttx1⌋1(n−2)!(−1)k(n−1k)(1−ttx1−k)n−1x1 dx1
From here it gets a little messy, but you should be able to interchange the integral and summation and then perform a substitution (e.g, u=tx11−t−k) to evaluate the integral and hence obtain an explicit formula for the pdf.

1

Assuming

"the N uniform distributions don't sum to 1."

This is how I started(it's incomplete):

Consider Y=∑i=1nXi and let X=Xi by a slight abuse of notation.

Consider, U=XY and V=Y:

X=UVY=V

Then following transformation of variables:

J=[VU01]

The joint probability function of (U,V) is given by:

fU,V(u,v)=fX,Y(uv,v)|J|

Where X∼U(0,1) and Y∼IrwinHall

fX(x)={10≤x≤10otherwise

And,

fY(y)=12(n−1)!∑k=0n(−1)k(nk)(x−k)n−1sign(x−k)

Thus,

fU,V(u,v)={12(n−1)!∑k=0n(−1)k(nk)(uv−k)n−1sign(uv−k)0≤uv≤10otherwise

and fU(u)=∫fU,V(u,v)dv


0

Suppose we already know sum of U(0,1) has a Irwin-Hall distribution. Now your question changes to find the pdf (or CDF) of XY when X had a U(0,1) distribution and Y has a Irwin-Hall distribution.

First we need to find he joint pdf of X and Y.

Let Y1=X1Y2=X1+X2Y3=X1+X2+X3

Then

X1=Y1X2=Y2−Y1X3=Y3−Y2−Y1

∴

J=|∂X1∂Y1∂X1∂Y2∂X1∂Y3∂X2∂Y1∂X2∂Y2∂X2∂Y3∂X3∂Y1∂X3∂Y2∂X3∂Y3|=−1

Since X1,X2,X3 are i.i.d with U(0,1), therefore, f(x1,x2,x3)=f(x1)f(x2)f(x3)=1

The joint distribution with y1,y2,y3 is

g(y1,y2,y3)=f(y1,y2,y3)|J|=1

Next let us integrate out the Y2 and we can get the joint distribution of Y1 and Y3 i.e the joint distribution of X1 and X1+X2+X3

As suggested by whuber now I changed the the limits

(1)h(y1,y3)=∫y1+1y3−1g(y1,y2,y3)dy2=∫y1+1y3−11dy2=y3−y1−2

Now, we know the joint pdf of X,Y i.e joint pdf X1 and X1+X2+X3 is y3−y1−2.

Next let find the pdf of XY

We need another transformation:

Let Y1=XY2=XY

Then X=Y1Y=Y1Y2

Then

J=|∂x∂y1∂x∂y2∂y∂y1∂y∂y2|=|101y2−y1y22|=−y1y22

we already the joint distribution of X,Y from above steps ref (1).

∴

g2(y1,y2)=h(y1,y3)|J|=(y3−y1−2)y1y22

Next, we integrate the y1 out we get the pdf of y2 then we get the pdf of XY

(2)h2(y2)=∫01(y3−y1−2)y1y22dy1=1y22(y32−13−1)

This is the pdf of X/Y i.e X1X1+X2+X3

We are not finish yet, what is y3 in (2) then?

We know that Y3=X1+X2+X3 from the first transformation.

So at least we know Y3 has a Irwin-Hall distribution.

I wonder can we plug the Irwin-Hall for n=3 pdf to (2) to get a explicit formula? or can we do some simulations from here as Glen suggested?


2
Simulation doesn't seem to agree with that pdf.
— Glen_b -Reinstate Monica

The logic and steps seem correct, but I feel uncomfortable about this solution.
— Deep North

2
Where you integrated out y2, you needed to account for the conditions y1≤y2≤y3 and y3−1≤y2≤y1+1.
— whuber
Korzystając z naszej strony potwierdzasz, że przeczytałeś(-aś) i rozumiesz nasze zasady używania plików cookie i zasady ochrony prywatności.
Licensed under cc by-sa 3.0 with attribution required.