数値例というお話が出ましたので,私が試行錯誤でやってみた経過を書いてみたいと
思います(長文).
有馬哲・石村貞夫1987多変量解析のはなし 東京図書 のP80表3.1の
数値例をためしに入力してみました.
サンプル 変量1 変量2
(スポーツ施設)(教育施設)
1 22.9 13.7
2 24.9 16.2
3 19.3 11.3
4 22.0 10.4
5 28.6 24.9
6 42.6 26.5
7 41.3 20.3
これをSASで
PROC FACTOR METHOD=P PRIORS=ONE COV;で分析すると以下のように出力されます
1 2
EIGENVALUE 119.3460 10.9421
DIFFERENCE 108.4040
PROPORTION 0.9160 0.0840
CUMULATIVE 0.9160 1.0000
FACTOR PATTERN
FACTOR1 FACTOR2
V1 0.98254 -0.18604
V2 0.90003 0.43582
また PROC PRINCONP COV;で分析すると次のようになります
EIGENVALUES OF THE COVARIANCE MATRIX
EIGENVALUE DIFFERENCE PROPORTION CUMULATIVE
PRIN1 119.346 108.404 0.916016 0.91602
PRIN2 10.942 . 0.083984 1.00000
EIGENVECTORS
PRIN1 PRIN2
V1 0.847878 -.530192
V2 0.530192 0.847878
HALBAU(Ver4)では以下のように出力されます.分散・共分散による主成分負荷量は
1を越えています.
【主成分分析:主成分負荷量】(相関係数による)
-----------------------------------
変数名 成分 1 成分 2
-----------------------------------
1)スポーツ施設 0.9495 0.3137
2)教育施設 0.9495 -0.3137
-----------------------------------
固有値 1.8032 0.1968
累積固有値 1.8032 2.0000
寄与率(%) 90.162 9.838
累積寄与率(%) 90.162 100.000
カイ2乗値 5.181 0.000
(自由度) ( 2) ( 0)
有意確率 0.07498 1.00000
-----------------------------------
【主成分分析:主成分負荷量】(分散・共分散による)
-------------------------------------------
変数名 成分 1 成分 2
-------------------------------------------
1)スポーツ施設 8.57559D+00-1.62371D+00
2)教育施設 5.36245D+00 2.59663D+00
-------------------------------------------
固有値 1.02297D+02 9.37892D+00
累積固有値 1.02297D+02 1.11676D+02
寄与率(%) 91.602 8.398
累積寄与率(%) 91.602 100.000
カイ2乗値 5.893 0.000
(自由度) ( 2) ( 0)
有意確率 0.05253 1.00000
-------------------------------------------
数値例を掲載していた本自身の解では(分散共分散行列の場合のみ)
固有ベクトルの値は以下のようにSASと一致しますが,
第1主成分 第2主成分
変量1 .8479 -.5301
変量2 .5301 .8479
固有値 119.3 10.94
(固有値の平方根*固有ベクトル)を求めても以下のようにHALBAUの結果
とは一致しません
第1主成分 第2主成分
変量1 9.261 -2.777
変量2 5.789 2.804
しかし 柴山(1992)欠測値を含む多変量データのための主成分分析的方法
教育心理学研究40,257- の方法をプログラム化したもの(MPCA.EXE 柴山氏自身の
作成)に,この(欠測値を含まない!)データを入力すると,「重み係数」として
出力される数値はHALBAUの分散・共分散による主成分負荷量と完全に一致します.
また同論文での「重み係数」の値も-1〜+1の範囲には収まっていないようです
(Table 4 P261).ということから,あながちHALBAUがオカシイとも言い切れない
ように素人考えでは思うのですが.
(とはいえ欠測値のためのプログラムに欠測値のないデータを入れるという無理を
しているので,同列の話として論じてよいものかどうか...)
数値例(欠損値なし)に関しては,他にも幾つかのものをやってみましたが,
どれも,一致するものと,一致しないもの のパターンは同じでした.
Keizo Hori <hori (at) ec.kagawa-u.ac.jp> wrote:
>>(唯一,一冊だけ「(固有値のルート)×固有ベクトルのことだ」とだけ,
>>書かれている本があっただけです.)
>
>どの教科書をあたったのでしょうか。教科書ではないですが、
この記述は次の本にみつけました.
菅民郎1993多変量解析の実践 上 現代数学社 Pp141-142
>-----------------------------------------------------------------------------<
Tsutomu Okada
E-mail : okada (at) ed.niigata-u.ac.jp
Nifty Serve : HCF00154 (at) niftyserve.or.jp
Niigata-University Department of Educational Psychology
<----------------------------------------------------------------------------->
ここは心理学研究の基礎メーリングリストに投稿された過去の記事を掲載しているページです。