2013-01-10 13 views
6

korelasyon s, biz zaman serisi ve aralık bantları açık mavi çizilen olacak güven ACF korelogram çalıştırabilir. Ama ACF nesnesinin yapısını çektiğimde, bu değerleri bulamıyorum. Güven aralığı gruplarının değerlerini nasıl çıkaracağını bilen var mı?Özü güven aralığı değerleri, R ise

örn.

List of 6 
$ acf : num [1:27, 1, 1] 1 0.06453 -0.06354 0.00213 -0.01324 ... 
$ type : chr "correlation" 
$ n.used: int 501 
$ lag : num [1:27, 1, 1] 0 1 2 3 4 5 6 7 8 9 ... 
$ series: chr "tser[i:(i + 500)]" 
$ snames: NULL 
- attr(*, "class")= chr "acf" 

enter image description here

+0

hakkında daha? hangi paket dava açıyor? – agstudy

+0

Bence istatistik paketinin içine yerleştirilmiş. http://stat.ethz.ch/R-manual/R-patched/library/stats/html/acf.html – pat

cevap

7

Ben fonksiyonuna bir göz yaşadım ve ben güven aralığını ayıklamak için kolay bir yol göremiyorum. Bölge, plot.acf işlevinde hesaplanır. Bu işlevi görmek için bu işlevde

getS3method("plot", "acf") 

kullanın bu peşindesiniz biridir değişken clim yoktur. En kolay yol, plot.acf'u myplot.acf'a kopyalamaktır, ancak clim değerini döndürür.

+0

Bunun için teşekkürler. Bir görev için çok kolay görünmüyor, ama başka bir yol bulamazsam ... kodu değiştirmem gerekiyor. Tam olarak nasıl değiştirileceğine dair hızlı mesajlar takdir edilecektir. – pat

+1

@pat Sadece 'return' ifadelerini yorumlayın ve bunun yerine "clim" değerini döndürün – csgillespie

+0

Bu cevap olarak kabul edilmelidir. İşte fonksiyon kodu: 'dump (" plot.acf ", dosya =" function_dump_plot_acf.txt ", envir = asNamespace (" stats ")). Ayrıca '' invisible() 'ile' return (clim) 'yerine – user1700890

1

bu soru süper eski olduğunu biliyorum ama herkes istiyor yaparsa güven aralığı, kullanıldığı gözlem sayısı sqrt bölü güven seviyesinin sadece z değerdir değer vermektedir. plot.acf fonksiyonda bu burada hesaplanır: with.ci mantıksal değeri (kullanıcı güven aralıkları çizmek istiyor ya da değil ve ci istenen güven düzeyi ise örneğin 0,95, 0,9, vb belirten bir

clim0 <- if (with.ci) 
    qnorm((1 + ci)/2)/sqrt(x$n.used) 

...)

DÜZENLEME: Bu bir düzeltme var durum değilse, beyaz gürültü gecikmeli değerlerdir varsaydık.Yaklaşık bazı okuyabilirsiniz

clim <- clim0 * sqrt(cumsum(c(1, 2 * x$acf[-1, i, j]^2))) 

uygulayabilirsiniz, bu güven aralığıdır Eğer acf çizmek nasıl bu here