社会調査データ分析のための汎用プログラム : HARA 007-800

全文

(1)社会調査データ分析のための汎用プログラム ⅡARA. -. 007-800. 原 A. -. 純. M:ultiobjective Computer. 輔*. Program Junsuke. for Survey. Data. Analysis. HARA. SUMMARY It is. a commonplace the socialscientists is desirable that everyone. For. or. OSIRIS The. in U.S.A. auther With. made this. to. processand are. who. analyze. not. proficient. data. survey. in computer. can. by. the. computer.. programming, as such. utilize easily program packages it is not yet realized in Japan. for elementary program short computer. it SPSS. But, a. (HARA. data. survey. program 007-800), compute and tables, (2) crossof correlation coe氏cients and (1) simple-tabulation types tabulation tables, or (3) double-cross-tabulation tables of various of data by This is in JIS Fortran. program only exebanging option Cards. written. analysis･ some. one. can. display. sorts. 目 1.ほ. じ. 次に. め. 2.プログラムの概要 2-1.演算の種類と作業の内容 2-2.入力. デ. -. タ. 2-3.結果の出力 3.デックの構成 3-1.プログラム部分 8-2.コントロール･カ-ド部分. 4.出力の例と見方 5.プログラムの内容. 1.. は. じ. め. に. 社会調査(統計調査)の集計や分析にコソビュータ(大型電子計算横)を用いることほ, *社会学教室(Dept.. of. Sociology).

(2) 43. 社会調査デ-タ分析のための汎用プpグラム. 今日でほ,既に常識であるといってよい｡しかし,データの集計･分析のために,誰もが自在にコソビュ-タを駆使できるという段階には,残念ながら至っておらず,計算会社に馬鹿らしい程の料金を支払って,計算を依頼しているのが実情である｡その点で望ましいのほ,種々の分析が可能なプログラム･パッケージが各地の計算壊セ SPSS, OSIRIS等のプログランターに備えられることである｡北米の社会科学界でほ, OSIRIS等を使用しム･パッケージが開発され,広汎に利用されている｡一度, たことのある老にとってほ,いかにこの種のプログラム･パッケージが便利であり,調査データの分析のレベルを向上させるものであるかは,身にしみて感ぜられるところである｡ SPSS,. 280-60)のため日本では,京都大学共同利用大型計算機セソターの使用機種(FACOM SPSSシステムを導入されるとともに,調査データ分析用のプログラに,三宅一郎氏が, (使 TAPを開発された1)o東京大学共同利用大型計算横センタム.パッケージSDAP,. 用機種HITAC. 8800/8700)でも,この種のプログラム･パッケージ導入の動きがあるよ. うだが,未だ実現していない｡筆者ほ,このプログラム･パッケージ導入の前段階に対応するものとして,社会調査デ 007-&00を作成し. -タの最も初歩的な集計･分析を行なうための汎用プログラムHA耳A たoもとより,この種のプログラムほ,多少の知識があれば容易に作成できるけれども, それにしても相当の時間を要するので,各人がやっていたのでほエネルギーの無駄でもあり,筆者なりの工夫も加えてあるので,その内容と使用法とを御紹介したい2'｡ 2.. プログラムの概要. 2-1.演算の種類と作業の内容コントロール･カードの指示によって. (1)単純集計(1変数) 作表〔実数,百分率〕日次の印刷. (2)クロス集計(2変数) x乞償およびX乞検定,クラマーの連関係数(￠係数),グッドマ1/-クラスカルPl暁位相関係数(r係数),積率相関係数作表〔実数,百分率,カテゴ1)一別のX2検定〕目次の印刷. (3)二重クロス集計(3変数) 作業の内容ほクロス集計の場合と同じのいずれかを行なう.但し,. 1回の演算で,. (1)-(3)のうちの2種額以上を行なうこと. ほできない｡また,クロス集計,二重クロス集計の作表ほ,省略することもできる｡.

(3) 44. 原. 2-2.入. 力デー. 純. 輔. タ数. 変 XI 1. 易品畢-････品 k[名草-････易. X3B---･･-----･･-･･････---X3,n. C.f. Xn3------･-------･･･--Xn,a. C]. 1. X28--･･---------.･･････-･X2n. C.7. 雅. X18--･--････---------･･･Xl仇. CI. 1. 1. 3･･-･･‥. ′レ. X8--･--･---･･-･･-･--･･--X,a. l. 2. サンプ. X2. (1)サンプル･変数 HARA. 007-800の入力データは,上図のようなマト1)ックスの形をしている.このマ. 1サソプルは, トリックスの行･列を,それぞれくサソプル〉およびく変数〉と呼ぶ｡調査単位に対応しているが,総サンプル数は99999個を, 1サンプルあたりの変数ほ800. 1. 個を超えてはならない｡. (2)変数番号･変数名変数ほ入力された順序に,. 1,2,.3,-･とく変数番号〉がつけられるo. この変数番号によって各変数を識別する｡なお,変数ほ,調査における質問項目にあたり,. プログラムでは,. 10字以内のく変数名〉をつけること. ができる(つけなくてもよい)0. (3)変数値･カテゴリー名各変数のとる値. Xi5,つまり. く変数値〉は, 0-9の1桁の整数でなくてはならない.. 10以上の整数がデ-タに含まれていてもよいが,計算からはまったく除外される. なお,変数値は,調査における各質問項目の回答カテゴリー(選択肢)にあたり,. 5字. 以内のくカテゴリー名〉をつけることができる(つけなくてもよい)0 外 (4)除値 D.K.,. N.A.,非該当など,作業には含めたいけれども,. X皇検定あるいほ相関係数の計算からは除外したいという場合のカテゴリーにほ,変数値ほ必ず9を用いなければならない(3-3の〔7〕を参照)o (5)入力のタイプデータの入力は,以下の二つの方法のいずれかによって行なわれる｡ (i)データ.･カードの最初から最後まで, 1コラムずつ,すべて変数として入力する方法(以下,くAタイプ･データ〉と呼ぶ). (ii)データ･カード上の特定のコラムのみを,. FORMATによって指定して入力する. 方法(以下,くⅠタイプ･データ〉と呼ぶ)0 2-3.結果の出力. 350表あるいほ200ページ以内.それをこえる計算の要求がでていても,何ら支障はないが,規定量以上の演算は行なわれない｡.

(4) 45. 社会調査デ-タ分析のための汎用プログラム. なお,このプログラムでは,各ページになるべく多くの表が印刷できるように,不用のカテゴリーを詰め,クロス表ほカテゴ1)-数の少ない方の変数を耗に印刷するようにしている｡. 3. 以下でほ,. デックの構成. Aタイプ･データの場合を中心に説明し,. Ⅰタイプ･データでやり方が異な. っているときにほ,その旨を注記するoまた,数字をカ-ドに′ミンチする(打つ)場合にほ,指定されたコラムの範囲の右端に寄せて.1ソチしなければならないことほ,いうまでもない｡. 3-1.プログラム部分〔1〕メイン･プログラム(HARAbo7-800) HENKAN). 〔2〕サブルーチソ(SUBROUTINE IA. 1サンプルあたりのデータほ,. (A)という形で入力されている(但し,. nほ変数番. 守)｡ここでは,二つ以上の変数を統合したり,変数を別の場所に移したりする作業を行なう8)o (同一変数内でのカテゴリーの統合等の変換作業ほ,ここでもよいが,後の〔7〕の. ⑧で行なう方が簡便である) 例1 .‥.★....1....書....2....斗.‥.3,..‥半....4....書....5. HEN監A甘( IA ) ェA(801 ). stJBROUT=NE DIMEINS. ION. ェⅩ‡9 =F =F. IA IA. ェA. 12. 5Q .EQ,7｡Am･=A ′≧0 LLノヽノ 0 i =IX. .苫Q.7.Am.=A. 壬i.I.丑.6i 51 51. .LE.2. =Ⅹ=2 ェⅩ=1. R丑TURN ENI). 上の例では,変数50の億が7であるもののうち,変数51の値が0,1,2の場合にほ 1を,. 3,4,5,6の場合にほ2を与え,新たに変数120を作っている(なお,上記のいず. れにもあてはまらないものについては,変数120にほ9という値が与えられている)o (注意)サブルーチソを用いる必要のない場合でも, ●.''*....1....*....2....首....ち.. IA) 宝冨器冨冨…諾E誓芸T8K8T() RE℡URN Em. だけは,デックに必ず入れておかねばならない｡ 3-2.コントロール･カード部分〔A〕 /*カード〔4〕演算の種類,作表の有無,百分率の計算法:カード1枚 ③. 演算の種類(コラム1-5).

(5) 原. 46. 純. 輔. ....#....1.. 単純集計A. S=肥A. 単純集計B. SI肝L. 指定された変数についてのみ単純集計する.. クロス集計A. CROSA. 指定された1変数と他の一連の変数とのクロ. 一連の変数すべてについて単純集計する｡. ス集計をする｡クロス集計B. CROSS. 指定された2変数間のクロス集計をする｡. クロス集計C. M皿=Ⅹ. 指定された複数の変数間のすべての組合せについてクロス集計をする｡. DOtlBIJ. 二重クロス集計 ⑧. 作表の有無(コラム7-ll) 作表を行なう場合のみTABLEと打つ｡. ⑧. 百分率の計算法(コラム1S-16) 従属変数の分布(百分率)だけでなく,タテ･ヨコ両方について百分率を求めたい場合にほ,. BOTHと打つ.. 〔5〕集計する変数の指定 (1)単純集計A. (SIMPA). (2)単純集計B. (SIMPL). (4)クロス集計B(CROSS). の場合(コラム1-,. 4コラムずつ):カード任意枚. (5)クロス集計C(MTRIX) 〈単純集計A〉の場合にほ,集計したい一連の変数の最初と最後の変数番号を交互に打っていく.0. く単純集計B〉の場合には,集計したい変数の変数番号を打ってい一く｡くクロス集計B〉の場合には.,集計したい独立変数,従属変数の順に,交互に変数番号を打っていくo. もし,独立変数,従属変数が,その一つ前の独立変数,従属変数. と同じときには,そこはブランクにしておいてもよい(→例3)0. くクロス集計C〉の場合にほ,組合せたい変数の変数番号を打っていく｡但し,変数の数は, 25を超えないことが望ましい｡いずれも,カードほ何枚になってもよいが,一番最後の変数番号の次の4コラムに9999 と打っておく｡例2. ････*･.･･･1-･･★････2････半･･･.う･ 1 5 7 129999.

(6) 社会調査データ分析のための汎用プログラム. 4?. 上のようにカードに打たれているとき, SIMPAならば,. 1,2,3,4,5,7,8,9,10,ll,12の単純集計が,. SIMPLならば,. 1,.5,7,12の単純集計が,. CROSSならば,. 1×5, 5×1,. 7×12のクロス集計が,. MTRIXならば,. 7×1,. 7×5,. 12×1,. 12×5,. 12×7のクロス集計が,. 行なわれるo 例3 ....'*.‥.1..‥★....2....羊....ち. 1 2 う9999 1×2,. CROSSの場合,上のようにカードが打たれていれば,. 1×3. のクロス集計が行. なわれる｡. (3)クロス集計A(CROSA)の場合:カード任意枚 1枚のカードに, ③. 4コラムずつ用いて,以下の内容を打っていく｡. 独立変数の変数番号(コラム1-4) 独立変数と組合せてクロス集計をしたい一連の従属変数の,最初と最後の変数番. ⑧. 号(コラム5-8,. 9-1琶). 例4 .…昔....1‥..甘‥..2. 2 1 4. ′. 上の例では,. 2×1,. 2×2,. 2×8,. 2×4のクロス集計が行なわれる.. この次に,別の独立変数あるいほ従属変数について集計する場合にほ, 同様に打っていけばよい｡. ,カードを改めて. 最後のカードの次に, ....*....1.. 9999. というカードを1枚つける｡ (6)二重クロス集計(DOUBL)の場合:カード任意枚 1枚のカードに, 4コラムずつ用いて,以下の内容を打っていく. ③. 第3変数の変数番号(コラム1-4). ⑧. 第8変数のカテゴ1)-の数(コラム5-8). ⑧. 第3変数の各カテゴ1)一に対応する変数値(コラム9-,. 4コラムずつ,不定個). ④クロス集計を行なうべき独立変数の変数番号(コラム不定,. 4コラム). ⑤クロス集計を行なうべき従属変数の変数番号(コラム不定,. 4コラム). 例5. 要義晋三表変数値 ←-㌔-1. ′ー----ヽ. 要義墓褒 .∼. ....*....1....★....2....半....ち. 1 2 11ううう 25. r∼. う1.

(7) 48. 原. 純. 輔. 上の例でほ,第3変数113の変数値1と2と3に対応する三つのカテゴ1) 独立変数25. -に関して,. と従属変数81のクロス集計を行なう｡. この次に,別の変数について集計する場合には,カードを改めて同様に打ってい桝fよい｡. 最後のカードの次に, ....*....1.. 9999. というカードを1枚つける｡ 3-3.データ部分〔6〕調査名,. 1サソプルあたりのカード枚数,入力のタイプ:カード1枚. ①. 調査名(コラム1-20). ⑧. 1サソプルあたりのカード枚数(コラム28-30) 1-10. (800変数以内という制限があるから) 但し, Ⅰタイプ･データの場合には1サンプルあたりの変数の数を打つ｡ (参入力のタイプ(コラム80) AあるいほⅠと打つo 〔7〕変数名,変数値の変換(変数値9の取扱い),カテゴリー名:カード任意枚 1枚のカードに,以下の内容を打っていく｡ ①. 変数番号(コラム1-4). ②. 変数名(コラム6-15). ⑧. 変数値の変換〔変数値9の取扱い〕,カテゴリー名(コラム21-,. 10. コラムずつ,不箆個) カテゴ7)-を統合したり,変数値を変換したりする場合,あるいほ,カテゴ1)I 名を指定する場合,または,その両方を行なう場合に用いるが,それらの必要がまったくないときには省略する.カテゴリー名を指定しない場合には,変数値がその t'-0-t'-b-8-b ままカテゴ.)一名として自動的に入る(. )｡. およびb-Ⅹ-b. 1コラム目:変数値 2コラム目:. 1コラム日に指定した変数値を別の値に変換したい場合に,. *(アス. テ7)スク)を打つ.. 3-4コラム日:新たに変換したいと思う変数値 5コラム目:. 1コラム日に指定した変数億(3-4コラム目で新たな変数値に変. 換したときには,その億)にカテゴ1)一名を指定したい場合に,. *(アステ. 1)スク)を打つo. 6･-10コラム日:カテゴリー名 ⑨のスペースが,カード1放では足りない場合にほ,次のカードのコラム21以下を用いて同様に打っていけばよい｡その場合,コラム1-芝Oは,. 1枚日とまっ.

(8) 49. 社会調査データ分析のための汎用プログラム. たく同じ内容を打つか,あるいは,まったくブランクにしておく｡ (変数値9の取扱い) 変数億9は,通常は,作表にほ含められるが,. X2検定あるいほ相関係数の計算. からは除外される｡これを除外したくないときには,ここで変数値を-9という値に変換してやればよいo. カテゴ1). -名ほ,特別に指定しなければ,自動的にb-9-b. が入る｡. (注意) ′二つ以上の変換を組合せたり,変数を別の場所-移したりする作業は, 〔2〕のSUBROUTINE. ⅡENKAN. め中で行なわなければならない｡ここでは同. 一変数内での変換作業に限られる. HENKANにおいても行なうことがでなお,変数値の変換ほ, SUBROUTINE SUBROUTINE HENKAN きる｡但し,同一の変数に関して両方に命令があると, →⑨という贋序で実行される(いずれか一方の命令のみ実行されるわけではない) のでo注意が必要であるo. (注意)カテゴ.)一名としてブランク(. )を指定することはできな. bbbbt･. い｡. (注意) 合には,. 9以外の変数値を,. ⑧あるいはSUBROUTINE. X2検定あるいほ相関係数の計算から除外したい場 HENKANにおいて,その変数値を9に変. 換すればよい｡. 特定の変数償を,作表からも除外したいときにほ,その変数値を10以上に変換してやればよい｡例6 ....1*....1‥‥*....2....★....う‥‥書..･‥ヰ..‥*.‥.5....*....6.. 12. NEmE=. 2半1. う★ 1*SEFGO4*. 上の例6ほ,調査対象者の年令(20才-59才)について, 代にほ5という変数値が与えられているデータを,戦前生れ(40, 戦後生れ(20,. 0. 5*. O*SENZ甘. 20才代には2,--,50才 50才代,変数値0)と. 30才代,変数億1)にまとめ,変数名とカテゴリー名を指定するための. ものである｡なお,二つ以上の変数値を同一の変数値に変換して,しかもカテゴ7)一名を指定したいときにほ,カテゴ1)一名ほ,そのうちの1個所で指定するだけでよい.. 2個所. 以上で同一のカテゴ1)一名を指定しても,勿論,支障ほないが,誤まってBljのカテゴ1)I 名を指定した場合には,最後に現われたものがその変数値のカテゴ1)一名となるo (変数番号に関係なく)どのような順序で入れてもよいが, このカードほ,何放でも, 最後にほ必ず(1枚も入れない場合でも), .‥.音‥‥1. 9999. というカードをつけるo. (注意)変数名,カテゴリー名の指定ほ,必ずしも必要はない｡省略した場合の変数名.

(9) 原. 50. はブランクになる.. 純. 輔. ㊨, ⑧の作業がまったく不要な場合にほ,. ③のみを打ったカードを,. わざわざ入れる必要ほない｡ (〔8〕はⅠタイプ･デ-タの場合のみ必要o. Aタイプ･データの場合は〔9〕へ). 〔8〕データ入力のFORMAT:カード3枚データ入力のためのFORMATの内容をパンチする.その際,両端のカtyコを忘れてはならない｡. カ-ドが3枚必要でない場合でも,ブランク･カードを添えて,全体として3枚入れる｡タ本体〔9〕デー. 最後に,. 1サンプルあたりのカード枚数だけブランク･カードをつけ,その冒頭に, ....★‥‥1｡. 9999. と打っておく｡ (注意) Ⅰタイプ･データの場合にほ,第1-4変数の変数億がそれぞれ9となるように打つ｡〔10〕 //EⅣDカード. 4.出力の例と見方付表1は,二重クロス集計の結果(部分)および日次の出力の例である｡各ページの1行目には,プログラム名(HARA. 007),演算の種類(DOUBL),調査名 CHOSA),サンプル数(787),ページ数が印刷される.. (A班-SM. 結果のうちで, ⅩSQほ, CR. X乏の値(DFは自由度)4). (SQRT)は,クラマーの連関係数の平方根. GKほ,グッドマソニクラスカルの順位相関係数 Rは,横率相関係数石) を,それぞれ示している｡. X乞検定の結果ほ, X皇値の頭に符号で示されるが, ≠は,危険率1%未満で有意差あり +ほ,危険率1%以上, NSほ,有意差なし. 5%未満で有意差あり. 杏,それぞれ意味している｡クロス表の中で,実数の頭にプラス(+)やマイナス(-)の符号のついている部分がある｡これほ,その桝に注目して,残りのカテゴリーは一つにまとめ,四分表を作成し, X乏検定を行なった結果を示したものである｡プラスは,実現値が期待値より大きいこと杏,マイナスは,実現値が期待値よりも小さいことを意味しているが,危険率の表示は, 上の場合と同じである(但し,. NSの表示は行なわない)o.

(10) 社会調査データ分析のための汎用プログラム. 5.. 付表2は,. HARA. 51. プログラムの内容. 007-800のプログラムである｡プpグラム中での作業の内容は,お. およそ,以下のとおりである｡作. カード番号より. 業. の. 内. 容. 配列宣言. 7｡｡｡7｡｡7認諾7演算の種類および初期値の指定〔4〕クロス集計B (CROSS)の変数の組合せ指定単純集計B (SIMPL)の変数指定二重クロス集計(DOUBL)の変数の組合せ指定クT2ス集計C (MTRIX)の変数指定 -67-. 単純集計A. 707--. クロス集計A (CROSA)の変数の組合せ指定調査名の読込み〔6〕変数名,カテゴリー名の読込み〔7〕. 〔5〕. (SIMPA)の変数指定. ｡芸芸耶闇78-774. データの読込みと集計〔8〕〔9〕集計表の整理(カテゴリーを詰める) x2償およびX2検定タラマ｣の連関係数環率相関係数グッドマソニクラスカルの慣位相関係数行数,ページ数のチェック. -6器機朋. 百分率四分表によるX皇検定集計結果の印刷日次の印刷. 註 19730 1)三宅一郎『社会科学のための統計パッケージ』束洋経済新報社, SOPHIA IIR 2)このプログラムほ,故佐々木伸夫氏(当時一橋大学大学院生)の努力により, OID という名称で,上智大学国際関係研究所のプログラム･ライブラ1) -に登録されているo このプログラムを用いた調査･研究で,結果が公刊(予定)されているものにほ,例えば, 以下のものがある｡青井和夫他｢主婦と就労の調査｣, 『第45回日本社会学会大会報告要旨』, 1972, 1972｡所『ジュリスト』499,有斐閲, 一彦･前田俊郎｢刑事事件処理と損害賠偉｣, pp. 純輔｢現代人と｢立身出世｣｣, 『現代社会学』2,講談社, 1974刊行予定｡ 28-36｡原 1973, 松村良之｢損害賠償に対する人々の意識｣,川島武宜編『法社会学講座8』岩波書店, pp.. 3). 257-720. Aタイプ･データの場合にほ,デ-タの中に,数字以外の文字がまじっていてもさしつかえない(Ⅰタイプ･データの場合には許されない)｡それらの文字のうち,特定のものに対してほ, 以下のように変数値が与えられる..

(11) これらは,いずれも値は10以上であり,このままでは,計算から除外されてしまうので, 9以下の数値に変換してやる心裏がある. 用いたい場合にほ,このサブル-チソの中で, ちなお, (プラソク)ほ,常に0と読まれることは,いうまでもない｡ 4)このプログラムでほ, X2検定の検定力を高めるために, 2変数が独立であれば,その尤度比ス 1ほ, に関して, -2loge1が近似的にX乞分布する,という定理を用いている(ここで, }=. Hut.i.1In.nj'] nnllnin,･ij. 196も p. 57を参照｡である)｡岩原信九郎『ノンパラメトリック法』日本文化科学社, (カテゴリーを詰めて整理した)クロス表の,第1列(第1行)目のカテゴリーから 5)ここでは, から順に1,. 2, 3,. --と整数値を与え,これを間隔尺度債と仮にみなして積率相関係教を求め Boyle, ている(与えられる整数値が,心ずしもデータの変数億とは一致しない点に注意せよ). Richard P., E`Path Analysis Ordinal Data", AJS, 75-4, 1970, and pp. 461-80を参照｡訂正〔付表2〕の中で･カード7156と7156の間に,. GOTO19というカードを追加する｡.

(12) 53. 社会調査データ分析のための汎用プログラA 〔付表1〕 ‖)uUaL). HA机007. AM-SM. 5. ヰ. TABLE. xミ)Q= 5. #. 我. 9. O. (. HIGH. 11.4. 7ち.4 49.a). (. 42.9). 12. 12. 8.5 (. 8.5. 25.5). (. 0.0. 'I. 12.8. (. 100.0 56.7) 1L+1. 47.5. b.0). 0. t. b7. ++. 5.0 (. 49.4). 0.a. 0/.O). (. 0.0). 51. T【汀AL. (. 91.8). 100.0 40'4). 1. 0. O.0). (. 87. 0.O 0.0). (. loo.0 1.4). 11/. /う1. b.0. 24.9. 33.5. 20.9. (loo.0). (loo.0). (1OO･0). (100･0). 14.も. AM-SM. (D□UBL). 0.0. (. 21. IIOO.0). 007. .5.1 5,5). (. 22･6). 1. .0. く. HARA. (. (. 4. --. 44.4). --. 30.5. 57.I) 0. -X-. (. 45. +. 1.3 1.Ll). 1OO.0. 40.b. 40.2). (. /9. 1. (. -0｡470. 701-^L. 一-. 52. +. 27.う. 7.a. >4.9〉. --. 55. GK=. -X-. '58. flo.きIー(. 9. '21.9 l. ++. 0.0. 28. ◆◆. MIDDL. 9. 0.0). i3.9. くど1.6)I. 0.402. BL-⊂L. --. ). (GENSY□KU. CR(SQRT)=. uH-⊂L. ll. .. a). SN-KN. LOU. 20S. )港'/0. (GAKUREKl. 89.303(DF. ++. GYOSY. '/0. 9. ‡e. ). (NENRヒl. 781[. C目口SA. PA. TABLE. Y n)0 K Y K. ヽ/ ヽ■′. 0 く E N 0 ■1らE N. ー. 0･∋) 549 100.0. ′u. ど一ュn 1 1 1 P 1. nヰ5も7. ▲■...ヽ. ヽ■一. ■I′. ー. 〓)‖)‖)〓〉〓-‖▼U. _■.′. 0 o EE∼S Y 00K 0 く N Y 氏 S55. A KK R E 氏 A 氏E 吃. ー l ー. G′q.D. (. R EK R EK. o く o. ∼S Y 0 吃 N らY C 紘 N S Y0K. E∈E ′りqG (.ー. ヽ■′. ′一. リ〓▼. I_一■. ー AK AK. 7777777. JIl. ′.ヽ. 〓u‖∨. }. 敬.v常溝沸落水沖. l. く. ′■ヽ. NE N R tn くN E Z 良. ヽ.′. ヽ一. L9G66. K NZ E. ー. ･′■ヽ. nう.Y)9q･99 777. ′.■ヽ. F 吃l一 R 紘王. 常北ハ瀕也ハ)祁)河じ河. ■■ヽ. A臣E. ′.ヽ. (N Z. FK. J)lt▲l▲III R=R JJヒじし〓u〓V〓}RR. A 吃 A 吃. G′u.D. n.q52tノ.I //′′/// 99q-Eノ仁一'Jt(A. 12うヰ5.07. ( (. ●l. (. (100･0.). 787. C‖OS^. 10O.0. ヽ一. ド. ー. 1. R=. -0.2t55b.

(13) 原. 64. 〔付表. 紙. 輔. 2〕. ⊂.‥#‥..1..‥韓‥‥2....戟.‥.3‥..紫‥ ,4.‥.鞍‥‥5..‥耗‥..8.‥.栄. ⊂ ⊂ROSS TABUしATION (-HAR人 007-800) DIMENSION IA(801)'tB(801)'lC(10I801)IN(11†11)IMD(801)1NF(ヨI3) DIMENS王ON P(ll,ll)'Q(11'11)'SG(2I30)IXX(llIll)▼蛸Z(801). DⅠMENSION lNTE6ER REAし C...i. ｢=JU. 4. .･.1....手壬‥‥2‥..手書.‥.3..‥准 READ(5†100) ZAME'ZTBL,之PER IF(ZAME.EO.4日IMPA) KAME=O lF(LAME.EQ.4H!MPい KAME=1 IF(ZAME.EQ.4日RDSA) KAME三2 lF(LAME.EG[.4HROSS) KAME=3 IF(ZAME.EQ.4HOUBL) KAPE=ヰ IF(LAME.EQ,4日TRIX) KAPE=5. IF(KAME.LT.0.OR.KAME.6T.5). KTBL三1 1F(ZTBし.NE.4日TABL) KPER=o IF(ZPER.EQ.4HBOTH) NTA巳ニ0 Lり( 1)ニ8H しIJ( 2)=8日 LIJ( 5)=8日 L王J( 4)=8日しIJ( 5)=8日 LZJ( 6)=8日 7)=8H LIJ( LIJく 8)=8日 9)=8日 LIJく. KTB｣王O KPER;1. ). (. ). 1. l. Lり(10)=8日. O12TrL4L^′D78. f. l. l. L. I. l. L. †. I. ]. 】. -. 一. -X-. D【I 2 K=1I8ol MD(K)=O N5ヱ(K)=O A(llK)=8H A(2IK)--8H DO 1 J=1IIO 1 MrJ(JIK)=LIJ(J) 2 ⊂ONTINUE IF(仁AME.【Q.1) IF(KAME:.EQ.0) SG(1▼ い=b.635 SG(1I 2)=9.210 S6(1I 3)≡ll.345 5G(1I 4)=13.277 SG(1▼ 5)=15.o8b SG(1I占)=16.812 SG(lI 7)=18.475 S6(1? 8】=2O,090 SG(1I 9)--21.666 0)=23.209 SG(1II 1)=24.725 SG(1II 2)=26,217 S6(111 3)=27.占88 SG(1ll SG(III 4)--29.1年1 5)=3C.578 SG(1II 6)=52.000 SG(1II 7)=33.409 S6(1II 8)=5/i.805 SG(1Il. Gロ TO GO TO. 507. leo. ｢. GO. TO. 1. LD ー一〇 0) l A ー '=hⅥ山. REAL兼8. IN(10,10,350)'MQ(2'350)'NQ(2I350))NPX(550). BGJBPIBR ND(a)IID(801),KAK(ll) LIJ(10)'MIJ(10 801)IPROJ. _. 72. 7. ‥載.‥.8 日ARA7000 HARA?001 日ARA7002 HARA7003 日ARA7004 HARA7005 日ARA700占 7‥‥♯.‥.8 HARA7007 =ARA7008 日ARA7009 日ARA7010 日ARA7011 日ARA7012 HARA7O13 HARA7O14 日ARA7015 日ARA7016 日ARA7017 HARA7018 日ARA7019 日ARA7020 日ARA7021 HARA7022 HARA7023 HARA7024 HARA7025 ト1ARA7026 日ARA7027 HARA7028 日ARA7029 HARA7050 日ARA7O31 HARA7032 日AFモA7035 HARA7Oヲ4 HAR^7035 日ARA7036 日APA7O37 HAR^7O38 HARA7059 日ARA7O40 HARA7041 HARA7042 HARA7043 日ARA7044 日ARA70与5 日ARA70与b. HARA7047 HARA7048 日ARA70年9 HARA7050 日ARA7051 HARA7052 日ARA7053 HARA705年 HARA7055 HARA705b HARA7￠57. ヲ孟??チ:;モ!!さ孟::ラブÅ!さ!:諾でワ亨.

(14) 55. 社会調査データ分析のための汎用プログラム. HARA7058 日ARA7059 日ARA70bO HARA7061 日ARA7062 日ARA7085 日ARA7064 日ARA70b5 HARA7O66 日ARA7O67 日ARA7088 =ARA7069 HARA7()70 日ARA7071 日ARA7()72 HARA7073 日ARA7074 日ARA7075 HARA707る =ARA7077 日ARA7()78 ト1ARA7079 HARA7080 =ARA7081 日ARA7082 HARA708ヨ HARA708年 HARA7O85. S与 (1ナ19)=36.191 SG (1I20)=37.56& S6 (1I21)芋38,9ヨ2 S(; (1I 22)芦4'0.289 S巧 (1I 23)叫1.638 5G (1I 24)=42.980 S6 (1I 25)=44.314 SG (王I 2b)=45.6年2 SG (1I 27)=46.9古3 S(; (1I 28)=冬8.278 S6 (1} 29)=4･9.588 5 0 892 SG (1I n ∧U s6(2I ヰ9 1)圭5 SG(2I 2)三5 1 SG(2I 3)= 7 SG(2I与.)=9.488 S与(2Iら)=11.070 SG(2I a)=12.592 SG(2I 7)=1年,067 8)=15.507 SG(21 9)=16,919 s6(2I S6(2IIO)=18.う07 SG(･2Ill)=19.675 S6(2†12)=21.026 SG(2I13)=22.う62 S6(2I14)=25.685 SG(2I15)=2年.996 SG(2I16)=2占.29b SG(2I17)=27.587 SG(2I18)=2日,869 5G(2I19)=3O,1(14 SG(2I28)=51.410 SG(2I21)=52.671 SG(2▼22)=33.924 sG(2I2声)=35.17Z SG(2I24)=36.415 SG(2I25)=37.占52 SG(2I2も)亡う8.885 SG(2)27)=40,113 sG(2I28)=41.337 SGく2129)=42.557 SG(2I30)=43.773 DO 5 Kニ1Il1 ( 3 KAK(K)=3H ヽ.∫〕l. 115. 89CO. LLt. rr). ≡-1. l‖DdD. 【川HU. 0∩こU. n)∩】∩U. G6t3G. pr1年LLI. Ln ▲U･q 0 T IF(仁AME.EQ T IF(KAME.EQ 10 7F(KAME.EQ lん8 T T O 5 1 IFくKAME.EQ ∼ 60 TO･72 2‥‥瀦....S‥‥-)f‥‥Lh..‥･)を.‥.5‥..与手....白‥..主を‥..7 C...#....1....#... 50ヰ NAME三8=⊂RO5S 4 READ(5IIOl) (lB(K)IK=1I20) Dロ 5 K=1I19I2 GO TO 19 IF(柑(K).GT.800i 7F(NTAB+1｡GT.350) GO TO a NTA8:NTABkl 【ヨL l●l■ KKl NG)(1INTAB) ● 田山 NQ(2iNTAB) lF(NQ(1INTAB).EQ.0) NQ(1INTAB)=NQ(1INTABI,1) 5 IF(NQ(2INTAB).EQ.0) NQ(2?NTAB)=NQ(2INTAB-1). HARA7086 日ARA7()87 日ARA7088 日ARA7089 HARA7090 HARA7091 日ARA7092 日ARA7093 日ARA709LL HARA7095 HARA709e) HARA7097 HARA7098 日ARA7099 HARA7100 HARA7101 HARA7102 HARA710う HARA710(+ HARA7105 HARA7106 ‥静.‥.8. 日ARA71O7 HARA7188 日ARA7109 =ARA711￠ HARA7111 日ARA7112 HARA711ヨ HARA711L+ HARA7115 日ARA7116.

(15) 56. 原. 6 GO TO ヰ ⊂‥.諌‥..1..‥寺壬‥. 507 NAME=8日SIMPL 7 READ(5●101) DO 8 K=ユI20. 純. 輔. HARA7117 2‥‥我.‥.3‥‥手蔓....4‥‥幹‥‥5‥‥Iモ‥..6.‥.手を‥‥7.. ...#....8 HARA7118 HARA7119. (I8(K)IK買1}20). GO 1F(lB(K).GT.800) TO 19 9 1F(NTAB◆1.GT.350) 60 TO NTAB=NTAB◆1 NQ(2事NTAB)=‡B(K) 8 NQ(1INTAB)=801 9 GO TO 7 ⊂‥.苓..‥1..‥鞍‥..2‥‥手e....!.. 510 NAME=8HDDUBL IO READ(5*101) (IB(K)▼K=1●20) IF(IB(1).GT.800) 60 T【】 19. 日ARA7120 HARA7121 日ARA7122 日ARA7123 日ARA7124 日ARA7125 日ARA712占 .井..‥4..‥寺壬‥‥5‥‥!モ....6.‥.書手.‥.7‥‥栽,.●‥8. 日ARA7127 HARA7128 日ARA712･9 HARA7130 HARA7L31 HARA7132 HARA7133 HARA713年 HARA7135 HARA7136 HARA71さ7 HARA7158. MA王ⅠB(2) Dロ11 K=1IMA 60 TO 12 IF(NTAB+1.GT.350) NTAB=NTAB+i MQ(llNTAB)ニー白く1) MQ(2INTAB)=lB(K+2)+1 NQ(1INTAB)=IB(MA+3) 11 NQ(2INTAB)=IB(MA+4) TO 12 60 10 c‥.輯‥..1....鮮....2‥‥炎....3‥‥†t....4‥‥鮮‥..5‥‥推.. NAME=8日〔ROSS 13 し1≡1 NMTRIX=O 14 L2=Ll+19 READ(5IIOl) (ⅠB(K)lK=LllL2) DO 15 K=LIIL2 IF(IB(K).GT.800) GO TO 16 NMTRZX:NMTRIX+1 15 Ll三Ll+20 GO TO 14 16 DO 18 l=2)NMTRZX 17 DO J=1Iト1 GO IF(NTAB◆1.GT.350) TO 19 NTAB=NTAB+1 NQ(1INTAB)=lB(I) 17 NGI(2)NTAB)=7B(J) 18 ⊂ONTINUE ⊂‥.准‥‥1....常‥..2‥‥米.‥.S‥‥M..‥4..‥米‥‥5.‥.拝‥..6. 180 NAME=8HSIMPL KAME=1 READ(51101) (IB(I)?!=1?20) 181 DO 183 Ⅰ=1I19I2 IF(IB(I).GT.800) GO TO 19 IHA=IB(I) ZHB=lB(I+1) DO 182 K=Ⅰ=A†IHB GO TO 184 IF(NTAB.+1.GT.350) NTAB=NTAB+i NQ(1INTAB)=881 NQ(2INTAB)=K 182 183 CONTZNUE GO TO 184 181 ⊂._.手t‥‥1....常‥‥2‥‥幸手.‥.さ‥‥井....4..‥簾‥‥5‥‥耗..‥8...'韓‥.･7 185 NAME=8HCRDSS -. KAME=5. ′○. #. 7. ..*....8 HARA713? HARA7140 HARA7141 HARA7142 HARA7143 日ARA7144 HAR^7145 HARA714古 HARA7147 日ARA7148 日ARA7149 日ARA7150 HARA7151 日ARA7152 HAF!A7153 HARA7154 日ARA7155. .舛‥‥7..‥溌..‥8 日ARA7156 HARA7157 HARA7158 HARA7159 HARA71占O HARA71bl HARA7162 HARA71b3 HARA71b4 日ARA71る5 日ARA7166 日ARA7167 HARA7168 HARA7169 ..書l....8 HARA7170 日ARA7171.

(16) 57. 社会調査デ-タ分析のための汎用プログラム. READ(5IIO1) lHCIIHA>IHB IF(IH⊂.GT.800〉 GO TO 19 DO K=IHA▼1日B 187 Gロ TO 188 ⅠF(NTAB+1.GT.･350) NTAB=NTAB+1 NQく1INTAB)==1⊂ 187 NO(21NTAB)=K TO 188 GO 186 ⊂‥.兼‥‥1..‥半....2‥..i....3.=.報=･･与･･‥脊‥‥5･･･t貨-･･6････きを-･17･ (PROJ(M)IM芸1I4)IN⊂ARDIZTY 19 READ(5'104) ITYヨO ) IF(ZTY.EQ.4日I lTy=1 NC8L=N⊂ARD相O NCロL=N⊂ARD IF(lTY.EQ.1) DO 20 K三1IN⊂OL 20 ID(K)=ヰH C,..A....1....# .潔....3.‥.栄..‥与....繋‥‥5..=it==8=.･亨さ‥‥7‥-翁････8 .2. M 〓ロ lh 1 '. ′○ H N ′t 21 LO(K)=4日ヨ9年 K=1I801 DO DO さ95 J=1IIO 393 1C(JIK)=J-i 39年 ⊂ロNTINUE 22 READ(5IIO5) MS,(A(K'801)'K=1'2)I(IA(M)'LD(M)'lB(M),NDく･u,LIJ(M),M=HARA7195 11I6) 0) MA=MS IF(MS.GT. IF(MS.GT.BOO) GO TO 23 IF(MS.EQ. 0) GO TO う95 A(1IMA)=A(1I801) A(2IMA)=A(2I801) 395 DO 598 M=1Ib MI三IA(M)◆1 TO 397 1F(LD(M),NE.4H書手 ) 60 MD(MA)ニ1 GO TO 39b 1F(lB(M).NE.p9) 186. (]. K･卑. 12K ･<OhU. NSZ(MA)=1 MIJ(10IMA)=8日. -9-. 王8(M)=9 王亡(MIIMA)=lB(M) MⅠ=Ⅰ8(Ml●1 TO 398 597 1F(ND(M).NE.4H米 ) GB MIJ(M7!MA)=LtJ(M) さ98 ⊂ONTINUE 60 TO 22 ⊂‥.鞍‥..1....凍,.‥2.‥.准‥‥3‥‥姥‥‥ヰ‥‥米‥‥5‥‥妊.‥.8‥.･米‥.･7･ 2b K=1INTAB 2さ DO 25 DO I31]10 24 DQ J=1†10 24 IN(IIJIK)=8 25 CONTINUE 26 CONTINUE NSAM-..o IF(iTY.EQ.1) READ(5IIO7〉 (LDくM)IM=1I占0) 399 27 GO TO IF(IVY.NE.1) READ(5ILD) (lA(M)lM=1IN〔OL) TO 481 60 399 READ(5IiOb) (ID(MIIM言1IN⊂OL) 00 400 M王1INCOL ヨ98. HARA7172 日ARA7173 日APA7174 日ARA7175 日ARA7176 日ARA7177 日ARA7178 日ARA7179 ..i....8 HARA7180 日ARA7181 HARA7182 日ARA7183 HARA7184 日ARA7185 =ARA718b HARA7187 日ARA7188 日ARA7189 日ARA719() HARA71?i 日ARA7192 日ARA719∋ HARA7194 HARA7196 日ARA7197 日ARA7198 日ARA7199 日ARA7200 日ARA7201 日ARA7202 日ARA720う HARA7204 HARA7205 日ARA7208 日ARA7207 HARA7208 日ARA7209 日ARA7210 HARA7211 HARA7212 日ARA7215 HARA721年 HARA7215 ‥雀....8. 日ARA7218 日ARA7217 HARA7218 日ARA7219 日ARA722() HARA7221 日ARA7222 HARA722さ HARA722年 HARA7225 日ARA7226 HARA7之27 HARA7228.

(17) 58. 原. !A(M)=99 IF(ZDくM).EQ.4日 IF(ⅠD(M〉.EQ.OHO. ). 王F(川(M),EQ.4日I. ) ) ) ) ) ). IA(M)=0 IA(M)三0 lA(M)=1 ZA(M)=2 lA(M)=3 lA(M)言4 IA(M)=5. 純. :l▼l. ZF(lD(M).EQ.4日2 1FりD(M).EQ.年目3 ZF(lD(M).EQ.4H4 !F(ID(M〉.EQ.4H5 IF(ID(M),EQ.4Hb ) IA(M)等6 IF(IDくM).EQ.4H7 )川(M)=7 1F ー_-D(M).EQ.4H8 ) IA(M)=8 IF ーH D(M).EQ.4H9 ) lA(M)言9 IF ー I D(M)..EQ.4Hx ) lA(M)=20 l▼ IF ー D(M).EQ.4HY ) ZA(M)=50 IF ( D(M).EQ.4日+ ) IA(M)=40 !F ー DD(M).EQ.4H) lA(M)=50 F ー 0O (M).EQ.4日. ) IA(M)=60 0 1 NSAMX:IA(1川1000+IA(2)#100+tA(3)*10+王A(4) IF(NSAMX.EQ.9999) 32 GO TO NSAM=NSAM●1 ⊂ALL HENKAN(IA) DO L=1I800 29 IF(MO(L).NE.1) GO TO 29 DO 28 M=1IIO 28 IF(IA(L).EQ.M-1) IA(L)=Ⅰ⊂くM'L) 29 ⊂ONTINUE DO 31 K=1INTAB IF(KAME.NE.4) GO TD 30 MA=MQ(1,吃) ZF(lA(MA).NE.MQ(2IK)-1) TO GO ヨ1 50 M【王NQ(1IK) MJ三NQ(2IK) l芦IA(MI)◆1 J=IA(MJ)◆1 ZF(I.GT.10.ロR.J.GT.10) 6口 T〔】 31 1N(lIJブK)=1N(lIJIK)+1 31 CONTINUE GO TO 27 C･･･祥--1････繋--2-･･壬モ.･･･3-･･与与･･-4･.‥きを‥‥5‥‥祥‥..6‥.｡拝‥..7 32 1A(801)=O NP=O NG=与4 DO b8 K=1INTAB DO 34 l=1IIO DO 53 J=1IIO N(lIJ)=1NくIIJIK) 53 ⊂ONT王NUE 3年 CONTINUE DO 36 J=1IIO Nく11I｣)=O DD 35 Ⅰ=1IIO 35 N(11ナJ);Nく11'J)◆N(IIJ) ∋b ⊂ONTINUE DO 38 Ⅰ=1Ill NくIIll)=0 DO さ7 J= lN一百E 0 …一日山◆ ∃ ー= I 37 N '1 l J ｢H一 N(lll 38 CONTINUE NSS=99 IF(N(llIll).EQ.O) NSS=0. 輔. HARA7229 HARA723O HARA7231 HARA7232 HARA723さ HARA7234 日ARA7235 HARA7236 日ARA72∋7 日ARA7258 HARA72ヨ9 HARA724O. H岬A7241. 年与. HARA7242 HARA72年3 HARA724年 HARA7245 HARA7246 日ARA7247 日ARA7248 HARA7249 HARA7250 HARA7251 HARA7252 HAFu7253 HAF!A7254 HARA7255 HARA7256 HARA7257 HARA7258 HARA7259 HARA72bO HARA72bl. H仰A72b2 HARA7263 HARA72b4 HARA7265 HARA7266 ..港....8. 口H一. nH一. 日ARA7267 日ARA7268 HARA72b9 HARA7270 日ARA7271 HARA7272 HARA7273 HARA7274 HARA7275 HARA7276 HARA7277 HARA7278 HARA7279 HARA7280 HARA7281 日ARA7282 日ARA728∋ HARA728年 HARA7285 HARA7286 HARA7287.

(18) 59. 社会調査データ分析のための汎用プログラム. う9. NSZI=1 NSZJ=1. 日ARA7299 HARA7300 日ARA7301 HARA7302 日ARA7503 日ARA750句･ HARA7505 HARA7う06 HARA73O7 日ARA73o8 HARA7う09 日ARA7310 HARA7う11 HARA7312 日ARA7315 日ARA7う1年. HARA7315 日ARA7316. ‡Y三Ⅰ×¶1. HAFiA7317. JY=JX-i. 日ARA7318. IZ=IY JZ=JY IF(NSZl.EQ.i) !F(NSZJ.EQ.1) ⊂‥.手モ. 妄1ARA7319. 日ARA7520 日ARA7521. [Z=lY-1 JZ=JY-1. ...1..‥与(∼.‥.2‥‥幹....3‥..与モ‥‥4‥‥-t･f..‥5‥‥き壬‥‥6...･寺号-･･7 R=0.0 Aし=0.0 Bし=0.0. ⊂し=0.0 ⊂R=0.0 GK=0,0 IDF=O. MXX=O NXX=O XXX=O.r'. lSG=2H !ST=エZ-1 IF(JZ,LT.Ⅰ7.) 1F(IST.LE,0). lST=JZt-i GO TO 55. MXT=O 45 1=1I王Z Dロ 42 J=1IJZ IF(N(7'J).E(a.0) MXX=M)(X◆N(l事J ALニAL+FLOAT(N CONTINUE IF(MXX.EQ.0). DO. 42. GO. To. 42. ). (IIJ))糊LOG(FLOAT(NE(.Ⅰ,J).).); GO. Lとう TO BL;8L+F'LOAT(MXX)粥ALロG(FLOAT(MXX) MXT=MXT+MXX与り NXX=Nxx◆MxX. HARA7288 HARA7289 HARA7290 日ARA7∼91 日ARA7292 日ARA729き HARA7294 日ARA7295 日ARA729b HARA7297 日ARA7298. l]..). l◆【U. L10 41. 〉ハX l--▲ヰ. GO TO 55 tF(NSS.EQ.0) M工芸NQ(1IK) MJ=NQ(2IK) 59 DO Ⅰ三1IIO ⊂lくⅠ)言MI｣くIIMⅠ) ⊂J(i)王MIJ(王IMJ) NSZ!=O NS之JニO IF(N(10†11).GT.0,AND.NSZ(Ml).EQ.O) IF(N(ll,10).GT.0.AND.NSZ.(MJ).EQ.0) CI(ll)=8HTOTAL ⊂J(ll)=8日TOTAL IF(KAME.EQ.1) CI(ll)=8日 IX-_0 41 DO !=1Ill JX=O '41 GロTO IF(N(IIll).EQ.0) ×◆ 1 IX= ⊂ U l くl ヽ/1 〓 D 0 1一1 IF(N(llIJ).EQ.0) GO TO 40 JX=JX+i N(Ⅰ×,JX)=N(IIJ) IF(I.NE.ll) GO TO 40 CJ(JX)=C.J(J) CO.NTINUE CONTINUE 60 TO 55 IF(KAME.E∈】.i). ). ..#....8. HARA7322 日ARA7323 日ARA7524 日ARA7325 HARA732占 HARA7327 日ARA7328 日ARA7329 日ARA753O HARA73う1 日ARA7332 日ARA73う3 日ARA7さ34 HARA7335 日ARA7き36 日ARA7337 HARA7338 日ARA7339 HARA7う40 日ARA7341 HARA7342 日ARA734き HARA734Lf HARA7345 HARA7348.

(19) 60. 原. ヰ3. 45. 輔. MXXsO CONTtNUE. DL芋FLDAT(NXX)祥ALOG(FLOAT(NXX). 44. 純. XMl7FLOAT(MXT)/FLOAT(NXX) MXTsO DO 45 J=1IJヱ D0 年4 T王1IIZ MXX=MXX◆N(IIJ) lF(MXX.EQ.0) GO TO 45 ⊂L=亡L+FLOAT(MXX)#ALOG(FLOAT(MXX) MXTsMxT◆MxX栄J MXX芦0 CONTINUE XMJ三FLOAT(MXT)/FLOAT(NXX) XXX言(AL...BL.･⊂L+DL)♯2.0 IDF買(IZ.{1)舛(JZ-1) lF(IDF.EQ.0) lSG=2H ZF.(lDF.6T.ヨ0) GO TO 46 Xl王XXX-S6(1●lDF). ). ). x5=XXX-SG(2?ZPF) GD TO 47 Xl言ABS(SQRT(XXX祥2.0)-SORT(FLOAT(IDF推2r,1) ))-2.576 ×5=xl◆0.616 47 lSG=2H++ IF(xl.GT.XO) tF(Xl.LE.XO.AND.X5.GT.×0) ● ･ISG=2日 lF(X5.LE.XO) lSG=2日NS C=.♯‥..1....井‥‥2‥..与i.‥.3....戟.‥.4.".-)t‥‥5.‥.推....6‥..井‥..7 DO 299 7=1IIZ LZJ=N(lIJY) IF(NSZJ.EQ.0) LZJ=O ZZJ=FLOAT(N(llJX)-LZJ) DO 298 J=1IJZ 舶. LZI毒NりY一J) IF(NSヱ!.EG).0) LZl王O ZヱⅠ=FLOAT(NくⅠ×I｣)-LZI) 298 29?. CR=〔Rヰ(FLOAT(NりIJ))/ZZl)#(FLOAT(N(I,J))/ZZJ) CONTINUE. ⊂R=〔R-1.0 CR:SQRT(CR/FLOAT( IST) ) ⊂=.井‥.=1....鞍‥‥2‥‥#....3‥‥萌..‥4‥‥う:I.‥.S‥‥)t..‥6....溌‥,.7 VI｣=0.8 VII三0.0 VJJ三0.0 DO 301 Ⅰ=1IⅠZ DO 300 JFlナJZ. HARA7347 HARA7!48 HARA7349 日ARA7!50 日ARA7351 HARA･7352 日ARA7355 HARA7354 日ARA7355 HARA7356 日ARA7さ57 HARA7558 HARA7!59 日ARA73bO HARA7361 日ARA73b2 日ARA73る3 日ARA736年 HARA73b5 日ARA7368 日APA7367 HARA7368 HARA7369 HARA7370 HARA7371 HARA7!72 ..#....8 日ARA7373 日ARA7374 HARA7375 日ARA737b HARA7ヨ77 日APA7さ78 HARA7379 HARA7380 日ARA7581 日ARA7382 日ARA7383 日ARA7584. XJ目貫FLOAT(J卜XMJ VIJsVIJ+FLOAT(NくIIJ))#XZH3!XJH VII=VIZ+FLOAT(N(王IJ))戟(XIH手ほ2) 3OO. ..1i..,.8 日ARA7585 日ARA7388 日ARA7587 HARA738B HARA7∋89 日ARA7390 HARA7591 HARA7592 日ARA739さ. VJJ三VJJ+FLOAT(∼(王IJ)川く央J日子捕2). HARA7394. 3Ol. CONTINUE. XIH=FLOAT(J卜XMI. HARA7395 HARA7う96. RsyIJ/SORT(yII戟VJJ) C=.准‥'.1..‥K‥...2 8P=￠ 8R=O 54 r=1IⅠZ-1 DO DO1 50 J=1IJ之-1 BG=O DO 49 L妄!+1IIZ DO 48 M=J◆1IJヱ. r<....手‥‥寺号‥‥4..‥与を.‥.5‥‥韓..‥6.‥.洋‥‥7....耕....8 HARA73g7 HARA7598 HARA7399 HARA7400 日ARA7401 HARA7402 HARA7403.

(20) 61. 社会調査データ分析のための汎用プ.pグラム. HARA7年04 日ARA7405 HARA7408 日ARA7407 HARA7408 HARA7409. ヰ8. BG-_BG+N(LIM) 49 CONTINUE 8P相P+B馴㈹(l'J) 50 53 J芦2IJヱ DO 8G=O DO 52 L=Ⅰ+1II之 DO 51 M=1IJ･,1 51も6三B6+N(LIM〉 52 CONTINUE 53 BR=aR+B6与さN(lIJ) 54 ⊂ONT!NUE. HARA74l_0 HARA7ヰ11 日ARAア412. 日ARA7叫3 日ARA7414. 日ARA7415. GK=FLロAT(BP-BR)/FLOAT(BP+BR) c‥.♯‥‥1..‥推‥‥2=..うモ..=3=‥栄･‥･4‥‥手書‥‥5‥‥iを････6･ 55 N6K$7◆(!X与的) IF(KAHE.EQ.5.OR.KPER.EQ.1) N6K王7◆(JX納) 王F(JX.LT.!X) IF(KAME.EQ.5.AND.JX.LT.ⅠX･ロR･KPER･EQ･1･ANO･JX･LT･IX) N6K=9 IF(KAME.LE.2) NGK;5 IF(KTBL.EQ.0) NGK=O tF(MSS.EQ.0) NG=NG+NGK GO TO 56 IF(NG.LE.も3) NG-_NGK NP三NP十1 GD TO b9 tF(NP.EQ.2O1) NAMEI(PRロJ(M)IM:1'4)?NSAM,NP uRITE(も,200) 5る NPX(K)=NP 60 TO 68 IF(MSS.EO.O) T〔】b4 GO IF(KTaL.EQ.0) c‥.井‥‥1..‥与を..1.2‥..栄...13==･;妄=･･4--Yc-･･5--韓-･･6･･･.黄-●･7･ DO 58 l=1▼!X 57 DO ｣=1IJX. p(I,J)=FLOAT(N(tIJ))･:壬100.O/FLOAT(∼(Ⅰ,JX)) 57 Q(I,J):FLOAT(N(l*J))*100.0/FLOAT(N(lX,J)) 58 CONTINUE c‥.井‥..1..‥韓‥‥2‥‥凱..'3==米==4･-･･::･--5--手書･.･･6-･.描-.･7･ I芸1I7X DO bO 59 DO J三1IJX 59 XX(!IJ)=4日 60 CONTINUE 6年 60 TO IF(KAME.EQ.1) GO TO 的 lFくIST.LE.0) Ⅰ=1IIZ Dロ 63 しZ｣=N(王IJY) LZJ=O IF(NSZJ.【Q.0) J=1IJヱロO b2 LZl;N(lYIJ) LZI王O lF(NSZ王.EQ.0) AL=0.0 BL=0.0 ⊂し王0.0 NFりIl)=NくIIJ) NF(1I2)=N(王IJX)IMF(1Il)-LZJ NF(2Il)…N(IXIJ卜NF(1Il卜LZI NFく212)=NXX廿NF(1)1)-NF(1I2卜NF(2Il) NF(1I3)=NF(1Il)◆NF(1.†2) NF(2I･さ)=NXXサNF(1I3) NF(3Il)qNF(lll)+NF(211) NF(さI2)=NXX･･NFくさIl) DO さ04 L=1I2. NGK=NGK+IX. .#....7....#....8 HARA7年1与 HARA7417 HARA7418 日ARA7419 NGK=NGK+JX 日ARA7420 HARA7421 HARA7年22 日ARA7423 HARA742年 HARA7425 HARA7年26. 日ARA7年27 HARA7428 日ARA7429 日ARA7430 HARA7431 ..#....8 HARA7432 日ARA74ささ HARA743年 HARA7ち35 日ARA7年さb ..#....8 日ARA7437 HARA7438 日ARA7439 日ARA7与40 日ARA7441 日ARA7年42 HARA744さ HARA744ヰ HARA74年5 日ARA7446 日ARA74ヰ7 HARA7448 日ARA7449 日ARA7450 HARA7451 HARA7452 日ARA745さ HARA7454 日ARA7455 HARA7年56 日ARA7457 日ARA7458 HARA7459 日ARA7帖○.

(21) 62. 原. 302. 303. 如4. 純. 輔. 502 DO M=1†2 1F(NF(L'M〉.LE.0) GO TO 302 AL王AL◆FLOATくNF(LIM))#ALOG(FLOAT(NF(LIM))) CONTINUE IF(NF(LIヨ〉.LE.0) 60 TO 303 BL王BL●FLロAT(NF(Lナ3))舶LOG(FLOAT(NFくL,3))). HARA74占1 日ARA7462 日ARA7465 HARA?4b4 HARA7485 日ARA7年66 HARA74b7 HARA7468 日ARA7469 HARA7470 日ARA7年71 HARA7472 日ARA747さ. IF(NF,(3IL).LE.0) 60 TO 304 亡L…⊂L.･.FLOAT(NF(3Iし))舶LOG(FLOAT(NF(5,L) ) ) 亡QNTINUE X=(ALqBし-⊂L+DL川2.O Xl=×-6.635 X5=Xr3.841 Xt)き0,0 ANIE;FLOAT(NF(111))-(FLOAT(NF(1I3))/FLOAT(NXX)〉韓FLDAT(NF(∋,1)) IF(ANZE.LT.0.0) GO To 61 1F(Xl.GT.XO) XX(一IIJ)=4日 ++. HARA7年7年 tlARA7475 日ARA747占 HARA7477 HARA7478 日ARA7479 日ARA7489 日ARA7481 HARA7482. + lF(Xl.LE.XO.AND.X5.GT.XO) XX(I,J)=4H TO GO 62 61 1F(×1.GT.XOJ XX(r7J)=4H IF(xl.LE.XO.AND.X5.GT.XO) XX(I,J)=礼H 62 CONT川UE b3 ⊂ONTINUE ⊂-･林-･･1････与さ････2‥=祥1..･3==手書.=.4.‥.I:‡‥..5‥‥詫‥..6‥..壬さ....7. 64 MA=MQ(1IK) --. -. MB=MQ(27K) TF(KAME･EQ･1) IF,(KAME･EQ･3･OR･KAME･EQ15). HRITE(も'201). KIMJI(A(LIMJ)IL=1●2) uRITE(6I2O2). ‥#.‥.8 HARA7483 HARA7484 日ARA7485 KIMtI(A(LIq!)●HARA748b. l盲F;吉真≡孟皇?ぎ喜;…;L'MJ''L喜喜:‡皇.'圭 ,呂芸≡芸………… ,. 1L;1I2)lMJl(A(LIMJ)IL=1I2)?ISGIXXX,IDF,CR,6K,R 1F(JX.LT.1X) 6n TO 66 MRITE(6?204) (CJ(J)▼J=1IJX) 1F(KAME.EQ,3.OR.KAME.EQ.5) HRITE(6I205) qlIMJ IF(KAME.EQ.年) bJRITE(占I20e)) MAIMI,MJ DO 65 l=1IIX ｢Hリ IF(KAME.EQ.1.AND.I N--L]} Gロ TO 65 E.I X 1L一L‖H l =U J ) l N WRITE(6I207) ⊂l( XQV∧ くl JR =汀Jlヒ JX) rH一 lF(KAME.EQ.4.AND 21占) MIJ(MBIM4】ヽ/ ド ▼l l J WRITE(6I208 l J 〓 1 l J 日ロ lF(KAME,EQ.5.OR.KPER Q 7 1 I R T E ー′0 (KAK(J)I9(TIJ)IJ=1IJX) 65 CロNTINUE 68 Gロ TO 66 MRITE(6I204) (CI(∫))l王1,ⅠX) IF(KAME.EQ.3.8R.KAME･.EQ.5) 叫RITE(bI205) MJIMI lF(KAME.EQ.4) URITE(b'206) MA,MJ,MI DO 与7 J=1IJX tJRITE(6'207) (J(J)I(XX(Ⅰ†J)'N(Ⅰ,J),Ⅰ=1,王X) IF(KAME,EQ･4.AND.J.EQ.1) MRZTE(6†216) MIJ(叫BI糾) lF(KAME.NE･5･AND･KPER.EG),O) URITE(6I208) (P(王}J)}Ⅰ=1IIX) lF(KAME･EQ･5･OR･KPER･EQ･1) 的RITE(a,2Oら=3(IIJ),l=1▼Ⅰ×)HARA7510 67 1F(KAME･EQ･5･OR･KPER･EQ･1) HRITE(6,217) (KAK(I),P(-?J),Ⅰ=1,tX) 68 〔ONT7NUE ●. ●■. ●. .▼Hリ. >∧. (ノー. ●. ｢‖り商山. LL). ●. lノO. '‥H. _. ｢一.=】L. I-▲. ●■l一. ●. C-･井-∴1･･･.1米‥‥2･‥･鞍････ヨ･=.洋=‥4‥..ラ(-‥‥5‥..米‥..6..,.予告..‥7. 69 70. IF(NP.EQ,201) NTAB=K,1 Ll=1 IF(Ll.GT..NTAB) 60 TO 72 uRzTE(bI209) NAMEl(PRロJ(M)IM=1I4)INSAM L2=Ll◆59 DO 71 Kl=し1'し2. HARA7490 日APA749) HARA'/492 HARA7493 日ARA7494 IIARA7495 日ARA7496 日ARA7497 日ARA749さ HARA7499 HARA7500. HARA7581 HARA7502 日ARA75o∋ HARA7504 日ARA7505 日ARA75o6 HARA75o7 HARA7508 HARA7509 HARA7511 HARA7512 ..!!....8. HARA7513 HARA751/tt HARA7515 HARA7516 日ARA7517 HARA7518.

(22) 63. 社会調査デ-タ分析のための汎用プログラム. GO TO 71 lF(Kl.GT.NTAB) MA=MQ(1IKl〉 MB=MQ(21Kl)甘1 MI=NQ(1IKl) MJ=NQ(2IKl)･ IF(把AME.EQ.1 IF(KAME.EQ 1,M=1I2)IMJI(A(LIMJ)IL=1I2))NPX(Kl) IF(KAME.EQ.4) uRITE(る1212). HARA7519 日ARA7528 日ARA7521 HARA7522 日ARA752さ HARA752年 HARA7525 日ARA752d. 主:oR.KAMET…去TE',6'210'Kl'"J'古紙T8!',芸:;!'…主:試写さ KIIMAIMBl(A(MIMA),M=1,2). 書MII(A(LIHARA7527 HARA7528 日ARA7529 日ARA7530 日ARA7531 HARA7532 HARA7533 HARA7534 日ARA75き5. 1Ml),L$1,2)IMJ)(A(KIMJ)?K=1事2)*NPX(Kl) K2=Kl◆占8 IF(K2.￠T.NTAB) MA=MQ(1IK2).. GO. TO. 71. MB:MQ(2,官2)亡1 Ml=NQ(1IK2). T≠謂宝島喜:ES:1, MRITE(6,21引K2,MJ,(A(M叫',q=1,2',NP=2'. uRI'E`は￠川2'"I'`A`…""HAi芸AA;5,…亨 1l,FH皇;?"2号:EMQJ･,子Å?RL;KMAJ?号こ≡;･,…',,NP…2, WRITF(b'215) K2IMA'MB'(A(M川A･),M;l･2)･M!,(A(L･HARA7538. IF(KAME.Ea,4) 1MT),L壬1,2)IMJI(A(KIMJ)IK=1I2)INPX(K2) 71 CONTINUE Ll=し1◆120 6【】 TO 70 72 WRⅠTE(6●218). HARA7539 HARA754O HARA7541 ‖ARA7542 日ARA7543 日ARA7544. ST9P. c‥.端‥..1‥..壬!‥‥2‥..溌....3.=.辞==4=..幸手==5==#･.･.18･･･.韓‥･･7--斗･-･8 ‥輔‥..1‥..壬!‥‥2‥..溌....3.=.辞==4=..幸手==5==#･･･.18･･･ 100 F【伯MAT(1XIA4IIXIA4I2XIA与) FORMAT(2014) 101 FORMAT(4A5I7XI13事49×IAl) 104 1O5 FORMAT(l与,1X,2A5I5X'6(IIIAl'12'AIIA5)) FロRMATく8()Al) 106 107 FORMAT(20A4) OO7 (IA5'1H)IIOXI4A5'‡5I50XI4日PAGE'Ⅰ4) 200 FORMAT(1日l,10日HARA 201 202. 203. HARA75年5 HARA754b HARA7547 日ARA7548 日ARA75年9 日ARA755() HARA7551 HARA7552. FORMAT(//1日I5日TABLEI‡ヰIZ6ク2日(I2A511日)) FロRMAT(//1日I5日TA8LEIIヰl王b'2=. # (,2A5,1日)/′1日,20日ARA7553 ('2A5'4日) 'Ⅰち'2H 1XI5HXSB= IA2IF8●3I3=(DFIIヨIl=),5X,9日⊂RくSQRT)=,Fb･3,5×,3日GK三,F7･3,5日ARA755年 HARA7555 2XI2HR=IF8.4) i ･ほ,2日ARA755占 FORMAT(/′1H ,5日TABLE,7年,Ⅰも'2日(I2A5'4H)推'13'2日(･2A5,4日〉 1日(,2A5,1日)′′1日,20X,5日X5Q-'A2'F8.3'5日(DF'!う'1日)'5x'9日CR(SORT)王,FHARA7557. ………二…≡…三…三==≡…: 呂芸≡芸…≡ …!_;1FF喜喜:MIAAI冒…E:';三m;冒喜て言…AA≦' _≡……… :H芸≡芸……… 007. (IA5ナ1=)IIOXI4A5'T5//1日,5日TAaLE,sox,5日PAGEHARA75帥. 209. FORMAT(1日IIIO日日ARA. 212. FORMAT(1日I21年IIH/IIIIIH(,2A5,2H)溌,Ⅰ5,1日(,2A5,2日)*,Ⅰ5,1日(,2A5,1日),HARA75る8. ,,X㈱H, ,,5.誉紙去?冒主A,, 1H ,. ,. …王;1E喜喜喜芸……王H:;……芸;≡;;;‡:;;…AA…;圭 215. 216 217 218. FロRMAT(1日◆IbOXI2t4IIH/IIl?1H(I2A512日川'Z5IIH('2A5'2H柑'!うIIH(-2A5,HARA7572 11日)I15IIHP) FORMAT.(1日◆IA5) FORMAT(1日I17･XIll(A3'F5.1!1H))) JIKKOQ OWARI榊兼/I 007 FORMAT(1HII29日柵韓日ARA END. HARA7573 HARA757ヰ HH･AARRAA 7, 7,. ;:. HARA7577.

(23)