2001 年度エコノメトリックスC.

XLISP-Stat による方法

経済学部の学生用計算機室にある PC では Excel や TSP が利用できるが、自分の PC でこれらのソフトを利用することは難しいかもしれない。そこで、フリーソフトで定評のある統計ソフトを紹介しよう。
XLISP-STAT は、ミネソタ大学の Luke Tierney 教授が開発した統計ソフトで、Lisp1というコンピュータ言語に統計処理機能を実装したものである。そのため、統計処理をしようとするには、TSP と同様にコード（プログラムと呼ぶべきか）を書かなければならない。
当然、プログラムの文法は Lisp に準拠している。また、Windows だけでなく、MacOS やUNIX など幅広いプラットフォームで利用することができる。
①XLISP-Stat の入手法
(a) 自分の PC の OS が Winodws であれば、参考文献に挙げた、垂水(1999)の付録 FDを使う。
(b) インターネットからダウンロードする。ダウンロード時間は Windows 版の場合、約１MB あるから、36.6K モデムで約 5 分程度。の２つがある。インターネットを使える環境にない人は(a)を勧める。(b)のダウンロード先はミネソタ大学（ftp://ftp.stat.umn.edu/pub/xlispstat/pub/xlispstat/）から行なうが、まず北海道大学工学部の山本義郎先生のホームページ(URL http://isw.main.eng.hokudai.ac.jp/~yama/lispstat/)でダウンロードとインストールの詳細についての日本語解説を読んでからの方が、間違いが少なくていいであろう。（注意：2001 年 7 月 10 日現在、上記 URL からのリンクが上手くいっていないようであるので、その場合は、山本先生の昔のホームページ(URL http://www.f7.ems.okayama-u.ac.jp/~yama/lispstat/)）からアクセスした方がよいかもしれない。
②XLISP-Stat のインストール（Windows95/98 の場合）
(a)

配布先から「WXLS32ZP.EXE」という名前のファイルを Xlisp-Stat をインストールするフォルダ（例えば「c:¥lispstat」）ダウンロードする。

(b)

スタート→プログラム(P)→MS-DOS プロンプトをクリックする。そののち、(a)のファイルが保存されているフォルダ（ディレクトリ）に移動する。
c:¥windows>c:¥windows> cd c:¥lispstatc:¥lispstat> wxls32zp –dと入力し、圧縮されているプログラムを展開する。

(c)1

インストールが完了したら、続けて設定を行なう。
c:¥lispstat> c:¥lispstat¥wxls32.exe config

UNIX に詳しい人は、Emacs エディタが Lisp で書かれていることを知っているであろう。

7/10/01©Y. Takeuchi, 2000, 2001.

C-1

2001 年度エコノメトリックスを実行する。
（「c:¥lispstat」は Xlisp-Stat をインストールしたフォルダ名）。
あとは設定用プログラムに従って、Configure and Quit を押す。（次図参照）
その後、スタートメニューに Xlisp-Stat が登録される。
（参考）XLisp-Stat の設定画面③Xlisp-Stat プログラムの書き方
プログラムの書き方の詳細については、垂水(1999)を参照のこと。また、前記の北海道大学山本義郎先生のホームページにも簡単な説明がある。
B の TSP のところで使ったデータの Xlisp-Stat プログラムは以下の通り。

7/10/01©Y. Takeuchi, 2000, 2001.

C-2

2001 年度エコノメトリックス④XLisp-Stat プログラムの文法
XLisp-Stat は、インタプリタ（逐次翻訳実行）処理系であるから、プロンプト（入力要求）の>に対してプログラムを入力すると、すぐにその内容が実行される。
Lisp という（コンピュータ）言語では、プログラムの 1 文は必ず括弧( )の中に記述されなければならない。すなわち( …記述… ) の形をとる。この点は注意する必要がある。

イ) データの入力
データに変数名をつけて入力するには、(def 変数名 (list データ … ))というように指定する。C-3 の例のばあい、> (def temp (list 27.1 … 中略 … 27.2))
で、（27.1, …… , 27.2）という 10 個のデータに temp という名前を付けていることになる。
ロ) 変数変換、新しい変数を定義するには既に入力した変数を使って新しい変数を定義（変数変換）するには、次のように指定する。 (def 新変数名

(数式))

。ただ、気をつけなくてはならないのは、Lisp では演算記号の置く位置が、通常の数式記法と異なる前置記法前置記法を採用していることである。下表
前置記法
は、通常の数式記法と Lisp での記法を対照したものである。
Lisp での記法(+ 5 1)(* 3 2 2)(/ 6 4)(- 8 3 4)

通常の数式記法
5+1
3×2×2
6÷4
8-3-4
例えば、 x, z から y = 0.5 +log xという式を作るには、(3 − z 2 )(def(log x) (-y(+ 0.5 (/3 (^z2)))))のように記述する必要がある。加減乗除の記号(+-*/)の他、べき乗(^)、自然対数 log、平方根 sqrt などの演算子・関数がある。
ハ) 基本統計量(mean

平均変数名)(median

中央値（メディアン）

変数名)(standard-deviation変数名)

標準偏差ニ) グラフ作成(plot-points(histogram

x 軸にとる変数名 y 軸にとる変数名) 散布図の作成

変数)

7/10/01©Y. Takeuchi, 2000, 2001.

ヒストグラムの作成
C-3

2001 年度エコノメトリックスホ) 回帰分析
回帰分析を行う関数 regression-model を用いる。
単回帰モデル重回帰モデル(regression-model(regression-model

説明変数 被説明変数 オプション)(list 説明変数 … ) 被説明変数オプション)

オプションには次のものがある。
:intercepet nil:weightww

定数項なしのモデル

重みつき最小 2 乗法（「重み」は”ww”という変数で与えている）

回帰分析の結果は C-2 ページの下の図の通り。図中、括弧()内は標準誤差（係数推定量の標準偏差）を表している。従って、係数推定値を標準誤差で割ったものが、係数=0 という帰無仮説を検定するためのｔ統計量になる。また、”R Squares”が決定係数 R 2 、”Sigmahat”が撹乱項の分散推定値の平方根 σ̂ をそれぞれ表している。
ヘ) その他
現在使われている変数（名）を確認するには、 (variables) と指定すればよい。
【参考文献】
垂水共之(1999)、『Lisp-Stat による統計解析入門』、共立出版。
Lierney, L., 垂水他訳 (1996)、『Lisp-Stat』、共立出版。

7/10/01©Y. Takeuchi, 2000, 2001.

C-4