Confronto tra AIC di un modello e la sua versione trasformata in registro


17

L'essenza della mia domanda è questa:

Sia YRn una variabile casuale normale multivariata con media μ e matrice di covarianza Σ . Sia Z:=log(Y) , cioè Zi=log(Yi),i{1,,n} . Come confrontare l'AIC di un modello adatto alle realizzazioni osservate di Y rispetto a un modello adatto alle realizzazioni osservate di Z ?



La mia domanda iniziale e leggermente più lunga:

Sia YN(μ,Σ) una variabile casuale normale multivariata. Se voglio confrontare un modello adatto a Y rispetto a un modello adatto al log(Y) , potrei esaminare le loro probabilità logaritmiche. Tuttavia, poiché questi modelli non sono nidificati, non posso confrontare direttamente le verosimiglianze (e cose come AIC, ecc.), Ma devo trasformarle.

So che se sono variabili casuali con pdf congiunto g ( x 1 , ... , x n ) e se Y i = t i ( X 1 , ... , X n ) per le trasformazioni one-to-one t io e i { 1 , ... , n } , quindi il pdf di Y 1 , ... ,X1,,Xng(x1,,xn)Yi=ti(X1,,Xn)tii{1,,n} è dato da f ( y 1 , , y n ) = g ( t - 1 1 ( y ) , , t - 1 n ( y ) ) det ( J ) dove J è il giacobino associato alla trasformazione.Y1,,Yn

f(y1,,yn)=g(t11(y),,tn1(y))det(J)
J

Devo semplicemente usare la regola di trasformazione per confrontare

to l ( log ( Y ) ) = log ( n i = 1 ϕ ( log ( y i ) ; μ , Σ ) )

l(Y)=log(i=1nϕ(yi;μ,Σ))
l(log(Y))=log(i=1nϕ(log(yi);μ,Σ))

o c'è qualcos'altro che posso fare?


[modifica] Hai dimenticato di posizionare i logaritmi nelle ultime due espressioni.


Sembra anche che tu abbia perso il giacobino nell'ultima espressione.
whuber

2
loglogYY

Risposte:


6

YZ

Gli investigatori dovrebbero essere sicuri che tutte le ipotesi siano modellate utilizzando la stessa variabile di risposta (ad esempio, se l'intero insieme di modelli fosse basato su log (y), non si creerebbe alcun problema; è la miscelazione delle variabili di risposta che non è corretta).

Inoltre, per utilizzare i criteri AIC o BIC, i tuoi modelli non devono necessariamente essere nidificati (stesso riferimento, pagina 88, sezione 2.12.4 Modelli non annidati), e in realtà questo è uno dei vantaggi dell'utilizzo di BIC.


5

log{y(n)+1}log{y(n)+1}n=1,2,...,N

Akaike, H. 1978. "Sulla probabilità di un modello di serie storiche", Journal of the Royal Statistical Society, Serie D (The Statistician), 27 (3/4), pp. 217–235.


1
capita che ci sia un approccio in R per farlo?
forestecologo
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.