目的

スピアマンの順位相関係数行列を計算する。

使用法

spearman(data.matrix)

引数

data.matrix　データ行列（列が変数）

ソース

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

spearman2 <- function(x, y)          # 2変数間のスピアマンの順位相関係数
{
        cor(cbind(x, y), use="complete.obs", method="spearman")[1, 2]
}

spearman <- function(data.matrix)    # スピアマンの順位相関係数行列
{
        cor(data.matrix, use="complete.obs", method="spearman")
}


使用例

d <- matrix(c(1, 3, 5, 2, 2, 4, 3, 1, 23, 2, 4, 4, 4, 3, 3, 5, 2, 1, 3, 5, 3), byrow=TRUE, ncol=3)	# 7行3列のデータ行列
spearman(d)

出力結果例

> d <- matrix(c(1, 3, 5, 2, 2, 4, 3, 1, 23, 2, 4, 4, 4, 3, 3, 5, 2, 1, 3, 5, 3), byrow=TRUE, ncol=3)

> d
     [,1] [,2] [,3]	# このようなデータ
[1,]    1    3    5
[2,]    2    2    4
[3,]    3    1   23
[4,]    2    4    4
[5,]    4    3    3
[6,]    5    2    1
[7,]    3    5    3

> spearman(d)
           [,1]       [,2]       [,3]
[1,]  1.0000000 -0.2129630 -0.6944444
[2,] -0.2129630  1.0000000 -0.3148148
[3,] -0.6944444 -0.3148148  1.0000000

> cor(d, method="spearman")	# cor 関数を使う場合
           [,1]       [,2]       [,3]
[1,]  1.0000000 -0.2129630 -0.6944444
[2,] -0.2129630  1.0000000 -0.3148148
[3,] -0.6944444 -0.3148148  1.0000000

 解説ページ