1. 程式人生 > >金融時間序列分析: 10. ARMA模型例項(R,Python)

金融時間序列分析: 10. ARMA模型例項(R,Python)

0. 目錄

1. 前言

建立ARMA模型的過程和前面提到的建立AR/MA模型基本一致,只是對ARMA模型的定階方法不一樣。

2. ARMA定階

2.1 R語言對ARMA模型定階

R語言的TSA包中包含EACF函式,可以用來ARMA定階
函式:eacf(data, max_p, max_q)

ths_pq = eacf(log_ret, 10, 10)

這個方法會直接輸出一個簡單的二維圖,選取左上角的“O”,其座標就是ARMA模型p,q階數,具體使用參考上一篇文章:
金融時間序列分析:9. ARMA自迴歸移動平均模型

結果如下圖所示:

這裡寫圖片描述

當然,也可以把具體的數值輸出

print(ths_pq$eacf)

這裡寫圖片描述

這些數字是對應ARMA模型的標準誤差 σ
XXOO的標準是:

|σ|<2/T 時,為O,
否者為X

其中T是資料大小。上面的資料量為1799,那麼閾值為:
2/1799=0.04715
用這個數值和上圖比較,你就能得到那張XXOO的圖了。

從數值看,應該選擇(p,q) = (5,2),但是這麼選的話模型的複雜度就太高了(7個引數),除此之外還有很多,等會我們對比下(0,2),(4,0),(1,2),(5,2)

再來看下ACF和PACF:
這裡寫圖片描述

這裡寫圖片描述

看到ACF和PACF我就有些蛋疼了………………
從ACF上看:AR(1)是一個比較合適的模型
從PACF上看:MA(1)是一個比較合適的模型

相關推薦

金融時間序列分析: 10. ARMA模型例項RPython

0. 目錄 1. 前言 建立ARMA模型的過程和前面提到的建立AR/MA模型基本一致,只是對ARMA模型的定階方法不一樣。 2. ARMA定階 2.1 R語言對ARMA模型定階 R語言的TSA包中包含EACF函式,可以用來ARM

2017.06.9 金融時間序列分析之Eview使用基礎

file 打開 操作 span bject 生成 com 工作表 時間序列分析 一.創建時間序列工作文件:首先將數據轉換為Eviews系統能夠分析的Eviews Workfile數據集 1.創建工作文件:工作文件結構類型:非結構/非日期型;日期-規則頻率型;平衡面板型;

金融時間序列分析:1. 基礎知識

0. 目錄 1. 金融時間序列 1.1什麼是時間序列 金融時間序列是屬於時間序列資料的一種,他們就是有很強的時間性,資料前後具有很強的依賴性,切無法調整順序,一般都是二維資料。 時間序列由於具有很強的序列行,而且資料前後一般存在依賴,週期等關

多層感知機MLP演算法原理及Spark MLlib呼叫例項Scala/Java/Python

多層感知機 演算法簡介:         多層感知機是基於反向人工神經網路(feedforwardartificial neural network)。多層感知機含有多層節點,每層節點與網路的下一層節點完全連線。輸入層的節點代表輸入資料,其他層的節點通過將輸入資料與層上節點

二十種特徵變換方法及Spark MLlib呼叫例項Scala/Java/python

Tokenizer(分詞器) 演算法介紹:         Tokenization將文字劃分為獨立個體(通常為單詞)。下面的例子展示瞭如何把句子劃分為單詞。         RegexTokenizer基於正則表示式提供更多的劃分選項。預設情況下,引數“pattern”為

MLlib--多層感知機MLP演算法原理及Spark MLlib呼叫例項Scala/Java/Python

來源:http://blog.csdn.net/liulingyuan6/article/details/53432429 多層感知機 演算法簡介:         多層感知機是基於反向人工神經網路(feedforwardartificial neural net

三種特徵選擇方法及Spark MLlib呼叫例項Scala/Java/python

VectorSlicer 演算法介紹:         VectorSlicer是一個轉換器輸入特徵向量,輸出原始特徵向量子集。VectorSlicer接收帶有特定索引的向量列,通過對這些索引的值進行篩選得到新的向量集。可接受如下兩種索引 1.整數索引,setIndice

隨機森林迴歸Random Forest演算法原理及Spark MLlib呼叫例項Scala/Java/python

隨機森林迴歸 演算法介紹:        隨機森林是決策樹的整合演算法。隨機森林包含多個決策樹來降低過擬合的風險。隨機森林同樣具有易解釋性、可處理類別特徵、易擴充套件到多分類問題、不需特徵縮放等性質。 隨機森林分別訓練一系列的決策樹,所以訓練過程是並行的。因演算法中加入隨機

二十種特徵變換方法及Spark MLlib呼叫例項Scala/Java/python

VectorIndexer 演算法介紹:         VectorIndexer解決資料集中的類別特徵Vector。它可以自動識別哪些特徵是類別型的,並且將原始值轉換為類別指標。它的處理流程如下: 1.獲得一個向量型別的輸入以及maxCategories引數。 2.基於

梯度迭代樹GBDT演算法原理及Spark MLlib呼叫例項Scala/Java/python

梯度迭代樹 演算法簡介:         梯度提升樹是一種決策樹的整合演算法。它通過反覆迭代訓練決策樹來最小化損失函式。決策樹類似,梯度提升樹具有可處理類別特徵、易擴充套件到多分類問題、不需特徵縮放等性質。Spark.ml通過使用現有decision tree工具來實現。

Spring Boot + Mybatis + 二級快取例項EhcacheRedis

使用Mybatis自帶二級快取 MyBatis 包含一個非常強大的查詢快取特性,它可以非常方便地配置和定製。MyBatis 3 中的快取實現的很多改進都已經實現了,使得它更加強大而且易於配置。 預設情況下是沒有開啟快取的,除了區域性的 session 快取,可以增強

ICMP協議及應用例項PingTraceroute

ICMP協議 在網際網路傳輸過程中,IP資料報難免會出現差錯,通常出現差錯,處理方法就是丟棄,但是一般,出現差錯後,會發送ICMP報文給主機,告訴它一些差錯資訊,以及對當前的網路狀態進行一個探尋。可以說,ICMP的主要目的是用於在TCP/IP網路中傳送出錯和控制

python時間序列分析ARIMA模型

原文地址:https://blog.csdn.net/u011596455/article/details/78650458 轉載請註明出處。 什麼是時間序列       時間序列簡單的說就是各時間點上形成的數值序列,時間序列分析就是通過觀察歷史資料預測未來的值。在這裡需要

ARIMA時間序列分析-----Python例項一週銷售營業額預測

以ARIMA模型為例介紹時間序列演算法在python中是如何實現的,一下是應用Python語言建模步驟: -- coding: utf-8 -- “”” Created on Mon Apr 2 16:45:36 2018 @author: hou

Python_Statsmodels包_時間序列分析_ARIMA模型

以前都是用SPSS/EVIEWS等來做的,現在用Python來搞搞看,前後弄了一個星期,總算基本走通了。 基礎庫: pandas,numpy,scipy,matplotlib,statsmodels : from __future__ import print_funct

統計思維例項11——時間序列分析

時間序列(time series)是來自隨時間變化的系統的一系列度量。本章使用的示例來自Zachary M. Jones。Jones的研究目的是調查像大麻合法化這樣的政策性決定會對市場產生何種影響。希望大家對本章內容感興趣,但藉此機會重申對資料分析保持專業性態度的重要性。藥品

時間序列分析——ARIMA模型預測(R)

讀取資料(scan()及read.csv()等函式,row.names引數可以用來指定索引列) 觀察檢驗時間序列是否平滑,對不平滑的時間序列要進行差分(diff()函式),差分的階數=arima(p,d,q)中d引數的值 acf()和 pacf()兩個函式分

R語言--時間序列分析步驟

align 如何 -- list arima test bsp nat 建立 大白。 (1)根據趨勢定差分 plot(lostjob,type="b") 查看圖像總體趨勢,確定如何差分 df1 = diff(lostjob) d=1階差分 s4_df1=diff(df1,

計量經濟與時間序列_時間序列分析的幾個基本概念(自相關函數,偏自相關函數等)

sig 永不 均值 blog 那種 屬於 class 觀察 自相關 1. 在時間序列分析中, 數學模型是什麽?數學公式又是什麽?數學推導過程又是什麽?... ...   一句話:用數學公式後者符號來表示現實存在的意義。數學是“萬金油”的科學,它是作為工作和分析方法運用到某

時間序列分析

數理統計 bsp 移動 分段 gpo 降雨量 基礎知識 -s data- 時間序列分析( time series analysis)方法,強調的是通過對一個區域進行一定時間段內的連續遙感觀測,提取圖像有關特征,並分析其變化過程與發展規模。當然,首先需要根據檢測對象的時相變化