【R言語】ARモデルの作成と検定

スポンサーリンク

 

 

今回はRでARモデルを作って検証していきます。

 

 

 一応ARモデルについておさらいしておくと、自己回帰モデルとはその名の通り、過去の自分のデータと回帰する分析手法のことを意味します。

 

 

これは色々ある時系列モデルの中で一番理屈が簡単なモデルで他のモデルの基礎にもなっています。

 

 

そして、1期前のデータだけから回帰したARモデルをAR(1)と表します。同様にp期前までの全てのデータを使って今日のデータを予測するモデルはAR(p)モデルと書きます。

 

 

・使うデータ

2005年から2016年までの任天堂の株価収益率のデータ

(※今回はr.nintenというオブジェクトにしています)

 

 

いざRで時系列モデルを検証していくわけですが、Rでは、ar( オブジェクト, aic=TRUE )で勝手にモデリングしてくれます。

 

>ar( r.ninten , aic=TRUE )      ##データから時系列モデルを検証する

 

aic=TRUE がデフォルトなので、AICで自動で最適なモデル選択してくれます。また次数は、Rであれば数値的に最適な次数を(AIC、Akaike’s Information Criteria等の比較によって)自動で判定してくれます。

 

 

ただ全部検証しているわけではないらしいので、一応当てはまりそうなモデルは自分で検証しておくというのもありかもしれません。

 

 

ただ今回は自己相関係数の記事でこのデータは1日前までの数値としか相関関係が保証されていないので次数は1、つまりAR(1)で計算します。

 

 

またaic=TRUEだと勝手に次数が決定されてしまうので、引数にaic=Fを加えAICを無効にしておく必要があります。

 

 

>ar(r.ninten, aic = F, order.max = 1)    ##AR(1)の算出

 

f:id:oruka199665:20170329022003j:plain

 

 

Coefficientsが次数で0.0522、そしてsigma^2は分散で6.956となりました。平均は出力されないので改めて計算する必要があります。

 

 

>(1-0.0522)*mean(ts.ninten)  ##平均の計算

[1] 0.02470081

 

 

これによって任天堂の株価収益率は下のような式でモデル化することができます。

 

 

rt= 0.02470081+0.0052×rt-1+et

 

 

これでモデル式は一応算出されましたが、これが本当に任天堂の株価収益率をちゃんと説明できているのかはまだ分かりません。

 

 

次は本当にこのモデルが正確なのかを検証するプロセスである残差分析をやっていきます。

 

 

www.tkstock.site

 

 

おまけ

 

・AIC による最適な AR モデルの選択

f:id:oruka199665:20170329023755j:plain

次数:8

モデル式: rt = 0.0554 rt-1-0.0019 rt-2 -0.0085 rt-3 +・・-0.0363 rt-8  + εt

誤差項の分散:6.929

 

 

ちなみに、時系列分析については、この本がとても分かりやすかったので、もし時系列分析でつまずいている人は是非一度読んでみてください。

 

プログラミングを挫折せず最短で習得する方法


プログラミングは小学校の義務教育にも導入され始めており、これから社会人として生きていく上にあたってはもはや出来て当たり前のスキルとなりつつあります。こういう話はテレビや雑誌でも最近増えてきたので、知らない人はいないでしょう。そしてそういう話を聞いて参考書なりを購入して独学でプログラミング勉強しようと思っている人も少なくないでしょう。しかしプログラミングを独学で勉強し始めようと思うものの



・「分からない箇所で詰まって挫折してしまった」

・「勉強する時間が足りない」

・「ネットの記事だと情報が断片的でよくわからない」

・「コードのエラーの原因が分からない」



という壁にぶち当たって、プログラミングの勉強を止めてしまう方が少なくありません。独学でプログラミングを勉強してる時間のほとんどはつまづいている時間です。実際僕も最初のころ独学でプログラミングを勉強していた頃はエラーの原因が分からず丸1日を不意にしてしまった・・・そんな苦い経験がありました。



それで僕は一度はプログラミングの学習を諦めてしまいましたが、就活で現実を知る中で「プログラミングを勉強して、いずれフリーランスとして自由な生き方がしたい」「エンジニアとして若いうちから高収入を得たい」という気持ちから一念発起して侍エンジニアのwebサービスコースに申し込み、プロのエンジニアの方に対面でマンツーマンでPythonによるWebサービス作り方とWeb技術の基本を教えてもらい、ようやくプログラミングが理解できました。


侍エンジニアでは様々なコースが用意されており、HTML・cssなどでWebサイトやWebアプリを作ってみたい人からPythonを勉強してデータサイエンティストやAIエンジニアになりたい人まで幅広いニーズにこたえています。


IT業界と言ってもエンジニアの仕事はプログラミング言語次第でサーバーから機械学習まで多種多様ですし、侍エンジニアの無料レッスン(カウンセリング)を受けてみて、自分のやりたいITの仕事は何なのか?を見つけてみましょう。ちなみに今「侍エンジニア」の無料レッスンを受けると1000円分のAmazonギフト券がもらえるので、試しに受けてみるだけもお得です。







自分は半端に独学やオンラインスクールで勉強して金と時間を無駄にするくらいなら、リアルのプログラミングスクールに通ってしっかりプログラミングを勉強した方がいいと思います。ちなみに今、侍エンジニアに申し込むと、25歳以下の学生の方であれば、受講料が20%OFFになるので超お得です。


そして、プログラミングは大勢で授業を受けたり漫然とオンライン学習をするよりも自分が分からない箇所をピンポイントでプロの講師に直接質問して、ちゃんと納得するというスタイルの方がお金は確かに少し掛かりますが、独学で学ぶよりも絶対にモノになります。
スポンサーリンク

シェアする

  • このエントリーをはてなブックマークに追加

フォローする