Bir modelin AIC'sini ve log dönüşümlü versiyonunu karşılaştırma


17

Sorumun özü şudur:

Let YRn ortalama bir çok değişkenli normal rastgele değişken olarak μ ve kovaryans matrisi Σ . Let Z:=log(Y) , yani Zi=log(Yi),i{1,,n} . Nasıl bir modelin uygunluğunun AIC gözlenen gerçekleşmelerine kıyasla nasıl bir Y gözlemlenen gerçekleşmelerine bir model uyum karşı Z ?



İlk ve biraz daha uzun sorum:

Let bir çok değişkenli normal rastgele değişken. ile uygun olan bir modeli uyumlu olan bir modelle karşılaştırmak istersem onların günlük olasılıklarına bakabilirim. Ancak, bu modeller iç içe olmadığından, günlük olasılıklarını (ve AIC vb. Gibi) doğrudan karşılaştıramıyorum, ancak bunları dönüştürmem gerekiyor.YN(μ,Σ)Ylog(Y)

ortak pdf ile rastgele değişkenler olduğunu biliyorumX1,,Xn ve eğer Y i = t ı ( x 1 , ... , x , n ) bire bir dönüşümler için t i ve i { 1 , ... , n } , daha sonra pdf Y 1 , ... , Y, n, ile verilmektedir fg(x1,,xn)Yben=tben(X1,...,Xn)tbenben{1,...,n}Y1,...,Yn burada J , dönüşümle ilişkili Jacobian'dır.

f(y1,...,yn)=g(t1-1(y),...,tn-1(y))det(J)
J

Karşılaştırmak için dönüşüm kuralını kullanmam gerekiyor mu?

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

l(Y)=günlük(Πben=1nφ(yben;μ,Σ))
l(günlük(Y))=günlük(Πben=1nφ(günlük(yben);μ,Σ))

veya yapabileceğim başka bir şey var mı?


Son iki ifadeye logaritma yerleştirmeyi unuttum.


Son ifadede de Jacobian'ı kaybettiniz.
whuber

2
günlükgünlükYY

Yanıtlar:


6

YZ

Araştırmacılar, tüm hipotezlerin aynı yanıt değişkeni kullanılarak modellenmiş olduğundan emin olmalıdır (örneğin, tüm model seti log (y) 'ye dayanıyorsa, hiçbir sorun yaratılmayacaktır; yanlış olan cevap değişkenlerinin karıştırılmasıdır).

Bu arada, AIC veya BIC ölçütlerini kullanmak için, modellerinizin mutlaka iç içe yerleştirilmesi gerekmez (aynı ref, sayfa 88, bölüm 2.12.4 Sınırsız Modeller) ve aslında BIC kullanmanın avantajlarından biri budur.


5

Akaike (1978, s. 224) model karşılaştırmasını sağlamak için AIC'nin dönüştürülmüş bir sonuç değişkeni varlığında nasıl ayarlanabileceğini açıklar. Şöyle diyor: “Değişkeni dönüştürmenin etkisi, karşılık gelen Jacobian'ın AIC'ye olasılığının çarpımı ile temsil edilir ...lÖg{y(n)+1}, −2 ΣlÖg{y(n)+1}, toplamın n=1,2,...,N-.”

Akaike, H. 1978. "Bir zaman serisi modeli olma olasılığı üzerine," Kraliyet İstatistik Kurumu Dergisi, D Serisi (İstatistikçi), 27 (3/4), s. 217-235.


1
bunu yapmak için R'de bir yaklaşım var mı?
theforestecologist
Sitemizi kullandığınızda şunları okuyup anladığınızı kabul etmiş olursunuz: Çerez Politikası ve Gizlilik Politikası.
Licensed under cc by-sa 3.0 with attribution required.