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.


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.


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.


1
apakah akan ada pendekatan dalam R untuk melakukan ini?
theforestecologist
Dengan menggunakan situs kami, Anda mengakui telah membaca dan memahami Kebijakan Cookie dan Kebijakan Privasi kami.
Licensed under cc by-sa 3.0 with attribution required.