bul gov. for. exp. sta. no. 29:3 1977 171~ 21(;月,の〆“e 4e ム...

46
Bul l .Gov.For.Exp.Sta. No. 29:3 , 1977 2 1(; 8 ) Yasuaki KUROKAWA: Reporton Computer (Researchnote) �J 0 ', 1;1: Fo r. Manage.-9

Upload: others

Post on 05-Oct-2020

0 views

Category:

Documents


0 download

TRANSCRIPT

Page 1: Bul Gov. For. Exp. Sta. No. 29:3 1977 171~ 21(;月,の〆“e 4E ム 林業試験場研究報告第293号 2. 主成分分析の考え方主成分分析は,理論的にはβ個の観測変量に対してある線型変換を施して企く新しいm個(ρ>m)

;H;~~\ìiJf fU Bull. Gov. For. Exp. Sta. No. 29:3 , 1977 171~ 21(;

(研究\tJ〆

15やi

成小

ログラミング十(8 )

主成 相刑V

黒川 A、主

Yasuaki KUROKAWA: Report on Computer Programmin立 (8)

一一-Princip乱l じomponent Analysis~一一

(Research note)

要 旨:主成分分析の新来から多くの有用な悩却を引き l.1 jすためには,変数の主11み換え, 什iのgz

数を&:燃したり,ザンプノレの集トJIの入れ誇えなど,いろいろ試みなければならない。

本報告では p 研究者の利用の使を図るために次の事項について配慮しているの①データ 1;): ,カート

から入)Jする場合は Free format を採用してデータカードの作成を容易にしている。②変数の11"'1主

は段~Iî 40 まで, サンプノレの倒数は 32, 767 まで計算処政が可能である。(3)人)J した変数の刊iを任志、

にとり出せる変数選択の機能を備えている。④ 20 将来日の変数変換ノレ{チンを備えているのくのサンプ

ノレ j刊の翌日択が可能である。⑥i11関係数行タIJ,分散共分散1 J 列のいずれからも主成分分析できる行ヂIJj通

訳の機能を備えているつ⑦ _l:成分スコアのf定数分布ならびに 2 つのスコアの剤lみ台わせによる散命位!

のllri\':が可能である σ

1 はしがき

ある対象について,その特徴を把握するには多律絹のデータを求めるのが普通であるの ζれらを多変景

データとよぷが,とのデータの的J肝 lこは , ]H~治は,多くの科 fl が同時総{'ì((j (と取り扱われ,ごく [ltQ られた

場介にのみ 1 穏怒iの特れが同居とされる。多特付のf占 Y~i{を評価iftfTすること:むその対象に関する 'nrllj'哀

の判断iこまつべきものであり,統計的 )ii去のなしうる仕事ではないとずる考えブ了は原則[内 lこけしいが,対

象 iこ関する rH!' り ((j知識が卜分でfよい段階においては歩なんらかの判断を下そうとしてもよるべき基準がH

もないという場合が生ずる。かかる場合に多変量のデータをM{j[lí,給。もする手iJ訟としてな場すろのが変

量併併である。その中でも p ここでとりあげる k成分分析は, TJ いに相関りある多1\11:均 i の特刊 11むのも~),出

減を立し、に知体!日司 tJ_少数個の総台特竹{ぬ:こ袋約する子 itであり,多変盤解併におけるんLノド的なカリ、の 1 つ

である。

本城fbは, 農林省林業試験場に設-iþ~j~ されている電子計算険 OKITAC4500 システムを対象にして開発

した二t]JX分分析のプ口グラムとそれの使用法を不IJRJ者に便利な形lごとりまとめ広く不IJ 片 Hと供ナるために研

究資料として発表するものである。また,木報告ーは竿ti可ーの東京工業大戸への国内儲γ'における研究成果の

一郎として作成したものである。|寸大学J1Jl学部情 w川学科鈴木光司教授ならびに i~1 村健二郎幼二f- Iこ;討さ;を

表するとともに,多くのJlJJê:を fった終?コ~ i~i~測定研究室Jll~MすijlYHl( 1ごお初 :11 し上げる次第でð5 �.J 0

1976"FIO月臼口叉胤

(1) 経宮部

',1;1: ~--9 For. Manage. -9

Page 2: Bul Gov. For. Exp. Sta. No. 29:3 1977 171~ 21(;月,の〆“e 4E ム 林業試験場研究報告第293号 2. 主成分分析の考え方主成分分析は,理論的にはβ個の観測変量に対してある線型変換を施して企く新しいm個(ρ>m)

の〆“

月,e

4Eム 林業試験場研究報告第 293 号

2. 主成分分析の考え方

主成分分析は,理論的には β 個の観測変量に対してある線型変換を施して企く新しい m 個 (ρ>m) の

相関のない蒸i'\f\化した変量を得る方法である。具体的問題からいうと,多数の要問についての観測データ

を少数の共通因子から成る指標にまとめる手法であり,方法論的 lとは相関係数行手1]あるいは分散共分散行

列の悶有値とそれに伴う l古[有ベクト Jレを求める問題に帰活される。との場合E とくに数少ない同有値で会

分散の大部分を抽出するととができ含しかも,その意味付けが的確にできるならば,もとのデータのもっ

ていた情報を要約するという了工場からきわめて有効である。主成分分析は,もともと心tm学において知能

テストおよび性特テストなどの結果を解釈するために考案され売援してきた統計技術であり,多数の観測

変量の中に合まれる共通院j内部因子の数を推定し,同時にその窓味付けを考えるものとして使用されてき

た。

ある対象に関する卸価の椋木が独立にとられ,その各々について戸伺の特性値が測られている場合多こ

の戸イ図の特怯 fnむの|習にはなんらかの相関があるのが通常である。相関が完全になければ各特性僚について

まったく独立に nj閣の標本を評価すればよい。しかし,相関がある場合は,その相関そもたらした要素を

考慮しながら解釈しなければならない。ことから,次の 2 つの婆:求が生ずる。 1 つは,戸{固の特性値が n

個の襟木に関して無杷関となる新たな統計量を作成するとと,他の l つは,乙うして作成した統計阜のう

ち,もとの情報をなるべく数/身ない統計量に要約したいというととであるのこうして作成されたのが主成

分である。主成分は数学的には h伺の特性債の玄み付き平均と考えられp その重み(係数)が以れと述べ

る条例ニを満足するように定められる。係数が定まれば (1) 式の左辺の主成分の値が決定できるの

X

も4i

J

PE付p

x

p

l

J

十十円“Z

2

1

J

十I

X --

J

一一

i z

... (1 ) drpy'B'R'

白目白日'SEa--ls-EE--111111''F

ノz

x

z

z

m

f

J

r-bPEュ

-

t

pZMe

一「

一一

.J

P

巴り

ー一ん

L叫・十

十一

十一向

2

E

X

白地

2

f

ず九日一十

+ぃ

'1

1

J

z

d

-

m

z

l

l

一一

一一

L

ゐ乙

ここで

p

Ik12 十九22十・ υ 十 Ikp2= Ehi2=1 (k=l , 2 ・ .., m)..................(2)

とする。

条件は,第 l 主成分 Zl の係数 {llú (i = 1, 2,…・, p) は (2) の条件の下でみの分散が長大にな

るように定める。第 2 主成分 Z2 の係数 {lzi} は (2) 式与を満足しかっ Zl と無相闘になる (j古交する)と

いう条件の下でゐの分散が最大になるように定める c 以下同様にして Z3, Z4, • …・.~ Zp;:;:で順次係数

を定める。換百すれば, 変数 Xh X2 , •.•••• xp に商交変般を施して得られる変数を Z1> Z2, …,

とする。この Zb Z2, ....・・ , Zp が互いに熊本目関となるように定め, このうちその分散が大きいものから

JIUむと Zh Z2 , ……, Zp とする。

Page 3: Bul Gov. For. Exp. Sta. No. 29:3 1977 171~ 21(;月,の〆“e 4E ム 林業試験場研究報告第293号 2. 主成分分析の考え方主成分分析は,理論的にはβ個の観測変量に対してある線型変換を施して企く新しいm個(ρ>m)

林業試験J若tE算機フ内ログラミング線作 (8) (研究資料) (黒川〉 qd

hi

令ha

3. プログラムの内容

( 1 ) プログラムの識別

プログラムお : PCA/45

作成者名;黒)I[泰 ~t

1'1 見交イ [ê n II : 1976 1[� 6 Yl 15 日

J[手 式・独立したプログラム

i七汗l 言語: FOR'1、RAN (JIS 71<司王 70(0)

(2 )白的

このプログラムば, l~;l準的な主成分分肝をい従i l:t':の解釈に必要な J虫の統計量在求め忍。相関係数ずj

ylJから一仁成分分析ヂるか予分散共分f&Uタljから主成分分析すおかは手IJ片j 設が指定する。

(3 ) プログラムの機能

とのアログラムは次に述べる計算処問機能をもっ。主成分分析するべき ρ似の変数をもっ nf同の総本の

デーダを読み込み, これらのもつ γJ{~ヲ t.\'!r手偏12 など )íl;1~1 ドj統計量ならびに分限共分散1 f711,相関係数行

列を計算する。使JlJ 苫の指定によって目白、ずれかの行列のもつ同有値とこれに対応する [iJi!?向ベクトルを

求め, これをもとにして討勺本?累積1寄与ネ,肉子負荷量を計算する。洋一日1 された各~l 成分のもとの変数

iと対する言I 与率ならびに'!i'([t'i寄与之容を百 I j立するとともに,各主成分のスコアをl~i /')ためのもとの変数に対

する係数そ~tlfìl~ 寸 00 そして, J没後に名}一成分のスコアを言i 釘するのなお,主成分分析の結束から右瓜な

械を引き i討すためには,変数のMlみ換え,ー部のさ足数の変換司 サンプル袋聞の入れ替え,スコア散布図

の作成などが必要となる。その穴めに,ィ:ての機能を備えている。

① 変数の選択:入力した変数の企部について主成分分析の計算をーするだけでなく,その一部奇任意に

j廷V~できる c

③ 変数の変換:主成分分析をする場介,もとの測に11立そのままでなく対数変換をしたり,あるいは,

Xi十 Jピj , Xì~Xh など適当な変換を施してからの方が良い結果が何られる ζ とがある c このために

20 孫類の変数変換ノいチンが組み込まれている。

③ サンツJ ノレ針。コ選択:入力データそあらかじめいくつかの群に分けておいて,自群ごとに,いくつか

の砕をまとめて,あるいは,全体について主成分分析できるようサンプノレ])yが係立に選択できるの

④ スコア散布閃の作成:各):_成分のスコアをもE用してサンプ)"の分対あるいはサンプノレの順位づけを

する場介,主成分スコアの散布図の作成が LF要となる p このため,名主成分の組み合わせによるスコア散

布悶の作成が 1111旨である。

⑤ データ入力の多総fb: ヂ…タ入力はカー!ごからだけでなく,磁気テ伊ーブ(11該宗〔ディスク)匂からも

トiJ能である υ また,印刷されるもωのうち l必要のない場介はその一部が街路できる。この場合l主計算時間

が短縮できるハカードから入j]ずる場介[土日 Ltl 斉式 Free format を採閉しているので, データの記入な

らびにパン一千が官易となり,また,カードの牧:数がiお約できる。

Page 4: Bul Gov. For. Exp. Sta. No. 29:3 1977 171~ 21(;月,の〆“e 4E ム 林業試験場研究報告第293号 2. 主成分分析の考え方主成分分析は,理論的にはβ個の観測変量に対してある線型変換を施して企く新しいm個(ρ>m)

4

斗晶

林業試験場研究報告第 293 号

4. 印刷結果の見方

このプログラムを使用して主成分分析をした場介の I':IJ )Iþl]結果と,それのもつ;意味をl阪を追って解説する

(凶 3 参照)0

( 1) PROBLEM NAME

とれば,とれから主成分分析しようとする問題の内容を表示するため lζ利用廷が任意に与えた標題であ

る c とれに続く 2 行lと変数の数ならびに標本の数が印刷される。

(2) ORIGINAL INPUT DATA

ことに原データが印刷される。まず,利用者が各変数!こ与えた変数名を印刷し,続いてこれに対応する

データを 1 標本ずつ印刷する。

(3) TRANSFORMED ANDjOR SELECTED DATA

原データあるいは原データに変数変換を施して新しく立成したデータのうち SELECT カードによって

選沢されたものを変数名とともに 1 標本ずつ印刷する。

(4) BASIC INFORMATlON ON DATA

( 3 )で選択したテ。ータのもつ基礎的統計ゐ量のうち、ド均 MEAN, 分散 VARIANCE, 標準偏去 S. D,

(standard deviation) , 3 x 標準偏差 3 キ SIGMA, 変動係数 C, V. (coe伍cient of variation) , 最大値

MAX. (maximum),最小11ll MIN. (minimum) を各変数について印刷する。

(5) VARIANCE COVARIANCE MATRIX

分散共分散行子1.]を印刷する。

( 6) CORRELATION MATRIX

相関係数行列を印刷する。

(7) RESULT OF PRINCIPAL COMPONENT ANALYSIS

これ以降に主成分分析の結果がl市刷されるのまず,主成分分析が相関係数行列から山手をしたときは (P.

C. A. ST ARTS FROlVl CORRELA TION MA TRIX) と印刷し, 分散共分散行列から i-!-\ 発したときは

w C. A. ST ARTS FROM V ARIANCE COV ARIANCE MATRIX) と印刷する。つづいて,固有値

EIGEN VALUE,寄与豆半( CONTRIBUTION RA TIO,累積寄与率 CUMULATIVE CONTRIBUTION

RATIO,回宥ベクトル EIGEN VECTOR,因子負荷量 FACTOR LOADING,各主成分のもとの変数

への寄与率 CONTRIBUTION OF COMPONENT TO ORIGINAL V ARIABLE,各主成分のもとの変

数への累積A好与卒 CUlVlULATIVE CONTRIBUTION OF COMPONENT 'I、o ORIGINAL VARIュ

ABLE,主成分スコアを得るためのもとの変数に士、j ずる係数 COMPONENT SCORE COEFFICIENT

を順次印刷する。

( 8) COMPONENT SCORE

主成分による標本スコアのfP !ijl]。スコアは相関係数行手IJから出発するとさは基準化されているが,分散

共分散行列から出先するときは説準化されない。

(9) COMPONENT SCORE SCATTER DIAGRAM

J主成分による標本スコアの散布図の印刷。 同一回~f)'il点 lこ 10 í町立|ーがプロッ卜される場合は X と印刷す

る。また2 散布図よりはみ1] 11 たドット数は欄外に NUMBER OF DOT5 FLO羽TED OUT OF DIAGRAM

Page 5: Bul Gov. For. Exp. Sta. No. 29:3 1977 171~ 21(;月,の〆“e 4E ム 林業試験場研究報告第293号 2. 主成分分析の考え方主成分分析は,理論的にはβ個の観測変量に対してある線型変換を施して企く新しいm個(ρ>m)

休業試験阪電機プログラミング f~告 (8) (研究官半りじょ 1) 』hd

門,tl

し C ilJ líI 1j する

P/_I の九統計量のjgnj ならびに数学的意味はx.献 (1) の pp. 160へ 175 を参照:されたいィまた, 3L休fr~

な芯i味については後で述べ 23 使用例を Jみられたい ο

5. データシートの書き

( 1 ) カードの種類

このフリログラムを {')UTl して主成分分析をする似合, 次の 10 極):fiのカードそむセッ卜する必裂がある c な

お q予のも ωl土イ\設のi弘子7はずTil!,!\できる。/こだし③は丸テ…プ, 磁気ディスク-;,I}!こ書き込まれt: デ

ータを用いる jぷùri不要である(表 3 左手照 c )

① PROBLM ガード

③ NAME カート

③ 変数名カ_- f

④ データカー 1.'

(与D TRANS カ ド

⑨ :交数変i負カ ドネ

⑦ SELECT カード

⑨~数i1~jJ日J ート‘

⑨ スコア散布ば|出燃H1;j~ カード*

⑩ STOP カード

( 2) 個々のカードの議き方

(I)は INTEGEH. の|終で整 J1X J\'J …て小数点を{~まない紋字, (R) は REAL の1ìYi\で尖数ニ小数点を含

む政4子を;意 uヰ:し,持政司,!は必す J~ íLさ才 l~ tニカラムの右詰めにさく。 )f îZ,A'i は括定されたカラムの純[:rHI付に

あればどこに浮いてもよい。以ドこの記 iJ を使用する。なお( )はブ 4 ーマットの形式をノJ'\すっ

① PROBLM カード CA6, I4 , 15, 212, lX , 60AlJ

データセり}び)先頭にくるカードはつねに PROBLM カードである。このカード l土佐川呂のデータセ

ット名をお;とするカードで,これによって名データセットをほ)Jlhr る。そして変数のお(ならびに燃本ω数

の持ノJし!京デーータの印刷の有却し以アータの入力デバイス-'I'.t伝亦するとともに問題の111闘を液JlJ すゐ。

1 ~ 6 カラム PROBLM とぎく υ

7 、 10 カラム:交Ji'zの紋そだく(1)。

11- 15 カラム a 隊本の放を !?f く(1)。

16 カラム ワど!日。

17 カラム: J息子 タを印刷寸るとき o :主 7こは作~ (11 印刷しなし、とき 1 とおく υ

18 カラム刷 ~t.: (-1 心

19 カラム: .Ij~~デ タ~~:' )刀 i、で人Jh るとさ 0 または~I/l ,的対テープ、 cm対ディスグ)から人

ブJ する_J-き 1 と;当く匂

20 カラム ツ日 Ho

21 、 80 カラム・ tW:mを-K く hhi:題は実数字,カタカナ,特殊又:子て、精l戎する。官:恕の最初JiζPCA!45

Page 6: Bul Gov. For. Exp. Sta. No. 29:3 1977 171~ 21(;月,の〆“e 4E ム 林業試験場研究報告第293号 2. 主成分分析の考え方主成分分析は,理論的にはβ個の観測変量に対してある線型変換を施して企く新しいm個(ρ>m)

-176. 林業試験場研究報告 第 293 号

と書く。

② NAME カード [A4, 9X , 12)

以下に変数名を与えるカードがくることを意味し,これらのカードの枚数を指示する。

1 C~ 4 カラム :NAME と書く。

5 ~13 カラム:空白。

14c~.15 カラム:変数名カードの枚数を害く(1)。

16~c80 カラム:室内。

③変数名カード[10X , 7( 12, A8))

原データの各変数に対して名前を与える。

1 ~10 カラム:空白,適当にコメントなどを書いてもよい G

11~12 カラム:後で使用する SELECT カードのために変数茶号を l から 1I民々に記入する(I)。

13~20 カラム:変数名をきく。変数名は英数字,カタカナ,特殊文字を用い 1 .立ご子以上 8 火字以ドで

!持iJ見する。

21~c22 カラム:変数番目白を記入する。

23~40 カラム;変数名を苫く。

以下 80 カラムまでくり返す。 2 枚組からもこの形式は同じである。

④ データカード [FREE FORMAT)

データカードは 1 標本ずつ小数点付きで記入する。燃本が変わるごと lと改行し左端から書:き始める。こ

れを最後の根本までくり返す。 白白書式 Free format を採泊しているのでとくにカラムの指定はない。

1 つの変数をあらわす数字の並びと他の変数をあらわす数r子の並びとの|閣に 1 カラム以七の空白そ置けば

よい。 したがって FORTRAN の欄記述子 10 F 8. 0 あるいは 8 F 10. 0 それ以外のいかなる形式で記入し

でもよい。 1 カラム lこ*叩をさーさ 2.~80 カラム lζ任意の英数字p カタカナ,特殊文字を高いたコメントカ

ードあるいは 1~80 カラム ~~t~l のフランクカードをデータカードの i習に適立何枚挿入してもよい。

⑤ TRANS カード [A5, 8X 歩 12)

以下 lこ変数変換カードがくることを怠l床し3 これらのカート‘の枚数を指ノ」くする。

1~5 カラム: TRANS とさく。

6~13 カラム:空白。

14~15 カラム:変数変換カードの枚数をき?く(I)。なお, 変数変換をしない助合\'立 O または空白と

ーヰーブヨ也。

16~c80 カラム:空白。

⑨変数変換カード* [10X , A5, 315, 2F10.0, 12, A8)

変数変換を行うもとの変数訴弓,新しく生成する変数の変数若手サ,変換式のコ…ド,必要な定数を勺え

る。

l'~lO カラム:空白。 適当にコメント古いてもよ L 、。

11'~15 カラム:変換式のコードをZ五詰めで書く。変数変換指定記号の1主主員とそれのもつ機能はぷ 1 の

とおりである。

16~20 カラム;変換式の L 変数の番号を議く(打。番号は,原データで Fでに使用した変数答号の

Page 7: Bul Gov. For. Exp. Sta. No. 29:3 1977 171~ 21(;月,の〆“e 4E ム 林業試験場研究報告第293号 2. 主成分分析の考え方主成分分析は,理論的にはβ個の観測変量に対してある線型変換を施して企く新しいm個(ρ>m)

体業試験場電~tic機プログ、ラミンク判日出 (8) (l!it'Jt'鶇D C引 1 1) 円Jウdl

1<1.玄数変換指定記サの稜匁l と機能

。、しP4

Jf

唱同J

Oむ

/。ば

すれ

ルホけ

出い心な

をしと定

b日

るず

β

な必

立とは

数〕ろ

店内出こ

VAと

D

ドる

,上あ

cNが

広沢一印

洛H仁

M一ち

cccc

ーか一ー十ネホ

NGPTSU川一段沢一川

ly

ちハし十

*hoxyrBf

一pu

/同

*Cf

しネド

CLLRWA日一N山則

M叫行珂

lilili--!11111111|

jE

JL第必

12345676901234067390

一〕

11111111112

一品上

寄り|コード

寸ー

:完投式

。。。

ωoocoo-J000000

X(L) X(M) 十 X(N)

X(L) ニ X(M)-X(N)

X(L) ニ X(M) ホ X(N)

X(L) ぷ X(M) / X(N)

X(L)~X(M) 十 C

X (L) ~ X (M) * C

X (L) = X (M) / C

X(L) =X(M)材C

X(L) ニニ X(M)料X(N)

X(L) ~ C -X(M)

X(L) ニ C ! X(M)

X(L)=X(M) ネじ +D

X(L) コ X(M) * C + X (N) * D

X(L) ~ーじ林X(M)

X(L) ~ALOG(X(M)+ C)

X(L) ~ALOG 1 1l (X(M) 十じ〉

X(L) ~EXP(X(M) + C)

X(L) ェー SQRT(X(M) + C)

X(L) ニ ABS(X(M) 十じ)

X(L)~SIGN(X(M) , X(N))

ccfuo

n Lノ

00000C000

ハ〉ト〉

(-, L ノ

次の{也から J邑し番弓で 11民々に与える 0

n~25 カラム交換式の M 変数の来日nを書く(I) 0

26'~"30 カラム・交換式の N 変欽の落弓を 33 く (!)α

31へ 40 カラム. .交換式のじの仰をー古く (R) 。

41.~50 カラム e 変換式のひの偵を勢:く (R) 0

51-v 52 カラム:変数変換によって新しく生成した変数の変数許可を古く。 ζ の波数喜子 ~t は限データで

すでに伏月j じた番りの次の他から通し番号でjijQ 々に',).える(1)れよって, 16 、 .20 カ

ラムて、 l l',ノjえした L と|司じ{直になるの

53-'-60 カラム回交;l~交 I史によ〉て;析しく正i:成した交をlUと与える変紋名をさく。 ðt JZi,g tま英数字予カタ

カナ, fJJ旅立_ J で 1 文字以上 8 文字以下で構成ヂる。

61. 80 カフム・ 2主r'1c

(主) SELECT カード (/\6 , 4X , 15J

j斗卜 Jと交政選れカードがくることを志1よし, これらのカート、の枚数をおJノト J る。

6 刀ラム: SELECI、とゴ子く ω

7 ヘ.10 カラム: ::'~(;。

ll'~ 15 カラム:変紋i哀れカ F ードの枚数をさく (1)。

Page 8: Bul Gov. For. Exp. Sta. No. 29:3 1977 171~ 21(;月,の〆“e 4E ム 林業試験場研究報告第293号 2. 主成分分析の考え方主成分分析は,理論的にはβ個の観測変量に対してある線型変換を施して企く新しいm個(ρ>m)

--178 ー 林業試験j誌研究液告\:fî 293 弓

16へ~80 カラム告白。

⑧ 変数選択カード (215, 412, 2X , 301 2)

サンブ}ノレ併の選択,選択したデータの印刷の有無,相関係数行列から IH発するか,分散共分散行列から

iUブ色グるかの指示,第fサ主成分スコアまで [:p刷するかラ怖かせようとするスコア散布図の枚数ならびに選

閃する変数を指示する。

1 ~ 5 カラム.サンブツレJl4'の始めの傑本番りを書く( 1 )。

6 ~lO カラム:サンプノレnrの終りの標本番りーを Zi く (1)。以上 2 つは空白にしておくと!ボデータの

すべてを l つの砕として選択する。

11 カラム:空白。

12 カラム:選択したデータを[=[]則するとき O またはさト'白,印刷しないとき l と者く。

13 カラム告白。

14 カラム:相関係数行列から主成分分析するとき 0 または空白 s 分散共分散行列から主成分分析

するとき 1 と省く。

15 カラムノ戸白。

16 カラム:第{Iリ J主成分のスコアまで求めるかをきく。 o または空白のときはスコアは計算しな

い( 1)0

n~18 カラム:泊かせようとする主成分スコア散布凶の枚数を指示する(1)。ただし, 16 カラムで

スコアの計算をしないと指示したときは散布閃は作成できないので 0 または祭白とす

る。

19~20 カラム:空白。

21~22 カラム:選択する変数の番号をきく(I)。

23~~24 カラム.選択する変数の番けを書く(1)。以下 80 カラムまで順次 2 カラムずつ必要なだけ伎

用する。ここを全部空主!とすると原データをすべて選択する。変数変換によ「て生成

したデータを選択する場合は,変換後,新しく生成した変数の番号を書く。

⑨ スコア散布凶膝根指示カードキ (5X , 31 60X)

スコア散布[2:1の照根と縮尺を指示する。このカードは③の 17~18 カラムで指示した枚数だけ必要とす

る。なお,縦軸と横軸 lと|苛じ主成分を指示するとそのスコアの度数分布を印刷する。

1 ~5 カラム:空白。

6 ~~10 カラム :"f'; M主成分を縦軸にするかを書く(I)。

1l-~15 カラム:第何主成分を横軸にするかを書く(1)。

16-~20 カラム:縮尺を自動的に決定したいとき 0 または空 (J とする。 1 , 2, 3,・ と指ぷすると,

縦軸,横紬とも最大値,最小値が土 2_ 5,土 5.0,土 7.5, • ーと 2.5 の倍数の単位が

とら il る(1)。

21u 80 カラム:主í' L.

⑮ STOP カード (A4, 76X)

計算の実行の終了を指/J<する。

1~4 カラム :STOP と吉く υ

Page 9: Bul Gov. For. Exp. Sta. No. 29:3 1977 171~ 21(;月,の〆“e 4E ム 林業試験場研究報告第293号 2. 主成分分析の考え方主成分分析は,理論的にはβ個の観測変量に対してある線型変換を施して企く新しいm個(ρ>m)

川、栄試験助活~:,t妓/。ログラミング縦汁 (8) (lりは〈資本ト) C,\U 11) 179

5 へ 80 カラム:空白。

以上でカードの準備はすべて完 fする。日IJのデータで計算

を続行する場企は STOP カードを入れずに PROBLM カ

ートずから lìili の)lil~iJ子にしたがってカードをセッ卜すればよ

い。また, 新たに変数変J免を f トコて計算1 を続行するときは

TRANS カードから後の;部分をセソトすればよい υ 根深1'1な

なデータデックの構IVì: は悶 1 のごとくなる。

6. 使用例

農山本、Iω計 1 1I1 1 を考える治合,地域が必主主な{,'[ ii汽を 1'-1めてい

ることはいうまでもないじ 地域および地域問題は iitllHiのより

どころとして重要(Chl、本条件ω1 つであり,地域農林業の日 i

i~IU立当該地域の是主体主主ω地域性にもとづいて検討されるの地

i或の誌は, )也Jft!学の \7J易からは 1 つまたはそれ_l;) 1 の村武

について間賞ないしj千十'kiJ ;fr するようなi克界内の|文主主であ

ると j るのが定説的である c

|ヨトデ タデックの格1&

さて?地域計i同とか地域 IJfl充など,いわゆる地域を対与とにした問題を考える陽子V, しばしば地域区分と

いう ζ とが行われる。この地域区分の←一般的な目的は,地域をいくつかに限うょすることによって,それぞ

れの JiliJ或の個性を発見すると正 にある。つまり,地域の l'iJ泌を考える場己、 lどは,まず,個々の地域の i凶 tl:

を比較して地域廷を発見し,これを泌識することから山発する。この地域区分の結果が地域;;t l!lr i とか地域

開発など地域を対象とする問題を:考える邸内、!こヲその地岐に関すゐ rttlNJlを設定する 11おj についての>f.袋

かつまま礎的な情慨を与えること!となるの

;也域をいくつかにじ五分する場f?には J也 )1人 気候, 土Pi などの li 然条件の他に農林業における経岱)1三

r!-~,経常上J~絞上限様式などの経済(i'-J条件や方、洛形態E 人し治皮などの社会的条V!:を同時に考 lむしなけれ

ばならないc しかし, これら多数ω弘~iHlを多)d(JにみfI関させて地域医分するととは容易でなし、。 iJt米か

ら,かかる場合には,単一J伝法にもとづいて透明紙lζ着色しp これを i~ ね余わせる )iiL つまり?オーパ

レイの方法が採用されてきた。しかし,このブ'i11;; は物犯(1'-)(とみてもせいぜいJ& íl~1の知的が京ねられるに寸

ぎない。そこで,この方法 l三代わるべきものとして主成分分析/ξ!11 し、た地域区分の方訟をとり上げるとと

にする。

地域区分ば防ノた県を対象とした。 地域区分のもっとも liki践となる地城市伏として農林業センサスの.1k計

単位の l つである IH市町村(自丹羽J 25 i十の [1,1.1持制査時点での l[i HIJ 十、J) をね:加するつ とこで 111 市町村を地域

単位として採用する盟山 ωtT~ 1 は, 昭和 30 咋前後からオすめられていお広域市lilr付心、併で",1← A 巾|町村 ω

中に純 111街地から山中j までが泌れしていて?まったく異質の地域が FT ま~している場 fTが多く, したが〉

て,現行の Ifj~1jト!を地域単位とすると広域的すぎて地域特引を把握オるのにj車引でない。しかし,これを

|日市町村!ととれば大きさも適ヴとなり,地域特性が比絞的よく J}~il',ずるものと/芯えられるつ活 2lふ長休

業センサスをはじめ多くの統計がし 1までも!日市町村単位に集汁されており,データ収集が容易である。第

3 は,り政の単位としていまなお|日市町付ごとに役湯の 11\必所などがある場合もあり,依然とじて Illtb判

Page 10: Bul Gov. For. Exp. Sta. No. 29:3 1977 171~ 21(;月,の〆“e 4E ム 林業試験場研究報告第293号 2. 主成分分析の考え方主成分分析は,理論的にはβ個の観測変量に対してある線型変換を施して企く新しいm個(ρ>m)

ハUO

O

寸1

晶 材業試験場研究報告 策 29:1 ';j

と城似の'に成成索。一メ'ず一不の松

川山地評くしよ門…合主製いル査存いを業は

じ小多多か恒例のて成なグ調依は在林

2

で各合もしの数し構ら

2

勢の切布る

p

部の場て

o

定少とのなの国業ーのめ市

会らのつる一く法そば際年就

hn

上門官

はれこいあらご方

E

れ指

mv

3

条に部

村こ'とでかしるずけす引へ均的家字

川町'が的と中いすまな一不

1

J体林は

市にる苅このな成'しをも部均主農

1

作次あ多る機一作合択件十林?のはけ守

す指単を尉選手行農たん

F4

るるけし断計な様るら叫蝋』ま林

xt

、わか判統う指すか業切は

o

襲。傑

でてるしでのよ成成中林持ねホな冷

o

出れす。度くる合作の農叶'一万う一一

47

理合ケい尺多すのをタと

ph

をよをな

な外凶な価す表

ζ

数一標川胎九肌る性く

主除をら-時一不代

o

変デ指

z;

す概う

dy

がお付工のをくる或刊すわれの不のり生ご

日域一町山一性よす州統一広明ふ?清合同た一のの

川崎刊に叫州出配付制悦樹市駅以蹴馴司齢れ

こ部のなの峡つパを既条問、ナ旬以の土担はは

切でい地も似性を境問お者日

A1績はの日お

るち

1

のし歩を耳特標環

f

とた印刷耐集ぬ業

X

いうてもまで椋のの旨く刊のし

1

〉口

λ

休仙のる

でのしる望

ζ

持活村削ま印欠事数都人ぬ農

x

夕、

き県成うもその共刊統りは比率率

4日本の

oは

-c

一て

生木作し長。らで市るとタ模様ギ系家家

J

山目白川叫す

KM

デ心

帆叩院畑山山叩い州州知予州出悦同時州川一川新日以制

L引い此相川

系でに述ノ\'を

GA

本のうり式〉〕一上九パ法

f'stしメト

x-

のと

出杓に江し様日ーと段られを

JUじ-ul

割農の七物腕ば古川ら名

ルャムの詳角断き択指の

wt'

れこ件以川市引の家者ば産

2

ぬのす叫れ村

の集る多判で川護見或週土ア

)O

F一P

出い丘

a

万引家自た叫バハ実符'会不況こ日

ι

収あをは断を仇出

ru引'る主等専旨

hh

〕農ちし出山ιけれぺ機を

U

山イいい小川市淵

ω成AUMU

岬間服州制引開制札川四川川峨刊開制蜘日川

アこ〕二日

fJ

人十件ま到来弄

r一面直完売完昔に\一一人へ}つ可ヲ

?l

Jい山跡什削Mいふ

Mが'hv

り官践産臥地林販販版事業的業入院川主持駅ぶ

域おら地あお

J一つ山いのの行かとは符次次率津山州湘山物物物従長

l

就収桁郎一山計指

VY

冶な明小がはくり則と択て計を人口

23

E

AE

有産院産業京地主卦の業家

UHるは

o

を各度合い作をも選け統業従人第第耕林業経保農林良農白耕民林ら応体此み

3

いたけは尺渇てを祈るの分ユ林::::::料:::ドド.トトじじじれ的良町を'

とつ特度価のつ桜分な数にシ良

hhhhv

川uh

良町内内川町川町川町向

ζ

市も

o

亙村

村なの凡評こ立指分と変プツ部れす比川

Page 11: Bul Gov. For. Exp. Sta. No. 29:3 1977 171~ 21(;月,の〆“e 4E ム 林業試験場研究報告第293号 2. 主成分分析の考え方主成分分析は,理論的にはβ個の観測変量に対してある線型変換を施して企く新しいm個(ρ>m)

川、、検伐ツログヲ ζ ンク局、Yti{' , (8) ((I)[ j'己資本 f) Uよ)11) 一 181

jて 2. ),~よテータ Jミ

7 I 7 日 4.

8 �J.971 32.

9 I l.6ól 乙¥0 I 1.

ハUGO

「/勺J

.• つりハ

υ

9

0

つ乙勺

d

lil �lJ 単位 I Jxl : 1 , 000 人) %2 尽人, χ 14 , X15

1:成分分れのS-'l.~~;t: (iJ 表 21:こ Jよされているデータの単位が li1] でないため村!関係数行列より L1J 党る。

表 2 にもとづき,これをデータシートピ記入ヂれば表 3 のようになる。計算対;長は凶 3 のとおりである。

計算時 1 :',11主力目→ドの沈み込み,ならびに l汁訴ふ'i -'1":ω r:jJ ir:1jをすべて合めて 13 分 40 秒であ〉たり

表 4 は百 131*山県のー却をまとめたものである。ここで凶行1な l土手千三Ir&;j}のもつ分散の六きさを/1三言。

1j ヰくは 15 主成分のもつ分以:ζ 対ずる会う}!訟の官IJ {j-を三し, これによって会打点分のもつ前I~Jノ:J-をあらわ

7。去 4 からも切らかなとおり, お1jニl点分がが~I 429$} :Gf~ 2 主成分が約 23~シゼある。よ寸て:2{~2 主成分

までで令体の約 6596 の,,~;lIIJj)Jをもっ。 1 <1 1 {Jベクトノレはも Jω交政の活みベクトルとしてj'l&分を形成す

るものであるいなお, こ乙では 130 の全サンブノレについての結::/~のみしかlt\げてないが, サンソ Jレをいく

っかに区分して主成分分肝した結J~~ も大休との去と近([J.した↑lうがj'} られp いわゆるキ}f!1~の蒋白川二ついて

も良好な減給与をぶしている ι 同子同守U',された主成分ともとのさ主主ý.との相関係以である。表 5 i i}.ミ

4 の凶 Fi'j荷量を分幻ノたものであるが, これか

らも I I) J らかなとおり,活 11 成分との制|苅を示

すもののうちその !n;( が1. 0 .1).8 である変数は

Xij , X7 , Xg , X15 , O. 8'~ ノ0.6 である交政は X12,O.S

XlSl Xl1 といずれも民主的特徴をがずふとの桁

O. 6~~ -0.8 0) ['J.の相聞をyjt LP こ l らはいずれ -I~O '0.5

日畑。吉

."

-!.O

関が出J (j -¥ u 反;Z;] iこ} Xsァ X6 , XI0 , Xlh X16 こはー

も休業 (1江村徴をぷ t交攻守あるにしたがって,

1 主成分はの占 IÎ'JIζ 農業度, Ho) 方 lí'j:こ主!文

~:. /]Z 9 ものと J::ll J1!f- ずることがでさる υ ノ弘前 2

足分につしp 、てみると Xj , X2、 Xs , Xl ど占i;

的特徴をン'j"ず急記数と強いJlの伺闘があ, 逆!こ :弐12 沼 1 分 2 ノ;ょ

!こ対する j;if 訟為、 Xs , Xlt1t XJ. 6 など JI 占:;11] (山特徴や示す&~お:と

Page 12: Bul Gov. For. Exp. Sta. No. 29:3 1977 171~ 21(;月,の〆“e 4E ム 林業試験場研究報告第293号 2. 主成分分析の考え方主成分分析は,理論的にはβ個の観測変量に対してある線型変換を施して企く新しいm個(ρ>m)

182

出コ司。炉

1合I u d o E e

百孟 b

出Q< 4

.L

!-< ふ与1

~ 対、

国 iE

〈ザ3 トト霊2凶Z ーー1ト

千吋

<Q

Cち

ヰ同

休業日点験協IJI究 i!i (1;- 3�1 293 リ

芝凶」国口広且

:11Jiiijjiliij111111jf!jJ1JJ11!::!:J1

jill111iJ1111jj1111:liJ11i:[11111i 11111111i1i111111j111:lull 1i11111;jillijilil:1111111111j Jillijillijjjjill;:111J111ijjjill 出iiiijill-ili:;自jilliijiiji111311iii1111112111-

Page 13: Bul Gov. For. Exp. Sta. No. 29:3 1977 171~ 21(;月,の〆“e 4E ム 林業試験場研究報告第293号 2. 主成分分析の考え方主成分分析は,理論的にはβ個の観測変量に対してある線型変換を施して企く新しいm個(ρ>m)

183 •

手:j 出

七il 之二足分 l 万う 3 :tf1文分

0.623:5

o. 5:591'

0.0841

CJ.4617

0.2166

- 0.1865

1 己0β

0.36/0

- 0 , 0139

cl 42'j5

0.1437

- 0. ~1961

�.3216

�.1838

0.CJ182

0.3969

竹、業叫ん主lß)J 1'IT主i 険 γc ログラミング組出ー (8) UíJト先手t1 り CiJjl!)

ワ. 637CJ

0.14白 1

0 吻 7084

0.7018

0 , 2594

一丸山/~

0.2100

- 0.5270

o 1799 - ワ 4404

0 , ,,1, 790

- o. ヌ917

0.3649

0.1 S94

0.10] ,)

o 4878

0.0394

0町。己 Oß

2701

l:)76

0.8193

0.74 て 7

()_ 89b3

o 4tl98 ヲ::::1 9

- 0 6.i3'!

()包 69;,:)

O. 6580

0./509

701:2

o 8806 -- 0 , óO~í 1

)成分分析の村長

γ; 31.1&分

1.8196

0.113γ

γ ó .33

CJ.I,621

0.4149

O.Oó之3

り 3419

0.1606

0.1383

0. 1117

_. O.

0.05 ぺ 8

0.31 ,;4

O. 106~

-‘ 29 ,16

- 0.23己4

。 1363

o 01 o. 2942

(). (_) 1 ~)ヌ

O.Ol96

O. 1040

() 060?

31 己b

0. :1 8己4

:5447

().1886

03Sb7

o 24/9 0 , 26 1' 8

(ん 20)33

0 , 2891

7700

O 司 3390

o 2330

同行悩

守 1] It:, 十〈

G

1

2

8

4

5

6

123

687891111111

Ix--1xxxxxxxxxxxz

討5 J オー• I!X; ヅト_iH :1 ~L)点介

Xth X7, X9~ X15

X12, X 13, X14.

X3, Xll , %15

X6, X7, X9~ X14

X8~ X12, X13~ X16

X2 , X4~ X10

Xl

Xh X2, X3~ X'l

.Xl1

λ\.';)

X9, X.14, X15

X7~ XUl, X18

X6 , X10 , X]曻

X1, X4

X2

X~

Xs, X10, X111 Xl�

[利子負副長

]旬。 O. 8

o 8 、日 6

O. モ _-,-, 0、 4

0 ペー りえ

0.2 0 0

O. () r" -0 ,

--O. 2 ゾ← O. 4

o 1, _v -0.6

0. 6 _'r~ -0.8

O. 8 _r' ~ 1.口

は(えの/卜IJ関が政い c したが-,て, ~~~ 2 � 成分は,1ーの方向に計]Çj!-j' ,1丸 山の )JiClj にゴi叩点目 dJìえを方、すものと解釈

ここでは私版的!とする ζ とができる。 3 ì 成分はすべての変数に対して相対的に悶了負荷量t ノう'/J 、さく

は避けた。図 2 (土川 l トー成分, ?.{~ 2 ])成分 iこ立J'J 変放のもつ D~ 了役符早ーを悶!としたもの解釈する

ちと lごして特引も明らかなとおりヲ判十111出[立大きく日つω グ _!l--~ --_/にけか;~ L 1 である〉このはカ、

!o'( fl 身をう〉引いj、ことも ilJí1伝である。

EÍl九分分析を i i うのはし 1;瓦分スコアを捌いてサンプノレのうあペパ 1 は吋ン/コノレ(ひげl{jJ:づ l' なするのがJþ:

オ乙の 111下;であろ o 130 0) サンプjレについて求めた L~わ 2 成分のスコアの -n'iz 市凶は電算機のアウト γ。

Page 14: Bul Gov. For. Exp. Sta. No. 29:3 1977 171~ 21(;月,の〆“e 4E ム 林業試験場研究報告第293号 2. 主成分分析の考え方主成分分析は,理論的にはβ個の観測変量に対してある線型変換を施して企く新しいm個(ρ>m)

184- 林業試験場研究報告俗 293 号

ットとして作成されている。との悶は 130 のサンプjレの 16 変数をJ~椋紬とする 16 次元空間における散布

の状態を第 1,策 2 主成分を酬とするτ]Z I出 tに射影したものである。 ζ の際,約 35%の情報のロスはある

ものの,大局的には 16 次元空間でもほぼこのように散布しているであろうと考えられる。 16 次空間での

散布図を可視的応分割することは不可能であるから,その代丹j としてとの 2 次元の悶を瓜いて会サンブ。ル

司会分割することになる。さて3 この散布図を眺めれば 130 のサンプルが大体 4 つのグループビ分かれて

いることがわかる。つまり,第 1 夕、、ノレーブ。は農業地域,第 2 グノレープは祁 rþ'地:lliふ箔 3 グループは都IÌÎ近

郊休業地域,第 4 グjレ,--プは 111村林業地域としてそれぞれのグループを特徴づけすることができる。もと

のサンプルがどのグループに属するか与を散布図のとで定め?これをもとにして同ーの属性をもつものをま

とめて地図川乙表示することによって地域医分が完了するのなお,紙幅の関係で地域区分をした地区iは苦'ì

略ずる。.Y:Jトーが主成分分所?とを用いた地域区分の考え方のi既日再であるが, j也l'Xl I だけでほじ属性をもっ地

域単位をグルーピングするだけでは各地域のもつ性格を明確にすることはできなし凡そこで全サンプルの

平均を 1 とする各該当地域の各変数の水準を多角形閃として示すと存地域のもつ性格をより鮮明に把揮す

るととができる。

かかる試みの主要な 11 的は,ある地域において農業なり林業に関する開発や~IT苅をすすめようとする場

公に,開発や計阿の対象地域のもつ性1与を計悶主体に判断の材料として提供しようとすることにあるが,

これには 2 , 3 の問題点を残しているのまず第 1 は, 統計データが 1ーとして農林業センサスから採られて

いるが,センサスは属人統計データであって属地的に収集されたものでないc したがって主成分として導

['1I した 2 つの総合特性{誌には二!匂地に関する物的情報が不足している。第 2 は, ここで使われているデータ

は現状のレベル与を示すものであって地域の将来を予測するものではない。市 3 は?ここセ;S1~づけした都

市皮,農業皮,林業皮などは対象地域内の相対 i'J!~t;強弱の関係を示すものであって,絶対的なものではな

い。以卜ーの点を卜分考廃すれば,かかる方法は地域区分の 1 つの)/i1 として有効に利用できるものと辺、わ

れる *10

7. プログラムの構造

このプログラムは 1 倒の主プログラムと 10個のサフツレーチン両Ijプログラムから構成されている。多くの

百Ijプログラムを使用する最大の理由はヰ記憶装聞の記憶容量;が小さく,ー→度にすべてのプログラムを主記

憶渓[~l の「十Itこ栴約できないためである。そとで,プログラムを多くの品目プ口グラムに分割し,オーパレイ

(over-Iay) 構造の下で実i jーする。 なお, オーバーレイの弱含は補助記憶装青からのアクセスが必要とな

るので計算速度が遅くなる。必要総メモリー認数はオーバレイ構造をとるとき 19 , 171 ~ïi (1語=16 bit) ,

オーパレイ構造をとらないとき 30, 115 訴となる。主プログラムはサブールーチン I;(IJ プログラムの呼び出し

のみを実行し,実際の計算はすべてサブールーチン副プログラムで行う。

各サブjレーチン副プログラムのオーパーレイ構遣と主主きな機能を解説すれば次のようになる。なお( )

内は必要メモリ一語数を示す。

*1 詳しくは,黒)11 泰亨 :j 成分分析による地域医分, 8711 1 1 円木林学会発表論文集,日木林学会, 1976 年 10JJ , pp.35-36 を参照されたい。

Page 15: Bul Gov. For. Exp. Sta. No. 29:3 1977 171~ 21(;月,の〆“e 4E ム 林業試験場研究報告第293号 2. 主成分分析の考え方主成分分析は,理論的にはβ個の観測変量に対してある線型変換を施して企く新しいm個(ρ>m)

林業l試験反;屯J兵q機フ。ログラミング減/; (8) (引先Yミキり(玉、 !II) A百

MAIN (10,

01 (1 , 240) SUBS 10 (1ヲ (62)

02 (1, 384) ~-SLJBS 20 (1, 211)

0日(1, 569)

SUBS 05 (1, 161) - SUBS 30 (1, 982)

SUBS 06 (1, 956) --SUおS 40 (6 , :,87)

( 1) SUBS 01 : PROBLM NAME,変数千;カドなど;;~.;山みi主ト須のトろ心

(2) SUおS 10 : n;i 1、 F タ涜み辺、み十iIl JlfJ ~i己似災市ヘ転j主雪る c

(3) 3UBS 02 :之絞変 Jí~カード ~t 読みみ変数変院の,γ(ii1J令

(4) 3UBS 20 :夜放変怜ぺヲ乙。

( 5) SlJ BS 03 :変 ')l:i選択カードをみp 選択した変数の hH1 ノ(/':]ti: どぞ計訴すわ

(6) 3UBS 04 :分以共分散行列p 相関係数行ダIJを弥j ろ。

(7) SUBS 05 分目立~~分散ヂIJ,相関係数行列,間半H仏国手J ベクトノレなどを印刷すろの

(8) SUBS30: ヤコビ-tJ;で四千JÚι 目白パクトノレを計算する判。

(9) SUBS 06 ‘凶了負荷 EI,許可ヰ成分スコアなどを口|灼し印刷する。

(10) SじBS 40: ):Ij立分スコアの散布灰|そイ午)!x; L, IW,:~!Jす乙 p

これらのソースゾログラムは凶 4 のとおりである。

8. J寄辺機器構成

とのプログラムの実行には次の}河川機慌を必要とする。

( 1 ) カード 1) ー夕、

(2 ) ライン.1 1) ンタ( 1 i� 1:i2 'ドHL)

( 3) tr!lDJJ三己可決(DlSじ 1 , DJSC 2, DISC 3 ())日つの論JiH. ~ν,)う ;í1: !fJ できるものれな;JccS,データ

補:日12「:ね法院から入力する場合のtG]腿オゴは必ず DISCO を{却ちするの原データ 7] ード,気子ープ (i涼

全i ディスク) ,Jこも 1 レコード 80.f{îのキャラグターモードごなけiLばならない c)

9. 注意事

(l) iこのプログラムのづ~1f rjJ必要に応じて次のメソ位ージを nJ!'I:Uして;, I'!"ì を 11'lL γ ろ。

① NAME CARD 13 INCORRECT 変数名カ… Fの枚数ヘ勺記ノ入が注ぺている。

。o DATA CARD IS INCORRECT : Jボデータが誤って疋入さがたりパンチさ t1 ている。

③ TRANS CARD IS INCORRECT 変数変検カードが~n コて記入されたりパンチされている。

④ SELECど CARD IS INCORRECT 変数選択カード、が誤って記入されパン子されている

( 2) 変数交換ノレーチンに希望の君主検式が見当らない均台は, 5j!j(こ変換プロタラムをfr:1点して /}~'iデ F 今夕

を変換してい〕たん俄気テーソって5へ書き iii し,政くに ζ1 しを読み込んで iil 京了すればよい c

(3) 1.0iデータを FORTlミAN の S ド 10.0 う:の tl版記述了を 1,、て作成寸ろ場介J :/~L~ [;"l をもって O とすース

的 1\1イT{I[[および[81Hベクトノレ奇求めるサブノレ…チンプ口グラム lZ,然J汁滞;11りと好「ノfL寝 11j]:!FIM \;,哉氏竹 fJ立のものを !'lJ}J

している。

Page 16: Bul Gov. For. Exp. Sta. No. 29:3 1977 171~ 21(;月,の〆“e 4E ム 林業試験場研究報告第293号 2. 主成分分析の考え方主成分分析は,理論的にはβ個の観測変量に対してある線型変換を施して企く新しいm個(ρ>m)

-186-"

ととは許されない。

乙のプログラムで処1河型

原テデ、一タの変数の数

J休業試験場研究桜子「 第 293 J%

10昭制限事項

主二 40

J)j(データの主主木の数 三三 32 , 767

変数変換によって新しく生成する変数の数十原データの変数の数三三 40

変数変換カードの枚数 三三 20

一度に選択できる変数の数

計算できる主成分スコアω数

文 高丈

1) 良野忠一他:多変最解析1];;, [j科技連日i版社, (1971)

2) 一一一一一一;続多変量解析法, 日科技連出版社, (1976)

三三 30

豆 7

3) COOLEY, W. W. and LOHNS , P匂R.:乱1:ul ti varia t巴 Data Analysis, John 羽Tiley & Sons, Inc.(1971)

4) 一一一一一一- and 一一一一一一一一: Multivariate Procedures for the Behavioral Scienc邑s , John

Wiley & Sons, Inc. (1962)

5) コンピコータ マネジメントサイエンス ハンドブック,オーム村ョ(1971)

Report on Computer Programming (8)

Principal Component Analysisュ

(Research note)

Yasuaki KUROKA W A (1)

Summary

Principal compon日nt analysis is a generally useful procedure whenever the task is to

determine the minimum number of independent dimensions needed to account for most of the

variance in the original set of variables. This analysis in practice is often conc邑rned with

the problem of reducing the original test space to a minimum number of dim印sions n巴巴ded

to describe the relevant information contained in the original observ旦tions.

The FORTRAN comput巴r program presented in this p昌per is proposed for solving th巴

principal component analysis. This program is divided into 10 sub-programs and reduces the

storage demand,乱nd permits use of the much small巴r computer (core cap乱city is 20 k wards ,

1 ward is 16 bits). This computer program can solv日 the problem which has fewer than 40

variables and 32,767 samples. In this paper a rural zoning problem which has 16 variables

and 130 samples is solved for example, and it takes about 13 minut日S.

Received October 9, 1976 (1) Forest Management Division

Page 17: Bul Gov. For. Exp. Sta. No. 29:3 1977 171~ 21(;月,の〆“e 4E ム 林業試験場研究報告第293号 2. 主成分分析の考え方主成分分析は,理論的にはβ個の観測変量に対してある線型変換を施して企く新しいm個(ρ>m)

林業試験場電算機プログラミング報告 (8) (研究資料) (黒jll) --187 -

図 3. 算: 結 果

-・"晴 DRINCI 戸 AL COMP口 NEN r A内 AlY515 ・噌・.,

~ð pRo日 LE 同 NAME •• PCA/-45 Z 白川 ING RY P~c~A* 事事 PC "145 寧.

MUMBER QF VAR! 向 BLES ヨ"NUN臼ER 口F SAMPLES

u ORIGJNAL IMPUT DATA 胡駒

1 ,-品(" 2 x-_( 2) ~ X--( 3) , x-唖(" 5 ,-由{旦】 , x・-(晶 3 7 X岨恥("

自 X~- ( 8) 唱,-恒(" 1O X申告(1 0) 11 x骨骨n lJ 12 X刷局(1 2) 13 x-・(1 3) 14 X曲-( 14' 15 x-ー(l 51 16 x--c 16 >

"噂'5'唱9 257.50口00 28.2899'ol 67.呂 9000 17464000 '.唖9000 ".h32999

19.34000 』自画67000 8 .42000 27.3499唖 42.5000口 53.00000 4'h57唖唖9

19.3299 号 3ω60000

2 哲也守699<) 25.2B999 36~26999 45.. 守 6 '1 99 3.99000 .3.吉申 9唖唖

26.71999 0.00000 52.98000 70$42明白骨 45ミ12999 25.吉Eヨ甲帝

22 ~Q999哩 31.42000 36.�l999 28.18999 39..28999 3. 司 9000 事~.8999唖

32.12999 39.4早口00 1.58000 lld2000 "‘039雪量 � .40999 54"065 唖申9

31φ10000 3.16000

3母也 2000 5.38000 19.78999 22.9唱 000 44.8299号 lO~99000 島,.吉0000

33 .3 1 唖" 3B.:nOQO 0.0000口 3.'16000 51.:31000 b2.48000 ".唖500025.84000 司 .81000

1 '..62弔99 37 ~20999 2 目 .1 白雪甲9 30 骨 789'<9 55. 守 2000 <.唖 9000 58~ ?3 000

11 ~06000 y,.自甲 000 O~OOOOO 10$0600口 3凸~73000 70 司 59000 46.15申唖申

30~5唖 .99 1Q.03999

4.94000 20.39000 22~ 76999 42.96骨9' 10ω99000 57 ., 45000

"排21999 o GOOOOO 13.68000 21~2ß999 65.569申昔 ヰ .h9Z(')QO

o. 口0000

1.<13000 <4.95000 i 自 .3官 000 21.6899号 2a.17000 47.宮."事 宮 7 ., 3400055.359甲 q 3s幽 68999 8&18000 lO~80000 48.1 自9申唖 6宮 .14000 ,書 .56事唱智

29.329守幸 5.1 品 000

33合%,ヨ甲 32000000 27.489Q9 33.59000 30~81000 36.98唱" "噂6099950.1499号 32 ~09000 1.09000 14 ,, 44000 43.21'<99 65 咽 189甲骨 <3.自1.唖甲

25.12骨'9 5 &80000

1. 品 600 口 2~46000 12.69000 1 \1 .18守" 26.39"99 50~OOOOO 畠1&92000

6!h569q9 43.17000 5骨 15000 a~23000 51.90唖99 '0.官"弔9 吉0.56喧 '9

28.939"9 5.60000

" 1.68000 L 色 46000 11.1 白 999 10.4 日守 99 24~61000 品 8.98999 E晶..1.4000

71..54999 Z自 .03999 10.25000 9.91000 40.23000 59.259甲唱 52~51 唖9唱

35 幽 20999 5.12自O口

II 32~939骨 9 36~459唖甲 31.09000 3 唖 .669争9 41.11000 31. 型自唖唖唖 62.23000

"皐06申9' 44.32999 2.16000 守.吉 0000 46.000QO 古l 曙 09000 50..2草00031.00000 ! .44000

12 ..~19000 3~87000 16.48999 2凸 .48999 42~43999 H..98999 6冒,哩9000

41 ~340QO 31a95999 0..47000 4G360QO 52 ,, 31000 70~0199事 41d8唖"

29.42999 0.9古 000

13 19 ~!)9甲99 120.21000 雪 4 惜 48999 411.18999 1O..HOOO 19~9699<) 唱."岳00

"‘ 45000 2.13000 8.23000 "噂7雪 000 25.45000 2fi.14唖99 31.唖串 000

II ~8品目。。 5.870口白

" 旦 .23000 且 .25000 32 ~189 'l 9 24.789甲唖 31.2唖9唖甲 32~989骨唖 11~ 75000

32.8499唖 9~04000 2.16000 13. 自甲000 31hlO弓" 51.1(>9唖事 43~l8唖"15. 自 7000 0.00000

15 5.01000 l'h46000 42.09000 ".晶8999 43.65999 26.98守唖甲 24.31000

2B φ00999 8.83000 1.95000 'iI ~OOOOO 33.唖 5唖甲9 38.00000 51.15号"

23 駒 599 '1'1 1.'l5000

16 9~82000 9.10000 :n.6B9'19 ! 9. 89000 1 吉咽 53000 6晶 .969唖9 昌 .19000

52.509弔q 3.38000 3.5坦000 27. 色 79 <;1 9 27 幽告 7000 28..81司甲事 35&14000

13.21000 2..98000

11 1..58000 1.800口。 "白 39000 15.39口 00 <;l.93000 80.98999 事.ー7000

44 ~ 19999 4.Z8000 8.69000 23.0499'1 32.20999 41085999 .H~34999

.10.88000 15~93000

" 0..86000 3.55000 .,.ち9000 11~4骨 000 11.21000 55. 司自甲唱弔 2.90000

26.12守" 1.93000 0.00000 3'.且4999 '6回 65号 99 1 唖 .209甲唖 32~12999

10~33000 0.00000

'9 5.00000 27.59000 75.09000 12亀且 9000 ".自 2999 28.989唖甲 5e91000

38.73999 3.93000 3.22000 35~20000 17. 自7唖" 33宥喧 0 '1 99 33~56999

15.940口0 0.00000

2' 12.00000 66.89000 30.98 司 99 ...吉 6000 4.99 口 00 23~549 <;1 9

"ω9自'" 17 唱 25000 0.0口 000 57.4699司 45.89000 47�92000 55.98000

18.31000 口 .00000

21 2~25000 i ヨ ~30000 2.2.39000 209 4 日甲" 20 暗 03000 43.9899<) 1.29000

ち 0.00000 口.自品 000 6.6òO口。 27.81 守" 13.2守口00 25.00000 :33.82 甲弔9

2ld;0骨甲事 0.00000

" 2.77000 2.1200口 .\3 .469骨q 11.69000 7~64000 自2..92骨 99 Z.llOOO

5l.0000Q 0.48000 呂.守口 000 51.31"99 13.18000 22 噛 20官 99 16.03999

早合 15000 6. ヲ20口。

H O~ 75000 1.04000 45.4 司 999 7..6900口 9479000 73 e 9B9唖申 5.04000

45.78000 3.020口。 <1.07000 3凸。 3499事 l 守 .32999 30 ,, 32唖" 31.16甲99

14.0130口口 0.0白 000

" ι.670白白 3"~87000 54.3')000 29. 司甲口 00 4 司 .39999 0.0000口 27.40999

66.65999 31 ~619<l 9 0.0口 00口 17.28999 45~81ヲ" 57 噌品499甲 --, 1. 'l 3 <J 99 ;;'7.12999 0.000口口

25 1.58000 13 ~25000 41.品目 999 21.18999 52.62"'99 O�OOOO 34 励 1 29'l<'i

"圃5也 O lJ O ".主99' 骨 0.0口 00司 ,.弓 8000 4 ヨ .019唖9 52 ,. 18999 10.65'199 27.54999 0.00000

1.48000 7.64000 23.7 且 999 65.20000 0.00000 3品 ~459司q

12.500司。 .H .53000 0.0口O 口口 7. 凸 00つo 41.509唖9 53.23000 自 O~200口o

0.00 口口 lコ

z7 ] .53000 17.7口官骨 9 52.18999 16.28司『号 59.l8唱 99 0.0口 00口 '1 ゅ 95999

12.50000 !1l .89ヲヲ喧 0.0(10')0 10.\主 100口 .39.51)000 4 唖 .56000 54 司 11000

23.4 <l守官官 0.00000

Page 18: Bul Gov. For. Exp. Sta. No. 29:3 1977 171~ 21(;月,の〆“e 4E ム 林業試験場研究報告第293号 2. 主成分分析の考え方主成分分析は,理論的にはβ個の観測変量に対してある線型変換を施して企く新しいm個(ρ>m)

。。 林業試験場研究報告第 293 号

図 3. (つづき〉

。。0

0

0

0

9

9

04

3

3

3

09

09

叩的

0

7

62

5

004

凹判

02

5

3

1

3

S

4

99

90

90

87

26

29

23

0

0

0

0

00

目印

2

0

6

。。。

090

090

i

Q

810

730

0

2

099

叩明判

385

699

g91

i

z

g

2

29 34.5<)000 95.8299<) ~1.28999 49.289<j9 3.3.0<)999 14.99000 3So?OO?? 37.5199') 2 1..う 999 'l 0 ・ 00000 15.94000 40.60999 46.18999 円 1. ち900022.31999 1.'.l9000

30 7.6雪 000 11.15000 31.5"000 20.39000 40.70999 36.98999 4 6-・ 20999,ち .98999 29.17999 0.00000 6. マ弓000 41.04999 60. 的 8999 53.29999 24 .42080 2~87000

1.E1400Q 5.30000 21 ・ 89000 13.49000 6O ~28999 0.99000 68 ・ 039995.25000 59 ・ 0699 '1 0.00000 5 ,51000 41)3-46999 611. 57999 ち 3.54999

34.20999 O~OOOOO

3~08000 6.!3000 22.39000 目 .789"9 60.26"99 6.99000 65.07999 14.55000 54.28000 0.00000 3~24000 S~~37999 70.<'3999 97~O!)OOO 33Q4BOOO 0.00000

9田明

74

13

3

8

23

q9

99

9

9

8s

q7

?

S

75

09

09

09

67

2

1

0

6

13

0

0

0

0

0

0

9

2

9

3

4

5

90

90

90

84

5

2

3

1

000

000

町印曲

624

285

000

000

000

147

210

495

R

A

4

3

9

0

0

9

9

0

0

9

9

0

0

9

0

ち?ち

2

4

9

4

9

q

5

0

3

5

1

4

0

9

9

9

0

9

9

9

099

9

5

8

1

9

3

9

7

1

6

1

7

1

5

5

3

9

9

9

9

9

9

9

9

9

9

9

9

41

2

5

3

3

4

5

?

5

6

4

3

2

2

0

0

9

0

0

0

9

0

0

0

9

0

9

6

8

7

0

5

4

8

5

9

0

7

2

2

A

0

0

9

0

0

0

9

0

0

0

9

0

9

0

8

4

8

4

7

4

bち

25

3

5

9

9

0

0

0

0

判円凹

mmm

5

9

9

7

7

3

2

0

6

0

8

6

6

6

2

3

6

3

3

2

l

9

9

0

0

9

0

明科的同判的

5

1

8

4

7

4

6

8

7

3

B

I

--7

8

0

2

0

0

5

E

5

1

4

2

5

6

3

3

37 1.40000 5.83000 26.89000 23.18999 65.56000 0..00000 47.51999 0.00000 39 , 90"99 0.00000 5.34000 40.07999 58.15000 8l.28000

42.45000 0.00000

38 l~ S6 000 13.:310::10 47d9000 2~.28Q9Ç 27_659Q9 27~ 'l 8Q99 43.289':;'9 45.89999 1'1.92000 fl ・ 77000 19.579'1'1 37~q2000 56.73000 44~40999

21.019'19 11.01000

9

9

9

9

2

7

22

7

5

9

9

9

2

0

4

0

0

四日

q2

3

0

90

0

0

0

0

ooo

m加叩

192

099

099

944

6l8

208

6

2

9

9

9

9

9

0

9

9

9

9

9

0

q

9

9

9

9

0

9

5

5

9

6

l

2

6

2

5

5

3

7

8

5

5

$

6

6

5

6

b

5

6

9

0

6

9

0

9

9

0

0

9

0

9

9

0

0

9

0

9

8

6

9

5

0

0

9

0

9

l

o

o

-

0

2

8

5

8

2

7

1

6

7

6

0

0

9

9

9

9

0

0

9

9

9

9

0

0

9

9

9

9

4

q

8

6

7

2

6

3

6

7

1

4

b

3

8

6

6

3

4

5

4

5

1

5

9

0

9

0

0

0

9

0

9

0

0

0

9

0

9

0

0

0

8

2

8

8

9

7

I

B

-

-

8

5

.

.

.

-..

-*切

1

6

5

6

3

5

2

2

1

9

0

0

0

0

0

9

0

0

0

0

0

9

0

0

0

0

0

8

9

9

1

9

7

1

1

B

O

B

-

8

1

4

?

0

8

3

2

2

2

0

9

0

0

0

0

0

9

9

0

9

0

0

0

0

0

9

9

0

9

0

0

0

0

0

9

9

6

5

0

Z

5

5

8

2

7

0

8

4

6

4

2

1

8

8

7027

5288

3

5

5

3

1

9

0

9

0

0

9

C

9

9

9

0

9

0

0

9

0

9

9

9

0

9

0

0

9

0

9

9

0

6

?

?

g

B

1

5

8

7

0

8

2

7

7

6

7

7

..••• 0

8

0

4

8

3

1

6

8

2

3

3

2

3

5

2

1

2

3

4

4

4

凶叩叩

96

90

05

5

9

9

99

99

9e3

87

9

9

99

5

l

027

24

0

0

0

0

0

0

1

6

9

0

9

0

S6

8

2

090

叫90

284

009

加開明

342 1.53000 0.39000 13.29000 16.98"'19 1.63000 94.<;18999 15.10000

ßO~34000 4.88000 3ち・ 8 ら999 7.55000 32. 雪 0999 47.90'199 46・ 42 '1 99U.86000 40.1 ち000

47 1.75000 1030000 20.89000 15.49000 9.60000 B4.9B'l9Q 40::'.18999 72.4899<;1 19 ・ 59000 24.12999 10&62000 33.84000 58..56000 <!l3.89999 18.95000 30.32999

q

9

0

8

4

0

8

5

0

9

0

ち0

6

0

6

0

9

0

9

0

2

3

0c

凹叩

9

0

4

0

5

0

0

0

別加

。。0

00

000

目的ω

lao

5

000

回目的

500

7Os

loa <49 25.64000 94 ・ 12000 ・~9.09000 39~18qq9 39.67999 2 ・ 99000 42 ・ 45 '1 99

iゐ. 53<;l99 23. ら 3999 12 ・ 50000 11~17000 37.37999 S3.9299Q 44.06999 23.73000 4.16000

雪o 1l ~93000 11.97000 34 ・ 489 <;1 9 16.48999 49.64999 1409 吟000 61 Q 18000 24 ,70000 44.000� 0 ・ 00000 6.92000 4 <;l ~71 'l'1 9 68~54q99 48 ・ 759 0;1 927" 7ヨ 000 O~OOOOO

51 8~70000 1<1.91000 <11.89000 30.18999 48.5199'1 12..99000 52.46句。。27.6099<;1 32.62000 0.00000 B~18000 43.-'9999 fJ4.aS9Q9 49.64000

"・OH99 O~OOOOO

52 2.14000 6.98000 )9..48999 15.29000 49.75000 0.00000 ~?o~0999 0.00000 38.06999 0 降 00000 5~2000Q 38.1099雪 (;0 .3 7000 50~34q99

28d5q9q O~OOOOO

99

9

6

5

叩円

08

06

9

9

9

9

3

3

S5

6

4

0

0

叩凶

O

Q

0

0

9

0

0

0

4

0

000

000

刊叫叫

5

日0

009

山間四

303

10i 9

0

90

90

62

51

52

6

6

叩刊

朝日

93

0

1T

5

0

9

09

09

02

8

ae5

90

90

90

85

78

46

2

叩問

。。日

曲B

ロ2

090

凹明叩

420

580

7bo

5

009

009

009

i60

6iQd

g35

3

4

4

5

o

q

9

9

9

0

0

9

9

9

9

0

0

9

9

9

9

0

2

2

5

5

9

8

9

1

3

5

0

7

56565

0

0

0

0

0

0

印加刊叩叩印

9

0

0

9

9

0

9

7

0

0

9

2

0

0

曲目帥

0

0

0

0

9

Q

m叩間四円叩

8

2

4

2

6

4

3

6

2

5

0

7

6373

4

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

0

3

1

i

9

0

0

0

0

0

明目印叩凹凹

柏町的加肝肝

3

1

Z

0

9

0

0

Q

0

0

9

0

mm側同制刑問問削

8

0

0

8

7

0

4

3

0

5

8

0

8

9

9

3

5

0

1

5

4

1

4

0

9

9

0

9

9

0

9

9

0

9

9

0

9

9

0

9

9

8

7

7

2

g

O

1

0

4

6b696

Y41

6

0

2

2

6

6

2

7

0

42312Z

5

6

7

5

s

s

~a 3~06000 5.97000 11. B9000 11.79000 52 ・ 40999 5.99000 44~67000 15.59000 2 '1 ~950 ヲ o s・ $7000 み之 e2QOOQ 25.5099<;1- 54.3799守 54 ・ 3900032.0 '1仁川 o O~OQO(}O

5守 nyQ

刊叩

80

10

L

0

2

990

明明∞

995

-50

676

2

3

099

099

Q90

85L

793

071A

123

51.37999 12.CJ9000 '3 7.,699 <1 45.5399守マ2031000. 51.53999

。。 ll. ,. 90の o 58.1599';¥ 16.98999 1;.守 ~40999';;;.;, 2000 :;4.06000 72~õ2999 50.81999

5.46000 34.849<;9 32 ~3 7999

Page 19: Bul Gov. For. Exp. Sta. No. 29:3 1977 171~ 21(;月,の〆“e 4E ム 林業試験場研究報告第293号 2. 主成分分析の考え方主成分分析は,理論的にはβ個の観測変量に対してある線型変換を施して企く新しいm個(ρ>m)

林業試験場電算機プログラミング報合 (8) (研先資料) (忠)11)

61

62

6> 6."6的01 r,, 0 1<)9') 29 ,'1 29 ワ 9

6'

ハ可d

O<U

図 3. (つづき〉

4.2flOOO 33~3 (J OOO ZO.S90()O <18.96999 lO~ <J'l OOO 61.50999 35.() C, 999 3~7600 <J 4.62000 '1 g.~.3 00a 70.4 日 999 49 ・ 81000l ~88000

?O.890QO lO~QQOOO 57Q31<J99 23.98999 55~色 20001.58000 f>..;37000 4B~62')9 <.l 69.95999 55.00999

>';・ 3 'l (!G() ち・ h8'h 'l 51.25999 6.0:)9000 60 ・ 610000.0';000 B.560DO 51~53999 る6.62000 告 4.64000

14.79000 22.39000 54.28999 5.99000 61.:34000 {).OQDOO 6u190JO 44.96999 61.87999 51.03000

もら

刊M

03

4

i

75

0

9

0

9

90

0

9

4m-0

0

明日

4

4

1

。。剛

山0

6

2e

-090

叩明凹

494

099

q79

40

~. 6 <, 0η0

~5 .'12999 2.21 ()円。

32. 明 90(]() 63.12999 8.99000 64~lO <'J 99 ;S.f.l800Q 4'J.09DOO 71 圃 84000 60.42000

67 3.匂 60 ,::0

".ω7 '1 99 _, <;.fl19 'i 9

t9

70

14~3匂000 ム 8 固 48 '1 9Q ち..ら 99<)9 1;?99000 64.35999 3.22QOQ <1 .480 乃 o 'D.34()QD !?6.12999 64 , 6099匂

4. 9~Oハ o 23.2 自 99 '1 22.!190rJO 51.20999 1~99000 11.;>3000 9.99000 SQ39(1)O 42.65999 50.93999 59~ Tl 999

しむ,加。

leHOOO 17~'7 B9 99 1S. 吃'1000 2So209'l9 61.989令') <;2.4�9<�9 33.18999 'h4200Q 8. ,jJ OQO 32.89000 45.07999 49.46999 9.42000

4~6?OOO 20~59000 21.390η46~34')9" 4.9<l000 62.269<1') (j 7 ~Oヮ。'N ü.ωOQO{) U.~lono 42.6~9ヨ 56.93999 54~J4000

0 , 00000

9

9

ヲ9

9

6

4

6

5

99

90

09

09

44

河刈

93

Q

0

3

00

00

2£

2

000

Mm的

I88

000

000

659

448

292

54

7

23.09000 13.09000 47.6.2999 23.93999 58.2300<l 2~81000 7.19000 6 ち .MOOO 69~20999 73.34000

1:¥.1899'1 8.1'1000 12.56000 83098999 43. 雪 199913.弓 1000 a 唱 61000 4q~OOOOO 72 白 09999 55.37000

70 コ .3'ちひ00

"'圃 1 .2 00039.56判 9

1 ‘5 ち 000 0.64000 9.89000 11.fl900Q ':¥.97000 88.98999 Jち .6299 'i

66.09000 72.65999 35.? ち 000 10 ・ 89000 38~4099 <:J �5.9200Q 62.10999 24 ・ 289'タ'} 40..310つつ

20.89000 17~18999 ;?�.}4999 67.98999 52.71999 2.34000 8 圃5(, 000 ち 9~50000 74.93999 96~59000

90.98999 21 ~12999 56.15999 39.90999

。。0

0

00

00

16

so

000

∞ω山

505

000

∞∞凹

,6

77 2.77000 73.04999 18.<)2000

マ8 6.81000 19.67明9

"・ 89999

明的

90

09

64

mW

09

94

6

99

99

00

0

0

7

8

開印

刊凶

3

0

090

叩刊∞

524

3

009

009

590

201

385

智伺74

19..28999 16.39000 63.78999 5.99000 60~01999 0 ,, 00000 6~ ら 7000 42.98999 70~51999 56 ・ 95000

10.4(JQOO 11.(,9000 4.05000 9l..9899CJ 24.069<;19 8.84000 41.98000 ヨ 8.09999 57010999

80 2.36000 Z. 叉 5000 13.29000 l7.39000 弓。 .57999 21.98<:J99 67~70000 50 ・ [>4999 51.51999 2.59000 2.73000 f. o.ど 3000 H .23000 63.5000C 34~48000 5~18000

0

0

9

0

0

0

9

0

0

0

9

0

6

1

5

2

0

8

9

9

3

1

0

3

4

5

S

6

9

9

9

0

9

9

9

0

9

9

9

0

8

8

8

0

9

7

9

5

-a

7

7

7

T

5

6

4

6

0

9

9

9

0

9

9

9

Q

9

9

9

?

5

9

9

1

3

0

8

4

7

3

8

7

3

3

4

0

0

0

0

0

0

0

0

0

0

0

0

9

h

9

9

835

2

3

0

3

∞的刊凹

0

9

9

0

9

3

8

6

8

9

2

2

@a

5

i

0

2

2

2

0

9

0

0

9

0

加問加田明叫

1

3

F

7

4

4

4

9

5

2

8

5

3

7

2

4

0

2

2

4

0

9

9

0

9

0

即刊羽田明田

2

I

8

5

2

4

3

2

9

0

b

6

5

2

6

4

2

9

A

g

5

2

1

2

8

8

83 4.08000 3.7BOOO 18.18999 12.5'1000 33.06000 43. '1 899ヲ(,0.4899949 , 40999 <11.98000 4.24000 4.06000 46.56000 67~5ó999 5ヲ・ 5099931.349切 9 .3.29000

B5

.134 B.l {泊 00 4.59000 2. 1 ・ 6~999 14.0QQOO 12.40000 79.9899<;1 2 ,'h42999 S.. , 8..99 り HI.989 '1 9 1~).3000 12.92000 45~3 ?9 ..g 63.18000 '52.54999 19.31000 .3. NOflO

2044.13'199

7.52000

7.[19000 12.65000 780989()9 42~98g99 ''1.4'-000 38 砂 59000 64.469'19 ;:'7.41099'き

2.3BOOO 67.349<;9 19.31000

8$ 4.60000 3.05000 19.18999 12.990勺 o J 7.73000 5 '5司 9 同 999 35.79999 58.92999 .3 Ú~42 (J OO 1.00000 5.14000 q3.0!999 69~95999 65.17000 20.~3999 <'.25000

17.43999 11~t89 'l 9 1ge3l000 ゐ 1 鋼 939<:)9 43.68999 2.46000 4~43:]OO 52.51.9'19 6(>~Q7999 50.75000

BS 13.05000 10.41000 15~4')OOO 16.4il <jヲ<;1 40.9699'1 31.98999 55.9899'1

3.;ぞ 10':'-0 49~ .2 5QOO 67.J9999 50~92999

1.\.4C,OOO 2"} .9299 ヲ S9~9Rq99 54.75999 Z 祖 28QQO 4.10COO S<;.84Qg9 71. ,1jヲ 99 55.42999

9よ

2.10000 3.GSQGO 19~ <l 8999 .I t~e9Cハ(,4.64')99 15.99000 ~4.1200() 32.19999 39.39<;1<;19 l~laOoo 5 回 43000 41.48999 63.17oQD ち 5.2300037.41 噌 99 .2.37000

1 1' ~78999 15.6'10CO 46.75000 26.98999 62.73000 1 ~4 2つつ,也41000 44.<:)8000 62.5Q999 46 回 43999

。。9

0

84

73

16

9

9

9

9

6

4

9

0

0

0

9

2

0

0

問問

60

090

090

190

200

目的判

388

9

93 8.22000 15.10000 37 園 8 ;1 000 28 ・ 18999 50...39"9 6.99000 58.14000 22.8<)000 4<?21999 0.00000 5. ち0000 45.14999 6 ち ~0 'l 999 5"~25000 7.7.2(>999 0.00000

Page 20: Bul Gov. For. Exp. Sta. No. 29:3 1977 171~ 21(;月,の〆“e 4E ム 林業試験場研究報告第293号 2. 主成分分析の考え方主成分分析は,理論的にはβ個の観測変量に対してある線型変換を施して企く新しいm個(ρ>m)

-190- 林業試験場研究報告第 293 号

図 3. Cつづき)

4 , 14000 20.759Y9 32..079"9

16.09000 53. z:', <;99 9.9<:)000 &3.329匂93~ 870勺 o 52.28<)99 色 7.4 (1<1 99 6S~70000

14.99000 0.00000

0

9

0

9

曲H

al

--06

75

mm

94

91

41

17

0

9

09

0

9

83

39

54

9

0

9

Q

9

0

8

8

64

9

0

9

0

9

0

Mm

•• 7

0

2

090

090

090

840

6B0

760

5

990

090

0qO

103

4f0

632

5

3

96 6.99000 32.34999 40.09ヨOO 35.890ηo 50 ・ 93999 圃 99000 56.28000 1 1I .~39 'l 9 47.00<;199 0.00000 13.2<1000 49.50<}<;I9 67.64999 58.12000

28.29999 0.00000

97 0.86000 1~58000 14..69000 22.78999 62.079<)'1 8.9<;000 68.82999 29.7 ,5<.99 58.17<J99 0.00000 5~400()O 55~7800D 61.03000 58.::nOOr. ;H~2a999 O~OOOOO

5.05000 7 ,1 0000 21.1899<'1 24.89000 63.92000 4.99000 61 ,, 6799。19..18999 44.25000 2.07000 7 ,<'1000 49.14999 66.3S999 S2~50000 27.90999 0.00回。

9 唱 01000 '1.58000 37 圃 50000 22.48999 47.65999 13.9'�OOO 57..989"9 n.ち '<000 35.12000 0.81000 7. ち 60 '1 0 38.70000 67.9299守 4 q~ 71999 27.50000 O~81000

100 1..33000 2.35000 15.3<;1000 21'0.89000 37 , -'5-'99 38.98999 5~&03qq9 46.31000 )Se ,HQÇ9 0 圃 00000 司 .45000 44.609号 11 48031000 ち 7 .1 8999286<)8000 1~13000

0

9

09

0

9

00

40

09

09

05

9

9

?

1

9

0

90

9o

mM

9

0

叩判的

S

$

0

009

m叩伺

452

24 凹

明1

7

0

0

叩帥

94

93

6

0

9

09

0

9

94

74

0

0

。き

69

0

80

60

9

0

oq

∞明開

250

640

009

009

462

683

0

9

0

9

3

0

9

9

3

3

0

9

98

92

9

9

08

0

0

0

0

96

0

0

0

0

00

0

0

0

0

090

m明叩

230

009

009

009

533

030

3

9

9

99

99

0

9

09

98

97

S

90

。o

mw凹

34

96

90

81

96

0

0

叫叩

伺叩

70

000

凹印∞

650

000

000

000

527

791

206

4

105 2.66000 9.21000 23.69000 ち .79000 46.50000 4.99000 43.25999 9.51000 24.40999 0.00000 5 圃 37000 43.50000 63~71 'l 99 4a~62999

22s21999 O~OOOOO

106 .3..67000 8.10000 22.4B999 _12.28Q99 号 0~43999 11.9899.. 49~S2999 21.28000 24 ,, 10Q99 0.00000 5.46000 31.7699'< 63.2099'< 42~ð700。21.31999 0.00000

107 8..96000 21~95000 36 ・ 28999 18.68'199 48.70000 23.98999 .3マ .9299。3...06000 21~069q9 0.00000 10*45000 31~3999CJ 48�13000 <19.21999

25.nooo ・ 60000

108 2.65000 2.55000 16.98999 11.79000 23.62999 61.98999 36~4899号43.78999 1l.<1S000 .3 ~57000 2.70000 30~ <l 5000 40.65999 4. 匂 6700025~45000 7e07000

9

9

2

2

1

0

刊伺

97

16

9

0

90

9

0

93

?

3

90

84

5

S

0

0

0

0

00

9

0

90

3

。。。

加凶叩

440

030

009

009

009

481

3

9

0

110 3.61000 9'・ 17000 2.. ・ 89000 26.39000 63.14999 6.99000 53.46999 10.82000 4"~09000 0000000 8.35000 49 ・ 93999 65.79999 7'hOQ999

3:h87999 0.00000

0

9

0

0

0

9

0

0

0

9

0

0

7

3

5

8

6

4

9

2

2

7

9

0

5

5

9

0

0

9

9

0

0

9

9

0

0

9

9

9

9

2

0

9

0

q

0

9

0

9

0

9

0

9

845

9

8

4

0

0

0

0

0

0

0

0

0

9

7

9

4

0

0

0

2

4

6

4

E

0

0

0

0

捌帥剛山

8

9

8

0

2

2

2

0

3

2

090090

凹明凹叩明叩

8L7520

978840

2

4

3

4

6

0

0

0

9

0

0

9

凹凹明叩凹刊

85077

7

2

7

E

L

5

6

8

3

4

8

5

1

3

2

1

3

113 4.20000 4 ・ 65000 22.69000 24&09000 0.00000 0.00000 33.59999 3.99000 10.65000 0.00000 16.37999 30.7 ち999 32.42000 ヨ 3.03000

17.ゐ 1000 O~OOOOO

114 2 ・ 21000 2.92000 14~79000 13 ・ B9000 72.12 999 3~ 99000 70.62000 47.35999 62.64999 0 圃 00000 2.79000 4ら .2 ち 000 5 円 .96999 <;4~6400C 40.43999 0.00000

115 2・ 21000 2 ・ 92000 母 .79000 13.89000 0.00000 O~OOOOO 78 ・ 1200019.98999 56.25000 0.00000 0.00000 45.25< 00 69 ・ 81999 49.14000 :53.81000

116 9.37000 6.02000 17.59000 17.39000 1>0.37q99 10.99000 ?3~~9?~~ 32~37000 50020000 ・'1 5000 4.01000 47.03000 52~10999 56~S9999 ;n~42000 l~95000

0

0

0

9

9

1

9

9

9

0

90

例目

。。凹

叫93

5

1

2

090

090

663

86a

4

099

叩判明

061

25

7

118 号 50000 8..66000 25.39000 28.18999 58.8ち 99 '1 20.<)I;199Q 58.25999 '1 ・ 87999 49.599"19 0.00000 lO~S4000 34.01'199 58.43999 ち 0.. 73000 31:..93999 -3e220DO

5.43000 3.68000 17.39000 13.89000 3C.S6999 48 ・ 9899 <;1 58. ち000063.06000 42~82999 3.01000 12 , 72000 41.23000 66.3明00 4<) ・白 00011.67000 2.51000

120 3.48000 2 , 19000 12.3'K)OO 13.190う o 35.&499'1 48.98999 77.620110 n.01999 50.1?0 つ o 3 , 01000 1.140 ,, 0 51.95999 73.18'199 48 ・ 100の032.81000 1~91000

0

0

0

0

9

9

9

9

8

6

9

0

q

O

D

O

Rマ5

0

9

0

0

凹判的

456

4

009

別的明

200 122: 6e5酬 o 3.64叩o 19.1fl999 lQ~5 t; 0~'0 :1 1. 84999 ち.叩刊 60.51>9<;9

ゐ 5 e 28999 39.59999 2.21000 5.:noηo "・ 18999 5ó.~8 '1 9今ちo. ヲ700026.S09'N 1 ~OOOOO

123 12 ・ 690つ o 16.59匂 99 41.39000 :W.390 つ o 32.75999 32 ・ 98999 15.?0000 .%~ü7999 9.05000 1.33000 弓 ~5BOOO 3:.:nOOo 33.9も匂 ~919&03削 o 6.410fJO

12<:> ;2・ 19000 3.56000 35.98999 i4 園 09000 19.50999 64."8'1'1'1 15. 】 '00の40.34999 ll.04000 l3 .~2000 9. ち 7000 40.10000 'j().11999 1Q.6l999 U~93000

12雪 l ・ 66000 O~Ó5000 ;>6 ・ 39000 9. 'j QO 勺 o 3.68000 91.'fB9<J9 10ol300C 104.50999 2 ・ 8..000 14.9り000 14~4 8 じり H ・ 87 輌<; 3~.η ", 000 38.67000

J 2. ,&000 23.98999

凹凹

9

0

9

0

2

1

84

0

0

2

0

0

0

0

0

9

6

0

0

0

0

Z

7

009

009

0nvo

田明凹

499

504

200

6

Page 21: Bul Gov. For. Exp. Sta. No. 29:3 1977 171~ 21(;月,の〆“e 4E ム 林業試験場研究報告第293号 2. 主成分分析の考え方主成分分析は,理論的にはβ個の観測変量に対してある線型変換を施して企く新しいm個(ρ>m)

林業試験lお電算機プログラミング組合 (8) (研究資料) (黒)11)

図 3. (つづき)

1.1 暗~'l(l{l n て o

70.'JQqサ., 14. 74 00 コ 1 ち。 I <., OO C)

10.49000

1.27 1.6 C, 口口円 2S 唱。 C,ûQü I ぅ.\')口門町

4.2., 0口口 1 りつつ ~ J .ちち 0つo

12S 自.,弓ウ00 40.28<)<) ワ Zl.fl9I)(I(l

'..2 .819守ヲ 口.'i 7Wl 口 29.799<; 官

7.13000

129 3.230ウ口 36.789守q /1.0<)0(10

15.1 口 000 15.720口0

¥}Q O.λ 14 , 4<)000 1/.390円。

4 ーも <<010 lO.円 t 口O口 』早岨 Ol 'ì 99

34.0h 口 ()O

'ψTR ANSFO円I1 l0 A/WICfl "F.'~lCTl:) llAiA .・ I ら'lE 仁Tl nl~ιARD N口. ~ 1 )

,-制( 8) x-鴨 (15)

..忌 AS] C ST AT 1ST !CS ••

26.20467

20.14756 宅呂 .11 fl38 31.1987 ち

4'1' .b792.3 \9. ち 0549

.32.i13517 4. ~6~2 6

ヨ 5044

42.18050 5 自 .1404 ¥ " l~ 78430 26.94411 ヨ薗1.11 , 2

•• vARJMiCE CO I' ÄH! A!'l仁ト門 r, iR � X ••

;0:--( 11 トω( 1) Ci l.78曲3

;0:--( 2) 之 ~5. 3Z 110

;0:--( :n 31.59653 )(--( 4) 59.05501 X 司白【 5) 8.78098

x司句( 6) -(, 1 圃 83923

X一司【 7 > l~26005

;0:--( B) -26 聞 14149

x_ー{号 -:'.98138

:':--(10) -6.0961:>0 x_ω(tl1 13.0ち '1'X 合会<1, > 2.7390" x--<i3 J -0.99E115

X ー・'"】 -35'H62& X匂ー 1).5) 令 2. ~290ß

x-ー (16)干11.2.HHI6

¥--( ¥)

x--( 2) ~--( 3)

" 51

" " Bl

" 10 ) 11)

I?' '" '" 15)

43q.~ , 450

-92.82945 381.58 之 80_90 , S 1 21 0 82.'0>7216 ).4.88�?3 -3.25824

-\-,.ら 7170

-'04.085γ 弓

司 ι2.436'>2

x-~ ( 2) ,柏崎.( 9)

H崎~ ( ! 1))

80.741 弓日

40 (, .27170 君 65.15760413.89140

391.5 日 2 自作294. 品凸 3 &0

56.54.B2 且 1.693ち 1

担 7.47964

180.83670 ム 46.26940

69. ち白且31

自 O. ち 6035

x 町田 1 2l

2l5oS2110 日 61.7\\)0

l 岳 Z 園 9B840

199.93830 9.91763

-235.15840 -6白色 40601帯 94.<ヲ白日口 2

-70.94742 争 14.032 \3

弔 0.67542

-11.7 ラ"白ーち 1.79747

耐 14 ~ 9 自 450

国.1 3 命自9萌 '3

ー J3 .O~4 ワ Z

x--( 守〉

唱副 1.0日 I 百申

明 70.9474.'

ー且 1.21 司%

7 咽白 9 7& 7

243.2ち4 口凸-2 内自 .19900

326 刷 33570

-90.51210

;>94.86.3t!O ち 9. 7:30S 1

-10~ .<f O 呂う O

l fJ ô.9 <'lち口

1 n .81640 12弔 .31390

126.4')740 ーら 7 .44 ;,~.2

町司静 131

" 百四司 (10 J x- 崎 (11)

9.58050 2<).35490

Z 白噂 956呂田

60.4ι8;' 88.27126 凸 1 ~032 ,, 4

5 E1 .602~ 0 51.':>141ι

22 ~ち『事『唖28.09343 28.05914 40~ 34 25'< 主唱 .28255

2 ヨ .0222自

,-輔【"

"園 59653162. 弔 8840

160.22唱6口

44.75145 -0.45688

司 53.94525

-113. 7<155口

-40.54747 -81 ~21996

-�4.26350 司 14.4718 4

64.8品ó2' 12.8647'l

-36.49'102 1. 日70M

-70.13516 同 6 令 :;QS 白

-2S.3A275 ,~.'!.?B 1.唱

-32 .1 41 号ら 品.,うち 1月自

開 24. &1 5 7. 0 -1 唱 .805 t 唱

百四-( 10) -6.09660

-1 4~032 1.] 明 \4.2古 350

時 14 幽 47184

+事4.81558

"日 .12070-66.3181:'> "司ち7<' 16

-5'1.73051 56.54332 l2.吉坦óJ O

-14.岳 29') (J

明23.fう 866B

-2司. 40170~

ó2.ち 48 4 7

:; .2 7 0口口 e?命 98999

n 剛 06 iJ OO 41.730口口

I 口 .41()n 口 71.98999 5 告.0;, 90口。 ~O .7 9999

9.68000 "噌"号号守32.7899申 19.53000

lÜ 司 37000 7~ .98999 -5 7.609 '1弔 , 5.9 弓 000

3.070口口 91~9 (1 999

.1/ 命 379'9 9 1\:10 75 唱 99

I 晶ー.( 5) う〈明四 1 " 日目抽 (12) x--( 13)

(~V • 1111.)1.

1.37弔 &5 "暢 759992.07202 2 ちZ 曽 50000

0.' 且 .305 75 , 0唱。。。

0.445号事 品 7.8900口

0~Sl87 自 72~64000 0.94.311 守 4.98999

0.4 2:品 69 ,..自白'"

0.<19447 80.34000 0.52296 62.64999 1.64712 35.85唱"O. 苗日 759 "晶4岳'99O. 221 7~ 65.&4000 0.23129 ßó場也5999

0.23355 98.59000

0.30956 43.92号'91.74'103 41.00000

X司咽( S 1 珂開予帽, .>8.71109 日 岨白 1.S39 l3

~.917ó3 胆 235.75 自 '0

-0.456且白 四 53.ヲ4ち"

43.24:223 -114.20300 406.271 {O 同465.38210

明4邑 5.38210 865.751品。

244.97480 働 327.10550

幽 263. 甲2. 770 '"個57450243.26400 -283.19900

-94.81 ち号事 ]58.7207口ー74.28929 ".ち4602

71.<16365 司 59.31395

10 ヲ .15790 ー 121.61430123.<f<f470 町 126.24 3'10

122.72'530 -142. 守7且 70

働 103.21940 183.17&50

,-嗣 (12 ) )(--(131 2.13907 -0.9961 弓

キ 11 伺 75436 -~ I .79747 制 36.49907 -10.13516

1.57064 ー‘白 .3う '56?l .9636~ 109.1579口

-5 'l .3139ち 伺'1 ?l .61430

121.73010 225.43510 4幽3 .25 81~ -17. ち 71 ,口

106.921">0 171.01640 ー 14.62 唱。口 司23.68668

輔ら 1 町 14622 明 91.7\ 657 87.47964 91.S 2" 7 ヨ l

'11.117151 180.8362"0 71.57460 ')6.77034 35 ,1 4 '1 7 口 白 7.005百 l

明 12. 守 1234 -21. 白 4374

10 島 53口 00

30.)4 口 00

10 , ;>4 口 00

3 面白 2ち守""

4.60000 19~ 119弔9

11.79000 3白骨唖3哩唖9

3~~;I OOO

x--t 71

河口-( 141

円 1 珂,

。砂 05000

0.29000 O~OOOOO

..唖甲 00口

口 ~OOOOO

Q~OOOOO

1 ~29000 o ~OOOOO Oa4800Q 0&00000 0.00000

1:h29000 1 日圃 75999H, .039弔唖

5 骨 67000

0.00000

宜崎骨-( 1)

2~26005 ~e.40601

-113.7<;1550 11.自"自4

24".974 坦O偽 327 Q I0550

41308予 140-92.fj2945

"岳骨 33 吉 70胃 66.31 目 13

-I3 1h52 品50121.13010 225.43510 125~280号。

137a68460 -72a6口 622

)(~-( 141 幽 3.97626

-・".弔8450

-,ち ~_~827'>}~92829

l?3命 99470

-1'古 a1437012'5.28050 齢笠';4.08575

129.313唖O

-22.47263 -60~67160 fl 057460 %命 71034

145.26 甲 .0&唱。65019

-27 陶 01SZ0

Page 22: Bul Gov. For. Exp. Sta. No. 29:3 1977 171~ 21(;月,の〆“e 4E ム 林業試験場研究報告第293号 2. 主成分分析の考え方主成分分析は,理論的にはβ個の観測変量に対してある線型変換を施して企く新しいm個(ρ>m)

-192 ー 林業試験場研究報告第 293 号

図 3. (つづき〕

恩.. CORREL向Tl ON HATRIX u

x-勧(" x-軸 r 2l X 骨骨( 3) x-世( ,> x-咽・(5) )(--{ 6} 買幽・(71

x伊助 (ll 1.00000 0.76563 。崎31002 0.68 ち" 0.04547 -O~ Zt 9 :S 1 ù~0 1l 60 R酔明(2) 0.76563 1.00000 0.43864 0.75800 0.01676 d・O~27Z哩S 軸O~ 1111 号

x--( 31 0.310口 Z 0.43864 1.00000 0.39345 -0.0017<) -0.144゚ 4 -O~44189

,旬._( 4) 。包6859甲 0.75800 0.39345 1.00000 。φ23875 -0.43195 0.064 自 Z

x__( S) 0.04547 O~0161凸 -0.00119 0.23875 1.00000 e・ 0.7自 470 O.雪唖 '41x柚働噌(6) 枠 0.21 ヲ 37 -0.21295 明 0.14484 -0.43195 -0.78410 1..00000 世 0.54645

x-四'71 0.01160 -0.11119 ー0.44189 0.064112 0.59741 -0.54645 1..00000 x叩曲(8) 山口 .139&8 -0.16562 -0 鳴 16398 -q0J..02唖4948879 -O~67032 O~ 76479 -0 圃 23359

x回世, " -0.0,'白" -0.14075 -0.31366 0.70285 画。~ 57041 O. 早34.14

X-~( 10) -0.084白 3 回 0.06357 -0.14985 -0.21418 -0.62558 O~ 71 731 -O~4335'1

x--( 11 } 0.14555 0.32986 O. う 4722 0" 15289 -0$3申'50 0.2晶画甲苫 ーO~72712耳目勘(1 2 ) O. つ 3057 -0.042日 1 -0.3口'" 0.02226 0.3611'3 情OÐ21 '5 53 0.63914 耳目-(13) -0.00775 情。畢 13122 -Oe41202 国 0.05261 O~40272 -0~30736 0..82402 x司P骨(1 41 -0 ~0343 2 骨Oe04221 明0.16580 0.03615 O~50a65 -0.35476 0~50..17

x-・(l S) 明。圃03040 圃0.138 "1 5 同 0.3C下'" 0.08878 0.13000 ー O. 吉 8260 0981140

x曲"06) -0 命 13070 ~O~12542 曲0.21716 .0 嘩 24557 帽 0.510 '5 5 。園 69360 由。 .39162

X--( 6) x-・, 9( X--( 10) ,.司.(11) 百園-(1 2 ) x-申 (13) x-副砂<1 4 }

,-岨(" ー0 .1.:3 96 且 曲。 ~0363古 画0.08463 0.145ち S OG03057 劇0.00175 ー0~03432

x--{ 2:) -0.16562 -O~14075 -0.0&357 0.32986 ー0.04281 岨0.13122 哩0~042n

x.ー伯母( 3) -0.16398 -0.::5 7366 -0 ゅ 1496雪 0~541 Z2 情 0.30829 -0.41202 -0.16580

x-・( 4' 醇 0.29481 0.0498哩 品0.2141 日 0.15289 0.02226 -0.05261 0.03串ts

x-且( " -0.ι7032 0.102閉 5 ー0~62 ヨ 58 -0.3号 35' 0.38173 O�40212 O. 晋0865

,-・, ., 。回 7641 甲 副 O. ち 7041 0.11131 0~26692 ー 0.21553 世0.30736 -白.35476

X--i 7J -Q~23J5号 。自 93414 -0.433Sl -Q~?Z712 0.63914 O@IIZ402 0.50917 R国情(" 1 ~OOOOO -0~269目4 0.51)215 0~0 1l 13 自 叫0.01163 ー・0.066自唖 帽。 ~228申亘

x--( 9) -0.26申自4 1.00000 世 0. 4 6259 国 Oe65ZotO 0.66574 O~"44107 0.62267

x--( 10) O~5621 5 副D 砂 46259 1..00000 。 $17814 幌 0.20800 圃0.23425 -O~2.7 1l

x-噌(11) O.口日 138 園。 $65240 O~17日 74 1.00000 -O~SB3弔 5 申0.72832 緬.O.S37.n

耳ー-( 12) -0.017 自 3 0.66574 個。 ω2080口 -O~58395 1 ~OOOOO O~ 73009 0.63275

x--( 13) -O�066S9 0.74407 咽0..23425 -0.72832 0~73009 1 ~OOOOO 0&59501 耳骨情(1 4 ) -0.22893 0.62267 明Oe24711 -0.53747 。 .63215 0..59501 1.00000

X曲回05 ) -0~32184 0.88293 -0445351 。 0.61931 0.45057 0.59740 0.69046

耳ー・(1 6) 0.52268 -0~43760 0.92自" 0.lH43 -0 哩 153131 -0~17932 圃O~2.8S7

宵ー-(lヨ} x--( 16) -0.03040 峰。 .U070

納O~l3S45 品O~ 12542 -0~304 <1 4 -0.21718 o.oa878 -0024557 0.73000 砕0~57055

圏、 0.58260 O~69360

0.81140 ~O~39762

-0.32184 0~52268

0.88293 暗D 圃 43760

ー0~45357 0~92676

恥0~61唖" 0.13143 。畠 45057 国0.15381

0..59740 -Od7932 O~6904晶 曲0.24887

1.00000 -0.45634 -0母 4S r;.34 1 ~OOOOO

u RE5ULT tlf PR I純CIPAl CO悦PONENT I¥NALY515 .・

{ρ.C~A. START5 FRO阿 CORRElATION MトTR!X)

1 ,!':I MP ~ "2 COMP. J C門門 P. .. CIJMP~ '5 COMP~ 6 C�P. 7 COMP~

E!GEN VALUE 6.74581 3 ~64 145 1 ~819ó3 O� 15186 0.H410 。 .49483

C白 NTRI6UTI0N RAT10 0.42163 0~ 2Z 799 。 .113 '1 3 0.04699 0.04651 0.030(13 0.02389

cu胤JlATJV正 CONTR i BUT 10M RAT ro 0.42163 O~ 76334 0~61 033 0.8568<1. 0.88777 。 .91166

g CQMP ~ 9 CCMP ~ 10 COMP� 11 CDMP. 12 C口 MP. 13 COMP~ 14 CO 匹 P.

EIGEN. VALUE 匂0.34344 0.25641 0.170司3 。‘ 13867 0.09472 0.056"9

COfHRI6UTHlN RATIO O.由 2147 0.01782 0.01603 0.010凸呂 0.00867 0600592 0.00356

CUI'¥ULATIVE CONTRIBUTION RATIO O. 守 331 J 0.95095 。島'16 697 0.97166 0.98632 0~992Z4 O.守甲 581

15 COMP ~ 16 COMP. EIGEfl VAしU在

0.04261 0.02448 CO品TRIBUTIOM RAT!O

0.00266 0.00153

CUMULATI'IE CONTRIBUTION RATIO 0.99且47 1.00000

Page 23: Bul Gov. For. Exp. Sta. No. 29:3 1977 171~ 21(;月,の〆“e 4E ム 林業試験場研究報告第293号 2. 主成分分析の考え方主成分分析は,理論的にはβ個の観測変量に対してある線型変換を施して企く新しいm個(ρ>m)

休業試験場電算機プログラミンク報告 (8) (研究資料) (黒)[1) -19:)-

図 3. (つづき〉

正 I (jEN V芭 CïCR

1 [C,~". , c口 MP. <1 C l1 M己. � "門P. , co 円P 岨 7 CヨMP.

)(--( 11 0 , 0151 自 -0."6206 0.10276 -0.1 唱'0 四O~04 0 12 。岨26177

X--( ;2) -0. 百 1956 国 O~414 95 O~01923 -0.01653 -0.02880 0.07 自由晶

賢司事, " w口 .10400 。 .06238 事。.,口 783 0.25る 12 0012060 曲 O.~4113

X骨酔, 4l 0.06067 -0 .3418~ 0.14907 0.00630 口 .15459 -0.0804 邑

胃椅-( 5) 0.31546 0.1 白 001 0.131323 0.2811;7 0..04713 相 0024140

x--( 6) -0.2 且537 明。企 13824 目。 .20345 -0.05478 0.08呼 58 0.1457"1

x.畠唱'" o ~34469 -0.11164 0.13698 -0.2627 日 0.10223 -0.26612

正明白'" 情0.18閉" -0.27210 -0 .3 613 白 -0.32998 0.4<;1653 -O~07396

x-由, 9l 0 .3弓5" 和0.05480 0.030口 q -0.06931 0.2晶 009 -口 ~24 191

;.::--( 10) -0.241 自 4 -0.31539 0.26081 0.37187 0 .1口m -o~ 18 468

x-畑 {11 ) "口 .26779 0.10654 -0.14587 口 .0ち'" 。咽 18156 -0.22733

,-・( 121 0.2 ち 132 -0.29353 冒 0.40552 0.04721 -0.42178 -0.1 目 005

x咽同(1.3】 0.2 日号[\ 開。坊 23 1:1 42 -0.11782 -0.11328 由。 .24685 -0.21941

x--( 14 J 司 .26官官 b 岨0 .1 3&22 -0.3<171" 0.56弓 32 0.0998<:) 0~50164

耳慣 -(15) 0.33ヲ 05 口.口 1343 。園 1121" 0.117 日 7 0.5717<:) 0.12'367

x司胆 (16) 値 0.23 2' 97 -0.29420 0.31191 0.376弔 E -0.05:340 司 0.21089

10 CÖ円 P.. p CCP"P. 13 COMP. 14 COMP.

l 、 0.49216 O.Of13()4 ~0.21299 -口 .024 '1 7

" -0.40050 -0.17甲.J3 口 .35088 0.11801

J) 0.06575 口 .17 弓 85 0.12'993 。固 1022 且

" 酔 0.030唱。 。圃 45CJ37 -0.03891 -0.01352 -口 .04728

当} 0 , 44:371 -0.196119 -0.62041 。 .03634 副 0.14430

" 0.270自 I 0*30545 -0.21530 0.10761 。闘 13839

" 悼 0014守自占 o .O~ 3 21 口..2 1942 口 .16159 。幽 127 吉q

思》 0.0490.{ 欄 o.rf8日 l 情0.24841 -0.37439 画 0.0曲目 3 弔

" 。砲'"凸。 -0.01594 0.26220 0.10329 。園 03905

'0' 曲。.口 1725 -0.013口 2 口 .20632 0.12134 -0.67488

'1' 'ロー l'白0; 0.33749 -0.04953 -0.1&守00 事 0.04676

'" 0.29397 -0.10263 。幽 249&5 -0.01050 唱。司02939

13> -0.368 目。 凸 .35789 -0.44913 甲口 .05261 帽。 .13122

,<) 勘当口.,唱 '70 0.10075 0.000口 2 -0.14882 0.02703

15' 。 .0ち 043 。回 035 月 4 0.071.33 0.1229t O~ II 755

)Ó, 0~O ,~9ち甲 担。 .00015 -0~05S09 -0.22126 0.64195

1; 仁白,0 園 " ιロ州,.x--( 1) 日 .04Q10 -0.0225 ,

I時四 , " 制 C 駒口 4043 0.06074

骨口 .00品" 目0.01052

0.0475.5 。 .02347

0.12697 踊口。 1'>571

口金口守由自由 品。 .025ワ 3

口 .187昌 2 ー口園 72507

0.06019 -0.06243

0.49 ,日 4 口 .58671

ー口肩口 2212 -0.0:'345 ー口 .06971 -O.1l45A

-o~ 31876 唱 0.02130

-0 園 1546 日 0.20171 0.30晶Z 弔 -0.1769 ¥

明 0.67白 5 自 0.08355

-0.00293 日。 0 .3 04 白

FA亡T門 R lQÄ円;ト'"

-} (コ同 P. , c刀 MP. 古同P. � U:'~P~ 7 (0 ドP 争

帽 0.16ゎ" 両日園口 2 且 n 001自『ち b

【巧 .01426 -0.つ2026 0.04864 引 0.27 275

-0.04975 ,口ーl<l 926

-0.74ll9 -0.5273.3 唱。 .18648 0.06302 0.0甲 0130.89526 四 0.209 '1ら ー口 .150!,口 0.11 ゚ 7f:1 0.071 喧 I -0.16 岨 54

-0.31314 ü.占 4928 -0.04 ちn

事0.0ち守 CA 0.1fl296 -0.14957 0.32078 O. 口 7159 時 O~ 11419 0.046 日目 0.131 ヲ4 -0.14056

事 0.2 '1 )6'> -0.351 巴 2 0.:)4072 -0';;: '1 67 口 品 0.11133

-0. .3白'" -0.10216 -0.09772 -0.17364 -0~13566

-0.15942 ー口骨 3012.? 口 .J ~76S 0.070':>7 0.31016 曲り .101d (, 0.097213 0.10168 0.40222 0.07646 -0.43774 日早 2704:; 0.32ら" 朝日骨03756 ー.O~16749

l 口 cc 円ロ帽 12 COMP. 13 COMP. 1 咽 COMP.メーー 1) -0.0955 自 。押 24921 0.03092 -0.06555 -0.00596 x__( 2) 041047 時 』口 .06678 0.1 口'" O~02 自 n珂 --c .3) 幽 0.29035 -0 聞 01041 口 .06ち" 0.03949 O. 口 2442,司-( 4) O. りm 日 0.1 日 '02 目。昏 01~ 49 幅 0.00416 刷。 .01129

司 E司'5> 口 .(}68:)8 0.224 白円 -O.081~0 -0.23103 0.0111 日 -O�03451 X輔-,岳 J 。明。4159 0.13713 口 .1':>6?8 -0.08017 。園 II 77自 0.03304 x--( 7) 幽 0.017"8 納 0.07ち" 0.017a9 つ .08171 0.04973 。也03046,-ー{聞} -0.01210 O. 口2A e ち -0.J15;>7 -0.09250 -0.11533 世 0.01 晶 33x 幽-( 9) 。色 l1542 O. 口 225 日 -0.0065甲 O. 口 97(, 4 0.03179 0.00932 〉ピ曹司{ 10) 朝日,り 18-'4 -0. 口)8 8 凸 -0.03019 0.076 日3 0.03735 幽 0 0 1 1:> 111

耳--{ 11) 晶 0.20345 -ー臼 .08408 0.13953 -0.01 日" 帽。 .0 ち 201 -O.Oll1� ば副-( 12) 0.09595 0.14 目目 b 甲 0.04243 。園口 9296 -0.00323 畠 0.0070lx-ー (13 ) 目 0. 1.\ 716 -0 園 18675 0.14797 -0.16747 -0.01619 -0.03 ~31 x-・ (14 ) -0.01805 -0.10065 0.04165 0.00001 -0.04580 D 咽 0064ヨ,-・(l S) -0. 口ら 1':>2 0.02553 0 , 01 <1 82 口 .02 目白口 口 .0:n83 0002806 x_帯 '1 白} 抽口 .01787 口園口 20口吉 ,口 .00口 06 ー口 .02051 白 0.06994 0.1546昌

iS こ口円". 16 CQ~;::>.

1) 2)

" X--( 4) ,-ー( '))

耳司阿南, "

x・-{ 71 X畑-( 8)

0.004'/7

Page 24: Bul Gov. For. Exp. Sta. No. 29:3 1977 171~ 21(;月,の〆“e 4E ム 林業試験場研究報告第293号 2. 主成分分析の考え方主成分分析は,理論的にはβ個の観測変量に対してある線型変換を施して企く新しいm個(ρ>m)

-194- 林業試験場研究報告第 293 号

COIHR 1 ElUT DN 口 F CO門ρO~~ENT T� QRIGI 川 ~l 'v'AR[トBLE

図 3. (つづき)

)(.~{ II )(--( 2)

x_町、"X-~t 4) )(--t 5) )(ma( t;)

X噌司【 7) )(m_{ a)

X司拘'",--< 10JX蜘旬 (11)

)(--{ 12) X伽 -(!3)

,-ー【 1 4J:<--{ 15) X‘-(1(;, ,

x--( 1)

,--<2JX争-( 3)

,ー-( 4)

x--( 5)

】ば一司.( '" J )(--( 7)

x一一'")(--( 9)

;<--( 10) )(--( 11) )(--( 12) X--( i31 ):--(14)

X--( 15 J

,-一ー( 16)

, ca例 P.0 ・ 001550.00258 O~07297

0.024 出 4

0.67127 0.54944 0.801 ら 8

0.23993

0.'ら "00.41439 0.48377 0.43299 0.56385 0 ・ 491680.17551 0.36611

EI COMP. 0.00914 0.01098 O.084J2 0.00047 0.00463 Q�00173

0.00029 0.00015 0.01332 ü.VOF, O 0.1 .3 9ι3 0.1)5315 0.02414 0.00030

-1708828974600155

P434421003

7442SD

H?24228722392865ψ

四位叫旧同町叩日叩回目町um伺出回

・・・

Eq------­

?0000000000000000

01053770831

54ql

P8428191041401974

007113S23578oa--1

000OL000230283060

C0000000200000010

60000000000000000

70Z5339371065446

P528002305926583?

70809211322I9705

62040505800000310

C0000000001000201

••••.•

5000000000000ocoo

-485t72t985064465

pqP9?31lib-064741

N7036414801630693

0009113iqo$121907

cooi0000000030000

・・

φ-a

・・

mom・曲

em--

40000000000000000

20853883625047SI

P-306976740684?35

83726424

l063307

celOA43q30823030

C330200otolo--oo-

30000000000000000

5992?19881-585201

P777a2lm63040l439

M

nvi写

7SEa-724

53507

005066747t392832γ3

C45r

50R020120iooZ

20000000000000000

-49039937962$8493

P076110920920907q

OOGo--00OS000003

田町田加別間関目凹回目印日開師同吋

.•••.•.•

40000000000000000

060233?119116059

P366014430376214B

mmumHmmmH旧日叫目白刊同日目別

C0000000000060000

30C00000000000000

1

-b691738630a45032

P94223465

936C084

044036689S088000

0000050C000002000

C0000000000000000

20000000000000000

901735294L7094

O

ps-20hq320qdR8720

2006650300ql1100

00

03030looi02000

C0000000000000000

100000000000000oo

--314816216767350

P11124876530181b4

肝何叫目的日日目叩叩叩肝日H目的凹

C0000000000000000

00000000000000000

-Z90

9058i5491I352

0477705oqQ23283b3

L218003752198020

C60020η0900401000

C0001000

00000000

90000000000

00000

CUMULAT IVE CONT 尺 18U Tl O ,", cr C¥'!MI'DN[NT TO C: RIGINt,L VARI閉店

x__( 1)

X抽担, ZJ

X均司'")(--( 41 x-ー《告 2)(--{ 6)

,-・(7)X~~{ Sl X--( 9) )(--( 10) )(--( 11 )

,-・【" J ,-ー( 13J x__( 141 ,.恥 (15J

x--( 16 1

)(--( l)

)(--{ 2.)

)(--{ .3) x__( 4)

x__( 5)

x__( 6)

,.目'",-ー【")(--! '9)

x--( 10) X--l ll) x--( 12) x--( 13) )(-~IH)

)(--{ 15)

,ー -(!ó)

1 C�P. 0.001 弓 S

O. 002 ~ 8 0.'J7297 0.�24/:14 0.67127 0.54<)44 0.801ち 80.239ヲ 3

0.65360 O~4143 9 0.48377 o ~4 3299 0~56 :l-8 5 0.49166 0.77551 0.36617

8 ca附 P.

O~86868

0.88991 0&<)9051 0.83443 O~88690 0.909ち 1

0.961;>39 0.9う 5800.96488 0.96314

0・ 92B370.95318 0.89312 0.98090 0.97375 O.'HOOI

-359?785627792550

P3156868574014495

M26498g3153200950

33721593394502704

C8a88889B89788889

50000000000000000

75713

6537837105

P7976844515850

弓667

35Z399320209801

4

00??19985

4491493

CB878788780n778678

40000000000000000

3?417ZE

ヲ452873739

P8689433544884815

mg

12528219540061

09780866598370586

C?85a788687?68

77

30000000000000000

0766357190260008

P33?2

56694200183

MF724067

7583877

4

のυ6793241

白eo--9183

C455

788586756576

20000000000000000

9 ca 円Pc 10 CDMP$ 11 COMP~ 12 COMP. 0.93009 0 恥 99220 0.99460 0.99555 0.6'1110 0.93283 0.98292 0.9a138 0.99129 0.<)9340 0 ,, <)93可 0.997800.96322 0 ,, <)6346 0.999ち 4 0~99915 0.88690 0 圃 <)3738 0.94401 0~99738 0.<)0986 0."12861> 0.94461 0.95104 0.91><:147 0.<)1523 0.91う 55 0.98222 0.96311 0.96433 0.97761 0.98617 0.96991 0.97042 0.<)1046 0.98000 0.96538 0.96574 0.95665 0.<:11256 0.96916 0.91683 0.99630 Oミ <:J9664 0~9629q 0~98515 0.<:J8695 0.99560 0.91193 0 ,, 94680 0.96870 0.99674 0.9BIZ:.'.o 0.9'<1:.'.06 0~99310 O.993l0 O,, 'H640 0 ,, 97705 0.<:1717.1 0.91810 0.97033 0.91074 0.91074 0.97116

仇脚脚mmm剛帥附剛山帥脚脚剛脚剛

印加凹叩向田叫凹凶凹凹叩叩伺同目印

6lil--22Iza--2131s

9209183073890338

榊刑制附叩叩刑肌刑制問州側叩泌隅山…

C9909999999999999

50020000000000000

123456T890122456

xxxxxxxxxxxxxxx

-430778ossa-44?077

P59292716S6269659

98632?6511908039

回目肝刊回開柑嶋田明拍可咽帥叫肝帖

.

.

.

.

. 。

.

.

.

.

. .

70000006000000000

369995269084?922

P15?4960526425379

2611999398980471

03733593E324685864

C888888q999788899

•..•••••••••••• 60000000000000

。。。

13 CC阿 p~ 1<1 CO I'IP~ O~99985 O~9q98e 0~9990S O~99984 0 因 99940 1 ~OOOOO 0.99976 O~99989 0 ,, 99750 0.99870 0.99848 0.99957 0 ,, <)8470 0.98562 0. ," 99411 0.99975 O~98101 t)e9Bl09 0~97395 0・ 999910~99935 O~99947 0~99S61 0.99566 0.99100 O�99i98 0.<)9519 O.99SZ4 O~91953 O~98032 0.97605 Os9999f.!

Page 25: Bul Gov. For. Exp. Sta. No. 29:3 1977 171~ 21(;月,の〆“e 4E ム 林業試験場研究報告第293号 2. 主成分分析の考え方主成分分析は,理論的にはβ個の観測変量に対してある線型変換を施して企く新しいm個(ρ>m)

本卒業試験j号電算機プログラミング~'[l~'(8)(研究資料) (黒J 11) -195

t6J 3. (つ dゴき)

,0何PONENT SC円 R[ (G L.~I';( !E:U

1 CClMP. 2 COM゙ . .3 C" ,~P. 4 UHP~ '" (('!I<!p$ ー〈俳押( 1) O.口 0584 。陶 17464 時 0.34255 口。 11851 申 O.2B13

耳目-( ;:) -0.007ヨ 3 0.20513 -0.307 自i 。眉 02218 品。 .01 ヲ 17

百四咽闘【 " -0.04004 0.19419 口 e 口4白" 司 0.5 且吉6岳 0.29691

z曲明 1 " 0.02336 口 .2061 ぢ -0.25341 O.i 71 中2 O~00962

耳 E・ 1 5> 0.12146 口 .07111 0.11 号oò 0.1 ち'" o ~32品 53

X.四 161 -0.10937 -0.141¥57 .0 骨 10248 -0.23463 ー 0.06350

z咽ー-( 1) 0.1.3271 市0.05756 情0.08276 0.1519 早 開0.304畠4

耳唱圃冒【 81 '0. 口 n61 曲口 .14447 隔世 .20171 司 0.4167 う -0.38253

K 開園1" 0.13695 -0.04934 伊 o ~04062 0.03470 働 0.0且 0 3<

,.帽(0) 由0.09542 -0.12口" -0.2.3381 0.30078 0.43109 ,.骨 (11 ) 帽。拘 10310 D 噌 13132 0.07自" -0.16822 0.06 .3口O

百科悼 112 】 O.守雪 75 .3 唱。 .07996 -0.21760 ー 0.'白 7古 7 0.05473

x--( 13 > 口 .11131 ー0 .1 0004 押口ー 17675 輔 0.13587 -Od3133 百ー-(H) 0.10394 -0.04371 -0.10司9 閉 -0.40063 0.65536

x--( 15) 。園 13054 -0.027日2 日 .00" '1 6 ロ.J2唱" Ool3�1>5

x--( 161 明日 .08'170 曲目 .133 '12 -0.2 ¥ 810 00359'11 0.43'702

@白 00荊 PONENT SClJflE ••

4 C:)f'P.

9 10 11 12 13

"‘ 15 16

" 18 曲 l ーヲ39i4

" -1.14795

20 -0.11842

1I 同 1.72243

22 同.2 .7 5916

23 -2.03188

2̂ 0.11116

25 0.37272

26 O. 阜 1182 0.23弓" l 市0< 1 呂う -0. 口 e3 ?7 2.52701 27 O.U0831 1.0284 岨 l 脅 31696 申。骨 76369 l 勾 55266

" -0.12351 l 舟 93g6B 崎 1.1 凸3吉 b 市 1. 守6353 0.25333

" 旬。 .267口口 2.43366 -2.03471 O. 口 0528 -0.44303

30 。喝 Og551 0.12262 0~3105l 咽 0.2577唱 O. 口事 1223J 1.14<l3.3 件。 010023 0.86035 0.' 日 703 0.26649

JZ 1.' 中>0, -0.2白O 日 3 0~036ラ 2 明 1.111 申 5 2.48444

" 団 0~ .5 7970 恥0.49511 0.40075 -0.731 口目 掛 0.8741 B

" .1 回 44791 曲。 .49589 且幽 0.1524 占 ー 0.5 <;l V1 噛0.14319

35 Od3白" 1. 品2249 ー1.77160 O. 岳2468 平 0.02073

" -1.26734 。自宅~5646 0.76276 ~0.657 雪晶 。 φ24919

31 1. l3 44 ヲ 口 .90071 。幽 34100 2' .69615 38 画。 e も 6287 0.0592') -0.655')5 0.38471

39 0. ., 14')4 戸。 .4 1'> 739 司 0.8471 自 -1.12040 4つ 0.14512 -0.72928 -0.220ち色 伺 O. 申占'"

" 口ー 77775 帯 1.OB21 -1.01150 世0.21063

" 0.94604 南口。 55553 -0.250唱 3 1.07 (:, 58

" -O~lóZ3 唖 間 1.33214 早口 .6守白聞3 1027038

" -1.25495 -1.5 胃 5'0 0.88125 O~94757

'5 ー-C.15433 -D.ó3ち" 田口 .40977 0.17537

" 開 '.1 白日目 b 開 1$75621 1.776-29 1 ,, 96602

" 母国 1.422/2 開 1.17736 0.'日'" 0.6自 775

" 0.51138 0.0 日 413 0.17010 咽民 2.3078自

49 -0. 占 4598 .,傘 55711 1.00中唱。 0 ,, 68525 50 0.67371 司 0.37088 酬。 ~50999

" 口 d ち 02' 0.7aóヨ l 0.052135 白-0.41875 -0.00609

52 0.60812 口ー 1 司 0'5 l 骨 18632 i 骨m つE 0~13663 53 0.9守 4 � ~ 幽 0.079つ9 1.23038 i 噸 Z 甲 671 。事告 102申

5' l.l 呂 01 雪 0.27242 0.62161 ), .39294 0.79316

55 1.:72')9 0.71380 0~07661 1.04743 。自 2587弔

56 O. 唱 347';' 0.1:3513 口 .97151 1..73442 唱。O.‘2341b1g54i

" 0.72884 0.15177 0.01109 -0.09 日 n

58 0.132白雪 0.20450 1.~ 026 5 1.46253 0.7坦'"

" 0.58754 。 .17369 制口会 1981 且 。令 54430 -0.17543

モO 0.94632 同 0.S8462 0.14987 0.04 ,318 聞 0.62602

" 0.544 昌弔 0.0¥035 0.21196 骨 0.' ち6 口早 0.0723哩

。z 0.57330 叫 0~5 219 2: 0.'型 908 -0.43352 -0.1429唱

" 口~15275 口 .42894 O~03112 -O~32~口 1 0.40331'

" 口 .61211 0.01245 0.522 百 5 。~ 788自7 -0.39503

65 0.98270 -0.486&0 0~2 8812 0.69422 “ 0.33361 66 1 .l 06~ 5 -0.53314 0.45190 口 Ð487 口 4 0.111977

67 0.'占 1729 守0.36238 0.14346 o.雫7守昌吉 0.8234也

" 0.88504 -0.24401 0.05592 0.8B7自6 0.64233 68 -0.4 岳 125 -0. 6.09弔3 0.10368 。骨ち ..946 位。 .39313

70 0.6 [1 0 日 5 酔 0.01412 0.50055 0.31076 戸。.2: 9643

71 0. ,, 1 目" 司口。ち 1730 口.,ち 536 0.029[17 -0 圃 53928

72 日,弔7487 -0.7351'f ーロ.3 15 \7 同 l 骨唱 4534 口 .856!l 0

73 ー口 .5岳807 -1.449甲4 -0.73ち" 持 1.34874 0.02419

" 事 l 園 3107 争 -2.23035 日 1. 甲骨'" 1.48317 2.700口7

75 O. 咽 4953 同 1.00958 曲 1.0~ 8う Z 自 3.06S l7 1.5190ヲ

7� 曲 1.<1153 日 0・ 1 .-1 ~'/43 -1.42861 1.2181 J 1.47020

77 掛 1.64019 四 2 .l B 口口3 -2.01320 1.1592 号 2-.41140

78 。抱 92499 0.054 う 2 0.45253 0.63622: 。舎 16779

79 0.61480 司恒 0.30546 口 .55275 0.93653 -0.25294 80 1.13441 -0.9367口 -0.50366 -0.5829邑 。命日百'"

81 僻 0.22709 輪。姐 36070 0.3666 時 時0.71002 -0.5B934

" 0.38029 明。争旦4品77 -0.0ち"申 .1 轟 24413 幽 0.23201

83 0.40113 -0 明 34晶 16 与・0.14029 -0.52447 園 0 .3 6835

" 帯 0.64264 -0.58056 -1. 岳唖062 -1.0646自

85 -0.703品 7 情 1.24995 -0.20764 輔の .84324 -1 ~l [1084 86 0.06弔ヲ 3 -0. 唖 7360 -0.63233 四之 .54693 -0.38..183 包ア 匹 0.28073 -0.90884 幽0.429弔 7 幽 1.' 耳目 05 -1.29675 88 0.40390 -司 .453阜 Z -0.46168 伸 0.22159 -1.27:832

" 0.27869 -1.7.1574 -0.60150 明 1.598 弓号 抽 1.22265

90 O. '4252 自 .14150 0.20702 口."口" 口 .37041

91 。骨 50084 -0.4068守 0.1 ~4 21 0.2号'" -口 .91462

" 0.91692 ー申告"日" O.40~ ♀ち 0.492唱。 掛 0 , 23 品 7Q

93 。眉 65372 。園凸 1281 0.04767 -O~27994 日.21'!7事4

9' l 傘 12003 冊。魯 36406 0.20458 0.13997 口 .56318

" 1.00-.31 -0.354阜6 。司 01614 O.口'"呂 -0.5 色!31.. O. ?i41 B O. 守 1855 -0 ‘ 336守 5 -0.58315 0.62710

" i 電 15249 明。 .32482 0.14676 0.07263 -0.05071

" 0" 81 !l~ 9 0.192 弔5 0.260昌7 0.58112 0.14725

" 。局 38380 O. 4~ 85 2 0.30816 咽 0.0216守 神口 .18225

100 0.18%4 -0.110691 0.2247Q 約。刊日 9601 -0. 3 <;i~ 97

Page 26: Bul Gov. For. Exp. Sta. No. 29:3 1977 171~ 21(;月,の〆“e 4E ム 林業試験場研究報告第293号 2. 主成分分析の考え方主成分分析は,理論的にはβ個の観測変量に対してある線型変換を施して企く新しいm個(ρ>m)

-196- 林業試験場研究報告第 293 号

ぼ1 3, (つづき)

101 。書甲骨082 0.3日 577 0 ,, 55160 Oelll92 1. ヨ:5 413

102 0.76424 0.41123 0.636昌吉 ー.0.27457 0.63003 103 開 0.3591 4 ~O .20918 1.115186 2.00683 -1.24194 104 0031495 O~17107 O.969Sa 1.J5zR5 -0.62016 105 0.2 串5也 5 0.22017 0.94163 0.444号。 -b。oe圃 b白O87b8468

'" Q.I744b 0.51151 0.472:67 1.037 且?

>107 冒0.16523 0.63757 0.61岳 13 -0.23545 0.139告唖

108 ~O.1Z5S4 -0.4守 551 o.aH21 0.77741 -0.55837 109 0.70706 0.51112 0.46841 -0.20892 0.74416 110 1.0565 自 0.2守'" 。園 26104 ー 0.23守守岳 1.10675 111 0.177自 3 0.29861 ー0.6 日 530 柵 0.65291 。晶口2420

112 0.62016 0.431'10 0.8�244 1.05339 0.13139

!l3 -0.132181 O. 守 2871 1.49754 1.55396 -0. '19983

'" 1.15.3'14 -0.48535 0.5閉"品 0.60217 酔 0.36205115 0.58076 -1.04267 明0.30965 -0.">6812 -2..71527 116 O.l)4611 曲 0.12593 0..25152 0.71.%2 0.14117 117 。ー 720 f>~ 0.< ち0口 q 日駒 .)3953 0.89814 0.31100 118 0.52087 0.31519 0- .2')[>4 口 O.B3Sf>9 -0.24208

'" 明0.69359 -0.03834 -0.411 ~ ~ -1.35391

1" 0.59959 恒 I ~ 25623 -0. 5389~ -0.64.H2 -1.98000 121 0.34171 -1. i ~O 自由 -0.03193 ー口.品品事 98 勝 1.42458

122 0~05 905 -0 ~55438 -0.21864 -0.ち 7764 -1.4551'1

123 -1.03415 1.07344 0.25329 0.20032 0.802口?

12< -0.99323 四 O~ 41460 0.01759 開。 .81 '>04 Iミ5di72

125 -0. 号 ~9 'l 6 -0.12016 -0.26372 0.46007

l2� -2.13113 司-0. 甲4 ち 20 白 ~21810 1.137自ら 骨口 .72138

127 -1.6405 古 -0.446'l8 0.2�136 -0.0ラ "0 ~0.17 !1 67

128 -2.2白~" 0.78107 。園 b号oeó -0.17452 -0.83055 129 -1.7 自4 日 Z ザ口。 l 町 261 合0 .1 1373 ー口 .19448 口。ち1315

130 -2.88616 -0.73 口 51 0.0"238 1 .4 7 3~6 O~ 234611

@・ COHÞON正 NT SCORE SCATTlR DIAGRAM ー・

10000000000000000000000000000

0864208642OBb

・-2086

・。208642024

--------B-a------aaa--ee

e

・-

544i

喝44333332222zi--210000000

1Il /"1

叩回目刊拍凹回目柏田抑制帥叫四回目的叫回∞回目叫担問目刊

544a

445333322Z22l12110000000

-0.60 1 / 1 11"'- (1 1 .........._ ........ 、-〆ケ 1 ・0.60

4・ 80 1 ! 1 ¥ : ¥ ""-... 一一-- .1. I -0・"~1.00 1 ¥ : 1 ""-.. I -1.00 ・ 1.20 1 11 i ¥: 1........... r ・ 1.20-1 ・ 40 1 ¥: \、'- I 叶.<0'1 ・ 60 1 \1 川\ 2 ----... 1 -1 ・ .0-1.80 I 1 1: \ リ 3 吋.,。~2&OO 1 い! ¥ I _.-/ 1 -2.00 ー 2.20 1 、'- I ¥ 1 ____ 1 -2~20

-2.40 J ..... ""-..・\ ---亡 I -2.40 -;:;� i J[ ~ / ¥ 1~ II � -2.60 4・ 60 1 '-.... 1/ - 1 ~2 ・"勘 3 ・ 00 I 、、・3.00

-3.20 t I -3.20

-3.40 1 1 -3 咽"・3660 1 1 -3.60 “ 3G80 i I -3~aO -<・ 00 1 -4.00 -4.20 1 噌 4*20

-,・" i 1 _" -4~60 1 ・4.60~4~80 1 1 -'*.80 -5 ・ 00 1 ・ 5.00

。白5 ・ 00 -4.00 -]園 00 -2.00 ‘ .00 0.00 目。 o 2~OO 3.00 4.00 5.00 Z CO!'lP.

〈制<< PC A,/45 END OF JOI:l >>沖〉

Page 27: Bul Gov. For. Exp. Sta. No. 29:3 1977 171~ 21(;月,の〆“e 4E ム 林業試験場研究報告第293号 2. 主成分分析の考え方主成分分析は,理論的にはβ個の観測変量に対してある線型変換を施して企く新しいm個(ρ>m)

リストず tユグラムソース図 4.

0001 PAGE SOURCE PR口 GRAM LlST FORTRAN QKITAC 4500

雲仙附弘一怒釘刑判落川判。ロ山、叫んド\、設咋(∞)〈翠持山政主〉のたこ一)

PCAOIOOl PCAOI002 PCAOI003 PCAOI004 PCAOI005 PCAOI006 PCAOI007 PCAOI00B PCAO 1009 PCAOIOIO PCAOI011 PCAOIOl2 PCAOI013 PCAOIOl4 PCAOI015 PCAOl 口 16

PCAOI017 PCAOIOl8 PCAOI019

(PR自己RAMED BY Y.KUROKAWA)

例 AIN PROGRAM

CO州問巴 N NO 曹 N 1 ,t>l2 , MO , Ml t M2 ,M3 , 1 T, 1 P t 1 R, 1 c , r G , 150 d 5 C口問問ON ISE~H , YH , X~DfV?YtA州, AN , A , 8

INTEGER J5[(30) REAL X(40} , D(40) , U(40) , Y(30) , A門( 30 ) ~ AN ( 30) 曹 A(30 , 30) , B(30.30)D口UBLE PRECISI 白 N H(40) ,HH(40) JT: 1 CALL 5U8501 ($30) CALL SU8502 ($30) CALL 5UBS03 ($301 CALL SUB504 ($20) CALL 5U白 505CALL 5U臼 506 ($10 ,$201 STOP END

STATEMENT

PRINCIPAL C由州PðNεNT MトLY515

10

20

30

rbFLFLPLPL

L!NE NO.

000] 0002 0003 0004 0005 0006 0007 0008 00口9

0010 0011 0012 0013 0014 0015 0016 0017 0018 0019

ト品

民主〉....,

Page 28: Bul Gov. For. Exp. Sta. No. 29:3 1977 171~ 21(;月,の〆“e 4E ム 林業試験場研究報告第293号 2. 主成分分析の考え方主成分分析は,理論的にはβ個の観測変量に対してある線型変換を施して企く新しいm個(ρ>m)

llHU∞il

図 4. (つづき)

斗山骨仙神山外調酔ぶ

gu出世間玲w明日明ωJ'

0002

PCA02001 PCAOZOOZ PCAOZ003 PCA02004 PCA02005 PCA02006 PCA02001 PCA0200S PCA020口9PCA02010

PCA02011 PCA02012 PCA02013 PCA02014 PCA02015 PCA02016 PCA02017 PCA02018 PCA02019 PCAOZ自20

PCA02021 PCA02022 PCA02023 POA02024 PCA020Z5 PCA02026 PCA02027 PCAOZ028 PCA02029 PCA02030 PCAOZ031 PCA02032 PCA02033 PCA02034 PCA020:SS PCAOZ036 PCA02037 PCA02038 PCA02039 PCA02040 PCA02041 PCA02042 PCA02043 PCA0204'早PCA02045 PCA02046 PCA02047 PCA0204B

PCA02049 PCA02050

PA品E50URCE PR 口GRA問Ll 5T

5UBROUT INE 5υB501 ( * )

CO判例 ON NO 事 N 1 ,N2 ,判 O. 何1>何 2. 例 3. !T , 1 P, 1 R d O,j G,j 50 , 15 COMMON ISE9HtHH , XtD , U , Y , A 円 tAN9A , 8INTEGぞ R 15E(30) REAL X< 40) ,D( 40 )tU( 40) , Y( 30) , AM( 30 ltAN( 30) ,A( 30030) ,5( 30 ,30) REAL DD(4) DQUBLE PRECISI口N -H(40) , HH(40) , H自

DATA DD/4HPR目白 t 4HNA 門E , 4HTRAN , 4HST百PhHB/BH5 REAO(5 , 100)HED.N. 問 9L , Y , (X(j) , j=1 , lS)100 FOR 阿 AT<A4 , 2X , 14f 15 , 212 , lX ,15A4)

00 10 K=hJT IF(HEo.EO.DD(K))6口 TO(1~213 , 4 }t K

CONτINUE 60 TO 5 NO=N 拘0= 門

lP=L MT=I IF(NO.GT.40)Gロ丁目 4

WRITE(6 ,110lCX(Jj.J=I , lS1 ,NO ,MO 110 .FORMAT(lHl , 44X , 42H. ・... PRINCIPAL CO例PONENτANALY5J5 ..・ 0..

lIII11HO , 18H.. PROBLE門的A何ε. .曹 4X , 15A4t31X , 12Hホ事 PCA/45 本キ2111HO , 5X , 19HNU伺 BER OF VARIAßLE5.4X.1H= , j5/1 村 .SX.9HNU例 BER 目F.3BH SA例 PLES , 6X , IH=.1511111HO.26H.. ORIGINAL INPUT DATA •• 111 !T=2 60 TO 5 I F<阿 .NE.((NO・-1 )/7.1) )GO TO 40 READ(5 ,120)(H(J).J=1.NOl

120 FOR 例 AT((10X.7(2X , A8111DO 15 Jぉ ltNOIF(H(Jl.EQ.H白 )6白丁目 40

E口 N Tl NUE凶RITE(6 , 130)(J , H(Jl , J=1 , NO)

130 F口R 問 AT((lH , 13x.7(SX , 12 , lX 予 A81) ) IF(MT.EQ 句。 )G口 T 口 18

REWINO DISCO REWIND DISCl 00 30 1=10 判。CALL SUB510(3 ,O,O.l IF( 阿 T.EO.01GO T口 6

READ(OISCO)(X(Jl ,J=I ,NOJ 品目 TO 35 D白 20.J=1.NOCALL SUB510(2 ,Q,X(J)) C口NT!NUE1Ft IP. 在日 .ll60 TO 25 凶 R!TE(6 , 140)j ,【 X(J) 守 J=l .NO) F口R向 ATCIH , 110 ,3X.7Fló.5/ClH ,13X.7flð.S))

FDRTRAN OK! T A� 4500

STATE 問 ENT

10

20 35

140

l 串

15

6

C

C

L1NE N口@

00.01 0002 0003 0004 0005 0006

0007 0008 0009

0010 0011 0012 0013 0014 0015

0016 0017 0018 0019 0020 0021 0022 0023 0024 0025 0026 0027 0028 0029 0030 0031 0032 0033 0034 0035 0036 0037 0038 0039 0040 0041 0042 0043 0044 0045 0046 0047 0046 0049 0050

Page 29: Bul Gov. For. Exp. Sta. No. 29:3 1977 171~ 21(;月,の〆“e 4E ム 林業試験場研究報告第293号 2. 主成分分析の考え方主成分分析は,理論的にはβ個の観測変量に対してある線型変換を施して企く新しいm個(ρ>m)

主総uh懇話同日制下落パ判。口、叫んに山正社日。(∞)(~号待成学)(川河)二)

0003

PCA02051 PCA0205Z PCA02053 PCA02054 PCA02055 PCA02056 PCA02057 PCA02058 PCA02059 PCA02060 PCA02061 PCA020�2 PCA020�3 PCA02064 PCA02065 PCA020�6 PCA02067 PCA02068

PAGE 50URCE PR口GRA門Ll ST

与TATEMENT

IF(NO.GE. 8IWR!TE(6d501 IF(MO口( 1,251. E 口 .OIWRlTE(6d501

150 F口 RMAT( lH ) よ 5 CONTINUE

WRITE(DI5C! )(X(JI ,J=1 ,NOI CONTINUE ENDFIL王 DISCl

lT=3 G白丁目 5

IT= 門

RETυRN WRITE(6.J5S1 FOR 問 AT(//l HO.SX , 22 i1 NA門 E CARD 15 lNC自 RRECTIG口 TO 99 WR lTE( 6 .160 1 FOR~AT(lHO , 9X , 29H<<<<< PCA/45 END OF JO同>>>>> )

RETURN 1 END

FC、 RTRAN。KITAC 4500

3 口

40 155

4

160 99

LHせに N白.

百 051

0052 0053 0054 0055 0056 0057 0058 0059 0060 0061 0062 0063 0064 0065 0066 0067 0068

にD<.D

Page 30: Bul Gov. For. Exp. Sta. No. 29:3 1977 171~ 21(;月,の〆“e 4E ム 林業試験場研究報告第293号 2. 主成分分析の考え方主成分分析は,理論的にはβ個の観測変量に対してある線型変換を施して企く新しいm個(ρ>m)

g 図 4固くつづき〉

勢鴻刊内調蹄節也市淀帯同時

皿叫出叩ωh叫e

0004

123456789912545678901Z345678901234567

901234567890

0000000ool--11111llZ22222Z222333333333344444444445

nunvnvnunv内υnvnunv内UAU

《υ内υnunu内ununuAUAυ内UAunvnv内uvnunu内UAV内uvnu

《unu内MV《ununU《unvnM《ununuvnuwnuAunuw

《unuw《U

をJwzd宅J官iw烹d官、“τd胃JV胃ザ宅d宅AW定d望d官d電JW烹d官d官、“定d胃ぽ官J驚4241d胃JW窓d電J'電iw宮AU電d怠dτdをiw理J1dzdzdzd胃J'電JvzdミIvtd電JV烹世定4zdtd官叫買J申

《υ

ハυnνハunuAUAUnu《unununuAU内νnunU《ununuハU内unv《U内υnuwnuAVAM《υnunu《υ内U内vnununuAU《Uハunununυnu内υAunvnvnU内υ

A門局内amaunAURA内良川AURaAAma舗内品開amA内&例A内凪幽na舗na内aHA

阿A岨naanaanA幽nA

伺a肉食mamamAMnA幽門&MnAMnaM門厄個nA開局斜身内-ma"局内AM門AMnamAnamAMna内aR

ECCCCECCECECCCCCCEECcccccccccccccccccccccccccczccc

pppppppppppppppppppppppppppppppppppppppppppppppppp

PAGE SOU且 CE PROGRAM L15T

5UBR口UTINE SUBSI0 (N.IN.RNI

蜘------ーーーーーー-ーーーー・・ーーーー・・・­LOGICAl INTEGR.MINU5 ,FRACTN DJMENSION JCHAR(13) ,INBUFF(80) DATA JCHAR/IH .IH. , l 付申 .IHO.IHl.lH2.1H3.1H4 , lH5.1H6.

llH7.1H8 , lH9/.ISTflH本/&0 T自( 192 , 41 事 42 )ON

INTEGR=.TRUE. G司 TO 3 INTEGR=.FAL5E. 何 INυ5 =.FALSE.

FRACTN=.FALSE. 5U問 1 =0. 5U何 2=0"15WTCH=1

4 IF( IC口l.EQ.80.AND.15凶TGH.E由 .2)60 10 120

5 IFI I C口L ・ lh 畠由 IG口 T由 7 6 READ(5 ・ 5010)IIN8UFF(JJ , J=1.80J5010 fORMATI80Al)

IF( lN8UFF( 1 loE自 .15T)6白 10 6

ICOL=O ICOL=ICOL+l DO 100 JC白L=l.13lF( lNBUFF( lCOL l. 正岳 .JCHARCJCOl) )GO T白 10 CONTINUE GO 丁目 42 lF(JCOL.GE.4lG口 T白 14IFIISWTCH.E白 .2)G口 T口 135口 TOI 5.110.20).JCOL

60 T口 1120.110.42 ), JC口しADD=JCOL-4 JFIFRACTNl 自由 TO 15 su 例 I=SU 関 l 事 lO@+ADD60 TO 21 su阿 2=5U附2ホ 10. 今 ADD

!5UM2=ISU例 2'1<>0 TO 4 IFlFRACTNIG窃 T口 42FRACTN=.TRUE. ISU例2 =0 GO T口 21

問 INU5 =.TRUE. ISWTCH=2 GO TO 4

IFI 阿 INUSISU例 1=・SUMl!N怠 SUMl

!FI FRACTN)G口 Tð 4Z

6口 T由 50IFC INTEGR)G口 T田 22

FORTRAN 。 K iTAC .4500

STAT 巨何 E f>l T

E

2 3

1 自0

110

120

15

20 21

10

!3 14

22

C

7

UNE 何回@

0001 0002 0003 0004

0005 口 006

0001 0008 0009

0010 0011 0012

0013 0014 0015

0016 0017

0018

0019 0020 0021 OQZ2 0023 0024

0025 0026 0027 0028 0029 0030

0031 0032 0033 0034

0035 0036

0037 0038 0039 0040 0041 0042

0043 0044 0045

0046

0041 0048

0049

0050

Page 31: Bul Gov. For. Exp. Sta. No. 29:3 1977 171~ 21(;月,の〆“e 4E ム 林業試験場研究報告第293号 2. 主成分分析の考え方主成分分析は,理論的にはβ個の観測変量に対してある線型変換を施して企く新しいm個(ρ>m)

葉山神科婦時革岡崎'蒸け

N。口、叫んて日可,,窓山村〈∞)(翠流域主)(知山}二)

0005

PCA03051

PCA03052 PCA03053 PCA03054 PCA03055 PCA03056 PCA03057 PCA03058 PCA03059

PCA03060 PCA03061 PC A03自62PCA03自 63PCA03自 64PCA0306S

PAGE SOURCE PR口GRA 阿L1 ST

IF(!SUM2l40 ,40 ,30 DO 130 I=J.jSUM2 su阿2=5UM2I IO.RN=5U阿lφ5UM2IF(MINUS)RN=明RN!FC .N口 T.FRACTN)GO 70 42 5自 TO 50 ICOL= 80 G口 TO 50

42 切 R lT E( 6 ・ 6010)(!拘BUFF(J6自 10 F自問抑 ATOHO. 雪 x ・ 9HDATA

113H 15 INC口RRE CT!l HO STOP RETURN END

FORTRAN OKlTAC 4500

STATE阿正問 T

30 !3Q 40

41

50

llN正純自省

0051 0052 0053 0054 0055 0056 0051 0058 0059 0060 0061 0062 臼063

0064 白065

"" o

Page 32: Bul Gov. For. Exp. Sta. No. 29:3 1977 171~ 21(;月,の〆“e 4E ム 林業試験場研究報告第293号 2. 主成分分析の考え方主成分分析は,理論的にはβ個の観測変量に対してある線型変換を施して企く新しいm個(ρ>m)

o t:Y

図 4. (つづき〉

業株民間掛釘宮河建叫国叫回申印町山

0006

-2345678901234567890123456789012345678901234567890

00000000011111111112222222222333333333344444444445

内U内U《unv内υ内υハU《U内U《υnvnυnv《υ内υ《unvnVAunu内unUAunu《UAVnυ内υAυnV内vnuAUAununu内U内U内υ《UAununuw《υnυ

《U内unυnu内υ

瓜苛aa噌a匂必崎Am可凋崎a品Ta句一必句aa守a匂必吟凋噌a句au司aa司a喝a匂a匂A拍va崎4崎必匂au守aqaa守a怜a句a崎必匂aa『a崎a句aa司au司A匂

a句aa可凋崎一a幽守S句a嗣可a邑守a岨守A時aa守a崎a匂d句aaマ

Avnu《ununu内Unu内UAUAU《ununu内ununu内ununu《υnu《vnununuAVnvnuw

《U《υ

仇ununununvnunvnvAununu

《υ

内υnυ内unu

《unuAυnu

AMnAma“nnma仙nAma内a邑na内anau内a品naanA例aunAMna

内a《a内魚肉a幽nauna幽na幽na内anA門店《俗的風内処幽門働側A紬na例局内ama晶na鍋円瓜岨na内店内a幽門A内乱州aumA飼HAMnAuna

列島内

FLPLPLP白wp目MPRMFLF」p-wpLPLW円LwpiwpしPLP也、rLPehniwpLnivppbpiwpup-MPLP

』pupLFL相phup」FL匹、PLPLnLVFLF回WPLPLvpivpLWFLPLFLPLwpiwph

句F、

p

p

p

p

p

n

r

p

p

p

p

p

p

p

P

P

D

-

D

r

p

p

p

p

o

r

o

r

p

p

p

p

o

r

p

p

p

p

p

o

r

p

p

p

p

p

p

p

p

p

p

p

p

p

p

p

p

FA'もE5DURCE PROGRA抑Ll ST

SUBROUT IN E SU8 5 02 ( * )

C白刊 MON NO , NltNZ , MO ,似 1~M2 , M3~IT'IP~IR'ICtIGjI50 , ISCOM 向。 N 'ISE , HtHH~X , O , UjY , AM , ANtA , B INTEGER lSE(30J , K(201 , L(20) ,阿 (20 hN(20) REAl X( 40 j,口 (401 ,U(40) , y( 30 J. A 阿 (3日 hAN(30hAC3白書 30) ,B( 30.30 l REAl 6(20) ,C(20) DOU 日 LE PRECIS]ON H(40l , HH(40l ,付 HH( 20) .H8 DATA G/4H+ 曹 4 H- , 4H本,4HI ~4H+C t4H本C ,4H/C

14H**C ,4H** 宮 4HC司 t4HCI ,4H*C+ t4Hホ令市 , 4HC* 本 .4HLN 24HL口 G ~4HEXP , 4HSQRT , 4 汁ABS , 4HSIGNf.HßI畠H I Nl=NO 15 =0 ]50=0 NOl=NO.l 白口 10 J=NOl ,40 H(J)=HB IF(IT.GT.20)GO T口 4lF( lT .EQ. 0 】 G口丁目 ZDO 20 1=1 r IT READ < 5dOO)U( 1 j,L( 1 ), M( I)oN( I), C( lh!H IJ 曾 HHH( !)

1 官 o F口R 問 AT(lOX9A4tlXt3IS , ZF10.0t2x , A8)20 C口 NTl NυE

D口 40 l=ldT 0030 J=I ,20 IF( U( [J.ε日 .G(JllGO 丁目 I C口問T1 NUE

GO TO 4 K( 1 l=J

J=L( 11 IF(J.LE.NO.OR.J.GT.401GO TO 4 íF(MAXOfM(!) , N(!)).GT.NO.口 R.M(!1.LE.01GO 10 4 JF < J.EQ.7.ANO.C(!1.EQ.O.1GO TO 4 lF(HHH(] ).EOφHB. 白R.H(Jl.NE.HB1GO T口 4陶 1=阿 AXO( NO ,L( 1)) H(JJ=HHH <])

E口 NT!NU正只王凶 !ND 015Cl REWINO O!SCZ D口 60 KK=!.MO READ(DI5C1)(X(Jl ,J=I ,NO) 臼045 J=NOl ,40 X(J )=0 聞lF( !T.王Q.O)GO TO 3

DO 50 1=I>1T CAlL SUBS2Q( X.K[ IltU [), M( 1) .N( I j, C( 1) , D( 1)) CONTINUE WRITE(DI5C2)(X(Jl ,J=1 ,NI1 CONTINUE

F口 RT 尺 ANOKlTAC 4500

5TATEMENT

10

45

nuumM

>

εJ?3AU

30

40 2

C

LlNE NO.

0001

0002 0003 0004 0005 0 口06

0007 0008 0009 0010

0011 0012 0013 0014 0015 0016 0017 0018 0019 0020 0021 0022 0023 0024

0025 0026 0021 0028 0029 0030 0031 0032 0033 0034 0035 0036 0037 0038 0039 0040 0041 0042 0043 0044 0045 0046 0047 0048

0049 0050

Page 33: Bul Gov. For. Exp. Sta. No. 29:3 1977 171~ 21(;月,の〆“e 4E ム 林業試験場研究報告第293号 2. 主成分分析の考え方主成分分析は,理論的にはβ個の観測変量に対してある線型変換を施して企く新しいm個(ρ>m)

募締民主謀邸州回戦藩沌口、中山父、設時(∞)(羽田同断悼む(州知」一一)

0001

PCA04051 PCA04052 PCA040S3 PCA04054 PCA040S5 PCA04056

PAGE 5白URCE PR口GRA併し1ST

ENOFILE DISC2 RETURN WRITE(6dIOJ FORMATIIIIHO ,5X ,23HTRANS CARS 15 INCOHRECTJ RETURN 1 正 ND

FORTRAN 白 K !TAC 4500

STATE問ENT

4

110

LlNE NO.

0051 0052 0053 0054 005'; 0056

INS-

Page 34: Bul Gov. For. Exp. Sta. No. 29:3 1977 171~ 21(;月,の〆“e 4E ム 林業試験場研究報告第293号 2. 主成分分析の考え方主成分分析は,理論的にはβ個の観測変量に対してある線型変換を施して企く新しいm個(ρ>m)

除。己〉事件

図 4. (つづき〉

0008

券株山一円懇話習浅惜同時連N宮市叩

PAGε

12345678901234567890123456789912345678901234567

DODODOGOOi--Ill--1122222222Z2333333333344444444

Aunvnu《υnunυハununu《uvAV内υ《UnυnunununvハU内U《U《ununu《U《U《υ内υnu《υnvnuAUnυnuハvnuAυ《uhunu《unununu《unU

5555555555555555555555-555555555555555555555555

《unu

ハυ

ハUハUAUnu

ハuハU

門unUハυnMAUnuハU《unuハυハ》内υハUハυハunυ

ハ》ハunuハunUハU《ununununvAununuハununuハununuAUnU

A州民“na晶nAMnA両府内倫凪na前向“na侃na内倫鮎na“nAMnAURAMHAMnh

“naanA供局内命的民飴HA品nhunA内A内向maunn-nhunA

内命“nAMna内AMnaMHa内a-na内a舗na内見組naMnRMnaMnaMn

fMFLPしF」rしPLPし併し

PしFLPLP』fiuFしF」

FLnLPL戸eMFLWFLPEWFLPLPLPLPLPLP

ルr」

FLPLPLPLFLPLPUnLPaupLPLFLPRhfνrhpa

、Fιw

pppppppppppppppppppppppphvpppppppppPPD-pppppppppp

5 口URCE PR口 GRAM LlST FORTRA伺口 K lT AC 4500

5TAfEMENT

SUBf~巴UTINE 5UBS20 (:ÚK , L , M 宮~~ 1 C, D )

REAl X(40) (,0 丁目 (1 , 2 , 3~4 ,ぅ, 6 , 7 tSf9 , lO d 1 ,12 , 13t14 t15t1 晶, 17.18 , 19 , 20) , KX( L) =X( 門) 0} X< N )

G� TO 100

z GO TLl 100

X( L) =X(M) 恒 X( N)

X (L ) =X( 問}ホ X( N) G 口 TO 100

X( L) =X( 円 )/X(IO

も口 10 100

6口 T 口 100X( L) =X( 抑 1+ζ

X( Ll =X{ M )*C

GO TO 100 X(LJ=X( 刊 J/C

G口 T 口 100

X( L) =X{MJ 存率C

GO TO 100

X( Ll =X( 門 J'*X( N) GO TO 100

X (Ll =C司 X( 阿〉10 G口 TO 10口

G口 TO 10口

自口 TO 100 X( L) =X( 門)キ己+口

X( L) =X(MJ* 亡+X( N) ホ白

12

13 GO TO 100

GO TO 1自 0

E口 TO 1 口口

16 ・GO TO 100

G口 T白 100

X(L J=SORH X( 何 )+C)

X( し}芯ABS(X( 問 )-+C)

G口 T口 100

GO TO 1 口0

己口 TO 100

X( L) =ζ IX (問}11

X(U怠C キ *X( 判}

X( Ll =AL口 G( X( M 1 令 CJ

X( Ll =AL口巴 lO( X( 門 ).C 1

X(LJ=EXP(X( 門 )+C)11

X( L) =S!GN(X{ 門) ,X(N))

RETURN EcND

14

15

18

19

100

20

。nv M川E

1234RJb789C1234567B90i23456f8QOI2345678901234567

NoonVOOOOODIlli--iliiZ222222222333333333344444444

10000000000000000000009000000DODO000000000000000

LO000000000000000000000000000000ODDOGoonuGO000000

Page 35: Bul Gov. For. Exp. Sta. No. 29:3 1977 171~ 21(;月,の〆“e 4E ム 林業試験場研究報告第293号 2. 主成分分析の考え方主成分分析は,理論的にはβ個の観測変量に対してある線型変換を施して企く新しいm個(ρ>m)

品券株抑止察統副将議u。口、w

んい\山司選時(∞)〈盟国間感輔さ(淵}ニ)iNC印|

0009

123456789012345678901

2345678901234567890l234567890

0090000Doll--Illi--2222222222333333333344444444445

ハV内υnunv《U《uhunUAU内υ《UnU《U内Unu《ununv《U

《υnunununununununu

内unu《VハUnununununvAυハunυnvnυハunUハν内υAUハunυハU

66666666666666666666666666666666666666666666666666

《vnununuハ)ハυハυハunU《UハυnunuハUハU

《ununu内U

ハUハυハυハ》ハunu内unυnu《unU《υ

内υハυハUハUハυnUAunυハununuハUハυAvnvハUAU

ハU

ハυ

ぬ偽府内amARAMn角的食州島内負"a“na鍾nhmaunaHAAHAMnAMnau

門ぬ内命una内食“naan凪MnAMnhunn

晶nAMnA“nAMna“nAMnaωna

臥na“na内A鍋nA内角州

AMnhnAun島内ANAMnaunAHA

RMnaA

CFLEECFLEECECCCCCCCFLEECECCFLEECECCFLEECEFLEECECECCC

戸LEErL

PpppppppPPD-ppppppppppppppppppppppPPD-ppρppppppppppp

PAGE SOURC正 PR口 GRA例1.' 15T

SU8ROUTINE SURS03 (場}

COM問 0 /11 NO , Nl f/'>J 2~ 例O. 祈 l' 阿2. 門 3. lT.IP.1R. IC. IG.ISO.IS COMMON 15 E.対 tHHtX , D , U1Y , AM , AN , A~BINT王 GER !5E(301 REAL X(40hD(40 ).U(40hY(30 ), A例 (30) , AN(30) , AC30 ・ 30) ., 8(30 ,30) DOU8LE PRECIS! 色N H(40)tHH(40 ), H8 DATA DD ,HB/4HSELE ,BH IF( !S.NE.O IG巴 TO 2

1 READ(S , lOOI HED , IS 100 FORMAT(A4 ,6X , 151

lF( HED.NE.DD)(,口 T口 lJF( 15.EQ. 0) (,白 Tn 4 恥 EAD(Sdl0)Ml ,阿 2 , !PdR , IC , IG .c ISE(JI.J=I.30)

FOR 似 AT(215t4I2 , 2X , 3012l

O� 12 J=lo3O ! F ( ! SE ( j 1. E 0.0 1 色白了目 5C口 NT!NUE

N2=_10

GO fO 6 N2=j- 1

lF( N2.GE.IIGO T 口 600 25 j=I ,NO 15E(j)=J

己口 TO 60 JF(N2 ・ GT.30. 口R.IC.GT.NZl GO TO 4 IF(Ml.GT ・阿 2.0R. 阿 2.GT.MOl GO TO 4

lF( !C.GT.8JGO TO 4 jF(Ml.EO.OIMl= 1 IF( 料 2.EQ.0) 阿 2=MO

150= ISO+l DO 10 J=I ,N2 K=ISE(J) JF 【 K.GT.40"OR* 村 {K}"EQ@H日}己 o 10 4 HH(JI=H(Kl

WR !TE(6 0120) 150 120 FORMAT( 1 Hl ,35H.. TRANSF口R例Eü ANDI口R S 1"し1" CTED DATA.

128H ・・ (5ELE CTI 口 N CARD NO.= , I4ilH}lll WRJTE(6 ,130)(HH(j) , J=1.N21 FORMAT({lH d3X ,7(8X ,A8))) D020!=I.N2

Y( J) = O. AN( 1)= 1.0E30 A州( I 1=-1 .O E3口

D� 20 J=I ,N2 A(J.JI=O. REWJND DJSC2

何 3 =0

D口 50 K=hMO

F 白内 TRANロ KITAC 4500

STATE 門ENT

2

110

60

130

10

12

25

20

6

lINE 制口.

0001

0002 0003 0004 00白う

0006

0007 0008 0009 001η

001J 0012 0013 0014

0015

0016 0017 0018

0019 0白 200021

0022

0023 0024

0025 0臼 2~

0027

0028

0029 0030 0031

0032 0033

0034 0035 0036 0037

0038 0039 0040 0041 0042

0043 0044

0045 0046

0047 0048

0049 005。

Page 36: Bul Gov. For. Exp. Sta. No. 29:3 1977 171~ 21(;月,の〆“e 4E ム 林業試験場研究報告第293号 2. 主成分分析の考え方主成分分析は,理論的にはβ個の観測変量に対してある線型変換を施して企く新しいm個(ρ>m)

t-;)

52 図 4. (つづき〕

立山総見簿一邸出山市河山叩雨時

PAGE 0010

PCA06051

PCA06052 PCA06053 PCA06054 PCA06055 PCA06056

PCA06057 PCA06臼 58PCAO�059

PCA06060

PCA06061 PCA06062 PCAO�063 PCA06064

PCA06065 PCA06066 PCA06067 ρCA0606日

PCA06069 PCA06凸 70

PCA06071

PCA06072 PCA06073 PCA06074

PCA06075

5自URCξPR口GRAM LJ5T

F

T也

sph

F仁恥

R

d

u

n

-

R

6

Tム「匂

r

r

M

m

吋''VA'

Vハ

CJ

3

I

l

nu,、

J

内υ

kdnυan

同川CJA

門u、,,ia-AF

'l

フh,、,

N

J

'

J

A

U

T

E

3

aV5

C

-irU司ム@、,rF』、,、,、,rah

N

}

S

6

0

T

}

}

J

L

9?ι1JMtaF

74TAT命,也、rL

IMn--FIR

<tX5

2

}79

XX

H

JT3J,6},,}4 ,ru{X{O}}1

、,@伺υvAマJCL'τAγ4f、事

4dkn78i,、9Tsnwrhr、VAV

(.,oIE}M

いNムマ5

いupn内UK1ム

D内a-iAmA丹、,、d'

r、内υ

ヲιru、dyhaM側、,?』,、,、,、ア』iJWIBM、,nV

7

N309}5}NX141N''DH

2

1

'

i

4

8

2

p

+

X

N

a

'

1

1

6

1

TC

I)

IH.,unl}AII{

1/

N5

=}Iae

・g-EklzIM

町内三

AAE9fl

cLTAγ111Ml

》【

ptFbfpLU,、r、マ4,、A門A円tuZ2HV'or、

MU内UJ』+,e‘Hυ

f、γlenυT'V92F句、,、,

N川削川,‘γEM内

E(

3OE--PEAZ

ADZ}}Oid--REAR

TDKM35}ITM

NM

似し4}114,,TUT

AA(=IJ

<

IR{{RI{(1JNTIRTD

TEF302tFR

FF

口口【

MUM河口{{口

ER

EN

5RiMHD-xiuwFIIFDYAADAACRWFRE

ハVAU内υ

4

5

6

I

1

3

4

1

F口 RTRAN。 KITAC 4500

30

40 50

lINE NO.

0051 0052 0053 0054

0055 0056

0057

0058 0059 0060

0061 0062 0063 0064

0065 0066 0067

0068 0064 口070

0071

0072 0073 0074 0臼 75

議NUω中

Page 37: Bul Gov. For. Exp. Sta. No. 29:3 1977 171~ 21(;月,の〆“e 4E ム 林業試験場研究報告第293号 2. 主成分分析の考え方主成分分析は,理論的にはβ個の観測変量に対してある線型変換を施して企く新しいm個(ρ>m)

菜料訟で添節副総一蒸UHO口、ゅんい\hw,,惑か(由)ハヨ同絞芝〉(円以」一一)

c ..._,

0011

-2345678901234567

901Z345618901234567890l

《434567890

000000cool--llil--lZ222222222333333333344444444445

内UAVnυ《vnυnuハununUハUAUハUハun》ハUV《U《UハununUAUAUnunununvAUハunUハununυnUAUハununuハunvハUnu川vnunνnununu《υ内U《U

777777777777777777777777777777777777777777777777?7

nUハununv《υハUハUAUハυnυハunUハunuハUWAυ内υハvnunu《U

ハuハυnuハunuhunvnununUハυハVnununununνnu《υハυnv門υ

ハνAUハU

ハunuハunυ

amaHhunAMnAm

瓜内角川島内角姐nA内角“na内典的A“n-AAuna円a“naハ泉州

amnMUA働na曲nA“nAunAM門倫内奥una“nAnAunAHaun奥州九州AmA州h“nAωnAHH《晶na“nA州A

内AMnA州民同

AU門AM門

ccccccccccFLEECEFLEEC

EDEEP

CCCCCCFLEEFLEFLEECECcccccccc

pppppppppppppppppppppppppppppppppppppppppppppppppp

su且ROUTl河 E 5U日 504 (キ}

COM問 ON NO.Nl , N2 ,向。,附 1. 伺2 ,問 3.IT , IP , IR , !C.IG.150 舎 ISC白河阿白N I5E9H , HH , X , D , U~Y , A 何 1ANtA , S

lNTEGER ISE( 30) REAL X{40J.D(401.U(40l , Vi30l , AM(30l , AN(30) , AC30 ・ 30 ,8(30 ,30) DaU9LE PRECISION H(40) ,HH{4Ql 。口 10 1 = 1 曹 N2

B四 10 j=1~N2 A(! , J)=(A(!.j) 倒y{ n'V( J lfFLOAH 何3) )!( Fしむ AT( 拘3 時1l)口百 20 1=1.N2 D(! J=SQRT< A( r.!)) D口 30 1=1.N2 DO 30 j=l 宥 N2呂 (1 曹 J)=A( hj) 〆〈口(l}宇D( j l) B( 1 d )=1. WR!TE(6.IOO)

100 FOR阿 AT(/ ノノ lHO , 30H.." 呂 A5IC STATl STlC5 白 N DAT゙. •• 1/IIHO , 25Xt4HMEAN~8XtBHVA火 IANCE.12X , 4HS.D. , 9X ,

Z7H3 事 SIGMA , 12x , 4HCQVe , 12X , 4HMAX~~12X~4H 何 1 N. l 由白 40 I::::ltN2 AI=Y( 1 )/FL百 AT( 阿 3 ) A2=3 歯*口(1)

A3=O( 1 )ノ AlWRITE{6 ,110

110 F口RMAT<!H

AM(! )=Al AN (1 )=D( 1 DO 50 1=1 lF ( IR. 正由 .O.AND. 日 (1 )冒し T.5.0E 幽 5 )品目 丁目 3CONTlNUE WRITE(6.J201 FOR川 AT{IIIIHO , 32H*~ VARIANCE C口 VAR l.ANCE 問 ATRIXNL=O NS=NL+l NL=MINO(NS+6 ,N2] WR!TE( 6 ,130) (HH( J) ,J=NS ,NL) F口 RMAT(/1HO , 13X , 7(8X , ASJl

D 口 60 1=I ,N2 I<RITE(6.JIO IF(NL.LT.N2 WR!T 王 (6 ,1 40 ) FORMAT I! I ノ lHO , Z4H.. CORRELATICN MATRIX ..1 NL=O NS コ Nl+l

NLニ州 INO( NSφ6 ,:'11 2)

WRJTE( 6 , 130 J( HH( J" J= ,~S , NL) む口 70 I=1 ,N2 WR ITξ(6 , 110)HH(!) ,(S(! , J) , J=NS , NU

PAGE 50URCE PR岳G再AM Ll5T FO詞TRANOK!T AC 4500

5TATE伺E外τ

10

20

3 む

1 ) ~ Al , A ( 1 J r ), D( 1 ), A2 , A3 , A 拘( 1 hAN{ 1) .7FI6.5J

r ) t ( A ( r , j ) 9 j::: N 5 , Nl ) 1

c

c

40

12つ

14�

70

130

50

60

@

向υ

創刊

EIZ34567890i2545678901Z3456789

12345678991234567890

NOODOOOOGo--illIll--22

内42222222333333333344444444445

7

みの

ununununuvnvハUAυハ》ハU八υハυnunvハununuハU

ハunvnununvnuハUハunuAvnvnuハunununuハυnuハunvnvnunυnνハMハUハυnvnuハununu

L00000090

OD0000050000000000000000000G909000DODDODO

Page 38: Bul Gov. For. Exp. Sta. No. 29:3 1977 171~ 21(;月,の〆“e 4E ム 林業試験場研究報告第293号 2. 主成分分析の考え方主成分分析は,理論的にはβ個の観測変量に対してある線型変換を施して企く新しいm個(ρ>m)

1Mgl

図 4. (つづき〉

PAGE 0012

PCA07051 PCA01052 PCA01053 PCA01054 PCA0705S PCA07056

PCA07051 PCA01058 PCA07059 PCA07060 PCA07061

S口 URCE PR 口6RA阿L1 5T

5TATE 例ENT

lF( NL.LT .N2)6口 TO 2 RETURN

3 WRITE(6.150) 150 FORMAT{IIIIHO.5X , 23HS ・ 0.:0 THl5 5正 LECT CARD.

1174 15 UNABlE TO RUNl IF( IG.EO.OlG口 TO 90 口o 80 !=ldG READ(5 ,160)!I.!2.!3 FORMAT(5X ,315) RETURI場 l

END

FORTRAN 白 K lT AC 4500

80 160

90

LlNE NO.

0051 0052 口0530054 0055 0056 0057 口口 580059 0060 0061

廿す糠訟伺蹄峠叫出羽凶難時

選出ω叩

Page 39: Bul Gov. For. Exp. Sta. No. 29:3 1977 171~ 21(;月,の〆“e 4E ム 林業試験場研究報告第293号 2. 主成分分析の考え方主成分分析は,理論的にはβ個の観測変量に対してある線型変換を施して企く新しいm個(ρ>m)

詳料叫比四冊焼酎戦議UH。口、,W川に山可議時(∞)(出早川市河湾笠)(一円(」一一)

11NO由li

0013

-刷駒唱酪置幡町四帽-畑噂勘町・暢暢 PCAoaOOl

SUBR百UTINE SU゚ S05 PCA08002 PCAO 8003

E口阿例目N NO , Nl 'l N29MOt 判 1 曹洞 2. 阿 3dT.JPdR.lC 智 !GdSO.1S PCA08004 E自 M 阿古内 ISEtH , HH , XtDtU , Y,AM9ANtA ,S PCA08005 INTEGER ISE(30) PCA08006 REAL X(40).D(40) , U(4 口), Y( 3 口 hÄ阿 (30hAN (3 0 , A (3 0 , 30l ,日 (30 , 30) PCAOB007 DOU8LE PRECISI 口N H(40) ,HH{40l PCAOBOOB lF( IR.E 臼 .1 )60τo I PCAOB009 D凸 10 lol.N2 PCA司自由 1000 10 J=1.N2 PCA08011 A(!.J)=B(! ,J) PCA白 8012CALL SUBS30( A. 再 .N2) PCAO 毘 013WRITE(&.1001 PCA08D14

100 FOR 阿 AT<ノ IIIHO , 44H ,. ê R王 SULT OF PRINCIPAL CDMP自 NENT ANALYSIS ..1 PCAQ8015 IF( JR.EO. 11G口 TO 2 PCA08訂 16

WRITE(�.llOj PCA08 口n!10 FORMAT(/IHü ・ 21X.39H(P.C.A. STARTS FROM C 自白RELA Tl口 N μATRIX)/) PCAOBOl 昌

GO T百 3 PCAD8019 2 IIRlTE(6d20) PCAOS02D ! 20 FOR例 AH ノ lHO , 21X , 39H(PtC 禽 A・ S T Aiミ TS F 罰百 M VARIANCE CDVARJANCE. PCAOBD21

18H MAτRIXI!) PCAOB02Z IF(A(N2fN2)~GT.()~)GC TO 17 PCA古田 0230口 16 1=!.N2 PCAOB024 AI!.!I=Alj , !)-AIN2 ,NZI PCA08025 CQNT! NUE PCAO 8026 TR エ O. PCA08027 口口 201=1 ,N2 PCA0802B Y(J)= AIl ,]) PCA08029 TR=TR+A( 1, I) PCA08030 0口 30 ! = 1. N2 PCA08031 U(] )=AI 1,] )!TR PCA08032 DI] 1 ヱ UIll PCAOB033 Da 40 I=2 ,N2 PCA08034 DIII ニ υ(1 )-t-臼( 1ω1) PCA08035 NL=O PCAOB036

NS=NL+j PCA08037 NL= 門 INO(N5-6.N2) PCAOB038 W司 ITE(6.130)(J , J=NS , NLI PCA08039

130 FORMAT(/IHO.13X ・ 7 (1)(, 13 号色対 CClM゙ .ll PζA08040 WRJT[(6.140)IYIJI.!=N5 ,NLl PCA08041

140 FORMAT(lH .5X.llHEIGEN VALUE/!IH 写 13X , ì Fl 6 ⑤ 5 1 PCA08042 WR! TEI 6 ,150)( U( 1 h l=NS ,NL) PCA08043

150 FORMATIIH .5X , 18HCONT科lBUTjON RIli!O!(!H d3X,7 Fl 6.5 ) PCA08044 レ~RITE(6 , 16 口)(百( ! ,,!零NS , NLI PCAOB045

160 F 白内 MATIIH .5X.29HCUMULATIVE C吉見lIUHUT! 巴 N RA Tl訂IUH 雪 13X , 7 Fl 6.5)) PCAOB046 lFINL ・ LT.N216母 T自 4 PCA08047 WRITE(6 , 17 口) PCAOB04B F 司 R~AT(IIIIHO , 5X , 1 嘉ド'.E IGEN VËCTO時 PCAOB049NL=O PCA08050

PAGE 5自U日 CE P読むGRAM LlST Fü悶 TRAN自K lT Aと 4500

5TAr E附 ENT

10 l

170

16 17

20

30

40

E

L1NE NO.

0001 0002 0003 口 0040005 0006 0007 0008 0009 口0100011 0012 0013 0014 0015 0016 0017 0018 0019 0020 0021 0022

0023 0つ 24

0025 0 凸 26

0027

002 自

0029 003口

0031 0032 0033 0034 0035 0036

0037 0038 0039

0040 0041 0042 0043 0044 0045

0046 0047 0048 0049 0050

Page 40: Bul Gov. For. Exp. Sta. No. 29:3 1977 171~ 21(;月,の〆“e 4E ム 林業試験場研究報告第293号 2. 主成分分析の考え方主成分分析は,理論的にはβ個の観測変量に対してある線型変換を施して企く新しいm個(ρ>m)

ト"Zコ

函 4. (つづき〉

FロRTRAN S白URCE PROGRAM LI5T PAGE 0014

PCA08051 PCA08052 PCA08053 PCAOB054 PCA08055 PCA08056 PCA08057 PCA08058 PCAOB059 PCAOB060

。KITAC 4500

STATE 何ENT

N5=NL+j NL=阿 INO(NSφ6 , N21WRITEI6.i30)(J.J.N5,NLl DO 5自 1= 1. N2 WR JTEI6 0180 )HH( n.I.B( hJ) ,J=N5.NLl

180 F口R 問AT(IH .5X.AB ,lFI6.S1 50 CONTlNUE

IF(NL.LT.N2)GO TO 5 RETURN END

llNE NCIo

0051 0052 0053 0054 0055 0056 0051 0058 0059 0060

業料開討察部酎期尚諮命瀦NSdl

Page 41: Bul Gov. For. Exp. Sta. No. 29:3 1977 171~ 21(;月,の〆“e 4E ム 林業試験場研究報告第293号 2. 主成分分析の考え方主成分分析は,理論的にはβ個の観測変量に対してある線型変換を施して企く新しいm個(ρ>m)

説明問料怒結成鴻一議パ

UO司、叫んい\札可。設命〈∞〉(唱七時均株

4)(況」ニ)

liMHMl|

0015

PCA09001

PCA09002 PCA09003 PCA09004 PCA09�05 PCA09006

PCA09007 P亡 A09008PCAO 9009

PCA09010

PCA09011 PCA09012

PCA09013 PCA09014 PCA09015 PCAO 9016

PCA09017 PCA09018

PCA09019 ρCA0 9020

PCA09021 PCA09022

PCA09023 PCAO 9024

PCAO 9025

PCA09�26 PCA09む 21PCA09028

PCA09029 PCA09030

PCA09031 PCA09032 PCA09033

PCA09034

PCA09035 PCA09036

PCA09037 PCA09038 PCA09039

PCA09040

PCA09041 PCAO 9042

PCA09043 PCA09044 PCA09045

PCA09046

PCA09047

PCA0904B PCA09049

PCA09050

PAGξ 5口URCE PR百GRA拘Ll 5T

5U自ROUTINE SUBS30 (A.V ,Nl

REAL A(30 , 30) , Ví30 , 30).BIP(30) ,自 IQ130 l

Nl=N-l 。純益 N

00 1 ! = 1 • N 00 2 J=l. 拘V( I.J)=O.

V(!o!)=l.

VNOR 拘 =0.00 10 l=l.Nl

Jξ=1+1

自由 10 J=JEoN VN口RM=VNOR 何φ2.キ A( j, j) 事 A (J .j)

JF(VNOR何 .LE.O.1G口 TO 79

VNORM=5臼R T( VN口尺州》FN口RM=VNOR何本 5.0E齢5

15WICH=O VNORM=VN口RM/DNB白 30 IQ ロ 2 , N

IPEND=j 白-1

DO 40 JP=I , IPEND !F� ABS( A( IP , IO)) .LT.VNO舟州 )60 .TO 40

15WICH=1 DLAMDA=-A(IP , IQ) DMU=.5 牢 (A(!P ,J P) 勘 A( 10 , 1 日))DO州ξGA=5IGN(I. , DMU)'日 LAM白A/ SQRT( DMU 辛口MU+DLAMDA.帯。LAMDA)

DW51=1. 四日口問 EGA*D 口 MEGAIF( DW51. Lt..0.)G 口 τロ 42

DWSl=l.+S白RT <D W51)

GO T口 44口凶 52= 1.

5IND=OOMEGA/50RTí2.*0 凶 S2 )

caSD=I. ・5IND*SJNDJFIC口50.6T.0.)C05口 =50トミ TíCOSDl

自白 50 J=!.N IFI I.EO.IP. 自 R.I.EQ .j白} 凸口 T 口 52 BJP(J)=A(J ,j Pl ホC自 50 ・ AI I ,J O )*5 IND 白 I 臼 IIJ=A< JdPl 物SINOφAI 1. 10) 本C口与DC口 NTJNUEVIiP =VIlolP)'COSû幽 Ví r.IC)*5INO Vlj.IO)=V� hIPJ*SINDφV II ,J Q) 本C口 5DV( I , !P )=Vl IP COSDD =c口SD 本ζ05D

SIl'll Dリ =SIND本 5lND

SINC05=SIND淑 ζ口 5 1)ßIμ( !P ) =A � 1 p.r P ) .COSOD.A 1I Q, 1 臼}ネSIND口司-2.*A( jP , 10) 牢 SINcas

BI 口 I IO)=A( !PdP) 傘SINDD 今 A( I 口,]口}牟COSDD+2.*A( IP , IO)'5INC口 551ρI!Q)=IAIIP , IP)-A(!O , IO)) キ 51NC口 S+A{IP , !Q) 牢 ICDSDD-5INDDl

FORTRAN OK lT AC 4500

STAγ定例ENT

c

Z I

10

50

20 25

42 44

52

C

L!伺 t NO.

0001 0002 0003 0004

0005 0006 0007

0008 0009 001� 0臼 II

0012

0013 0014 0015 0016 0017 0018

0019 0020 0021 0022

0023 0024 0025

0026 0027

0028

0029 0030

0031

0032 0033 0034 0035 0036 0037 0038

0039

004日

自由41

0042 (l043 0044 0045

0 巴46

00.41 0048

0049 005 。

Page 42: Bul Gov. For. Exp. Sta. No. 29:3 1977 171~ 21(;月,の〆“e 4E ム 林業試験場研究報告第293号 2. 主成分分析の考え方主成分分析は,理論的にはβ個の観測変量に対してある線型変換を施して企く新しいm個(ρ>m)

r、2

E、2

図 4. (つづき)

#骨川神叫舛嗣冊邸、問中日業時

同明日甲山岨

0016

PCA09051

PCA09052 PCA09053 PCA09054 PCA09055 P0A09056 PCAQ9057 PCA09058 PCA09059 PCA09060 PCA09061 PCA09062 PCA09063 PCAOQ064 PCA09065 PCA09066 PCA09067 PCA09068 PCA090�9 PCA09070 PC'A09071 PCA09072 PCA09073 PCA09074 PCA09075 PCA09076 PCA09077 PCA09078 PCA09019 PCA090B。

PAGE 5自 URCE PROGRAM L!5T

0

2

吋ζ

。o

向U

向U

マ'TE'

門UρURU

T

G

G

}

、,,胃

d、,

.d

向UM

門。。、

JVA

守A

G

R

J

A

5

3

})口口,

Mmum-d

Qhe})}INγlde'

XI

V-'ivaTA

,a@

戸』r

,、

1JZAm

《d

tftttE

15AX}muw

pNPGNPaNTNNG-

]

3AJN]'D

T'II

11

G'')EJIMFX'IJZ

TBIBaulBEN--1Iru''{Ala

,{)

NE

司三

2Z2--zEEC05MH4-=

・・

JEIAMH

luvJE

Et'1311}}UUI=ZRIJGX{U

EF'Jr

=xu

M円

nypsnU7aTANNUWMn

門UF

AハA,

JM内

AM円、JSJuv

、a,

AMHNM同

F』

7ふ

nuyι-Y-令'h''τaAyaFbfu

内νN川ハν2awMH=27&=TAVA4Z7

MHTAQH

T{rb

,96POTTIITV88JAXJTi

,A8l

,,TU

A門口

ITA--NN{

{{{AVAN51MHSJJNTU

T-ot(

口{{自由

F5

FDDFFmADuw((

uwf

〈百

EN

SBDAADAAEC--G10

IIAMnCDAADDVVCRE

FORTRAN 口KITAC 4500

60

70 ・79

84 80 999

ヲιハU内υ

643

83

82

LlNE N 自@

0051 0052 00 ち3

0054 0055 0056 0057 0058 臼 059

0060 0061 0062 0063 0064 0065

0066 0067 0068 0069 0070 0071 0072 0073 0074 0075 0076 口077

0078 0079 0080

Page 43: Bul Gov. For. Exp. Sta. No. 29:3 1977 171~ 21(;月,の〆“e 4E ム 林業試験場研究報告第293号 2. 主成分分析の考え方主成分分析は,理論的にはβ個の観測変量に対してある線型変換を施して企く新しいm個(ρ>m)

++日附刊

H円岬蹄蕊同組一一市開パ可。口、,

JW川て、起昨〈∞)(~半減邸玄)(~吉一一)t、コト~

乙心

PAGE

PCAIOOOl

PCAI0002 PCAI0003 PCAI0004

PCAl 由 005PCAI0006

PCAI0007

PCAIOOOB PCAI0009

PCAIOOIO PCAIOOll PCAI0012 PCAI0013

PCAI0014 PCA10015 PCAI0016 PCAI0017 PCA100l8

PCAI0019 PCAI0020 PCAI0021 PCA10022

PCA10023 PCAI0024

PCAI0025 PCAI0026 PCAI0027

PCAI0028

PCAI0029 PCAI003由

PCAI0031 PCAI0032 PCAI0033 PCAI0034 PCAI0035 PCAl 自036

PCAI0031 PCA10038 PCAI0039 PCA10040

PCA10041 PCA10口 42PCA10043

PCAI0044 PCA10045 ρCAI0046

PCAI0047 PCAI0048

PCAI0049 PCA1D050

5UBROUTINE SU日 506 (半,本}

E口問問白内 側口.fn , N2 , 問。,向 1 .拘 2. 問 3 ,!T dPdR 曹 !C を IG. 150 , 15 (0月刊白 N ISE , H , HH , X , O , u , y , A阿, AN , A 'l ß

INTEGER ISE( 30) REAl X(40) , DI40) , U(40J , Y(30) , A科 (30 )・ AN(3司 )~A(301 五0)'8(30 , 30)

DOUBlE PRECISJ 口 N H(40J.HH(40) WRITE(6.J05 )

105 FOH 開 AT 【 11 ノ 1 HO ・ 5X.14HFAC丁目R L白 AD l NG)

対 L=ONS=Nl+l Nl= 何 I 制O!N5+6 , N2lWRJTE(6 , lOO)(j , j=NS , NL)

100 FORMAT(/1HO , 13X , 1(1X , j3 , 6H CO例 P. ))

00 20 1=l>N2 U( 1 ) =0. S!!=AN(!)

!FI IR.E日 .0)511=1.D口 10 j=I.NZ AJJ コ Y( J J

JF(AJJ.LT.O.lAJJ=O.

D(J)=50R T( AJJl'白(j, JJ/SIIA(!.JJ=OIJJ 意DI J)

U( 1 )=U(! ).ト( 1 ,J) WRITE(6 ,1l 0JHH( 1 J ・{臼 (J).J=NS.NLF自R 例 AT<lH , SX , AAt7F16@ 宮》

IF(NL.lhN2JGロ TO 1 WRITEI6.200J

200 F口RMATIIIIIHO.5X.2BドE口 NTRIBUTION OF C 口内 PONENT T白 B

llBH 口R!GINÄL VARIABLEI DO 21 1=I.N2 白白 21 J=t.N2 A( 1 , J > =A < I.J J/U( J)

NL=O N5=NLφL

NL= 阿 INOCN5+6 , N2)

WRITE(6.100)(J.J=NS.NLI 口百 13 1=I ,N2 凶 RITE(6dlO)HH( 1 ),( A( I , j) , j=N5 , NL) IF(Nl.LT.N2) 己百丁目とWR lTE< 6 , 210)

210 F口R 門 AT(//IIHO , 5X , 26HCU 阿υLATIVE CONTRIBUTION OF.

131H COMPONENT T口百R 1 凸 INAL VAR!A邑 L ElNL=O NS=NLφl

NL=MINO(NS+6 ,N2) WRITE(6 , lOO)(J.J=N5.NL) DO 24 I=19N2 υ(1 I=A( 1 d)

。011SOURCE PROGRAM UST F 口RTRANOKlTAC 4500

STATEMENT

10

20 110

13

21

3

C

C

LINE 陶口.

0001

0002 0003 0004

0005 0006

0007

0008 0009

0010 0011 0012 0013 0014 0015 0016 0011 0018 0019 0020 0021

0022 0023 0024

0025

0026 0027

0028 0029 0030

0031 0032 0033

0034 0035 0036 0037

0038 0039 0040 0041 0042

0043

0044 0045 OQ4b

0047 0048 0049

0050

Page 44: Bul Gov. For. Exp. Sta. No. 29:3 1977 171~ 21(;月,の〆“e 4E ム 林業試験場研究報告第293号 2. 主成分分析の考え方主成分分析は,理論的にはβ個の観測変量に対してある線型変換を施して企く新しいm個(ρ>m)

"" トーゐ長h図 4. (つづき)

議仙川民間時節~ヨ凶噺間滞

目前凶器叩

。 OIB

ρCAI0051

PCAI0052 PCA10053 PCAI0054 PCAI0055 PCAI0056 PCトI0057 PCAI0058 PCAI0059 PCAI0060 PCAIOQ61 PCAIOO�2 PCAI0063 PCAI0064 PCAI0065 PCAI0066 PむAI0067PCAI0068 PCAI0069 PCAI0070 PCAI0071 PCAI0072 PCAI0073 PCAI0074 PCAI0075 PCAI む 076PCAI0077 PCAI0078 PむAI0079PCAI0080 PCAI0081 PCAI0082 PCAI0083 PCAI0084 PCAI0085 PCAI0086 PCAI0087 PCAI0088 PCAI0089 PCA1009口PCAI0091 PCA10092 PCAI0093 PCAI0094 PCAI0095 PCAI0096 PCAI0097 PCAI0098

口口 23 J=2d'a UlJl=Al j ,J lφU l J-I 1 WRjTEl6.II01HHl 11 ,( U(Jl.J=NS ・ N LlIF(NLcLT.N21GO T 口 3lF( Ic.EQ.01G口丁目 45

WRlfc(6ol15l 115 F白RMAτ{ ノノIl HO.5x.27HCO 例PQNENT SC口RE COEFFICIENT)

00 11 J=I.IC lFIYIJl.GT.5.0E-41GO T 口 11lC=J- 1 GO TO 12 CONTINUE WRITC(6.1001(J ,J=I , ICl 00 35 1=I.N2 00 25 J=I ,J C 日 (1 , J)= 日 II.J)/SQRTlY(Jll

凶 RITEl6tl101HHlIld 白 (loJl".!=I ,J CwR ITE l 6.120 )

120 FORMATlII11HO.21H.. COMPONENT SCORE ..1 WRITEl6.100)(J ,J =1.ICl RE凶 IND Dl5C2 REWIND OI5C3 DO 40 K:::19 例。HEAD lD I5C21(O(J) ,J=1.Nl) IF( K.L T ・阿 I.OR.K.GT.M2)Gむ T 由 40 D口 3 の J:::l , NZ

1= 15El ,) 1 XlJ l=DlIl jF( IR.EO.I lGO TO 30 X(J l=l X( J)-AM( J 1 l/ANl J) ylJl=O. 00 50 J=ldC 00 50 I=1.N2 YlJl=Y(Jl+X(!l 本 tl( l ,j) IF( IG.EO.01GO TO 55 WRITE(OI5C31(YlJl.J=1.!C) WRITE(6t130 >K ,( Y(J) ,J=l tIC) F口R 門 ATlIH , II0.3x.7F16.5) IF( 門口口 (K , Z5).EQ.O)悼RiTElóol40)

F口R 阿 AT< IH ) C白 NT!NUEENDFILE DISC3 E白NTINUElT=4 IFC IG.GT. 臼 )CALL SU8S40 IF( ISO.餡.IS) RETU問N 1 RETURN 2 巨ND

PAGE 50URCE PR口 GRA 例 LI5TFQ終 TRAN口 KITAC 4500

STA rE例 cNT

23 24

11 12

25 35

Ll NE NQ ・

30

5S 130

45

140 40

50

0051 0052 0053 0054 0055 0056 0 日 57

0058 0059 0060 0061 0062 0063 0064 0065 0066 0067 0068 0069 0070 0071 0072 0073 0074 0075 0076 0077 0078 0079 0080 0081 0082 自 083

0084 0085 0086 0087 0088 0089 0090 0091 0092 G093 0094 0095 0096 0097 0098

Page 45: Bul Gov. For. Exp. Sta. No. 29:3 1977 171~ 21(;月,の〆“e 4E ム 林業試験場研究報告第293号 2. 主成分分析の考え方主成分分析は,理論的にはβ個の観測変量に対してある線型変換を施して企く新しいm個(ρ>m)

亨治uh削ゆ」出減努惑い

wo白川叫んい\凶Y

JJ設な〈∞)ハ草山閃成金)(一河川」ニ)

N511

0019

-23456789012345678901234567890123456748901234567890

08000OBool--lllil--2222222222333333333344444444445

000006000

GOOBCGOODODDDODD000

00000000000000008060

11111111111111lllIll-1111111111Illi--lIIIll--llliI

lli--Ill--lIlli--Ill--11lIll-lil--illili--liIlli-ュ

A内auna内角川

AmAωnAmRMnAURamaunA

邑naAAωmAunaHAunaunAMnamA

且門前例AU【aunaMnaunA

柿良ωnAHAHAM門AU門会例AM門病院AMnA“nAM門A門AωnA“円九州氏州内向内“nAMnhmAHA同

AH

FLEECccccccccccFLEECECCCCCCCCFLEECEccccccccccccFLEECcc

ppppppoppppppppppppppppppppppppppppppppppppppppPPD

PAGE 5口URCE PROGRAM LIST

SU臼 R 口U Tl制 E SU日 540

CO 阿例 ON NO ,Nl ~N2 ,向。,問 10 何2~M3dTdP t1 R ,I C 宮 IG ,J SOdSCOMMON ISE9HtHH , X , D , U , Y~A 門 9ANtA , 8INTfGER IS E< 301 ,DO T< SI ,91) ,5!(131 REAL X{40lo0(40 J. U(40) , V <3 0hA問 (30 ), AN(30) 予 A(30dOI 手自 (30 , 30 ) REAL TT(5! l>YY{91 1 口口 UBLE PRECl510N H(401 宮村村( 4日}DATA 5I/IH , lHl~lH2 宮 l ト-1 3tlH4dHち dH6tlH7 , lH8 容 lHq , lHXfIH:~lH-11 M:.:O

5 READ(5 , !10Ijl.j2.13 110 FORMAT(SX ,315)

1F( 向 INO(II , I21.LE.0. 四R ・ IC. 正 O.O)GO T 口 901 F ( M AX 0 ( 11 , J 2 ) • G T • 1 c)巴 o TO 90 (00.

Coι >1.

JF( 13.GE.l )C=13 TT(!1=2.5 キ仁

YY( 1 l=-TT( 1 1 DO 20 j 0 2 ,51 TT(j)oTT(J 曲 1)-O. 1* C

D凸.3 0 j=2 ,91 YY(Jl=yy(j-! )ゃ 0..5*仁 19.D口 10 1 = 1 ,51 口口 10 j=I.91 D白T(j, j)=O

K3=0 REWJND DISC3 D� 100 K=l ,何 OIF(K.LT.Ml.OR.K.GT.MZIG口 T 也 100REAQ<D ISC})(U( j I ,J=l , !c) Ul=u(!11 U2=U(J2) M=1 lF(Ll.lT.TT(251) 刊 =25

00 40 1 =問, 50

JF(Ul.LT.fT(]).AND.Ul.GE 旬rT(]令 1) )GO H1 1 C口 NTINUE

60 T 白 25

例= 1 1F(U2.GE.YY(45)) 円 =45

D口 50 J=M ,90 IF (L'2 ・ GE.YY(jl.AND・ U2 ・ LT ., YYCJ+l) )(J 口 TO Z C 口N 了 INUE

K3=K3+1 1F( J3.EO.O.AND.K3.GT.O IGO T(1 3 G自 TO 100 DOT( j, jl=DO T( 1 ,j )'1 t口 N Tl NUE

FORTRAN 。 κ lTAC 4500

STATE何 ENT

C

20

Z 10il

30

10

40

50 25

C

-内υ

創刊E

123456789012345b789012345678901234567S901234567890

NGOOOOOOGO--1l1111112222222222333333333344444444445

7A

内Uハυ

《U

ハU《υハU

門》内υ内υAUAV《》ハV

ハUハ》ハνハununUハUAυハununuハU内unリハリハununuハUハυハυnunυハUAυ

凡Uハυハυハununu

ハυ

ハυ

ハリハυ

ハυ

ハu

?」

nu《U《υ凡UハUハVハUnuハυハVハU内unυnυ内υnunυnvハunuハvnununυハUハUハυ

ハUAυハUハU

ハvnuハυAUnuハU

ハυハυハUハunυハυ門》ハυnuハリハunvハU

Page 46: Bul Gov. For. Exp. Sta. No. 29:3 1977 171~ 21(;月,の〆“e 4E ム 林業試験場研究報告第293号 2. 主成分分析の考え方主成分分析は,理論的にはβ個の観測変量に対してある線型変換を施して企く新しいm個(ρ>m)

'" σコ

図 4. (つづき)

設株mw付添邸嵩明開設吟地州N甲山中

PA凸E 0020

PCAII051 PCAII052 PCAII053 PCAl1054 PCAII055 PCA! 1056 PCAl1057 PCAl105B PCAII059 PCAII060 PCAII061 PCAII062 PCAl1063 PCAII064 PCAII065 PCAII06� PCAII067 PCAl106B PCAII069 PCAII010 Pむ AllOllPCAII072 PCAII073 PCAII074 PCAII015 PCAII076 PCAII077

ae aw M円

A円

n門rnv

a

内-

}

3

nu

今Fh

3 .

.

,

onvA『,E

nu

戸tF

白よト

MHF』

TT

,zw

T

T

I

'

D

A偶n'uu目、,t』

C

}

1

1

F

319

画}

GFVA幽@《J

Fぃ,,A

nF7・

R

l

'

E

M

E

Dzl,DD FLldaM円働FL

qd

l

ZF

、目pnヲe・74H円内U

2

3

T

J

'

H

'

'

b

ilN'X9},3R 3

Z

E

-

-

4

1

4

K

E

J

E-ht』州内r、.Env'4Y&、,au

胃帽

、,、,内UTiy&E

''nvMnr‘

l

l

p

m

u

'

1

z

z

u

a

E

--MHnuTιaz・nd

削川Z

09

{H骨J99UH5

EEC'3

'F61MH

@-、ae,M

Fι、、s'ag、ヲ'ham門u

l」!L.,i

弓4

弓ザ事

3

》1・広L-wpnT'

I・e・}{e'5《'TaVAG

1-nυnu屯g

可』品川MH・V-zivV《@VZ

》A守AmUAMnmuv

ZM

川M刊BL74守gnr,hE1・Rdnvg、唱'bpnヲ,yapb

T&2

τAqLAMnA内,、、s'aMH'A

dc・、,,ae'a骨+、,1eU

'nu、,

RdadAV、eea-v‘nu,内URJnv'nV9PL'目働、,内uav、,nuru

,,、,TA&AMWPO円、dF334rupn》

'A

av--94nvuHE'7・

T11J14zzlH12HZH5202HeI

口-

M閃

az,

.

.

. 、ae'τ

MHZ,t‘,官晶-ZEM'1

ゐ明!J'T1

EIJIEq

J6{616{6

2Jf6{G{T1L

{FUEE,tT,《TfTmyfr・-UA?-N

仁』nvnvT'e@@rAF-AMnn幽ynvza匂AnF』府内=nuzCLAm'Ja品目門UMHMHpn

T66

LJI{TMnI772

TM円383TMHKM

11U

a門nu'也,‘,

Rも

T'VApnyE・V-

ロ門TARN-A,uyapn,‘onM叶怠,‘?,nu

?smU

向UZ戸pw十戸「円uon向U,F向υnwn向Unwn向ν,‘内u,、nk向U'十mUマaMnprFt創刊

QJwnunνna』T47・6YA

九日WM鴨E'J'nνM開F・Mmzsvsnuv'M腕CYYAF0・-TATAD円EI

1

l

AVAVnunν

ハυ

F3

0123

&78

今494ヲ'h9g-

SOURCE PR口GRA門 llSTFORIRAN OKJT AC 4500

品目

LlNE NO.

0051

80

90

0054 0055 9自56

0051 0058 0059 0060 0061 0062 0063 0064 0065 0066 0067 0068 0069 0010 0071 0072 0073 0014 0075 0076 0077