Membandingkan AIC dari model dan versi log-transformed-nya

17

Inti dari pertanyaan saya adalah ini:

Misalkan YRn menjadi variabel acak normal multivariat dengan mean μ dan matriks kovarian Σ . Misalkan Z:=log(Y) , yaitu Zi=log(Yi),i{1,,n} . Bagaimana cara membandingkan AIC model yang sesuai dengan realisasi yang diamati dari Y versus model yang sesuai dengan realisasi yang diamati dari Z ?



Pertanyaan awal dan sedikit lebih lama:

Misalkan menjadi variabel acak normal multivariat. Jika saya ingin membandingkan model yang cocok dengan Y dengan model yang cocok untuk log ( Y ) , saya bisa melihat kemungkinan log-nya. Namun, karena model ini tidak bersarang, saya tidak dapat membandingkan kemungkinan log (dan hal-hal seperti AIC, dll.) Secara langsung, tetapi saya harus mengubahnya.YN(μ,Σ)Ylog(Y)

Saya tahu bahwa jika adalah variabel acak dengan pdf g bersama ( x 1 , ... , x n ) dan jika Y i = t i ( X 1 , ... , X n ) untuk transformasi satu-ke-satu t i dan i { 1 , , n } , lalu pdf dari Y 1 , ,X1,,Xng(x1,,xn)Yi=ti(X1,,Xn)tii{1,,n} diberikan oleh f ( y 1 , ... , y n ) = g ( t - 1 1 ( y ) , ... , t - 1 n ( y ) ) det ( J ) di mana J adalah Jacobian yang terkait dengan transformasi.Y1,,Yn

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

Apakah saya hanya harus menggunakan aturan transformasi untuk membandingkan

ke l ( log (

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

atau ada hal lain yang bisa saya lakukan?


[Sunting] Lupa untuk menempatkan logaritma dalam dua ekspresi terakhir.

Stijn
sumber
Anda juga tampaknya telah kehilangan Jacobian dalam ekspresi terakhir.
whuber
2
Saya tidak mengerti loglogYY

Jawaban:

6

YZ

Peneliti harus yakin bahwa semua hipotesis dimodelkan menggunakan variabel respon yang sama (misalnya, jika seluruh rangkaian model didasarkan pada log (y), tidak akan ada masalah yang akan terjadi; itu adalah pencampuran variabel respon yang tidak benar).

Omong-omong, untuk menggunakan kriteria AIC atau BIC, model Anda tidak perlu bersarang (referensi yang sama, halaman 88, bagian 2.12.4 Model Tidak Berikat), dan sebenarnya itulah salah satu keunggulan menggunakan BIC.

Stat
sumber
5

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

Akaike, H. 1978. "Tentang kemungkinan model deret waktu," Jurnal Masyarakat Statistik Kerajaan, Seri D (Ahli Statistik), 27 (3/4), hlm. 217–235.

Gord B
sumber
1
apakah akan ada pendekatan dalam R untuk melakukan ini?
theforestecologist