La stazionarietà viene preservata sotto una combinazione lineare?


12

Immagina di avere due processi di serie temporali, che sono fissi, producendo: .xt,yt

È , stazionaria anche?α , β Rzt=αxt+βytα,βR

Qualsiasi aiuto sarebbe apprezzato.

Direi di sì, dal momento che ha una rappresentazione MA.


1
Perché è garantito che sia MA? ci sono processi AR stabili. Ad ogni modo, se stai parlando di stabilità BIBO, allora la somma è banalmente stabile perché puoi calcolare i nuovi limiti. La stabilità asintotica è valida anche perché limtzt=αlimtxt+βlimtyt
Steve Cox

Relativamente ad alcune estensioni: Nota nell'analisi numerica, usi quello che viene chiamato un precondizionatore (una particolare trasformazione lineare) per guadagnare stabilità, quindi dubito che la risposta sia sì.
Surb

Risposte:


26

Forse sorprendentemente, questo non è vero. (L'indipendenza delle due serie storiche lo renderà vero, tuttavia.)

Capisco che "stabile" significa stazionario, perché quelle parole sembrano essere usate in modo intercambiabile in milioni di risultati della ricerca, incluso almeno uno sul nostro sito .

Per un controesempio, sia una serie temporale stazionaria non costante per la quale ogni è indipendente da , e le cui distribuzioni marginali sono simmetriche attorno a . DefinireX t X s s t , 0XXtXsst,0

Yt=(1)tXt.

! [Figura 1: grafici di X, Y e (X + Y) / 2 nel tempo

Queste trame mostrano parti delle tre serie storiche discusse in questo post. stato simulato come una serie di disegni indipendenti da una distribuzione normale standard.X

Per dimostrare che è stazionario, dobbiamo dimostrare che la distribuzione congiunta di per qualsiasi non dipende da . Ma questo deriva direttamente dalla simmetria e dall'indipendenza . ( Y s + t 1 , Y s + t 2 , , Y s + t n ) t 1 < t 2 < < t n s X tY(Ys+t1,Ys+t2,,Ys+tn)t1<t2<<tnsXt

Figura che mostra alcuni grafici a dispersione incrociata di Y

Questi grafici a dispersione ritardati (per una sequenza di 512 valori di ) illustrano l'affermazione secondo cui le distribuzioni bivariate congiunte di sono come attese: indipendenti e simmetriche. (Un "diagramma di dispersione ritardato" visualizza i valori di rispetto a ; vengono visualizzati i valori di .)Y Y t + s Y t s = 0 , 1 , 2YYYt+sYts=0,1,2

Tuttavia, scegliendo , abbiamoα=β=1/2

αXt+βYt=Xt

anche per ed altrimentit

αXt+βYt=0.

Poiché non è costante, ovviamente questi due espressioni hanno distribuzioni differenti per ogni e , dove la serie non è fermo. I colori nella prima figura evidenziano questa non stazionarietà in distinguendo i valori zero dal resto.t t + 1 ( X + Y ) / 2Xtt+1(X+Y)/2(X+Y)/2


1
L'indipendenza delle due serie storiche è ovviamente una condizione sufficiente. Ma non sarebbe sufficiente anche il requisito più debole di stazionarietà congiunta?
Dilip Sarwate,

1
Sì, è vero @Dilip. Grazie per quell'osservazione.
whuber

5

Considera il processo bidimensionale

wt=(xt,yt)

Se è strettamente stazionario, o in alternativa, se i processi e sono congiuntamente strettamente stazionario , allora un processo costituite da una funzione misurabile sarà anche strettamente stazionario.( y t ) f : = f ( x t , y t ) , f : R 2R(xt)(yt)f:=f(xt,yt),f:R2R

Nell'esempio di @ whuber abbiamo

wt=(xt,(1)txt)

Per esaminare se questo è strettamente stazionario, dobbiamo prima ottenere la sua distribuzione di probabilità. Supponiamo che le variabili siano assolutamente continue. Per alcuni , abbiamo c RwtcR

Prob(Xtc,(1)tXtc)={Prob(Xtc,Xtc)t is evenProb(Xtc,Xtc)t is odd

={Prob(Xtc)t is evenProb(cXtc)t is odd

Prob(Xtc,(1)tXtc)={Prob(Xtc)t is evenProb(|Xt|c)t is odd

Seguendo l'esempio di whuber, i due rami sono distribuzioni di probabilità diverse perché ha una distribuzione simmetrica attorno allo zero. xt

Ora per esaminare la rigidità stazionaria, sposta l'indice di un numero intero . abbiamok>0

Prob(Xt+kc,(1)tXt+kc)={Prob(Xt+kc)t+k is evenProb(|Xt+k|c)t+k is odd

Per rigorosa stazionarietà, dobbiamo avere

Prob(Xtc,(1)tXtc)=Prob(Xt+kc,(1)tXt+kc),t,k

E non abbiamo questa uguaglianza , perché, diciamo, se è pari e è dispari, allora è dispari, nel qual casot,ktkt+k

Prob(Xtc,(1)tXtc)=Prob(Xtc)

mentre

Prob(Xt+kc,(1)tXt+kc)=Prob(|Xt+k|c)=Prob(|Xt|c)

Quindi non abbiamo una stretta comune , e quindi non abbiamo garanzie su cosa accadrà a una funzione di . f(xt,yt)

Devo sottolineare che la dipendenza tra e è una condizione necessaria ma non sufficiente per la perdita di una rigorosa stazionarietà. È il presupposto aggiuntivo della dipendenza di dall'indice che fa il lavoro.xtytyt

Ritenere

qt=(xt,θxt),θR

Se si esegue il lavoro precedente per si troverà qui la stretta stazionarietà rigorosa.(qt)

Questa è una buona notizia perché per un processo dipendere dall'indice ed essere rigorosamente stazionari non è tra i presupposti di modellizzazione che dobbiamo fare molto spesso. In pratica, quindi, se abbiamo una stazionarietà rigorosa marginale, ci aspettiamo anche una stazionarietà rigorosa anche in presenza di dipendenza (anche se ovviamente dovremmo verificare).


2

Direi di sì, dal momento che ha una rappresentazione MA.

Un'osservazione Penso che avere una rappresentazione MA implichi una debole stazionarietà, non sono sicuro che implichi una forte stazionarietà.


1
Ri "Non posso immaginare": vedi la mia risposta per un controesempio.
whuber

oneloop, rimuovere la parte relativa alla stazionarietà rigorosa e lasciare quella relativa alla stazionarietà debole. Ti darò un +1, poiché mi ha anche aiutato. ;)
Un vecchio nel mare.

@Anoldmaninthesea. Come questo?
oneloop

si Così. La rappresentazione MA implica una debole stazionarietà, anzi.
Un vecchio nel mare.

1
Questo viene automaticamente contrassegnato come di bassa qualità, probabilmente perché è così corto. Al momento è più un commento che una risposta dai nostri standard. Puoi espanderci? Puoi anche trasformarlo in un commento.
gung - Ripristina Monica
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.