Jeśli ,

9

Załóżmy następującą konfigurację:
Niech $Z_i = \min\{k_i, X_i\}, i=1,...,n$ . Również $X_i \sim U[a_i, b_i], \; a_i, b_i >0$ . Ponadto $k_i = ca_i + (1-c)b_i,\;\; 0<c<1$ tj. $k_i$ jest wypukłą kombinacją granic odpowiednich podpór. $c$ jest wspólne dla wszystkich $i$ .

Myślę, że mam prawidłowy rozkład $Z_i$ : jest to rozkład mieszany .
Ma część ciągłą,

X_{i} \in [a_{i}, k_{i}), Z_{i} = X_{i} \Rightarrow Pr (Z_{i} \leq z_{i}) = \frac{z_{i} - a_{i}}{b_{i} - a_{i}}

$X_i \in [a_i, k_i), Z_i=X_i \Rightarrow \Pr(Z_i \le z_i) = \frac {z_i-a_i}{b_i-a_i}$ a następnie nieciągłość i dyskretna część gdzie prawdopodobieństwo prawdopodobieństwa koncentruje się:

Pr (Z_{i} = k_{i}) = Pr (X_{i} > k_{i}) = 1 - Pr (X_{i} \leq k_{i})

$\Pr(Z_i=k_i) = \Pr(X_i > k_i) = 1- \Pr(X_i \le k_i)$

= 1 - \frac{k_{i} - a_{i}}{b_{i} - a_{i}} = 1 - \frac{(1 - c) (b_{i} - a_{i})}{b_{i} - a_{i}} = c

$= 1- \frac {k_i - a_i}{b_i-a_i} = 1-\frac {(1-c)(b_i-a_i)}{b_i-a_i} =c$

Więc we wszystkich

{fa}_{Z_{ja}} (z_{ja}) = {\begin{cases} 0 z_{ja} < {za}_{ja} \\ \frac{z_{ja} - {za}_{ja}}{b_{ja} - {za}_{ja}} {za}_{ja} \leq z_{ja} < k_{ja} \\ 1 k_{ja} \leq z_{ja} \end{cases}

$F_{Z_i}(z_i) = \begin{cases} 0\qquad z_i<a_i\\ \\ \frac {z_i-a_i}{b_i-a_i}\qquad a_i\le z_i<k_i \\ \\1\qquad k_i\le z_i\end{cases}$

podczas gdy dla mieszanej funkcji „dyskretna / ciągła” masa / gęstość wynosi $0$ poza przedziałem $[a_i, k_i]$ , ma ona część ciągłą, która jest gęstością jednolitego $U(a_i, b_i)$ , $\frac {1}{b_i-a_i}$ ale dla $a_i\le z_i<k_i$ , i koncentruje masę prawdopodobieństwa dodatniego $c >0$ przy $z_i = k_i$ .

Podsumowując, sumuje się do jedności nad rzeczywistością.

Chciałbym móc wyprowadzić lub powiedzieć coś o rozkładzie i / lub momentach zmiennej losowej $S_n \equiv \sum_{i=1}^n Z_i$ , jako $n\rightarrow \infty$ .

Powiedzmy, że jeśli są niezależne, wygląda na to, że jako . Czy mogę „zignorować” tę część, nawet jako przybliżenie? Wtedy zostałbym ze zmienną losową, która zawiera się w przedziale , wyglądając jak suma cenzurowanych mundurów, w drodze do zostania "nieocenzurowanymi", a więc może jakieś centralne twierdzenie graniczne ... ale raczej rozbieżę się, zamiast się tutaj zbiegać, więc jakieś sugestie? $X_i$ $\Pr(S_n = \sum_i^nk_i) = c^n \rightarrow 0$ $n\rightarrow \infty$ $[\sum_{i=1}^na_i,\; \sum_{i=1}^nk_i)$

PS: To pytanie jest istotne, wyprowadzając rozkład sumy zmiennych ocenzurowanych , ale odpowiedź @Glen_b nie jest tym, czego potrzebuję - muszę pracować nad tym analitycznie, nawet przy użyciu przybliżeń. To są badania, więc proszę traktować je jak pracę domową - ogólne sugestie lub odniesienia do literatury są wystarczająco dobre.

distributions central-limit-theorem censoring minimum Alecos Papadopoulos
źródło

Jeśli potrzebujesz, napisz rozkład jako , z odpowiednim , w którym jest zestawem Borela.

Z_{i}

$Z_i$

μ_{Z_{i}} (B) = P (Z_{i} \in B) = \int_{B} g (t) d t + c I_{B} (k_{i})

$\mu_{Z_i}(B)=P(Z_i\in B)=\int_B g(t)\,dt +c\,I_B(k_i)$

g

$g$

B

$B$

Zen

@ Zen Już napisałem w pytaniu, że dystrybucja jest nieciągła. Również RHS sprawia, że oczywiste jest, że oznacza gęstość w , ale prawdopodobieństwo prawdopodobieństwa - i wolę notację zwartą.

f

$f$

f

$f$

[a_{i}, k_{i})

$[a_i,k_i)$

k_{i}

$k_i$

Alecos Papadopoulos

O ile mi wiadomo, ten zapis był pdf, a pmf nie istnieje; i mamy odpowiedni język matematyczny, aby precyzyjnie opisywać mieszane rozkłady. Wątpię, aby ta notacja została zaakceptowana, kiedy opublikujesz swoje badania. Oczywiście tylko moja opinia. Zawsze powinieneś robić to tak, jak lubisz.

f

$f$

Zen

@Zen Publishing jest daleko przed nami - i rzeczywiście, recenzenci marszczą brwi, widząc nieokreślony zapis. Ten jest tylko skrótem, gdy chce się opisać rozkład stopniowy w wielu wierszach. Nie ma „argumentu za” i przeciw ustalonemu zapisowi, jak na przykład ten, którego użyłeś w poprzednim komentarzu.

Alecos Papadopoulos

5

Podążę za wskazówkami Henry'ego i sprawdzę Lyapunov z . Fakt, że rozkłady są mieszane, nie powinien stanowić problemu, pod warunkiem, że i zachowują się poprawnie. Symulacja konkretnego przypadku, w którym , , dla każdego pokazuje, że normalność jest w porządku. $\delta=1$ $a_i$ $b_i$ $a_i=0$ $b_i=1$ $k_i=2/3$ $i\geq 1$

xbar <- replicate(10^4, mean(pmin(runif(10^4), 2/3)))
hist((xbar - mean(xbar)) / sd(xbar), breaks = "FD", freq = FALSE)
curve(dnorm, col = "blue", lwd = 2, add = TRUE)

CLT

Zen
źródło

Rzeczywiście całkiem normalne. Dobrze wiedzieć. Zwykłe warunki dla CLT nigdy nie były tutaj problemem, moje pytanie brzmiało, czy istnieją inne, być może subtelne problemy, które przekręcają asymptotyczny wynik i wymagają zmodyfikowanego CLT. Twoja symulacja pokazuje, że w rzeczywistości prawdopodobieństwo nieciągłości dyskretnej staje się nieistotne, gdy do zmiennej sumuje się więcej zmiennych.

Alecos Papadopoulos

Nic konkretnego, ale nie stwarzają żadnych problemów. Pomyśl o nich również o zachowanych liczbach skończonych, niezależnych od indeksu . Mogą one zwiększać lub zmniejszać, jak rośnie (nie reguła konkretnego), a nie każdy z nich jest nieproporcjonalnie większy niż inni ... reprezentują różnice w wielkości Niemniej „porównywalny” podmiotów. Zatem stan Lindeberga z pewnością się utrzymuje

i

$i$

i

$i$

Alecos Papadopoulos

Miły. Powodzenia w kolejnych krokach. Wygląda na ciekawy problem.

Zen

3

Poradnik:

Zakładając, że jest ustalone, a są niezależne, możesz obliczyć średnią i wariancję każdego : na przykład i wiesz, że . $c$ $X_i$ $\mu_i$ $\sigma_i^2$ $Z_i$ $\mu_i=E[ Z_i] = c\frac{a_i+k_i}{2} + (1-c)k_i$ $k_i = ca_i + (1-c)b_i$

Następnie, pod warunkiem, że i nie rosną zbyt szybko, możesz użyć warunków Lapunowa lub Lindeberga, aby zastosować twierdzenie o limicie centralnym, z wnioskiem, że zbiega się w rozkładzie do standardowej normalnej lub w sensie machania ręką jest w przybliżeniu normalnie dystrybuowane ze średnią i wariancja . $a_i$ $b_i$ $\displaystyle\frac{1}{\sqrt{\sum_1^n \sigma_i^2}}\left(\sum_1^n Z_i - \sum_1^n \mu_i\right)$ $\sum_1^n Z_i$ $\sum_1^n \mu_i$ $\sum_1^n \sigma_i^2$

Henz
źródło

Dzięki. Nie ma problemu z i , nie rosną wraz z indeksem, po prostu się zmieniają. Mówisz więc zasadniczo, że CLT może obejmować również zmienne losowe o rozkładach mieszanych?

a_{i}

$a_i$

b_{i}

$b_i$

Alecos Papadopoulos

Gdyby np. i zostały ustalone, wtedy mielibyśmy niezależne identyczne zmienne losowe o skończonej wariancji, więc zastosowanie miałoby centralne twierdzenie graniczne. To, czy jest to rozkład mieszanki, czy nie, nie wpływa na ten wynik. Mówię o tym, że można to rozszerzyć na przypadki, w których zmienne losowe są niezależne, ale nie identycznie rozmieszczone, pod warunkiem, że średnie i wariancje pozostają rozsądne.

a_{i}

$a_i$

b_{i}

$b_i$

Henry

2

Moje główne zmartwienie w tym pytaniu dotyczyło tego, czy można zastosować CLT „jak zwykle” w przypadku, który badam. Użytkownik @Henry stwierdził, że można, użytkownik @Zen pokazał to poprzez symulację. Tak zachęcony, teraz udowodnię to analitycznie.

Najpierw sprawdzę, czy ta zmienna z rozkładem mieszanym ma „zwykłą” funkcję generowania momentu. Oznaczają wartość oczekiwana , jego odchylenie standardowe, a na środku i skalowany wersja o . Stosując wzór zmiany zmiennej, okazuje się, że część ciągła to Funkcja generowania momentu powinna być $\mu_i$ $Z_i$ $\sigma_i$ $Z_i$ $\tilde Z_i = \frac {Z_i-\mu_i}{\sigma_i}$

{fa}_{\tilde{Z}} ({\tilde{z}}_{ja}) = σ_{ja} {fa}_{Z} (z_{ja}) = \frac{σ_{ja}}{b_{ja} - {za}_{ja}}

$f_{\tilde Z}(\tilde z_i) = \sigma_if_Z(z_i) = \frac {\sigma_i}{b_i-a_i}$

{\tilde{Z}}_{i}

$\tilde Z_i$

{\tilde{M.}}_{ja} (t) = mi ({mi}^{{\tilde{z}}_{ja} t}) = \int_{- \infty}^{\infty} {mi}^{{\tilde{z}}_{ja} t} re {fa}_{\tilde{Z}} ({\tilde{z}}_{ja}) = \int_{{\tilde{za}}_{ja}}^{{\tilde{k}}_{ja}} \frac{σ_{ja} {mi}^{{\tilde{z}}_{ja} t}}{b_{ja} - {za}_{ja}} re z_{ja} + do {mi}^{{\tilde{k}}_{ja} t}

$\tilde M_i(t) = E(e^{\tilde z_it}) = \int_{-\infty}^{\infty}e^{\tilde z_it}dF_{\tilde Z}(\tilde z_i) = \int_{\tilde a_i}^{\tilde k_i}\frac{\sigma_ie^{\tilde z_it}}{b_i-a_i}dz_i + ce^{\tilde k_it}$

\Rightarrow {\tilde{M.}}_{ja} (t) = \frac{σ_{ja}}{b_{ja} - {za}_{ja}} \frac{{mi}^{{\tilde{k}}_{ja} t} - {mi}^{{\tilde{za}}_{ja} t}}{t} + do {mi}^{{\tilde{k}}_{ja} t}

$\Rightarrow \tilde M_i(t)=\frac {\sigma_i}{b_i-a_i}\frac{e^{\tilde k_it}-e^{\tilde a_it}}{t} +ce^{\tilde k_it}$ z

{\tilde{k}}_{ja} = \frac{k_{ja} - μ_{ja}}{σ_{ja}}, {\tilde{za}}_{ja} = \frac{{za}_{ja} - μ_{ja}}{σ_{ja}}

$\tilde k_i = \frac {k_i-\mu_i}{\sigma_i},\;\; \tilde a_i = \frac {a_i-\mu_i}{\sigma_i}$

Używając liczb pierwszych do oznaczania pochodnych, jeśli poprawnie podaliśmy funkcję generowania momentu, powinniśmy uzyskać od tego czasu jest wyśrodkowaną i skalowaną zmienną losową. I rzeczywiście, obliczając pochodne, stosując wielokrotnie regułę L'Hopital (ponieważ wartość MGF na zero musi być obliczona przez granice) i wykonując manipulacje algebraiczne, zweryfikowałem dwie pierwsze równości. Trzecia równość okazała się zbyt męcząca, ale wierzę, że tak jest.

{\tilde{M.}}_{ja} (0) = 1, {\tilde{M.}}_{ja}^{'} (0) = mi (\tilde{Z}) = 0 \Rightarrow {\tilde{M.}}_{ja}^{″} (0) = mi ({\tilde{Z}}_{ja}^{2)}) = Var ({\tilde{Z}}_{ja}) = 1

$\tilde M_i(0) = 1, \;\; \tilde M_i'(0) = E(\tilde Z) = 0 \Rightarrow \tilde M_i''(0) = E(\tilde Z_i^2) = \operatorname {Var}(\tilde Z_i)=1$

Mamy więc odpowiedni MGF. Jeśli weźmiemy rozszerzenie Taylora drugiego rzędu około zera, to mamy

\tilde{M.} (t) = \tilde{M.} (0) + {\tilde{M.}}^{'} (0) t + \frac{1}{2)} {\tilde{M.}}^{″} (0) t^{2)} + o (t^{2)})

$\tilde M(t) = \tilde M(0) + \tilde M'(0)t +\frac 12\tilde M''(0)t^2 + o(t^2)$

\Rightarrow \tilde{M.} (t) = 1 + \frac{1}{2)} t^{2)} + o (t^{2)})

$\Rightarrow \tilde M(t) = 1 + \frac 12t^2+ o(t^2)$

Oznacza to, że funkcją charakterystyczną jest (tutaj oznacza jednostkę urojoną) . $i$

\tilde{ϕ} (t) = 1 + \frac{1}{2)} (ja t)^{2)} + o (t^{2)}) = 1 - \frac{1}{2)} t^{2)} + o (t^{2)})

$\tilde \phi(t) = 1 + \frac 12 (it)^2 + o(t^2)= 1 - \frac 12 t^2 + o(t^2)$

Z właściwości funkcji charakterystycznej wynika, że funkcja charakterystyczna jest równa $\tilde Z/\sqrt n$

{\tilde{ϕ}}_{\tilde{Z} / \sqrt{n}} (t) = {\tilde{ϕ}}_{\tilde{Z}} (t / \sqrt{n}) = 1 - \frac{t^{2)}}{2) n} + o (t^{2)} / n)

$\tilde \phi_{\tilde Z/\sqrt n}(t)=\tilde \phi_{\tilde Z}(t/\sqrt n) = 1 - \frac {t^2}{2n} + o(t^2/n)$

a ponieważ mamy niezależne zmienne losowe, charakterystyczną funkcją jest $\frac 1{\sqrt n}\sum_i^n\tilde Z_i$

{\tilde{ϕ}}_{\frac{1}{\sqrt{n}} \sum_{ja}^{n} {\tilde{Z}}_{ja}} (t) = \prod_{ja = 1}^{n} {\tilde{ϕ}}_{\tilde{Z}} (t / \sqrt{n}) = \prod_{ja = 1}^{n} (1 - \frac{t^{2)}}{2) n} + o (t^{2)} / n))

$\tilde \phi_{\frac 1{\sqrt n}\sum_i^n\tilde Z_i}(t)= \prod_{i=1}^n\tilde \phi_{\tilde Z}(t/\sqrt n)=\prod_{i=1}^n\left(1 - \frac {t^2}{2n} + o(t^2/n)\right)$

Następnie

lim_{n \to \infty} {\tilde{ϕ}}_{\frac{1}{\sqrt{n}} \sum_{ja}^{n} {\tilde{Z}}_{ja}} (t) = lim_{n \to \infty} {(1 - \frac{t^{2)}}{2) n})}^{n} = {mi}^{- t^{2)} / 2)}

$\lim_{n\rightarrow \infty}\tilde \phi_{\frac 1{\sqrt n}\sum_i^n\tilde Z_i}(t) = \lim_{n\rightarrow \infty}\left(1 - \frac {t^2}{2n}\right)^n = e^{-t^2/2}$

przez jak liczba jest reprezentowany $e$ . Zdarza się, że ostatni człon jest funkcją charakterystyczną standardowego rozkładu normalnego i według twierdzenia Levy'ego o ciągłości mamy

\frac{1}{\sqrt{n}} \sum_{ja}^{n} {\tilde{Z}}_{ja} \overset{re}{\to} N. (0, 1)

$\frac 1{\sqrt n}\sum_i^n\tilde Z_i \xrightarrow{d} N(0,1)$

którym jest CLT. Zauważ, że fakt, że zmienne - nie są identycznie rozmieszczone, „zniknął” z widoku, kiedy rozważymy ich wyśrodkowane i skalowane wersje i rozważymy rozszerzenie Taylora drugiego rzędu ich MGF / CHF: na tym przybliżeniu funkcje te są identyczne, a wszystkie różnice są zagęszczone w pozostałych kategoriach, które zanikają asymptotycznie. $Z$

Fakt, że idiosynkratyczne zachowanie na poziomie indywidualnym, ze wszystkich pojedynczych elementów, znika jednak, gdy weźmiemy pod uwagę średnie zachowanie, uważam, że bardzo dobrze jest to pokazane przy użyciu paskudnego stworzenia, takiego jak zmienna losowa o mieszanym rozkładzie.

Alecos Papadopoulos
źródło

Naprawdę fajnie, Alecos. Wydaje mi się, że argument powinien zależeć od bardziej szczegółowych warunków i . Na przykład: czy dowód psuje się, jeśli szybko ? (Wiem, że w twojej aplikacji tak się nie dzieje.) Co myślisz?

a_{i}

$a_i$

b_{i}

$b_i$

(b_{i} - a_{i}) ↓ 0

$(b_i-a_i)\downarrow 0$

Zen

@Zen Kwestia wariancji niezależnych, ale nie identycznie rozłożonych RV jest bardzo subtelna, nie sądzę, że nadal rozumiem ją jasno. Znane warunki Lyapunova lub Lindeberga są wystarczające tylko do utrzymania CLT. Są przypadki, w których CLT utrzymuje się, mimo że te warunki nie. Myślę więc, że jeśli nie ograniczymy wariancji, nie ma jednej odpowiedzi, a problem staje się całkowicie specyficzny dla każdego przypadku. Nawet książka Billingsleya nie jest jasna w tej sprawie. Pytanie brzmi, jak będzie wyglądać reszta i co możemy o tym powiedzieć.

Alecos Papadopoulos

Jeśli ,

Odpowiedzi: