2013-01-10 5 views
6

से आत्मविश्वास अंतराल मूल्य निकालें आर में, हम समय श्रृंखला के एक एसीएफ correlogram चला सकते हैं और आत्मविश्वास अंतराल बैंड हल्के नीले रंग में प्लॉट किया जाएगा। लेकिन जब मैं एसीएफ ऑब्जेक्ट की संरचना खींचता हूं, तो मुझे ये मान नहीं मिलते हैं। क्या किसी को विश्वास है कि आत्मविश्वास अंतराल बैंड के मूल्यों को कैसे निकाला जाए?एसीएफ correlogram

उदा।

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

के बारे में अधिक जानकारी आप अपने एसीएफ कैसे प्लॉट करते हैं? आप किस पैकेज पर मुकदमा कर रहे हैं? – agstudy

+0

मैं इसे आँकड़े पैकेज में बनाया गया है लगता है। http://stat.ethz.ch/R-manual/R-patched/library/stats/html/acf.html – pat

उत्तर

7

मैं समारोह पर एक नज़र किया है और मैं एक आसान तरीका विश्वास अंतराल को निकालने के लिए नहीं देख सकता। क्षेत्र की गणना plot.acf फ़ंक्शन में की जाती है। इस समारोह देखने के लिए,

getS3method("plot", "acf") 

का उपयोग इस समारोह में, वहाँ एक चर clim है, यह एक के बाद हो रहा है। सबसे आसान तरीका plot.acf को myplot.acf पर कॉपी करना है, लेकिन clim मान वापस करें।

+0

इसके लिए धन्यवाद। यह एक काम की भी आसान की तरह नहीं दिखता, लेकिन अगर मैं किसी अन्य तरीके से नहीं मिल सकता है ... मैं कोड को संशोधित करने की कोशिश करनी होगी। वास्तव में संशोधित करने के तरीके पर कोई भी त्वरित पोस्ट की सराहना की जाएगी। – pat

+1

@pat बस 'return' बयान टिप्पणी और बदले वापसी' clim' – csgillespie

+0

इस जवाब के रूप में स्वीकार किया जाना चाहिए। 'डंप (" plot.acf ", फ़ाइल =" function_dump_plot_acf.txt ", माहौल = asNamespace (" आँकड़े "))': यहाँ समारोह को खींचने के लिए कोड है। इसके अलावा 'अदृश्य() 'को' रिटर्न (क्लाइम) 'के साथ बदलने की आवश्यकता है – user1700890

1

मैं जानता हूँ कि इस सवाल का सुपर पुराना है, लेकिन अगर किसी को चाहते हैं करता है विश्वास अंतराल में इसका इस्तेमाल किया टिप्पणियों की संख्या के sqrt से विभाजित आत्मविश्वास का स्तर का सिर्फ z-मूल्य है महत्व देता है। plot.acf समारोह में यह यहाँ की जाती है:

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

जहां with.ci किसी तार्किक मान दर्शाते हैं कि उपयोगकर्ता या विश्वास के अंतराल साजिश करना चाहता है नहीं और ci वांछित आत्मविश्वास का स्तर है (उदाहरण के लिए .95, .9, आदि है ...)

संपादित करें: यह विश्वास अंतराल है अगर आप, सफेद शोर मान लेग्ड मान हैं आप

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

आवेदन कर सकते हैं आप कुछ पढ़ सकते हैं कि कहीं ऐसा तो एक सुधार है नहीं है उस here