Disuguaglianze di probabilità


37

Sto cercando alcune disuguaglianze di probabilità per somme di variabili casuali illimitate. Lo apprezzerei davvero se qualcuno potesse darmi qualche pensiero.

Il mio problema è trovare un limite esponenziale superiore alla probabilità che la somma delle variabili casuali iid illimitate, che sono in realtà la moltiplicazione di due iid gaussiane, superi un certo valore, ad esempio , dove , e sono generati da .Pr[Xϵσ2N]exp(?)X=i=1NwiviwiviN(0,σ)

Ho provato ad usare il limite di Chernoff usando la funzione di generazione del momento (MGF), il limite derivato è dato da:

Pr[Xϵσ2N]minsexp(sϵσ2N)gX(s)=exp(N2(1+4ϵ21+log(1+4ϵ21)log(2ϵ2)))

dove gX(s)=(11σ4s2)N2 è la MGF di X . Ma il limite non è così stretto. Il problema principale nel mio problema è che le variabili casuali sono illimitate e, sfortunatamente, non posso usare il limite della disuguaglianza di Hoeffding.

Sarò felice se mi aiuti a trovare qualche limite esponenziale stretto.


3
Sembra un problema correlato al rilevamento compresso. Guarda le note di R. Vershynin sulla teoria della matrice casuale non asintotica, in particolare i limiti di ciò che chiama variabili casuali subsponenziali . Questo ti farà iniziare. Se hai bisogno di più puntatori, faccelo sapere e proverò a pubblicare qualche informazione in più.
cardinale

1
Ci sono almeno un paio di domande e risposte relative a questo argomento su math.SE (disclaimer: incluso uno a cui ho partecipato).
cardinale

1
Il prodotto wivi ha come distribuzione 'prodotto normale'. Credo che la media di questo prodotto sia zero e la varianza sia σ4 dove σ2 è la varianza di wi e vi . Per N largeish, è possibile utilizzare il teorema del limite centrale per ottenere norality approssimativo di X . Se riesci a calcolare l'inclinazione della normale distribuzione del prodotto, credo che tu possa applicare il teorema di Berry-Esseen per limitare il tasso di convergenza del CDF.
shabbychef,

1
@shabbychef, Berry-Esseen ha piuttosto lento convergenza, dal momento che è una divisa legata sulla classe di tutte le funzioni di distribuzione . F
cardinale

4
@DilipSarwate: mi dispiace che sto vedendo il tuo commento da poco tempo fa. Penso che potresti essere interessato al seguente piccolo articolo, che ho collegato un paio di volte anche su Math.SE: TK Phillips e R. Nelson (1995), Il momento legato è più stretto di quello di Chernoff per la coda positiva probabilità , The American Statistician , vol 42, no. 2., 175-178.
cardinale

Risposte:


1

Usando il limite di Chernoff che hai suggerito per alcuni che verranno specificati in seguito, dove vale la seconda disuguaglianza grazie a per qualsiasi . Ora prendi e , il lato destro diventa che produce per qualsiasi .s1/(2σ2)- log ( 1 - x )

P[X>t]exp(st)exp((N/2)log(1σ4s2))exp(st+σ4s2N)
log(1x)2xx(0,1/2)t=ϵσ2Ns=t/(2σ4N)exp(t2/(4σ4N)=exp(ϵ2N/4)
P[X>ϵσ2N]exp(ϵ2N/4).
ϵ(0,1)

Un'altra strada è applicare direttamente le disuguaglianze di concentrazione come la disuguaglianza di Hanson-Wright o le disuguaglianze di concentrazione per il caos gaussiano di ordine 2 che comprende la variabile casuale a cui sei interessato.

Approccio più semplice senza usare la funzione di generazione del momento

Prendi per semplicità (altrimenti, si può ridimensionare dividendo per ).σ=1σ2

Scrivi e . Stai chiedendo limiti superiori su .v=(v1,...,vn)Tw=(w1,...,wn)TP(vTw>ϵN)

Sia. Quindi per indipendenza di e è indipendente da con la distribuzione con gradi di libertà.Z=wTv/vZN(0,1)v,wv2Zχ2n

Per limiti standard su variabili normali normali e casuali, Combinando con il limite di unione si ottiene un limite superiore di della forma .χ2

P(|Z|>ϵn/2)2exp(ϵ2n/4),P(v>2n)exp(n(21)2/2).
P(vTw>ϵN)2exp(ϵ2n/4)+exp(n(21)2/2)


0

Il limite che ottieni è dell'ordine come . Non penso che tu possa fare molto meglio per il generale . Dalla pagina Wikipedia sulle variabili di prodotto la distribuzione di è dove è una funzione di Bessel modificata. Da (10.25.3) nell'elenco delle funzioni DLMF , modo che per sufficientemente grande che non ti darà un limite sub-gaussiano.eϵϵϵwiviK0(z)/πK0K0(t)et/txP(wivi>x)xet/tdt

Utilizzando il nostro sito, riconosci di aver letto e compreso le nostre Informativa sui cookie e Informativa sulla privacy.
Licensed under cc by-sa 3.0 with attribution required.