【EViews】置換変数、コントロール変数、サブルーチン|プログラムで都道府県別に推計

大量な処理に必要なこと

EViewsは、エクセルのように、いろいろな変数を入れたり、推計方法を変えたりできるのが魅力です。

分析には試行錯誤が必要なので、その点では非常に便利なソフトです。

しかし、いったん分析手法が決まった後、データだけを変えることがありますが、推定を最初からやり直すのが面倒なことがあります。

また、47都道府県について、それぞれ推計する場合もあります。こういう場合はいちいちクリックして行うのは非効率です。

そういう場合に有用なのがプログラムです。この記事ではプログラムに必要な基本的な知識について解説します。

経済統計の使い方では、統計データの入手法から分析法まで解説しています。

https://officekaisuiyoku.com/category/software/eviews/

Contents

置換変数　%x
コントロール変数　!counter
サブルーチン　subroutine
実例

置換変数　%x

置換変数は文字に関する変数です。％をつけることでできます。%Xや％Aなどさまざまに作れますが、％Xを基本にしていると良いと思います。

％Xに文字列を代入するには以下の式を使います。文字列なので” “が必要です。

%X=”HK”

その変数を変数名などに使うときは、以下のように｛　｝で囲んで使います。

GDP｛％X｝

これはGDPHKを表します。

コントロール変数　!counter

コントロール変数は数字に関する変数です。X！やA！などどのような数字でもいいですが、!counterを基本にするといいと思います。!counterに１を入れるのは次の式です。

!counter=1

これを使うと、サブルーチン（後述）を実行するごとに数字を一つずつ増やすプログラムを書くことができます。

!counter=!counter+1

サブルーチン　subroutine

ルーチンはひとまとまりの処理のことです。サブルーチンはその一部です。プログラムの中で変数を変えて何度も繰り返すような作業がありますが、それをサブルーチンと呼びます。

基本的な構造は以下の通りです。

subroutine サブルーチン名
コマンド
endsub
call サブルーチン名

まず、サブルーチン名を宣言します。その後処理を書きます。処理の最後には、endsubを書きます。これで処理の内容が決まります。

その処理を呼び出すはcallというコマンドを使います。

実例

これらの仕組みを使ってプログラムを使った例を示します。

今回の例は、都道府県別の実質GDP（gdp{%X})を人流（職場）(wp{%X})に回帰した例です。

都道府県別のIDは、北海道がhk、青森がaoなど47都道府県別あります。

サブルーチンとして、人流（階差）を実質GDP（対数階差）に回帰します。

equation gdpelast{%x}.ls dlog(gdp{%x}) c d(wp{%x}/100)

この推計結果の係数と自由度修正済み決定係数を表(atab)に書き込むようにしています。

counterは都道府県が変わるごとに一つずつ増えていくので、行が増えていき、結果が記録されます。


%x="hk"
table atab 
!counter=1

subroutine gdpelast

equation  gdpelast{%x}.ls dlog(gdp{%x})  c d(wp{%x}/100)

atab(!counter,1)=gdpelast{%x}.@coef(2)
atab(!counter,2)=gdpelast{%x}.@rbar2

!counter=!counter+1

endsub

call gdpelast
%x="AO"
call gdpelast
%x="IW"
call gdpelast
%x="MG"
call gdpelast
%x="AT"
call gdpelast
%x="YG"
call gdpelast
%x="FS"
call gdpelast
%x="IG"
call gdpelast
%x="TG"
call gdpelast
%x="GM"
call gdpelast
%x="ST"
call gdpelast
%x="CB"
call gdpelast
%x="TK"
call gdpelast
%x="KG"
call gdpelast
%x="NG"
call gdpelast
%x="TY"
call gdpelast
%x="IK"
call gdpelast
%x="FI"
call gdpelast
%x="YN"
call gdpelast
%x="NN"
call gdpelast
%x="GF"
call gdpelast
%x="SO"
call gdpelast
%x="AC"
call gdpelast
%x="ME"
call gdpelast
%x="SI"
call gdpelast
%x="KT"
call gdpelast
%x="OS"
call gdpelast
%x="HG"
call gdpelast
%x="NR"
call gdpelast
%x="WY"
call gdpelast
%x="TT"
call gdpelast
%x="SN"
call gdpelast
%x="OY"
call gdpelast
%x="HS"
call gdpelast
%x="YI"
call gdpelast
%x="TS"
call gdpelast
%x="KW"
call gdpelast
%x="EM"
call gdpelast
%x="KC"
call gdpelast
%x="FO"
call gdpelast
%x="SG"
call gdpelast
%x="NS"
call gdpelast
%x="KM"
call gdpelast
%x="OT"
call gdpelast
%x="MZ"
call gdpelast
%x="KS"
call gdpelast
%x="ON"
call gdpelast

この結果として以下の表ができます。都道府県別に弾性値と自由度修正済み決定係数が記録されたものです。