Gaussiani indipendenti a coppie


12

Dato (tra gaussiani con media e varianza ), è possibile (come? ) (per ) tale che siano gaussiani indipendenti a coppie con media e varianza . 0 1 m = k 2 Y 1 , , Y m Y i 0 1X1,,Xk01m=k2Y1,,YmYi01


1
@Suresh, quindi non sembra funzionare. E[(Xi+Xj)(Xi+Xk)]=E[Xi2]=1
Kaveh,

4
Non so perché, ma trovo la risposta MO a questa domanda abbastanza divertente (a parte il puntatore a stats.SE): mathoverflow.net/questions/46180/…
Suresh Venkat,

2
Quello che stavo cercando era qualcosa come prendere combinazioni lineari (che ovviamente non funzionano) o polinomi ecc. (Che non funzionano immediatamente) ma non riesco davvero a pensare a nessuna nozione ragionevole che la risposta di Shai su mathoverflow non soddisfa.

2
forse dovresti aggiornare la domanda sottolineando la risposta su MO?
Suresh Venkat,

2
Hai bisogno di una distribuzione congiunta gaussiana? In tal caso, ciò di cui hai bisogno sembra impossibile poiché tale distribuzione è determinata dalla sua matrice di covarianza e quindi l'indipendenza a coppie e la piena indipendenza sarebbero le stesse.
MCH

Risposte:


4

La pubblicazione su MathOverflow spiega come passare da un piccolo numero di variabili casuali Uniform [0,1] indipendenti a un numero maggiore di variabili casuali Uniform [0,1] indipendenti. Ovviamente puoi andare avanti e indietro tra Uniform [0,1] e Gaussian invertendo il CDF. Ma ciò richiede un'analisi numerica poiché il CDF non è in forma chiusa.

Tuttavia, esiste un modo più semplice di passare dall'uniforme gaussiana all'uniforme. Dati due Gaussiani indipendenti , l'angolo è uniforme nell'intervallo . arctan ( X 1 / X 2 ) [ 0 , 2 π ]X1,X2arctan(X1/X2)[0,2π]

Allo stesso modo, il metodo Box-Muller trasforma due variabili Uniform [0,1] indipendenti in due variabili casuali gaussiane indipendenti.

Usando queste due trasformazioni, consumi due gaussiani per produrre un'uniforme o due uniformi per produrre un gaussiano. Quindi c'è solo un fattore di nell'efficienza di campionamento. Inoltre, non è richiesta alcuna inversione del normale cdf.O(1)


-2

Questa costruzione NON fornisce variabili indipendenti a coppie (anzi, sotto) come richiesto da Anindya, ma fornisce variabili non correlate a coppie che sono sufficienti per ottenere buoni limiti di concentrazione per la somma attraverso la disuguaglianza di Chebyshev (e questo è molte volte l'obiettivo finale).|Yi,j|=|Yi,j|

Per ogni coppia distinta , lascia , dove è la funzione del segno. È chiaro che ogni è una variabile normale con media 0 e varianza 1. Per vedere che sono ortogonali, per , nota che che può essere facilmente verificato per uguagliare 0 osservando i vari casi di possibili uguaglianze tra .(i,j)([k]2)Yi,j=|Xi|σ(XiXj)σ()Yi,j(i,j)(i,j)

E[Yi,jYi,j]=E[|XiXi|σ(XiXiXjXj)]
i,i,j,j

PS: una versione precedente rivendicava falsamente l'indipendenza a coppie.


Non riesco a capire perché la media del prodotto pari a zero implicherebbe l'indipendenza.
Tsuyoshi Ito,

@TsuyoshiIto: le tue critiche erano corrette, ovviamente. Ho ancora lasciato questa risposta, poiché penso sia interessante.
Arnab

2
Se desideri mantenere il tuo post, utilizza le precauzioni necessarie per evitare di confondere i lettori. Puoi sostenere che la versione corrente (revisione 3) del tuo post non indica nulla di errato. Vero, ma la domanda pone qualcosa e il tuo post risponde a qualcos'altro senza dichiararlo. Per favore, comprendi che è estremamente confuso per i lettori.
Tsuyoshi Ito,
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.