Python Statsmodels的時(shí)間序列Ljung_Box檢驗(yàn)

Ljung-Box檢驗(yàn)即LB檢驗(yàn),是時(shí)間序列分析中檢驗(yàn)序列自相關(guān)性的方法。LB檢驗(yàn)的Q統(tǒng)計(jì)量為:


image.png

用來(lái)檢驗(yàn)m階滯后范圍內(nèi)序列的自相關(guān)性是否顯著,或序列是否為白噪聲,Q統(tǒng)計(jì)量服從自由度為m的卡方分布。

LB檢驗(yàn)可同時(shí)用于時(shí)間序列以及時(shí)序模型的殘差是否存在自相關(guān)性(是否為白噪聲)。Python的statsmodels包提供了該檢驗(yàn)的函數(shù):

from statsmodels.stats.diagnostic import acorr_ljungbox as lb_test

函數(shù)輸入 lb_test(x,lags=None,boxpierce=False):
x:檢驗(yàn)的時(shí)間序列

lags(int,list or None):檢驗(yàn)的延遲數(shù),

若為None則輸出min((nobs // 2 - 2), 40),其中nobs為觀測(cè)樣本數(shù)量,樣本較大的情況下輸出40

boxpierce:若為True,則同時(shí)輸出boxpierce統(tǒng)計(jì)量的檢驗(yàn)結(jié)果

(Box-Pierce檢驗(yàn)為白噪聲檢驗(yàn)的另一個(gè)版本,是LB檢驗(yàn)的前身)

函數(shù)輸出:

LB統(tǒng)計(jì)量值(array)

LB-p值(array)

若boxpierce=True,則繼續(xù)輸出BP統(tǒng)計(jì)量的值和相應(yīng)p值

作為樣例,選取了一段市場(chǎng)指數(shù)的回報(bào)率進(jìn)行檢驗(yàn),輸出p值序列:

plt.plot(lb_test(df_list[index_list[0]]['Return'])[1])
plt.show()
image.png

在延遲超過(guò)5時(shí),p值下降到0.05置信度以下,可以認(rèn)為出現(xiàn)顯著的自回歸關(guān)系,且序列并非白噪聲

?著作權(quán)歸作者所有,轉(zhuǎn)載或內(nèi)容合作請(qǐng)聯(lián)系作者
【社區(qū)內(nèi)容提示】社區(qū)部分內(nèi)容疑似由AI輔助生成,瀏覽時(shí)請(qǐng)結(jié)合常識(shí)與多方信息審慎甄別。
平臺(tái)聲明:文章內(nèi)容(如有圖片或視頻亦包括在內(nèi))由作者上傳并發(fā)布,文章內(nèi)容僅代表作者本人觀點(diǎn),簡(jiǎn)書系信息發(fā)布平臺(tái),僅提供信息存儲(chǔ)服務(wù)。

相關(guān)閱讀更多精彩內(nèi)容

  • 我這里補(bǔ)充一篇之前我做的時(shí)間序列部分的工作記錄,留個(gè)紀(jì)念怕忘。 這里用的數(shù)據(jù)說(shuō)我們公司自己的數(shù)據(jù),分析過(guò)程是隨著統(tǒng)...
    曦寶閱讀 13,292評(píng)論 2 2
  • 公式只能傳圖片?,好多公式都沒(méi)寫上,但是,不管有沒(méi)有公式,寫得一樣無(wú)趣至極,在中國(guó)人民大學(xué)出版社出版的王燕編著的<...
    今天芒種閱讀 13,072評(píng)論 0 10
  • from statsmodels.tsa.stattools import adfuller 可以看到里面的參數(shù)x...
    曦寶閱讀 83,052評(píng)論 5 15
  • 01 引言 上一篇推文【Python量化基礎(chǔ)】時(shí)間序列的自相關(guān)性與平穩(wěn)性著重介紹了時(shí)間序列的一些基礎(chǔ)概念,包括自相...
    CuteHand閱讀 5,584評(píng)論 0 7
  • 爸,咱今天回來(lái)老家了。 我們熱熱鬧鬧的送你走,您可一定要照顧好自己。 我們大放悲聲,也找不回您的回頭一顧,您是有多...
    微風(fēng)1991閱讀 182評(píng)論 0 0

友情鏈接更多精彩內(nèi)容