• 検索結果がありません。

大阪府立大学工学部狩野裕 1.はじめに

N/A
N/A
Protected

Academic year: 2021

シェア "大阪府立大学工学部狩野裕 1.はじめに"

Copied!
3
0
0

読み込み中.... (全文を見る)

全文

(1)

106 統計数理 第39巻 第1号 1991

Rubin(1956)).因子数を后より大きくしたとき,特殊因子に対応する負荷が追加されるだけで 共通因子行列に関しては一意であるための十分条件(Tumura and Sato(1980)).

 新しい結果として以下の3点を報告した.

 (1)Anderson and Rubinによる必要条件を拡張し一意であるかどうかを調べやすくした.

 (2)行列の次数が大きいと部分行列の階数がおちることがある.そこで因子負荷行列が

(二1:㌻)

という形に対して,AndersonandRubin,T㎜uraandSatoの十分条件をみたすための必要/

十分条件を与えた.

 (3)「大部分の要素が一意である行列」を提案し,そのための十分条件を与え利用例を示し た.行列が一意でないときには,何を推定しているのか一般に不明とたる.しかし,行列の一 部の要素のみが不定で他の多くの要素が一意であることがある.両者を区別することにより,一 意てたい行列のすべての要素の推定を無意味とせずに一意である要素に対する推定を有効とし

うるので,有意義である.

 本報告はその後の進展を含めて掲載予定(Sato(1992))である.

      参考文献

Anders㎝,T.W.and Rubin,H、(1956).Statistica1inference in factor ana1ysis,P肌r励〃B励e妙     ∫ツmヵ.om Mα砺.∫勉左ゐ左Pmろ.,Vol.5,111−150,Univ.of Caiifomia Press,Berkeley.

Konishi,S.(1979).Asymptotic expansions for the distributions of statistics based on the samp1e     corre1ation matrix in principa1component ana1ysis,〃κo∫〃mαMα肋.∫,9,647−700.

Sato,M.(1989).Some comments on Shapiro s paper:identiiability of factor analysis,Tech.Report,

    No.249,Statistica1Research Group,Hiroshima University,Hiroshima.

Sato,M.(1992).A study of an identi丘。ation problem and a substitute use of principa1component     ana1ysis in factor analysis,肋ms〃mαMα励、∫,22(to appear).

Shapiro,A.(1985).Identiiability of factor ana1ysis:some resu1ts and open prob1ems,〃m〃λ像eろm

    /1力ψム,70,1−7.

Tumura,Y.and Sato,M、(1980).On the identiication in factor ana1ysis,皿σMα肋.,16(2),121−131.

因子分析モデルにおける不適解の発生構造       大阪電気通信大学工学部猪原正守       大阪府立大学工学部狩野裕

 1.はじめに

 因子分析モデルによって多次元データを解析する中で,独自分散推定値の一部がゼロになる 問題が出現することがある.この問題は不適解問題と呼ばれ,その発生メカニズムについては,

Drie1(1978)・によるモンテカルロ実験を始めとして多くの実験研究が行われている.一方,

Law1ey and Maxwe11(1971)は,最尤因子分析モデルを用いた解析において,力十1(=がと

おく)次元観測変量κ=(κ、,...,伽。、)τに対する尾因子分析モデルの適合によって,初めのm

変量κ、,...,κmに対する独自分散ψ、,...,ψmがゼロと推定されたとき,モデルを

(2)

(1.1)

多変量解析における潜在変数モデルの理論と応用

1一

i:1)一(二1:t)(ll)・(二)

107

と変更することを考えた.ただし,X、=(κ、,_,κm)τ,κ。=(κm。、,_,加。、)Tであり,λ。へ,λ。、,

λ・・のサイズはそれぞれm×m,(が一m)×m,(が一m)×(尾一m)であって,誤差変量e・の分 散ψ。は正値定符号行列である.そして,彼等はMaxwe11(1961)のグ=10変量に対する々=

4因子をもつモデルによる解析によって,第8変量に対する独自分散ψ。の推定値がゼロとたる が,モデル(1.1)は適合していると主張した.しかし,このデータに対してはκ。を除去した周 辺データがに対してが=λ8.十eの因子分析モデルを当てはめると,后=3因子分析モデルが 適合し,モデル(1.1)は必ずしも識別可能でたいことを明らかにする.即ち,彼等の方法論と Maxwe11(1961)データに対する結果は重大た疑問が残っていることを指摘する.

 2.モデルにおける幾つかの命題

 ここでの議論は因子分析モデルにおける識別問題と直接関わる.この領域の研究はAnderson and Rubin(1956)を始めとしてTumuraandSato(1980)など多くの人々によって行われて

きている.

 因子分析において,が次元観測変量κの分散共分散行列Σは  (2.1)       Σ=λλT+Ψ

と表現される.ここで,行列λ(が×后)はランク尾,が次元の対角行列Ψは正値定符号である.

 さて,我々の議論の大前提として次の仮定を導入する.

 仮定.が次元観測変量κを第1変量κ。と残りの力次元ベクトル兀。に分割するとき,条件

付き変量κ・1κ・の分散共分散行列Σ・・.1に対して,

       τ

 (2.2)         Σ・・.・=λ…λ…十Ψ…

の分解が存在する.ただし,行列λ。.1は(力X后)であり,Ψ。.。は正値定符号対角行列である.ま た,(2.2)における行列λ。.、は,任意の1行を削除しても互いに素なランク尾の2つの部分行列 が存在する.

 このとき,次の3つの命題を証明することができる.たお,以降においてm=1とし,σ、、=

亙(κ書),σ。、=亙(κ。κ、),Σ。。=亙(エ。㌶)とする.

 命題1.仮定の下σ。、がただ1つのnon−zero要素をもったらばx。に対する分散共分散行 列Σ。。は識別可能た后因子分解をもつが,Σに対して存在する因子分解はあらゆる因子数に対 して識別不可能である.

命題2.仮定の下で,次の3条件は互いに同値である.

(i) rank(λ2.1,621)=后

(ii)行列Σが后因子分解をもつ:Σ=λ五丁十Ψ

(iii)行列Σ。。が尾因子分解をもつ Σ。。=λ。。五£十Ψ。.、

このとき,(ii)と(iii)の分解は識別可能である.

(3)

ユ08

統計数理 第39巻 第1号 1991

 命題3.仮定およびΣ。。が后因子分解        Σ。。=λ。。λ品十Ψ。

をもち,行列λ。。(力x后)は,任意の2行を削除しても互いに素たランク后の2つの部分行列が 存在するとする.このとき,rank(λ。.。,σ。、)=尾十1とたるのは適当たスカラーa(>0)と添字 プ(2≦ブ≦力十1)が存在して,

       ψ。=伊。.、十あeタ

が成り立つときであり,かっこのときに限る.ここに,e5は第ブ要素のみ1の力次元単位ベク

トルである.

系.命題2の条件の下で,rank(λ。.、,σ。1)=々十1ならば,行列Σのあらゆる因子数に対す る分解は識別可能てたい.

 3.データヘの適用と討論

 以上,我々はが次元観測変量κの特定の1変量κ。に関する残りの力変量兀。の条件付き変 量κ。1κ、が后因子分析モデルに従うという前提条件の下に,観測変量κに対する分散共分散行 列がどのような構造をしているかを明らかにしてきた.とくに,Law1ey andMaxwe11(1971)

によるモデル(ユ.1)に対して,変量 。自身が同じ因子数をもつ因子分析モデルに従う場合の疑 問点を明らかにした.実際,得られた命題をMaxwe11(1961)データおよびDavis(1944)デー タに適用するたらば,それぞれにおいて,モデル(1.1)による解析には重大な疑問のあることが 明らかにたる.

 ここでの議論は,多次元観測変量xの分散構造が因子分析モデルに従うかどうかをエの2つ の部分要素ベクトルκ、とκ。の条件付き変量κ。1κ、が因子分析モデルに従うという条件のも

とで考えたものであった.しかし,得られた命題はκ。がスカラーである場合に限定されている.

応用面からの要求は満たされるかもしれたいが,理論的には一般のベクトルにまで拡張できる かどうかを検討するべきであり,今後の課題として残っている.

       参考文献

Anderson,T.W.and Rubin,H.(1956).Statistica1inference in factor analysis,P肌r肋〃Be伽妙     ∫ヅmφ.om Mα肋.∫械赦Pmろ.,Vol.5,11ユー150,Univ.of Ca1ifomia Press,Berke1ey.

Davis,F.B.(1944).Fundamental factors of comprehension in reading,Psツ。ゐ。me切肋,9,185−197.

Drie1,O.P.von(1978).On various causes of improper so1utions in maximum1ike1ihood factor ana1ysis,

    P砂。ゐ。me切肋,43,225−243.

Law1ey,D.N.and Maxwe11,A.E.(1971).地。C07λmめsタs㏄α∫肋κ∫枕αZ Me肋。a,2nd ed.,Butterworth,

    London.

Maxwe11,A.E.(1961)、Recent trends in factor analysis,∫Roツ∫勉桃た∫oα∫mλ,124,49−59.

Tumura,Y.and Sato,M.(1980).On the identiication in factor ana1ysis,TRσMα肋.,16,121−131.

参照

関連したドキュメント

岡山県立大学保健福祉学部紀要 第25巻1号2018年

患者への治療機会を最大限確保するため、重症度分類が中等症以上、または、65 歳以上の高齢者及

今回の研究では産科受診者を対象にすることで、妊娠を契機に退職した者を

[r]

16 ABCX モデルは Hill によって開発された.その後,McCubbin と

温州ミカンの生理落果は,一般に5月上旬の開花磯潮   から6月下旬まで起こる 早期落果

可逆水素電極の合金材料として最もよく取り上げられているのはLaNi5である。初期の

 各項目に対するプロマックス回転後の因子負荷量は Table 3 に示すとおりである.因子の命 名について,第 1 因子は, “皆が友人を大切にすること” ,