rotatefactors
因子負荷量の回転
構文
B =旋转因子(A)
B = rotatefactors(A,'Method','orthomax','Coeff',gamma)
B = rotatefactors(A,'Method','procrustes','Target', Target)
B = rotatefactors(A,'Method','pattern','Target', Target)
B = rotatefactors(A,'Method','promax')
[B,T] =旋转因子(A,…)
説明
B =旋转因子(A)
は,d行米列の負荷配列一个
を回転して,バリマックス基準を最大化し,B
に結果を返します。一个
とB
の行は変数に対応し,列は因子に対応します。たとえば,一个
の(i, j)番目の要素は,j番目の因子のi番目の変数の係数です。配列一个
は通常主成分分析
またはpcacov
で作成された主要要素の係数,またはfactoran
で推定された因子負荷量を含みます。
B = rotatefactors(A,'Method','orthomax','Coeff',gamma)
は,一个
を回転して,係数γ
をもオソマックス基準を最大化します。たとえば,B
は以下を最大化する一个
の直交回転です。
sum(D*sum(B.^4,1) - GAMMA*sum(B.^2,1).^2)
γ
の既定値1は,バリマックス回転に対応しています。この他にγ
に設定する可能性がある値には0 m / 2およびd (m - 1) / (d + m - 2)があり,それぞれクォーティマックス,エカマックス,パーシマックスに対応します。また,“方法”
パラメ,タ,に“最大方差法”
、“quartimax”
、“equamax”
または“parsimax”
を指定し,多项式系数的
パラメ,タ,を省略できます。
“方法”
が“orthomax”
、“最大方差法”
、“quartimax”
、“equamax”
,または“parsimax”
の場合,追加パラメ,タ,は以下のとおりです。
“正常化”
-回転のために負荷行列を行正規化するかどうかを指定するフラグ。“上”
(既定の設定)の場合,一个
の行は、ユ、クリッドノルムをも、ために回転前に正規化され、回転後に非正規化されます。“关闭”
の場合,生の負荷は回転されて返されます。“Reltol”
- - - - - -T
を見けるために使用される反復アルゴリズムの収束の相対許容誤差。既定の設定はsqrt (eps)
です。“麦克斯特”
- - - - - -T
を見けるために使用される反復アルゴリズムの反復制限。既定の設定は250
です。
B = rotatefactors(A,'Method','procrustes','Target', Target)
は,一个
からd行米列のタ,ゲット負荷行列目标
の斜交プロクラステス回転を行います。
B = rotatefactors(A,'Method','pattern','Target', Target)
は,負荷行列一个
からd行米列のタ,ゲットパタ,ン行列目标
の斜交回転を実行し,結果をB
に返します。目标
は,B
の“制限された”要素を定義します。すなわ,目标
の0要素に対応するB
の要素は,小さく制限され,目标
の非ゼロ要素に対応するB
の要素は,任意の大きさを取ることが可能です。
“方法”
が“普罗克汝斯忒斯之
または“模式”
の場合,追加パラメ,タ,は回転の型である“类型”
です。“类型”
が“正交”
である場合,回転は直交であり,因子は無相関な状態から変化しません。“类型”
が“斜”
(既定の設定)の場合,回転は斜交で回転した因子は相関となる可能性があります。
“方法”
が“模式”
の場合,目标
にいて制限があります。一个
が米個の列をも場合,直交回転では,目标
のj番目の列には少なくとも米-j個のゼロが含まれていなければなりません。斜交回転では,目标
の各列は,少なくとも米- 1個のゼロが含まれている必要があります。
B = rotatefactors(A,'Method','promax')
は,一个を回転して,一个プロマックス基準を最大化します。これは,オソマックス回転によって作成されたタゲットをも斜交プロクラステス回転に相等します。プロマックスにより内部で使われるオーソマックス回転を制御するために,4つのオーソマックスパラメーターを利用します。
“プロマックス”の追加パラメ,タ,は,プロマックスタ,ゲット行列を作成するための指数,“权力”
です。“权力”
は1
以上でなければなりません。既定の設定は4
です。
[B,T] =旋转因子(A,…)
は,B = a * t
でB
の作成に使用される回転行列T
を返します。发票(T ' * T)
を使って,回転係数の相関行列を算出できます。直交回転に対しては,これは単位行列で,斜交回転に対しては,単位対角要素をもちますが,これは非ゼロの非対角要素です。
例
rng('default') % for再现性X = randn(100,10);默认(归一化变量)旋转:%前三个主成分。LPC = pca(X);[L1,T] = rotatefactors(LPC(:,1:3));平均旋转:%前三个主成分。[L2,T] = rotatefactors(LPC(:,1:3),…“方法”、“equamax”);Promax旋转:%前三个因子。LFA = factoran(X,3,'Rotate','none');[L3,T] = rotatefactors(LFA(:,1:3),… 'method','promax',... 'power',2); % Pattern rotation: % first three factors. Tgt = [1 1 1 1 1 0 1 0 1 1; ... 0 0 0 1 1 1 0 0 0 0; ... 1 0 0 1 0 1 1 1 1 0]'; [L4,T] = rotatefactors(LFA(:,1:3),... 'method','pattern',... 'target',Tgt); inv(T'*T) % Correlation matrix of the rotated factors ans = 1.0000 -0.9593 -0.7098 -0.9593 1.0000 0.5938 -0.7098 0.5938 1.0000
参考文献
《现代因素分析》第三版。芝加哥:芝加哥大学出版社,1976年。
[2]劳利、d.n.和a.e.麦克斯韦。《作为统计方法的因素分析》第二版,纽约:美国爱思唯尔出版社,1971年。
バ,ジョン履歴
R2006aより前に導入