偏相関係数     Last modified: Jun 28, 2004

目的

偏相関係数を計算する

使用法

partial.cor(x)

引数

x	データ行列(行がケース,列が変数)

ソース

インストールは,以下の 1 行をコピーし,R コンソールにペーストする
source("http://aoki2.si.gunma-u.ac.jp/R/src/partial_cor.R", encoding="euc-jp")

# 偏相関係数行列
partial.cor <- function(x)                                   # データ行列
{
        x <- subset(x, complete.cases(x))                    # 欠損値を持つケースを除く
        i <- solve(cor(x))                                   # 相関係数行列の逆行列
        d <- diag(i)                                         # 対角成分
        i <- -i/sqrt(outer(d, d))                            # 偏相関係数行列
        diag(i) <- NA                                                # 対角成分は未定義
        rownames(i) <- colnames(i) <- paste("Var", 1:ncol(x))
        return(i)
}


使用例

x <- matrix(c(	# 5ケース,4変数のデータ行列例(ファイルから読んでも良い)
	1,  5, 6, 4,
	2, 14, 5, 3,
	3,  3, 4, 2,
	4,  2, 6, 6,
	3,  4, 3, 5
), ncol=4, byrow=TRUE)

partial.cor(x)

出力結果例

           Var 1       Var 2       Var 3      Var 4
Var 1         NA -0.35532521 -0.34405087  0.3732220
Var 2 -0.3553252          NA  0.03407991 -0.2332674
Var 3 -0.3440509  0.03407991          NA  0.3823565
Var 4  0.3732220 -0.23326740  0.38235645         NA

・ 解説ページ


・ 直前のページへ戻る  ・ E-mail to Shigenobu AOKI

Made with Macintosh