偏相関係数を一度に計算する方法

[質問]

CORRプロシジャでPARTIALステートメントを使用すると偏相関係数を計算できますが、他の変数による影響を除去した個々の2つの変数間の偏相関係数を、単相関係数と同じように行列形式で出力することはできないようです。何かよい方法はありますか？

[回答]

FACTORプロシジャを用いることにより実現できます。PROC FACTORステートメントでMSAオプションを指定してください。
出力の「Partial Correlations Controlling all other Variables」の箇所に行列形式で出力されます。

プログラム例

proc factor data=test msa;     /* MSAオプション*/
  var x1-x20;
run;