マ イ ク ロ コ ン ピ ュ ー タ に よ る統 計 解 析
―
地
方
自
治
体
の
社
会
福
祉
財
政
関
係
デ
ー
タ
を
用
い
て
―
坂
田
周
一
本資料 は,昭和56
年度文部省科学研究費補助 金 (奨励研 究 A)に よる資金援助 を受 けて現在続 行 中の研究 「地方 自治体 の社会福祉予算における 意思決定の計量的研究」 のための分析 ツール とし て作成 した ものの一部であ る。 (なお,本 資料 に掲載 した プログラムの無断転載 を禁 じます。) 日 次 は じめに 0 マイコソBASICの限界 Ⅰデータ ・ファイルの作成 (プログラム 1) ⅠⅠデータ と分析 プ ログラムとの リンク (プ ログ ラム2) IIⅠ分布統 計量 と度数分布 とグラフ (プログラム 3)Ⅰ
Ⅴ相関 マ トリクス とⅩ-
Y
プ ロッターと回帰分 析 (プ ログラム 4) Ⅴ質的 デ ータの ク ロス集 計 とカ イ2
乗 ・ク ラ マ-係数 (プ ログラム 5) おあ りに は じめ に マ イ ク ロ コン ピュータの普 及 に よって, コン ピュータに よる統計解析が個人 の机上 で 日常的 に 行 えるよ うになった。 これ は, データとプ ログラ ムを持参 して計算機 セ ンタ-に依頼す るバ ッチ処 理方式 の面倒 な手続 きか ら解放 され ることを意味 す るので,サ ンプル数が50,100,200とい った程 度 の大 きさのデータを 日常的に処理す る-必要 のあ る研究者に とって大 きな救 いであ る。 しか し, マ イコン利用の場合, ユーザ ーがみず か らプログ ラマ となって プログラムを組 まなけれ ばマイ コンは走 らない。 プログラム言語が障壁 に な ってマイコン利用が広が らないの も大 きな損失 であろ う。 そ こで, プ ログラム言語 について何の 知識 ももた ない人で もマイコンによる統計解析が で きるよ うにす るために, アプ リケーシ ョン ・フ ログラムを発表す ることにした。 なお, プログラ ムの実行 にあた ってほ,地方 自治体 の社会福祉財 政関係 の実際のデータを用いている。 アプ リケシ ョソ ・プ ログラムであ る以上,特定 の機種 と言語仕様 を予定 してい る0枚種 はシャー プ ・パ ー ソナル コンピュータMZl
8
0B
,言語仕様 は倍精度BASICインタプ リターSB-5620に基づ く。テープ・ベースのBASICであ ることをお こと わ りしてお く。 システムの拡張 は基本的 に必要 と しない よ うに各 プ ログラムを作成 した。 しか し、 第1
Ⅴ節でX-
Y
プロッタ-を作成す る関係か ら,ど うして もグ ラ フィックRAM
ページ1
を拡 張 す ることになった。 なお,本機種でユ-ザ ーが 自由 に使用で きる記憶容量 は6
4
キ ロ ・バ イ トであ る。0
マ イ コンBASI
C
の限界
プ ログ ラム言語 に は事 務処 理 向 のCOBOLと か科学計算向のFORTRAN
をは じめ とした各種 の言語があ る。BASICもその一つであ る。Begi n-ner'sAlトpurposeSymbolicInstructionCodeの 略称であるが, この言語 は事務処理 も科学計算 も で きる うえに入門者用 として習得 も容 易であ るこ と,結果の出力がFORTRAN
の よ うに面倒 で な い こと,実数 と整数 を区別 しな くて もよい ことな どのいろいろの利点 を もつ ことか ら,現在 ほ とん どのマイコンで使用で きる言語 である。 しか し,マイ コンで使用で きるBASICは,大型 コンピュータで使用で きる トップ ・レベルの BA-- 43-SICに くらべ て機能 に限界があ ることも事実であ る。現在,統計学 の教科書 のなかには統計量 の数 理的説 明を行 うはか りでな く同時 にBASICによ る計算 プ ログ ラムを併記 して解説 を加 えてい るも のがある。 これ は,実際家 に とって極めて有益 な 書物であ る。 しか しなが ら, これ らの書物 で用い られ て い るBASIC言 語 は大 型 コ ン ピュータの TSSター ミナルでの利用 を予定 した ものが多 い。 この種 のBASICはBASICの なかで も最 も高 い 機能 を もった ものであ るため, マイ コン ・レベル のBASICで は対応 で きず,せ っか くの良書 が無 益 になる場合 がでて くる。具体例 を示 そ う。 BASLCの最 も大 きな魅力 は行列演 算機 能 を も つ ことだ と思 う。少 な くとも多変量解析に とって この較能 の役割 は大 きい。次 に示す プログラムは 3次の正方行列A A= 1 0 1 2 3 2 2 0 1 r L の逆行列 A 1を求め るための もので あ る。 こん な小 さな行列 で も実際に逆行列 を計算す るとなる と実 に面倒 な ことはご承知の と うりであ る。 しか し, プ ログラムは極めて簡単 であ る。 0 0 0 0 0 0 0
1
2
rf.)
4
5
.8
7
REM 守 7・ワ千"ヨつレ':J DATA 212,1,0,3,0.11I2,1 DIM A(吉
事3),B(3,3)MAT READ A(3,3) MAT B=INV(A) PRINT B この秘密 は,MAT文であ る。上記の プログラムの 40行,50行,60行でMAT文 を使用 しているが, 40行 は20行 のDATA文 に書 かれ た行 列Aの要 素 を読み込み,50行でその逆行列 を計算 し,60行 で結果 を出力す る。しか し,マイコンのBASICで は, この便利 なMAT文が使 えない。 したが って, マイ コンでは同 じことを実行す るために別 の プロ グラムを工夫 しなけれ ばならない。 そ こで, これ と同 じ内容の ものをマイ コン向に MAT文 を使 わず に書いた のが次 に示 す プ ログラ ムであ る。 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 5 0 0 0 0 C . 占 0 0 0 0 0 0 C ・ l つふ rjlJ 4 5 も 7 8 9 C . 1 2 TZU 4 5 も 7 8 9 0 0 1 2 1 2 rrJ ,I) 4 tJ=) 占 7 日 9 1 1 1 1 1 1 1 1 1 1 2 2 2 2 (・4 2 2 2 2 2 P:] M rJ.I r,ll ∧H D T T M T M 3 R R A X X R E A : ll O ︻U E E E hU R D D N F F R N N F 守"T・ウキ"ヨウLJ'.I) 27271709370?1,2,1 A(3▼3)IB(3) Ⅰ=l TO J=l TO A(
Ⅰ
,∫)
∫ I N N t(=l TO N P=1′A(ド,ド)
FOR Jニl TO N a(∫)=AH=:,∫)半P A(K事
J)=B(∫) NEXT∫
FOR Iェl Tロ N IF I=K THEN 24O Q=A(Il
に〉
IF Q==O THEN 240 FOR J=l TO N IF J=K THEN 220 A(Ⅰ
IJ)=A(Ⅰ
,.I)-GはB(J) NEXT ∫ AくIlに)ニーP*Q NEXT I A(K,ド)こP NEXT r: FOR I=l TO N FロR J=l TO N PRINT A(IIJ)芋 NEXT ∫ PR王NT NEXT I END 行 列 要 素 の 読 み 込 み 逆 行 列 の 計 算 マ イ コンではMAT文 が使 えないので,その代 わ りにFOR文 とい うものを使 っている。このため プログ ラムも計算時間 も長 くなる。先 の プ ログラ ムとこの プ ログ ラムの対応 は次の とうりであ る。 先 の プ ログラムの40行 に相当す るのが ここで は io行か ら80行 まで,50行に相当す る逆行列の計 算が90行 か ら260行 まで,60行 の出力 に相 当す るのが270行か ら310行 までであ る。 マ イ コンBASICは この よ うに機 能 が 劣 るた め, プ ログラムが複雑 になる。 しか し,計算がで きない とい うことではない とす ると, マイコンの ための細かなプ ログラ ミングのテ クニ ックを積み 重ねてい くことが肝要である。価格面で極めて有 利 な コ'/ビュ-タをいかに上手 に利用すれ ばよい 一44-か とい うこ とで あ る。 筆 者 が この資 料 の発表 を思 いた った理 由 は,釈 計 解 析 関 係 のBASICに よる ア プ リケーシ ョン ・ プ ログ ラムが ほ とん ど発 売 され てい ない とい うこ とが ひ とつ あ るが, い まひ とつ は, プ ログ ラ ミン グ ・テ クニ ックの情報 交 換 に資す るためで あ る。 大 方 の ご助 言 ・ご批判 を乞 う次 第 であ る。 な お,以下 で は5つ の プ ログ ラムを発表 した が, こ こで は, 基 本統計量 の計算 と
2
変量 の解 析, し か も記 述統 計 の レベル に限定 した。 この他 に, 多 変量 解 析, 時 系列分析 お よび推測統計 の レベル を 計 画 して い るが, これ は次回 に譲 るこ とに した。Ⅰ
デ ー タ ・フ ァ イ ル の 作 成(プログラム 1)
コン ピ ュー タにデ ー タを イ ソプ ッ トす る方法 に は3
通 りの ものがあ る。 1つ は,プ ログ ラムの中 にDATA文 を使 ってあ らか じめ書 き込 んで お く方法 。先 に示 した逆 行 列 の プ ログ ラムが これ に当た る。 2つ め は,tNPUT文 を使 って,キーボ ー ドか ら 直接 入力す る方 法 であ る。 3つ め は, デ ー タをあ らか じめ カー ドや テ ー プ や デ ィス ク な どの外部記憶 装置 を用 いて これ に記 図1 データファイル と分析 プログラムの関係「
ラー
タ(n)ト
ファ
イル
作 成 プ ロ グラム
憶 させ フ ァイルを作成 ・保存 してお いて, 必要 に 応 じて必要 なデ ー タを読 み込 ませ る方法 で あ る。 通常 は,3
通 りの方法 の うち どれ か一 つ を選 ん で用 い るので ほな く, 用途 に合せ て適切 に組 み合 せ て用 い る。 統計解 析 の場合, 分析 され るべ きデ ー タを プ ロ グ ラムの中 に書 き込 む ので はな く, プ ログ ラムか ら切 り離 して フ ァイル化 して お くの が有 利 で あ る。 プ ログ ラムの中 に書 き込 んで しま うと, 同 じ プ ログ ラ ムで 異 な るデ ー タを処 理 した り, 同 じ デ ー タを異 な るプ ログ ラムで処理 す る必要 が あ る た び毎 に, デ ー タや プ ログ ラムを書 き変 えなけれ ば な らない。 これ は実 に "しん どい''ことで あ り, 実 際,全 く ``や る気''が起 らない。 キ ーボ ー ドか ら入力す る場 合 も同 じで あ る。 この場 合, プ ログ ラムを書 き変 え る必要 はない けれ ど も,処理 のた び毎 に同 じデ ー タを入力 しなけれ ば な らず,全 ぐ しん どい"0"これ で は何 のため の コン ピ ュー タが ' と思 いた くな る。した が って,最 も有 利 な方 法 は, デ - タを分 析 プ ログラ ムか ら切 り離 して テ -プな どに フ ァイル と して保 存 した うえで,任 意 に必要 な デ ー タを必要 な分析 プ ログ ラムに リン クす るや り方 となろ う。 この関係 は図1に示 す と う りであ る。 図1
に よれ ば,分析 プ ログ ラムに よって結果 が 出 され るまで に,(1)ファイル作 成 プ ログ ラム, (2) リンク ・プ ログ ラム, (3)分析 プ ログ ラムの3種類 の プ ログ ラムが必要 で あ る。 本節 で は,(1)の フ ァ イル作成 プ ログ ラムを紹 介す る。 プ ログ ラム1が それ で あ る。 14 5-プログラム1の説明 プ ロ グ ラ ム1 (デ ー タ フ ァイル の作 成 ) 1REH【7●Drラム1
テ
■
-
977仙J叩tイ Tが シユつイチ ヤカ■Jサーィ.8.11 36RAPHC:CONSOLESOI21lPR川T'7-7-ヲテ'一拍ウニhITTイ':PRINT一m力Ⅰ粥 日】ヲ才シ〃イ1 46ETR:IFR=OTHEN4 占PR川TCHRItもI:INPUTI I】ウ /ヒ,'サJ
l
l;ZI 86RAPHI
II
CIOl!L川EOI40I
OI128I300I128I300140I
OI40 王CONSOLESもl19 10INPUT'77イIi
-
ム
J;FfりNPUTl
†
:
I71Zウ N=';∼:lNPUT■1:tスつ.J力ztM:-;M IOPRINT'lI78Zt乃ヨ:l'tPRINTI21肘ム
シ
ソ
一
㌢
■
;:lNPUT.r朽テ'zJ)一紳 PRINT bOJFA=JTHEN川PUT■M'ウフケ:)1■†デスカ YeS=J No:2 .;8!YI=■†シフ-1■ ム2IF柑=lII(Nく)47160T88 701FA=2THENINPUT'テt-タJ151‡叩 1:17ン折 -;C:B:2 iY事三一7Jitt1 71IFB=lTHEND=1:DIMB引付):687880 74PRINT:PR川T'lIl;YI;"「ウニユウlJ3r :PRINT1211ンZll'†二ユウlJ37 .;:lNPUT'ド事ラテ'Z力 I;D 80mPIA川I
r
l
hAHPlHCONSOLES
OI24 糾 6MPHIII
CIDl:LINE4
0I
OI40I199-152,199I15270王LINE72I
OI72I199:LINE40I
8,152,8 鮎 PRINTTABH0HI竹27メイ?l:COJISDLE52,24 88FORJ=lTO【:PR川TTAB柑-LENtSTRI(J日日J;!lNPUT-I;A引J):NEXTJ 95CURS8R20I20:川PUTリr!1イテtzA?];I:lFX=1THEN100 9占60SU別00 100CONSOLES
Ol24:6--1 102IFtBく)日暮柑=lITHENAf='†:lフ
I
l
■
1041FD=2THENA事=■竹コr loblFfA=2日(D=2JTHENC暮
=l子
ン
I
108IF川=lI暮LD=2)THENCI=.N0. -1140NDGDTOll
占l150 118FORl=lTON:OND60TO120I12I 120READA暮:BI(I)=A事壬PRINT一 No.
I;AI:687012占 121PFuNTI NE
I
.
I;AfH 12もCONSDLES
2124!FDRJ--1TOPl:PRINTTAB柑-LEN(STRI(J日日J;I-1NPUT一
一;A(
ⅠTJ):NEXTJ 130CURSOR20I20:INPUT■叩1514fl加?I;
XHFX=lTHEN140 13460SU8400 IloCt)NSOLESOI24:NEXTl:GOTOlq0 154FDRJ=l
T
D
M
!PR川Tl I;C事;■ ●;At;I:CONSOL【S
2I24:FORl=lTDN 155IFJl=2THENE三Ⅰ◆C-I 15もIFA=lTHENE=1 158PRINTTA別8-LENtSTRI(EIII;E;HNPUT.-;At
lIJ):NEXTI 159CURSOR20I20:INPUT一門'シ作●加71;XHFX
=lTHEN170 1舶 60SU糾00 170CDNSOL【SOI24:NEXTJ lワoREM -【Nrit亡】 1926月APNCiPFuNTCHRfはI
;'CoIPUteriSWitin9thedatatotape.
I
210-8PEN王PR川T/TZIIFII
NI
P
l
l
AI
8I
CID:FロRJ=lTO【王PRINT/TAftJ):NEXTJ 21‡OND60TO214I254 2111FD=lTHENFロRl=lTDN:PRINT/TB引日:NfXTI-46-2
2
0F
O
R
I
=
l
T
O
N
王
F
O
R
J
:
l
T
O
r
l
:
P
R
I
N
T
/
T
朋Ⅰ
,
J
I
:
N
E
X
T
J
,
Ⅰ
:
6
0
T
O
2
5
占
2
5
4F
O
R
J
=
l
T
O
M
!
F
O
R
I
=
l
T
D
N
:
P
R
I
N
T
/
T
A
‖
,
J
HN
E
X
T
I
7
J
2
5
もC
L
O
S
E
:
R
E
V
:
E
N
D
4
0
0C
U
R
S
O
R
2
0
I
2
0
:
I
N
P
U
T
一群77■イJJ)rI
l?
I
;
I
I
O
IF
8
R
K
:
l
T
O
X
:
C
U
R
S
P
R
2
0
I
2
1
:
t
N
P
u
T
lマ和■ィ./E
l
t
:
l311I
;
Y
王
C
U
R
S
O
R
2
0
I
2
2
:
P
R
I
N
T
Y
;
4
0
21
F6
=
l
T
H
E
N4
0
5
1
0
3川P
U
T
一
三
一
;
A
引Y
I:6
8
T
O4
0
7
4
0
II
FA
=
2T
H
E
NY
-
-
Y
-
C
+
1
1
0
5I
FD
=
lT
H
E
NI
N
P
u
T
■
=
I
;
郎I
I
Y
)
t
C
BI
FD
-
-
2T
H
E
N川P
u
T
.
=
-
;
A
t
Y
,
J
I
I
O
TR
E
T
U
R
N
5
0
0D
A
T
A
tI州 -1
,
7
7t
IJIイワテI三円t
,
7
千1,f胡`†I
7
7
!l77イJl.ラ‡ll桝 ■Irン7'†イ17'事が'糊 ;つ'
耶
功tワ5
0
5D
A
T
A
:
棚・
サI桝I
棚ワ
I
7
7イI†7T・yIl力・JIr7I:'Z・榊I7Lf
I
三Ⅰ
I
i,n・I‡州 Im力
Iヒ31jt
,T5'ワ・附 ,hM
5
1
0D
A
T
A:}マ717mI
ヒ
防
マI†7r事IhrJ?,
約
・ワ
,
Ⅰ
ヒ3,3ウ事I77柑I
†
力IIlh・叩Iワ
脚 I抑 タI三fll‡'氾■シ7'柑 ワ図2 くデークマ トリ.クス) 1 2-All----jI-I---- m l ・・ -・・ ・・ -・ ・・J ・・ ・・ ・・ -花 〟 A .. .. --." " A ." -. " .. A A12--I-AlJ--一一--Aim A,2- --A,)---Azm An2---AnJ 図3 データ入力順序の選択 (l) クロスセクションデータ この プ ログ ラ ム は,図
2
の よ うな サ ン プ ルn 個,変数m個 か らな る多変量 デー タマ ト1)クスを カセ ッ トテープに記録 してデ -タフアイルを作成 す るものであ る。プ ログラムは同時 に≪ファイル作 成 の 日付, ファイルの名前,サ ソプルの個数 ,変 数 の個数,変数 の名前,データの性質,入力順序≫ をテ-プに記録す る. また時系列 デ ータの場合 に はデー タの開始年,都道府県別 デー タの場 合 には 都道府県名が テ ープに記録 され る。 この うちデ ー タの性質 とはクロスセクシ ョン ・デー タか タイムシ / (1'圃 道鵬 如± 旦ト
ー 一 一 一 都道府県別^h ・2)唾 二二
二
二二:; … 芸 喜喜月芸 ≡ (2) タイムシ T)-ズデータ (1)年
次別
入 力 (2) 変 数 別 入 力 リーズ ・データか,ク ロスセ クシ ョンの場合都道府 られ る。 これが デ ータの性質 と組み合 わ さって図 県別 か否 か とい う意味 であ る。 これ は, 図3に示 した よ うに分岐点 ごとに選択 を求めて くるので1 又 は2の数字 キ ーを押せ ば よい。入力順序 は大 別 してサ ンプル (年次)別 入力 と変数別入力 に分 け-4
7-3
の最 右 欄 に示 す5
通 りの選 択 肢 に至 る。 コン ピュータでは, デ ータマ トリクスの入 力 は行 ベ ク トル又 は列 ベ ク トルを一組 として行わ れ る。つ ま り図2でいえば行 ベ ク トル とは図2の横 の行,列ベ ク トル とは縦 の列 で あ る。 この両者 は入 力 の順 序 が完 全 に逆 に な るわ けであ るが, この プ ログ ラ ムで はいず れ の順 序 に も対応 で きるよ うに した。 とい うの は,社 会 調査 な どの よ うに調 査票 が サ ン プル ご とに作 成 され てい る場 合 に は行 ベ ク トル ご と (サ ンプル別 ) の入 力 が便利 であ る し, 表
1
の よ うに デ ー タの一 覧表 が作成 され てい る場 合 で も 変数 の個数 が多 い場 合 に は列 ベ ク トル ご と (変 数 別) の入力が便 利 であ るか ら,何 れか一方 だ けを 採 用す るの は フ ァイル作 成 プ ログ ラム と して汎 用 性 に欠 け る と判 断 したた めであ る。これ に よって, あ らゆ るタイ プの数値 デ ー タの フ ァイル作 成 がで きる。 なお, 1つ の変数 につ いて ク ロスセ クシ ョ ナルに時系列 デ ー タが作 成 され る ことがあ るが, 表1 (チ-タ 1)都道府県の民生費 と人 口 この場 合 は年 次 又 は月次 な どの系 列 を縦 の列 に置 いて変数 とみ なせ ば よい。 なおサ ンプル の個数 は最 大2
5
5
まで可 能 。 変数 の個 数 も最 大2
5
5
まで可能 であ るが, メモ リーの使用状態 に よって制 限 され る。 メモ リー許 容量 を超 える と プ ログ ラ ムの実行 が中止 され る。 サ ンプル数 を最 大2
5
5
に設 定 した場 合,変 数 は1
0
個 ぐらい に とどめ ておい た方 が安全 であ る。 サ ン プル数 を少 な くす れ ば変 数 を多 くす る こ とがで き る。 変数 の名前 は,1
0
文字 以 内 に と どめ て お いた方 があ とあ と便利 で あ る。 プログラム1の 実行 (単位 100万円 :100人) 歳 出 民 生 費 児 童 老 人 人 口1
9
歳以下6
5
歳以 上 総 額 福祉費 福祉費 児童人 口 老齢人 口事
I
J
イ
r
つ 122830
87
0
1
4
0
13恥2 13515 555糾 17072 13707
才
t
l
J
3占97
1
82
3
も
3
4
ム522 33糾1
5185 一808 1357 ーワテ 355
L
4
1
2
1
9
7
7
5556 33871
一ObB 4289 1374 ミ†‡'3
糾
1
7
22
2
3
0
1 8
9
04 3225 20719 占2仙 1707 7I,3
4
1
4
2
l
1
9
5
d
2 4
9
9
3
3占一3 12も94 3570 1351 †7rサ3
2
仙
1
71
2
9
3
8 4
2
4
7
2017 12510 3520 川59 刀5I7一
3
9
日
51
8
3
9
3 5
ー
4
4
3217 202も2 占15占 2074 LJrラ
‡
'
一
2
3
0
5
22
2
5
3
2 占
2
7
占
3822 255恥 7981 23占8けf
'
3
I
恥
0
3
1
7
〇
一
l
&
2
1
5 2
4
4I 179鮎 5514 l占52rン
7
3
2
2
5
2
8
1
7
㍍
0引
2
8 3
2
7
5 I別 13 5占B7 1833 †イ175
8
2
4
1
4
3
0
5
8
9
1
0
7
7
4
8
1
0
2
54002 18231 3327 tJr5
7
1
5
3
8
2
一
5
9
09
3
7
7
3
8
7
9
47231 15379 3301 トつ事3ウ2
占
2
0
2
1
3 1
9
7
5
3
8
5
5
4
3
7
5
1
8
2
5
1
1589も 32597 8舶7力
1
力lワ占
8
5
0
8
8
3
1
3
3
5
?
?
4
7ム
7
砧ふ
9210 2201占 一291=
イ
カ'サ5
4
8
4
7
8
2
0
2
7
4
7
5
9
24
4
5
5
2
4
3
丘も 7098 2ム52n
72
5
0
3
1
3
l
o
川ュ
仙1
3
2
3
5
2
1
0
9
0
9 3157 1215イ
シ
脚
2
5
4
相一
1
1
9
一
7
4
8
1
9
2
8
4
7
l
i
O
9
2
5292 1137刀イ
2
0
3
3
8
3
1
0
3
4
1
3
4
7
8
2
0
6
27
9
0
2
2307 930 †71シ1
9
0
9
5
9
1
0
7
3
t
3
1
0
0
2
3
7
98
0
0
5
2478 9 1 5 1rJ4
1
2
5
?
7
2
2
9
1
占
ム
0
ム
†
5
2
4
1
2
0
7
8
8
5
9 d Ll 2 5 0 0 I-7
3
5
3
1
0
9
1
5
2
3
3
5
引払
3
2
7
0
l
叫
7
4
ム
0
7 7 1 8 9 3 5TT一
九
4
?
4
0
1
0
2
0
2
3
2
7
7
4
0
4
1
5
1
3
4
5
3
0
l
O
も
1 0 3 1 4 4 7イ事8
5
3
0
占
1
4
1
4
1
4
1
2
9
9
8
9
9
一
2
ふ
2
4
9
9
2
0
1
8 5 1 5 4 4 三Ⅰ
3
0
砧
5
0
1
7
5
3
14
1
2
9
2
7
3
1
1
ム
9
0
7
5
0
7 8 1 8 5 6シ
カ
'
2
2
9
ム
2
2
l
J
3
8
53
古
さ
占
2
4
0
0
J
O
8
2
5
3
小
吉 1 0 7 7I
ヨ
ウ
ー
3
2
一
2
7
3
1
9
3
7
3仙
1
93
5
8
d
2
5
0
9
7
7
5
9 3 2 5 0 87
才
ー
力
1
0
3
3
8
7
ム
ム
2
3
0
3
1
8
0
ふ
5
2
4
4
2
8
8
5
0
3
7
2
7
1
5
一 ム 2 1 0t
一
7
ラ
つ
コ
'
2
7
3
0
1
0
0
9
も
1
7
1
3
1
5
4
6
5
5
?
5言
7
1
0
古
5
5
5
0
0
1
0
1
1
5
5
1
5
?
3
1
5
9
5
7 t も 7 13
0
1
l
1
2
0
8
5
3
7
7
9 1 0 8 8ワ
h
カ
ー
M
マ
2
1
9
3
4
1
4
0
2
2
0
5
34
5
7
4
3
1
2
ふ
1
0
8
8
3
3
1
2
5 1 3 0 42
3
3
9
1
0
7
一
5
3
7
4
7
2
2
5
2
ム
0
8
1
1
ム
9
2 J 7 1 7 シマi2
一
8
0
2
5
l
帖
2
9
3
8
1
0
3
ム
9
5
7
8
8
9
2
1
9
8 l o 仙 才肘73
5
1
2
J
l
2
0
0
9
占
5
8
5
1
3
9
1
7
J
8
7
1
8
5
5
2
9
2 2 2 9 tE)5lマ4
4
8
2
5
02
7
5
1
0
1
0
q
b
ふ
0
3
8
2
7
4
1
7
8
3
仙
2 8 8 0†
ト
マ
ワ
r
シ
マ 2
事 3
1
2
7
7
別
2
7
8
8
1
7
8
8
7
5
1
8
0
仙5
5
1
5
7
3
5
I
一
伯
1 8 5 91
8
3
8
2
5
2
1
1
3
2
ム
5
8
3
2
4
2
3
7
2
9 占 7力
.r
ワ 2
0
2
3
1
8
1
5
0
5
3
3
5
3
4
3
5
1
2
9
8
7
5
2
7
8
9
1 1 8 1 Ⅰヒメ3
0
1
0
占
9
1
8
3
2
も
5
1
5
5
一
2
5
3
1
5
0
5
5
4
1
5
3
1 7 ふ 3 コつ事2
7
も
3
8
3
1
7
9
も
1
4
2
0
7
3
3
7
7
8
3
5
4
2
2
8
3
1 1 2 17
1
7
力
b
t
9
5
も
2
ふ
7
1
3
2
9
2
4
9
7
9
9
2
4
5
ふ
占
5
1
3
7
5
9
4
3 5 2-4
8-出所)地方財政統計年報 お よび昭和55年 国勢調査 で は表
1
のデータをカセ ッ トテープに記録す る こ とに し よ う。操 作 は極 め て 簡 単 で あ りRUN 匝司 とキー入力すればモニタ (テ レビ)画面が画 面 I- 1の よ うに写 し出 され る回 とは キ ャ リッジ ・リターン ・キーの ことで, このキーを押 す ことによって コンピュータは命令 を理解す る。 RUNはプ ログラムの実行命令であ る。あ とは画面 の指示 に従 って必要 な情報 を入力すれば よい。 画面 Ⅰ-1に従 ってテープを交換 した ら[司 を 押す。 次の瞬 間,画面 は 日付等 をたずねて くるの で答 えてやれ ば よい。画面I- 2はその結果であ る。 次 に変数 名 を要求す るので1
0
文字 以 内程度 で 変数名 を入力す る。画面Ⅰ-3
がその結果であ る。 入力が正 しけれ I重 ]を押す,誤 入力があれ ば訂正 のため かこ且 を押す.次に数値 デ ータの入力 を求め 画面Ⅰ- 1 テープ交換チ
-
7
'ヲテ'-サヨウニ
カⅠ叩イ M如
タラ 【l
】ヲねMイ 画面Ⅰ-2 初期設定f
3
1J
ヒ
T
ケ
I
l1
9
8
1
.
1
0
.
2
3
7
7
仙7
-
ム
三ンt
イヒ
1
9
7
8
IrJコ
1
1
9
8
0
†
ン
7
●
l祁N
=
1
7
㌔
J
Z
ウJ
加I
r
l
=
7
I
)
7
D
Z
t
ワ
如ン
2)1LL的-
T
r事ラテ■
Z
) 1
H■?乃ン11†デZ
力
Y
e
5
-
-
1 N
0
-
-
2 1
画面Ⅰ-3 変数名の入力 画面Ⅰ-4 数値データの入力 て くるので 「ホ ッカイ ドウ」か ら順 に入力す る。 入力が終れば訂正 の有無 を求めて くる。画面 Ⅰ-3では変数1が誤入力になってい るので訂正 して いる。訂正が終 ると同時 に 「7オモ リ」 と表示す るので青森のデータを入力す る。 同 じ要 領で 「オー4
9-キナ ワ」の入力が終 るとテープが 自動的に回転 し てデータが記録 され る。 以上が プログ ラム
1
の実行例であ るo プ ログラ ム言語 に知識 のない人の利用を考 えて,画面の指 示 に従 って進行すれば 目的が達成 され るよ うに し たつ も りであ る。I
I
デ ー タ と分 析 プ ログ ラ ムの リン ク(
プログラム
2)
ファイルに納 め られたデータを分析す るために は, チ-プか らデータを読み出 して コンビュ-メ 内部の記憶装置 に移 さなけれ ばな らない。 プ ログ プログラム2 (データのリンク) ラム2
はそのための ものであ る。文番号1
0
でデー タの 日付, ファイルの名前, サ ンプルの個数,変 数 の個数, デ-タの性質,入力順序 な どを読み出 す。文番号1
6
と2
0
で入力順序の判断を行 ない, それ に応 じて, サ ンプル (都道府県,年 次)別 に 入力 され た データは文番 号3
0
に よって データを 読 み取 り,変数別 に入力 された ものにつ●いては文 番号4
0
によって読み取 る。データの読み取 りが終 る と文番号5
0
に よって ファイルを閉 じて テープ は巻 き戻 され る。 データの読 み 取 りが 終 る と文 番 号5
2
以 降 に よって データの説 明が画 面 に表 示 され るO申 面1R
E
Pl【フ
I
D
rラ
ム2
テ'-チIl
J
)
ワ
I
J
)
1知
り付 1A■/rイhl?]26
朋P
H
C
王
C
O
N
S
O
L
E
C
仰
3P
R
I
N
T
C
H
R
別
も
)
,
■
テ
ープ●ヲテ'-乃ウニ如
Mイ
'
:
P
R
I
N
T
-ト
l
J
hI関 目】ヲ恥
Mィ.46
E
T
R
HF
R
-
-
O
T
H
E
N
4
占P
R
I
N
T
C
H
R
I
t
も
)
;
T
bA
P
U
t
e
ri
sr
e
a
d
i
n
gt
h
ed
a
t
a
.
:
P
R
I
N
T
7R
O
P
E
N
1
0I
N
P
U
T
/
T
Z
事
I
F
f
I
N
I
P
I
I
A
,
B
I
C
I
D
:
D
l
【
A
t
N
I
P
HI
A
f
川1
1
4F
O
R
J
=
l
T
8
【
:
l
N
P
U
T
/
T
A
f
t
J
HN
E
‡
T
J
l
占I
FB
く
)
I6
0
7
02
0
1
8D
I
N
B
I
t
N
l
:
F
O
R
l
=
l
T
D
N
HN
P
U
T
/
T
B
引日:
N
E
X
T
2
00
ND印7
0
3
0
I
1
0
3
OF
O
R
I
=
l
T
O
H
:
F
O
R
J
=
l
T
O
r
(
:
川F
U
T
/
T
肘Ⅰ
,
J
l
:
M
E
X
T
J
I
I
:
6
r
J
T
O5
0
4
0F
O
R
J
-
-
1
T
D
【
:
F
O
R
I
=
l
T
O
N
!
l
N
P
U
T
/
T
A
uI
J
HN
E
X
T
I
I
J
5
0C
L
O
S
E
:
R
川
占
OF
R
I
N
T
C
H
R
f
t
ふ
HC
O
N
S
t
)
L
E
C
8
0
占
21
FA
=
JT
H
E
NY
f
=
t
ワ
E
I
Z
門お
ン
ー
占
3I
FA
=
2T
H
E
NY
f=
■サ
イL
SIIJ-Z-IL
41
FD
=
lT
H
E
NY
f=
'
H.
ウ7
竹 11γも
5C
I
=
S
T
R
I
N
6
引
1
-
I
I
3
8
日F
R
I
N
T
C
f
7
0P
R
I
N
T
-
テ
'
-
9It
Z
l
ケ '
;
Z
f
:
P
f
uN
T
'
7
7
仙 J17
I
■
;
F
I
:
P
R
川T
-
†
ンフ●▲Ihl■'
;
N
7
1P
R
I
N
T
IlンZウJ如t'
抑P
R
川T
.
コ
J
テ一一号Jl【
-
;Y
事;']テ
■ス'7
2I
F
A
=
l
T
H
E
N
P
R
川T
:
P
R
I
N
T
C
I
:
6
0
T
D
8
0
7
3P
R
I
N
T
'テー-1I Iわ^51377I
;
C
;I71的一
;
C
+
∼ll;.ill7
㌢-:
P
f
uN
T
C
I
8
0P
R
I
N
T
' 1ンコウイ事ラl
I
:
P
R
I
N
T
C
f
王
F
O
R
J
=
l
T
O
M8
21
F
t
J
)
=
1
1
川 J
く
=
2
1
)
T
H
E
N
C
u
R
S
O
R
l
l
I
J
-
1
8
3I
F
J
)
=
2
1
T
H
E
N
C
U
R
S
D
R
2
7
I
J
-
I
1
8
4P
R
I
N
T
J
;
-
1;A
I
f
J
)
:
N
E
r
T
J
:
P
R
Z
N
T
C
I
9
06
0
S
L
J
B5
0
0
5
0
0P
R
川T
■
カワ:)沖ラ【
D
K
】ヲね1
7
イ
'
5
0
46
E
T6
l
:
l
FG
I
く
)
■
D
■T
H
E
N6
D
T
O5
0
4
5
0
8R
E
T
U
R
N
ー5
0-ⅠⅠ-1がそれ であ る。 これ は,大量 の ファイルを 作成 した と き,読 み込 まれたデ ータが 目的 のデ ー タであ るか ど うか を確認す るための ものであ る。 この プ ログ ラムはRUN[垂 ]命令 だけです べ て 自動的 に実 行 され る。 ところで, この プ ログラムは分析 プログ ラムに 必須 の補助 プ ログ ラムであ るか ら,分析 プ ログ ラ ムは, この プ ログ ラムに続 けて例 えば文番 号100 番 か ら書 き始 めれ ば よい。 以下 では,統計解析 の 中 で も最 も基 本 的 と考 え られ る分 布 統 計 量,相 関 ・回帰分 析,質的 デ ータの ク ロス集計 とカイ2 乗 ・クラマ -係数 の分析 プ ログ ラムを紹介す る。 画面ⅠⅠ- 1 データの説明
Ⅰ
Ⅰ
Ⅰ 分布統計量,度数 分布 お よび グラ
フ (
プログラム 3)
プログラムの説明
(1)分布統 計量 分布統計 量 には表2
に示す よ うにいろい ろの も のがあ る。 測度の役割 に よって分類す れば,代表 値 ,散布度 ,歪度,尖度 とな り,算出方法 で分類 す れ ばモ ー メン ト系 と分位数系 お よびその他 とな る。 モーメ ン ト系 とは平均値 か らの偏差 の積和 を 中心 に求 め る もので最 も一般的 に使用 され てい る ものであ るo分位数系 とはデ -クを大 きさの順 に 表2 主な分布統計量 モーメント 系 統 分位数系統 その他 ・代 表 値(詣 ) 算術平均 メジ7 ソ 幾何平均調和平均モー ド 散分 度 分整 準鰻重散 四分位 レンジレンジ (ひろが り)平均偏差 十分位 レンジ 変動係数 四分位偏差四分位変動係数 歪 度 (積率)歪度 四分位歪f受 ピアソン (ゆがみ ) 十分位歪度 の歪度 並べ て求 め るものであ る。 ここで は表2
の うち最 もポ ピュラ リテ ィーを もつ と思われ る算術平均, 標準偏差,変動 係数, メジアソを求 め これ にデ ー タの最大値 と最 小値 を追加 した。 歪度や尖度 を求 め るのは容易であ るが現実 にそ の数値 の大 きさを評価す るの は難 しい。 (積 率) 歪 度 につ いて は 0よ り大 きいか小 さいか,(積率) 尖度 につ いては3よ り大 きいか小 さいか とい った 程度 の数字 の読み方 しか しないのが一 般 であろ うノ。 もしそれだ けの ことであ るな らば, 算術平均 とメ ジ7 ソを比較 した万 が分布 の形 を判 断す るのに実 際的 な意味 があ る よ うに思 う。す なわ ち,算術平 均 よ りもメジアンの方 が大 きけれ ば分 布 は負に歪 んでい るので あ り,算術平均 よ りもメジア ンの方 が小 さけれ ば分布 は正 に歪 んでい る こ とがわか る。 これ に最大 値 と最 小値 が分かれ ば分布 の形 ほほは 判 断で きるだろ う。 散布度 は標準偏差 と変動係数 だけに とどめた。 分布 の広が りを知 るには標準偏差 は絶 対必要で あ る。多変量問 で分布 の広が りを比較 す るに はや は り変動 係数 が必要 であ る。 プ pグラム3で は, まず メジア ン と最大値 と最 小値 を求 め るため に,文番号102か ら109でデ ー タを′l\さい順 に並べ か えてい るo最大値 と最 ′j\値 は後 に度数分布 の分析で度数 を カウン トす る時 に も必要 にな るので, ここで は思 い きって もとの配 列 の内容 を こわ して小 さい順 に並べ変 えてい る。-5
1
-プ ロ グ ラ ム3 く分布統計量,度数分布)
1
0
0R
E
M【
フ●
8rラム3
7
■
ンフ●Mサイlhつ Tカタシコウイ事
T
h
-
Irl
'
力
'
7】
1
0
1D
川 S
H
I
I
占
)
,
C
(
2
0
1
1
0
2R
E
E
【テt一号/丁ラAr.7)Ⅰ】1
0
3F
O
RZ
Z
ZT
ON
-
1:
F
8
RJ
=
lT
D【
l
o
tI
FA
tl
l
J
)
)
=
A
H+
l
l
J
IT
H
E
NK
=
A
HI
J
)
!
A
fI
I
J
I
=
A
H
◆l
l
J
HA
H
◆l
I
J
)
=
K
I
O
5舵X
TJ
I
1
l
o
bF
8
Rl
=
lT
ON
-
2;
F
O
RJ
=
lT
O
M
I
O
7汀 A
t
l
]
J
)
)
A
H+
l
I
J
IT
H
E
N1
0
3
l
o
o1
はX
TJ
I
1
1
0
91
FA
H
I
J
-
1
I
)
A
t
N
-
l
l
J
l
l
IT
H
E
N1
0
3
1
1
0R
EIH -【T'Jフ●トつ14I)ヨつJサイヤシ】1
1
2F
O
RJ
=
lT
OH
!
S
l
=
0
:
S
2
=
0
!
S
3
=
〇
三
S
4
=
0:
S
tJ
1
4
)
=
A
fl
I
J
):S
tJ
I
ム
)
=
A
(
N
-
J
)
1
1
1I
FN
=
I
N
T
t
N
I
2
)i
iT
H
E
NS
(I
,
5
)
=
t
A
t
N
/
2
,
J
)
+
A
E
N
/
2
1
1
I
J
I
I
/
2
:
6
D
T
O
1
2
0
1
1
8S
f
J
,
5
)
=
A
l
川+
J
)
/
2
I
J
l
1
2
0F
O
Rl
=
lT
ON:
S
l
=
S
l◆
耶Ⅰ
I
J
)!
S
2
;
S
2
+
A
H
,
J
)
I
A
H7
J
I:
N
E
X
T
I
:
S
l
=
S
l
/
N
1
2
2S
2
=
S
e
R
f
S
2
/
∼
-
S
l
暮
S
l
HS
3
-
-
S
2
/
S
l
暮
1
0
0!
S
(
J
I1
)
=
S
l
:
S
(
J
I
2
1
=
5
2
:
S
t
J
I
3
1
=
S
3川E
X
T
J
1
4
5R
E
.NF 【加力/ilコ†L
J
3
7
】
1
4
占P
R
I
N
T
C
H
R
I
t
も
I
;
`【
l
;
F
f;']I7
'
ン
7
日:D
I
=
S
T
R
I
N
6
日JJI
1
0
1
1
4
7F
D
R
J
:
J
T
D
ム
:
P
R
川T
D
事
;■
T一
;
;
N
E
X
T
J
:
P
R
I
N
T
D
事
1
4
8F
O
R
l
=
l
T
8
2
:
P
R
川T
T
A
8
日0
)
;●l
■
;
:
F
O
R
J
:
l
T
O
も
:
R
E
A
DA
f
1
4
9P
R
I
N
T
T
A
B
H1
1
◆
1
1
暮
J
)
-
L
E
N
tAI
H;
A
I
;
:
N
E
X
T
J
:
P
R
I
N
T
:
N
E
X
T
I
1
5
0F
O
R
J
三
l
T
O
i
:
P
R
I
N
T
D
I;.十
I
i
:
N
E
X
T
J
:
P
m T
D
I
:
F
O
R
J
:
l
TO
M:
P
R
I
N
T
A
M )
;;
F
O
RK
:
1
10も
1
5
4P
R
I
N
TT
A
B
(
1
0
I
K
)
;
'
l
.
;
:
P
R
川TU
S
川6
.
1
対日l
l
l
.
書
J
;
S
t
d
,
K
)
;
:
N
E
X
TK王
P
R
川T
1
5
51
FJ
=
HT
H
E
N1
5
8
1
5
8F
O
R
K
-
-
I
T
O
も
:
P
R
川T
D
f;I十`
;
:
N
E
X
T
r
:
P
R
川T
D
I
1
5
8N
E
X
T
J
: F
O
R
J
-
-
1
T
O
d
:
P
R
川T
D
事
;
■
⊥
一
;
:
N
E
!
T
J
:
P
R
川T
D
I:
6
0
S
U
8
5
0
0
1
6
OD
A
T
Al・作ン
事I
I
ヒョウiltユ
JJハンドつI
I†付ヨつ
事I
l乃
■7
:llI
I
イヴ
.
J事I
'l
')■
竹
I
I
'1'ケ
イコつl
■1
占
lD
A
T
A
'l■
I
一I
.I.I1
7
0R
E
.L'.- Erzウノ.Wン日1
7
2P
R
I
N
T
C
N
R
I
は
);
'「.
7
つ
7
㌧
フ●ノ
ブ`
:
t
t
‡
-
:
P
R
I
N
T
:
P
R
I
N
T
■J
ll
ン
コ■ウヲシテ仙川イ
■
王
P
R
I
N
T
1
7
4F
O
RJ
=
lT
OH!
P
R
I
N
T∫
;
'
)I
;
A
f
t
J
HN
E
X
T
J
:
P
R
I
N
T:
l
N
P
U
T
■1:(ll':I:51マスカL
;
I
.
i:
P
R
川T
1
7
81
N
P
U
T
一?ラzJ舶-J
l
イワ†二シ7
Z
カー
;
C
:
I
F
C
)
2
0
T
H
E
N
1
781
8
2D
:
(
S
t
J
T
占
卜S
f
J
I
4
日/
C
:
氏
-
-
0
:
F
O
RK
=
lT
OC
:
C
は)
=
0
!
N
E
X
TK:
F
O
R
K
=
l
T
D
C:
F
O
R
I
=
l
T
O
N
】
鮎 1
F
I
A
Ll
,
J
I
)
-
-
S
LJ
,
4
)
+
D
I
K
-
D
川 A
LI
I
J
I
く
S
LJ
I
4
)
+
D
I
K
JT
H
E
NE
L
K
.
)
=
C
はI
r
l
1
8
7N
E
X
T
I
1
8
8l
FK
=
CT
H
E
NC
f
C
)
=
C
(
C
)
+
i
1
8
9N
E
X
TK:
P
R
I
N
T
C
H
R
I
t
ム
)
l
ワ
oR
E
.L,-【
ドZ77
㌧
フ
●
ヒョウ】1
9
ZP
R
川T
■【
-
;
A
I
L
J
)
I'L】ノ1.
才
17
-
ンフ
●t
ヨ1'・
.
L
l
=
S
T
R
I
N
SIL)-I
,
J
7
1
1
'
u
I
-
-
S
T
R
I
N
G
H暮-1
1
占
)
1
9
2P
R
川T
L
事
;
{T
l
;
:
F
O
RI
=
l
T
O
3
:
P
R
I
N
TU
t;-Tl
;
:
N
HTI:
P
R
Z
N
TU
事
1
9
3P
R
I
N
T
l I rZ?l 王 JI
L
t
手I
王.
1
9
4P
R
J
N
T
' 伯'ヨつ ミマ
ンJ
I
I
r.7つl
-1
9
8P
R
川T
L
l;一十■
;
:
F
O
RI
-
-
1
T
D
3
!
P
R
I
N
TU
I;一十
■
;
:
N
E
X
TI:
P
R
川TU
I
2
0
0F
O
RK
-
-
JT
DC:
R
=
R
+
E
L
K
)
:
6
8
5
日
83
0
0
2
0
1P
R
I
N
TU
S
I
N
G
-
1
1
1
日l
■
;
C
t
K
l;
.(
-
;
!
P
R
I
N
TU
S
I
N
6
1
日I
I
.
l
'
抑 K
)
/
N
I
I
O
O
+
,
0
5
;
Il一;2
0
2P
R
I
N
TU
S
I
N
8
1
日
Jl
f
l
l
;
R;
Ll■;:
P
R
I
N
TU
S
I
N
G
■
日日.
l
-
;
R
/
N
I
1
0
0
◆
.
0
5
!
N
E
X
T
K
2
0
もP
R
I
N
T
L
事;一十
■
;
;
F
O
RI
=
l
T
D
3
王
P
R
川TU
t;
一十
一
;
:
N
E
l
TI王
P
R
I
N
TU
事
2
0
7P
R
I
N
T
' コ一つケイJ
-
;
:
P
R
I
N
TU
S
I
N
6
一
書
榊J
I
l
■
;
N
;
'
Ⅰ
'
;
Il
o
o
.
O
l
l' -52-2
0
9P
R
川T
L
I
;
.
⊥
■
;
!
F
O
RI
=
l
T
O
3
:
P
R
I
N
TU
I
;
'
上
.
;
!
N
E
X
TI:
P
R
I
N
TU
I:
6
8
S
U
D5
0
0
2
1
0R
E
㌃・一一一【とかγラ
ム
】
2
2
0P
R
I
N
TC
H
R
f
f
占
)
;
-【
.;A
f
f
J
)
;
一
】/とかrラム●:U
2
事
-
-
S
T
R
I
N
6
8(
'JI
5
0
)
2
2
2P
R
川T
L
I;'T'
川2
事
!
F
O
RK
-
-
1T
OC
:
6
0
S
U
83
0
0川事
=
S
T
R
川6
引.--
I
C
I
K
‖:
P
R
川T
N
t
2
2
bN
E
Z
T
K
:P
R
川T
Lt'I■.
L
liU2事:6
0
S
U
B
5
0
0
2
3
0R
【ト- 【
l
L
t
か zl7
17
'rラフ】2
3
1P
R
I
N
T
C
H
R
I
t
も
I
;-【●;A
I
L
J
Hl
】II
イ
t
I
rz
17
㌧7
●rラ7
.ミ
P
F
uN
T
L
I;■Tl
川2
事
2
3
IR
=
0
!F
O
RK
:
IT
OC
:
R
=
R
+
C
t
K
)
!
6
0
S
U
B3
0
0
:
H
事
=
S
T
R
I
N
6
引7-1
I
R
HP
R
I
N
TH
I
2
3
もN
E
X
T
X
!P
R
川T
L
I
-
,
'
⊥
l
;
U
2
事!
6
0
S
U
B5
0
0
!
G
D
T
O1
7
0
3
0
0P
R
I
N
TU
S
I
N
6
■
l
日I
l
榊T
;
S
fJ
,
4
)
+
D
I
K
-
a;I+.;3
0
1P
R
川TU
S
I
N
G
l
I
日I
I
I
暮
]
;
S
L
J
,
4
)
+
D
I
K
;
■
1
'
;
:
R
E
T
U
R
N
5
0
0P
R
I
NT一㍍:1㍍ラl
O
K
)ヲ祖叩イ
1
5
0
4G
E
TG
I
:
I
F6
事
く
)
-
0
'T
H
印 5
0
4
5
0
8R
E
T
U
R
N
データを並べ か えたあ と,文番号1
1
0
か ら1
2
2
で分布統計量 を求めて2次元配列S
(J,K)に格納 してい る。 それぞれの測度 は次の とお りであ る。 算術平均X
-
請
,
?.
標 準偏差 S= 変動係数 cu-
且 ×100(%) x -メジアン (Me)刀が偶数 の時 Me-(xn/
2+
xn'2÷1)/2 nが奇数 の時 Me-x(n+1)/2 ただ しxl≦x2≦・・・・・・≦x. 画面m- 1 分布統計量一覧(
2
)
度数分布表 とグラフ 文 番号1
7
0
か ら1
8
9
で度 数 を カ ウン トして い る。 カウン トの方法 は,最大値 と最 小値 の差 (範 囲, レソジ)を クラスの数で割 って間隔を求め, 各 クラスの下限以上,上限未満 に属す るサ ンプル を IF文 によって判断 させ るや り方 である。クラス の数 は2
0
を最大 に した。2
0
までな らい くつで も よい し,最 も適切 な度数分布が得 られ るまで何度 で もクラスの数 を変 えて度数分布表 を作成 で きる ように した。 ここでは度数分布,相対度数分布,累積度数分 布,相対累積度数分布 を同時 に表示す るよ うに し た 。 文番号2
1
0
か ら2
2
6
で ヒス トグラムを作成。 文番号2
3
0
か ら2
3
6
で累積度数分布 グラフを作 成 してい る。 【ミン
t
L
t
1
9
7
8
,シ
ー
ン
コ
,
1
9
8
0
】J
フ■
ン7
● W事ンfb?
ヽ
ン
ー
yJ
l ヽケ
LンZー●?サ TLおっ
事
J
r7
)
ーィ
,'付†
イ
シ
ユ
ワIイ4
4
5
7
3
8
.
I 3
8
5
9
1
3
.
i
8
ふ
.
5 1
9
0
9
5
9
.
0_3
1
1
4
2
1
.
0 2
ふ
2
0
2
1
3
.
0
-
rJtLt2
占
7
0
0
.
3 2
8
5
1
L占
l
o
b
.
7
1
0
1
4
2
.
0
1
9
‡
7
3
.
0 1
9
7
5
3
8
.
0
yr1
7
7
9t7
5
0
8
.
2
7
占
7
8
.
7
1
0
2
.
2
3
1
0
0
.
0
5
5
弘.
0 5
5
4
3
7
.
0
□ウ
シ
ー
ン
刀yt5
7
一
0
.
A
773
3
.
I
1
3
1
.
7
2
0
1
7
.
0
3
5
8
8
.
0
5
1
8
2
5
.
0
シ
'
ココつ2
1
8
7
5
,
8 2
2
1
5
5
.
A
9
0
.
2
帥8
1
.
○
托9
0
7
.
0 1
1
5
8
%.
0
1
9
-
5
1
.
ユ
コ
ウ
7
ム
0
8
.
7
占
8
9
8
.
9
9
0
.
占
l
占
9
2
.
0
5
0
7
8
.
0
3
2
5
9
7
.
0
占
5
+
>
■
ココウ2
2
4
7
.
7
1
5
1
2
,
0
占
8
.
5
7
4
7
.
0
1
7
8
7
.
0
8
砧7
.
0
-53-画面III- 2 変数 とクラスの指定