• 検索結果がありません。

7 ( 7 ( Workfile Excel hatuden 1000kWh kion_average kion_max kion_min date holiday *1 obon 7.1 Workfile 1. Workfile File - New -

N/A
N/A
Protected

Academic year: 2021

シェア "7 ( 7 ( Workfile Excel hatuden 1000kWh kion_average kion_max kion_min date holiday *1 obon 7.1 Workfile 1. Workfile File - New -"

Copied!
19
0
0

読み込み中.... (全文を見る)

全文

(1)

1

EViews

の使い方 第

4

小西葉子

 伊藤有希

初版

2007

7

4

目次

7 基本的な使い方(時系列データ) 2 7.1 Workfileの作成 . . . 2 7.2 データの読み込み . . . 4 8 ダミー変数(休日効果) 6 8.1 単回帰 . . . 6 8.2 回帰直線 . . . 8 8.3 土曜ダミーと日曜ダミーの作成 . . . 10 8.4 休日を考慮したモデル . . . 12 9 系列相関 14 9.1 誤差項に系列相関があるかどうか検定. . . 14 9.2 AR(1)モデル コクラン-オーカット法 . . . 16 9.3 ラグ付き内生変数 . . . 18 一橋大学経済研究所 [email protected] † 一橋大学経済学研究科博士課程 2年 [email protected]

(2)

7 基本的な使い方(時系列データ)

7

基本的な使い方

(

時系列データ)

本節では時系列データを分析する際のWorkfileの作成方法とデータの読み込み方法について述べる。用い るデータは以下の通り、2003年7月1日から2003年9月19日までの日次データである。 表1 時系列データの系列の意味 Excelの系列名 意味 単位 hatuden 東京電力の日次の発電量 1000kWh kion_average 東京の日次の平均気温(気象庁方式) ℃ kion_max 東京の日次の最高気温 ℃ kion_min 東京の日次の最低気温 ℃ date 日付番号 holiday 国民の休日*1のダミー変数 obon お盆期間のダミー変数

7.1

Workfile

の作成

1. Workfileを作る。「File」-「New」-「Workfile」をクリック

(3)

7 基本的な使い方(時系列データ) 7.1 Workfileの作成

2. 日次の時系列データ(例:2003年7月1月から2003年9月19日まで)のWorkfileを作る。以下のよ うに入力する。

・Workfile structure type Dated-regular frequency(デフォルト設定) ・Frequency Daily - 7 day week

・Start date 2003/7/1

・End data 2003/9/19

3. 以下の画面が出れば OK。Workfileの上側の表示が Range:7/01/2003 9/19/2003、Sample:7/01/2003

(4)

7 基本的な使い方(時系列データ) 7.2 データの読み込み

7.2

データの読み込み

7.2.1 コピー&ペースト方式

1. Workfile画面で「Quick」-「Empty Group (Edit Series)」をクリック

2. Excelデータファイルを開き、読み込みたいセルを範囲指定してコピー(例:C1からI82まで)。この

際日付は選択する必要はない。*2

3. EViews のワークシートに貼り付ける。Excel の 系列名をコピーしている場合は obs の右のセル

(7/01/2003の右上のセル)を右クリックして「Paste」をクリックする。

(5)

7 基本的な使い方(時系列データ) 7.2 データの読み込み

7.2.2 直接読み込む

1. Workfile画面で「Proc」-「Import」-「Read Text-Lotus-Excel」をクリックし、読み込みたい任意のExcel

ファイルを指定(例:denki_789.xls)すると以下のような画面が出る。

2. 今回はExcelファイルの系列hatuden, kion_average, kion_max, kion_min, date, holiday, obonの7つを 読み込みたいので以下のように入力する。

・Data order By Observation - series in columns*3

・Upper-left data cell C2*4

・Name of series… 7

・Excel 5+ sheet name 空欄のまま*5

*3系列を縦方向に読み込むという指定。

*4データの読み込みをはじめるセル (この例では Excel の spreadsheet で数字が入っているセルの中で一番左上) の場所を指定する。1

行目の系列名は自動的に読み込んでくれる。

(6)

8 ダミー変数(休日効果)

8

ダミー変数(休日効果)

8.1

単回帰

被説明変数に hatuden、説明変数に kion_averageを用いて以下のモデルで回帰分析を行う。

hatudent= c +βkion_averagett (1)

1. OLSを行う。Workfileを開いた状態で、「Quick」-「Estimate Equation」をクリックする。以下のよう に入力しOKをクリックする。その他の欄はデフォルトのままとする。

・Equation specification hatuden c kion_average

(7)

8 ダミー変数(休日効果) 8.1 単回帰

3. 推定式に名前をつける。推定結果を開いている状態で「Name」をクリックし、任意の名前(例:eq01) を入力しOKをクリック。Workfile上に eq01というObjectが表示される。

4. 残差の状態とfitをグラフで見てみる。推定結果を開いた状態で「View」-「Actual, Fitted, Residual」

-「Actual, Fitted, Residual Graph」をクリックする。

5. 残差のJarque-Bera Test(正規性の検定)を行う。推定結果を開いた状態で「View」-「Residual Tests」

(8)

8 ダミー変数(休日効果) 8.2 回帰直線

8.2

回帰直線

1. 回帰直線を書いてみる。Workfile画面を開く。Ctrキーを押しながら hatudenと kion_averageを選 択し、右クリック-「Open」-「as Group」をクリックしてGroup化を行う。

2. グループ化されたSeriesのspreadsheetが表示される。

3. Groupを開いている状態で、「Name」をクリックし、任意のGroup名(例:group01)を付ける。Workfile

(9)

8 ダミー変数(休日効果) 8.2 回帰直線

4. Group(例: group01)を開いている状態で、「View」-「Graph」-「Scatter」-「Scatter with Regression」 をクリックする。ダイアログが出るがそのままOKをクリック。

(10)

8 ダミー変数(休日効果) 8.3 土曜ダミーと日曜ダミーの作成

8.3

土曜ダミーと日曜ダミーの作成

1. 土曜ダミーを作る。Workfileを開いている状態で「genr」ボタンをクリックし、以下のように記入す る。*6

・Enter equation sat=@weekday=6

2. 日曜ダミーを作る。Workfileを開いている状態で「genr」ボタンをクリックし、以下のように記入する。

・Enter equation sun=@weekday=7

*6@weekday=x は曜日を表す関数である。x は1から 7 の値をとり、以下のような対応である。1=月曜日、2=火曜日、3=水曜日、4=

(11)

8 ダミー変数(休日効果) 8.3 土曜ダミーと日曜ダミーの作成

(12)

8 ダミー変数(休日効果) 8.4 休日を考慮したモデル

8.4

休日を考慮したモデル

被説明変数に hatuden、説明変数に kion_average、 sat、 sun、 holiday、 obonを用いて以下の 回帰分析を行う。

hatudent= c +β1kion_averaget+β2satt+β3sunt+β4holidayt+β5obontt (2)

1. 上記のモデルで回帰分析を行う。

・Equation specification hatuden c kion_average sat sun holiday obon

(13)

8 ダミー変数(休日効果) 8.4 休日を考慮したモデル

3. 推定式に名前をつける。推定結果を開いている状態で「Name」をクリックし、任意の名前(例:eq02) を入力しOKをクリック。Workfile上に eq02というObjectが表示される。

4. 残差をグラフで見てみる。推定結果を開いた状態で「View」-「Actual, Fitted, Residual」-「Actual, Fitted, Residual Graph」をクリックする。

5. 残差のJarque-Bera Test(正規性の検定)を行う。推定結果を開いた状態で「View」-「Residual Tests」

(14)

9 系列相関

9

系列相関

9.1

誤差項に系列相関があるかどうか検定

9.1.1 Durbin-Watson Test 先ほどの、休日を考慮したモデルの回帰分析の推定結果を見るとDurbin-Watson統計量が1.29と2から離 れているので、誤差項に以下のような1階の系列相関が存在する可能性がある。

hatudent= c +β1kion_averaget+β2satt+β3sunt+β4holidayt+β5obont+ ut (3)

utut−1t,|ρ| < 1 (4) Durbin-Watson統計量を用いて有意水準5%で両側検定(H0:ρ= 0, H1:ρ, 0)を行う。説明変数(定数項 を除く)が5で標本数が81のときのDurbin-Watson統計量の上方の限界の分布に関する臨界値(d81,5,0.025U ) が約1.71、下方の情報の限界の分布に関する臨界値(d81,5,0.025L )が約1.45なので、Durbin-Watson統計量 < dL 81,5,0.025となり系列相関がないという帰無仮説H0が棄却され、系列相関が存在する可能性がある。*7 9.1.2 コレログラムの作成

1. コレログラムを作成する。 eq02を開いている状態で「View」-「Residual Tests」-「Correlogram -Q-statistics」をクリックする。ダイアログが出るが36のまま「OK」をクリックする。

2. 結果を見る。ACはAutocorrelation(自己相関係数)、PACはPartial Autocorrelation(偏自己相関係数)、

Q-StatはLjung-Box統計量、ProbはLjung-Box統計量に対するP値をそれぞれ表している。

(15)

9 系列相関 9.1 誤差項に系列相関があるかどうか検定

9.1.3 Breusch-Godfrey LM test

1. Breusch-GodfreyのLM testを行う。 eq02を開いている状態で「View」-「Residual Tests」-「Serial

Correlation LM Test…」をクリックする。ダイアログが出て、ラグの長さを尋ねられるので今回は2と

(16)

9 系列相関 9.2 AR(1)モデル コクラン-オーカット法

9.2

AR(1)

モデル コクラン

-

オーカット法

以下のような1階の自己回帰(autoregressive)モデル(AR(1)モデル)を考える。

hatudent= c +β1kion_averaget+β2satt+β3sunt+β4holidayt+β5obont+ ut (5)

utut−1t,|ρ| < 1 (6)

1. AR(1)モデルの推定を行う。Workfileを開いた状態で、「Quick」-「Estimate Equation」をクリックす る。以下のように入力しOKをクリックする。その他の欄はデフォルトのままとする。

・Equation specification hatuden c kion_average sat sun holiday obon ar(1) *8

2. 推定結果を見る。推定結果のAR(1)の行が(6)式のρ に関する推定結果である。

*8最 後 の AR(1) と い う 指 定 は 、コ ク ラ ン-オ ー カ ッ ト 法 で AR(1) モ デ ル を 推 定 す る と い う 指 定 で あ る 。HATUDEN C

KION_AVERAGE SAT SUN HOLIDAY OBON AR(1) AR(2) と書けば 2 階の自己回帰モデル(AR(2) モデル)を推定する指 定。

(17)

9 系列相関 9.2 AR(1)モデル コクラン-オーカット法

3. 残差をグラフで見てみる。推定結果を開いた状態で「View」-「Actual, Fitted, Residual」-「Actual, Fitted, Residual Graph」をクリックする。

4. 残差のJarque-Bera Test(正規性の検定)を行う。推定結果を開いた状態で「View」-「Residual Tests」

(18)

9 系列相関 9.3 ラグ付き内生変数

9.3

ラグ付き内生変数

以下のようなラグ付き内生変数が含まれるモデルを考える。

hatudent= c +β1kion_averaget+β2satt+β3sunt+β4holidayt+β5obont+β6hatudent−1+ ut (7)

1. ラグ付き内生変数が含まれるモデルの推定を行う。Workfileを開いた状態で、「Quick」-「Estimate

Equation」をクリックする。以下のように入力しOKをクリックする。その他の欄はデフォルトのまま

とする。

・Equation specification hatuden c kion_average sat sun holiday obon hatuden(-1)*9

2. 推定結果を見る。

(19)

9 系列相関 9.3 ラグ付き内生変数

3. 残差をグラフで見てみる。推定結果を開いた状態で「View」-「Actual, Fitted, Residual」-「Actual, Fitted, Residual Graph」をクリックする。

4. 残差のJarque-Bera Test(正規性の検定)を行う。推定結果を開いた状態で「View」-「Residual Tests」

参照

関連したドキュメント

或はBifidobacteriumとして3)1つのnew genus

を塗っている。大粒の顔料の成分を SEM-EDS で調 査した結果、水銀 (Hg) と硫黄 (S) を検出したこと からみて水銀朱 (HgS)

を派遣しており、同任期終了後も継続して技術面での支援等を行う予定である。今年 7 月 30 日~8 月

タップします。 6通知設定が「ON」になっ ているのを確認して「た めしに実行する」ボタン をタップします。.

③  「ぽちゃん」の表記を、 「ぽっちゃん」と読んだ者が2 0名(「ぼちゃん」について何か記入 した者 7 4 名の内、 2 7

パキロビッドパックを処方入力の上、 F8特殊指示 →「(治)」 の列に 「1:する」 を入力して F9更新 を押下してください。.. 備考欄に「治」と登録されます。

等に出資を行っているか? ・株式の保有については、公開株式については5%以上、未公開株

・電源投入直後の MPIO は出力状態に設定されているため全ての S/PDIF 信号を入力する前に MPSEL レジスタで MPIO を入力状態に設定する必要がある。MPSEL