Come possiamo limitare la probabilità che una variabile casuale sia massima?


21

Supponiamo di avere variabili casuali indipendenti , , con mezzi finiti e varianze , \ ldots , \ sigma_N ^ 2 . Sto cercando limiti senza distribuzione sulla probabilità che qualsiasi X_i \ neq X_N sia più grande di tutti gli altri X_j , j \ neq i .X 1X n μ 1μ N σ 2 1NX1Xnμ1μNσ12σ 2 N X iX N X j j iσN2XiXNXjji

In altre parole, se per semplicità assumiamo che le distribuzioni di Xi siano continue (tale che P(Xi=Xj)=0 ), sto cercando limiti su:

P(Xi=maxjXj).
Se N=2 , possiamo usare la disuguaglianza di Chebyshev per ottenere:
P(X1=maxjXj)=P(X1>X2)σ12+σ22σ12+σ22+(μ1μ2)2.
Vorrei trovare alcuni limiti semplici (non necessariamente stretti) per N generale N, ma non sono stato in grado di trovare (esteticamente) risultati piacevoli per N generale N.

Si noti che non si presume che le variabili siano iid. Eventuali suggerimenti o riferimenti a lavori correlati sono i benvenuti.


Aggiornamento: ricorda che per ipotesi, μjμi . Possiamo quindi usare il limite sopra per arrivare a:

P(Xi=maxjXj)minj>iσi2+σj2σi2+σj2+(μjμi)2σi2+σN2σi2+σN2+(μNμi)2.
Ciò implica:
(μNμi)P(Xi=maxjXj)(μNμi)σi2+σN2σi2+σN2+(μNμi)212σi2+σN2.
Questo, a sua volta, implica:
i=1NμiP(Xi=maxjXj)μNN2i=1N1(σi2+σN2).
Ora sto chiedendo se questo legato può essere migliorata a qualcosa che non dipende linearmente N . Ad esempio, vale quanto segue: \ sum_ {i = 1} ^ N \ mu_i \ P (X_i = \ max_j X_j) \ geq \ mu_N - \ sqrt {\ sum_ {i = 1} ^ N \ sigma_i ^ 2} \ enspace?
i=1NμiP(Xi=maxjXj)μNi=1Nσi2?
E se no, quale potrebbe essere un controesempio?

3
Questo limite può essere più stretto se si utilizza l'indice che ti dà la più piccola limite superiore al posto di . Si noti che questo valore dipende sia dalla media che dalla varianza. NjN

5
@MichaelChernick: non credo sia corretto. Supponiamo ad esempio di avere tre distribuzioni uniformi su . Quindi, se non sbaglio, , mentre . Non so se intendevi scrivere , ma lo stesso esempio mostra che non è ancora valido. [0,1]P(X1<maxjXj)=2/3P(X1<X2)=P(X1<X3)=1/2P(Xi>maxjXj)
MLS,

2
@Michael: Purtroppo non è ancora vero. Gli eventi per fissi non sono indipendenti. Aj={Xi>Xj} i
cardinale il

2
@cardinal: tra le altre cose, è legato ai banditi multi-armati. Se scegli un braccio in base alle ricompense precedenti, quanto è grande la probabilità che tu abbia scelto il braccio migliore (che sarebbe nella notazione sopra), e possiamo limitare la perdita attesa per scegliere un sub braccio ottimale? P(XN=maxjXj)
MLS,

2
Crossposted to MathOverflow: mathoverflow.net/questions/99313
cardinale

Risposte:


1

Puoi usare la disuguaglianza multivariata di Chebyshev.

Caso di due variabili

Per una singola situazione, vs X 2 , arrivo alla stessa situazione del commento di Jochen del 4 novembre 2016X1X2

1) Se allora P ( X 1 > X 2 ) ( σ 2 1 + σ 2 2 ) / ( μ 1 - μ 2 ) 2μ1<μ2P(X1>X2)(σ12+σ22)/(μ1μ2)2

(e mi chiedo anche della tua derivazione)

Derivazione dell'equazione 1

  • usando la nuova variabile X1X2
  • trasformandolo in modo tale che abbia la media a zero
  • prendendo il valore assoluto
  • applicare la disuguaglianza di Chebyshev

P(X1>X2)=P(X1X2>0)=P(X1X2(μ1μ2)>(μ1μ2))P(|X1X2(μ1μ2)|>μ2μ1)σ(X1X2(μ1μ2))2(μ2μ1)2=σX12+σX22(μ2μ1)2

Caso multivariato

La disuguaglianza nell'equazione (1) può essere cambiata in un caso multivariato applicandolo a più variabili trasformate per ogni i < n (si noti che sono correlati).(XnXi)i<n

Una soluzione a questo problema (multivariata e correlata) è stata descritta da I. Olkin e JW Pratt. "Una disuguaglianza multivariata di Tchebycheff" negli Annali delle statistiche matematiche, volume 29 pagine 226-234 http://projecteuclid.org/euclid.aoms/1177706720

Nota teorema 2.3

P(|yi|kiσi for some i)=P(|xi|1 for some i)(u+(ptu)(p1))2p2

in cui il numero di variabili, t = k - 2 i e u = ρ i j / ( k i k j ) .pt=ki2u=ρij/(kikj)

Il teorema 3.6 fornisce un limite più stretto, ma è meno facile da calcolare.

modificare

Un limite più netto può essere trovato usando la disuguaglianza multivariata di Cantelli . Tale disuguaglianza è il tipo che hai usato prima e ti ha fornito il limite che è più nitido di ( σ 2 1 + σ 2 2 ) / ( μ 1 - μ 2(σ12+σ22)/(σ12+σ22+(μ1μ2)2) .(σ12+σ22)/(μ1μ2)2

Non ho avuto il tempo di studiare l'intero articolo, ma comunque, puoi trovare una soluzione qui:

AW Marshall e I. Olkin "Una disparità unilaterale del tipo Chebyshev" nel volume Annali di statistiche matematiche 31 pagg. 488-491 https://projecteuclid.org/euclid.aoms/1177705913

(nota successiva: questa disuguaglianza è per uguali correlazioni e non un aiuto sufficiente. Ma comunque il tuo problema, per trovare il limite più netto, è uguale alla, più generale, disuguaglianza multivariata di Cantelli. Sarei sorpreso se la soluzione non esistesse)


Potresti fornire una chiara dichiarazione della disuguaglianza multivariata di Chebyshev?
whuber

1
Ho modificato la soluzione fornendo l'intero teorema.
Sesto Empirico

-1

Ho trovato un teorema che potrebbe aiutarti e proverò ad adattarlo alle tue esigenze. Supponiamo di avere:

exp(tE(max1inXi))

Quindi dalla disuguaglianza di Jensen (poiché exp (.) È una funzione convessa), otteniamo:

exp(tE(max1inXi))E(exp(tmax1inXi))=E(max1in exp(tXi))i=1nE(exp(tXi)

exp(tXi you have to plug in whatever the moment generating function of your random variable Xi is (since it is just the definition of the mgf). Then, after doing so (and potentially simplifying your term), you take this term and take the log and divide by it by t so that you get a statement about the term E(max1inXi). Then you can choose t with some arbitrary value (best so that the term is small so that the bound is tight).

Then, you have a statement about the expected value of the maximum over n rvs. To get now the a statement about the probabilty that the maximum of those rv's deviates from this expected value, you can just use Markov's inequality (assuming that your rv is non-negative) or another, more specific rv, applying to your particular rv.

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.