時間序列的ADF檢驗(yàn)(單位根檢驗(yàn))

from statsmodels.tsa.stattools import adfuller

我使用的是python的statsmodels哥攘。
1532680234(1).jpg

可以看到里面的參數(shù)
x:array_like裆赵,1維拥诡,時間序列

maxlag:int型,最大延遲階數(shù)

regression:Constant and trend order to include in regression艇潭,回歸中包含的常數(shù)和趨勢階數(shù)坑填。解釋一下:ADF檢驗(yàn)和DF檢驗(yàn)一樣可以用于三種類型的單位根檢驗(yàn)。(王燕虹蒋,應(yīng)用時間序列分析糜芳,6.3.2ADF檢驗(yàn)飒货,分了三種類型,代碼里給了四種類型)

    * 'c' : constant only (default)峭竣,默認(rèn)塘辅,僅有常數(shù)均值。
    * 'ct' : constant and trend皆撩,有常數(shù)均值扣墩,有趨勢。
    * 'ctt' : constant, and linear and quadratic trend扛吞,有常數(shù)均值有線性和二次趨勢呻惕。
    * 'nc' : no constant, no trend,無常數(shù)均值滥比,無趨勢亚脆。

當(dāng)選擇的回歸類型不同時,檢驗(yàn)結(jié)果也會不同盲泛。我在這里深深地踩了一個坑型酥。可以看一下時序圖查乒,你的數(shù)據(jù)是否有趨勢弥喉,是否有常數(shù)均值。
autolag : {'AIC', 'BIC', 't-stat', None}玛迄,自動選擇延遲階數(shù)
* if None, then maxlag lags are used由境,如果選擇none,則使用最大延遲階數(shù)
* if 'AIC' (default) or 'BIC', then the number of lags is chosen
to minimize the corresponding information criterion蓖议,如果選擇AIC或者BIC虏杰,則延遲階數(shù)是根據(jù)相應(yīng)的最小信息準(zhǔn)則決定的。(BIC就是SBC準(zhǔn)則)(王燕勒虾,應(yīng)用時間序列分析纺阔,3.3.6模型優(yōu)化)
* 't-stat' based choice of maxlag. Starts with maxlag and drops a
lag until the t-statistic on the last lag length is significant
using a 5%-sized test,這個選擇要基于最大延遲maxlag修然,從最大延遲階數(shù)開始笛钝,每次減少一個延遲,直到某一階延遲的t統(tǒng)計(jì)量對5%來說是顯著的愕宋。(也就是說p_value值小于0.05玻靡,顯著拒絕原假設(shè),證明無單位根中贝,也就是說序列是平穩(wěn)的)


image.png

這是我的數(shù)據(jù)的檢測結(jié)果囤捻。
第一部分是τ(tao)統(tǒng)計(jì)量的值。
第二部分是p_value的值邻寿。
第三部分是結(jié)果使用的延遲階數(shù)蝎土。
第四部分是ADF回歸和計(jì)算臨界值所使用的觀察次數(shù)视哑。
第五部分是臨界值。
第六部分是最大的信息準(zhǔn)則的值(如果autolag 非空)誊涯,也就是AIC或者BIC的值黎炉。
當(dāng)我們看序列是否平穩(wěn)的結(jié)果時,一般首先看第二部分的p_value值醋拧。如果p_value值比0.05小慷嗜,證明有單位根,也就是說序列平穩(wěn)丹壕。如果p_value比0.05大則證明非平穩(wěn)庆械。
源碼里有一句note,如果p_value接近于0.05時菌赖,則要通過臨界值進(jìn)行判斷缭乘。也就是說如果p_value接近于0.05就要通過第一部分τ(tao)統(tǒng)計(jì)量的值和第五部分的臨界值進(jìn)行對比。τ(tao)統(tǒng)計(jì)量的值比臨界值小琉用,就證明平穩(wěn)堕绩,反之就是非平穩(wěn)。我的檢測結(jié)果τ(tao)統(tǒng)計(jì)量的值再臨界值5%-10%之間邑时。比5%的臨界值大奴紧。這里的1%,5%晶丘,10%對應(yīng)的是99%黍氮,95%,90%置信區(qū)間浅浮。
ADF檢驗(yàn)只適合AR(P)模型沫浆。
且對方差齊性效果好,對異方差性效果不佳滚秩。異方差可用PP檢驗(yàn)专执。

最后編輯于
?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請聯(lián)系作者
  • 序言:七十年代末,一起剝皮案震驚了整個濱河市郁油,隨后出現(xiàn)的幾起案子本股,更是在濱河造成了極大的恐慌,老刑警劉巖已艰,帶你破解...
    沈念sama閱讀 206,482評論 6 481
  • 序言:濱河連續(xù)發(fā)生了三起死亡事件痊末,死亡現(xiàn)場離奇詭異,居然都是意外死亡哩掺,警方通過查閱死者的電腦和手機(jī),發(fā)現(xiàn)死者居然都...
    沈念sama閱讀 88,377評論 2 382
  • 文/潘曉璐 我一進(jìn)店門涩笤,熙熙樓的掌柜王于貴愁眉苦臉地迎上來嚼吞,“玉大人盒件,你說我怎么就攤上這事〔涨荩” “怎么了炒刁?”我有些...
    開封第一講書人閱讀 152,762評論 0 342
  • 文/不壞的土叔 我叫張陵,是天一觀的道長誊稚。 經(jīng)常有香客問我翔始,道長,這世上最難降的妖魔是什么里伯? 我笑而不...
    開封第一講書人閱讀 55,273評論 1 279
  • 正文 為了忘掉前任城瞎,我火速辦了婚禮,結(jié)果婚禮上疾瓮,老公的妹妹穿的比我還像新娘脖镀。我一直安慰自己,他們只是感情好狼电,可當(dāng)我...
    茶點(diǎn)故事閱讀 64,289評論 5 373
  • 文/花漫 我一把揭開白布蜒灰。 她就那樣靜靜地躺著,像睡著了一般肩碟。 火紅的嫁衣襯著肌膚如雪强窖。 梳的紋絲不亂的頭發(fā)上,一...
    開封第一講書人閱讀 49,046評論 1 285
  • 那天削祈,我揣著相機(jī)與錄音毕骡,去河邊找鬼。 笑死岩瘦,一個胖子當(dāng)著我的面吹牛未巫,可吹牛的內(nèi)容都是我干的。 我是一名探鬼主播启昧,決...
    沈念sama閱讀 38,351評論 3 400
  • 文/蒼蘭香墨 我猛地睜開眼叙凡,長吁一口氣:“原來是場噩夢啊……” “哼!你這毒婦竟也來了密末?” 一聲冷哼從身側(cè)響起握爷,我...
    開封第一講書人閱讀 36,988評論 0 259
  • 序言:老撾萬榮一對情侶失蹤,失蹤者是張志新(化名)和其女友劉穎严里,沒想到半個月后新啼,有當(dāng)?shù)厝嗽跇淞掷锇l(fā)現(xiàn)了一具尸體,經(jīng)...
    沈念sama閱讀 43,476評論 1 300
  • 正文 獨(dú)居荒郊野嶺守林人離奇死亡刹碾,尸身上長有42處帶血的膿包…… 初始之章·張勛 以下內(nèi)容為張勛視角 年9月15日...
    茶點(diǎn)故事閱讀 35,948評論 2 324
  • 正文 我和宋清朗相戀三年燥撞,在試婚紗的時候發(fā)現(xiàn)自己被綠了。 大學(xué)時的朋友給我發(fā)了我未婚夫和他白月光在一起吃飯的照片。...
    茶點(diǎn)故事閱讀 38,064評論 1 333
  • 序言:一個原本活蹦亂跳的男人離奇死亡物舒,死狀恐怖色洞,靈堂內(nèi)的尸體忽然破棺而出,到底是詐尸還是另有隱情冠胯,我是刑警寧澤火诸,帶...
    沈念sama閱讀 33,712評論 4 323
  • 正文 年R本政府宣布,位于F島的核電站荠察,受9級特大地震影響置蜀,放射性物質(zhì)發(fā)生泄漏。R本人自食惡果不足惜悉盆,卻給世界環(huán)境...
    茶點(diǎn)故事閱讀 39,261評論 3 307
  • 文/蒙蒙 一盯荤、第九天 我趴在偏房一處隱蔽的房頂上張望。 院中可真熱鬧舀瓢,春花似錦廷雅、人聲如沸。這莊子的主人今日做“春日...
    開封第一講書人閱讀 30,264評論 0 19
  • 文/蒼蘭香墨 我抬頭看了看天上的太陽。三九已至堰怨,卻和暖如春芥玉,著一層夾襖步出監(jiān)牢的瞬間,已是汗流浹背备图。 一陣腳步聲響...
    開封第一講書人閱讀 31,486評論 1 262
  • 我被黑心中介騙來泰國打工灿巧, 沒想到剛下飛機(jī)就差點(diǎn)兒被人妖公主榨干…… 1. 我叫王不留,地道東北人揽涮。 一個月前我還...
    沈念sama閱讀 45,511評論 2 354
  • 正文 我出身青樓抠藕,卻偏偏與公主長得像,于是被迫代替她去往敵國和親蒋困。 傳聞我的和親對象是個殘疾皇子盾似,可洞房花燭夜當(dāng)晚...
    茶點(diǎn)故事閱讀 42,802評論 2 345

推薦閱讀更多精彩內(nèi)容

  • 1 概念 ARIMA模型,全稱為自回歸積分滑動平均模型(Autoregressive Integrated ...
    風(fēng)逝流沙閱讀 41,138評論 1 48
  • 公式只能傳圖片?,好多公式都沒寫上,但是,不管有沒有公式,寫得一樣無趣至極,在中國人民大學(xué)出版社出版的王燕編著的<...
    今天芒種閱讀 12,597評論 0 10
  • 時間序列簡介 時間序列 是指將同一統(tǒng)計(jì)指標(biāo)的數(shù)值按其先后發(fā)生的時間順序排列而成的數(shù)列雪标。時間序列分析的主要目的是根據(jù)...
    熙淺閱讀 177,033評論 15 55
  • 本文結(jié)構(gòu): 時間序列分析零院? 什么是ARIMA? ARIMA數(shù)學(xué)模型村刨? input告抄,output 是什么? 怎么用嵌牺?...
    不會停的蝸牛閱讀 65,206評論 6 97
  • 銳眼視點(diǎn): 以色列初創(chuàng)企業(yè) Utilis 在2017年 IH2O 比賽中獲得第一打洼; Insurify 推出智能保險(xiǎn)...
    銳眼看世界閱讀 262評論 0 0