目的

相関係数のジャックナイフ推定値を求める。

使用法

jack.knife（x, y)

引数

x, y  二変数のベクトル

ソース

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

# 相関係数のジャックナイフ推定
jack.knife <- function(x, y)                                 # 2 変数のベクトル
{
        OK <- complete.cases(x, y)                           # 欠損値を持つケースを除く
        x <- x[OK]
        y <- y[OK]
        n <- length(x)                                               # サンプルサイズ
        est <- cor(x, y)                                     # 標本相関係数
        cor.est <- sapply(1:n, function(i) cor(x[-i], y[-i]))        # i 組目のデータを除いたときの相関係数
        jk <- n*est-(n-1)*cor.est                            # ジャックナイフ推定量
        plot(1:n, cor.est)                                      # 推定値の図を描いてみる
        list(sample.cor=est, jk.est=mean(jk))
}


使用例

x <- 1:10*10
y <- c(100, 150, 154, 178, 40, 204, 211, 220, 279, 290)
jack.knife(x, y)

出力結果例

> x <- 1:10*10
> y <- c(100, 150, 154, 178, 40, 204, 211, 220, 279, 290)
> cbind(x, y)
        x   y
 [1,]  10 100
 [2,]  20 150
 [3,]  30 154
 [4,]  40 178
 [5,]  50  40
 [6,]  60 204
 [7,]  70 211
 [8,]  80 220
 [9,]  90 279
[10,] 100 290
> jack.knife(x, y)
$sample.cor
[1] 0.769872	# ピアソンの積率相関係数

$jk.est
[1] 0.7342372	# ジャックナイフ推定値



 参考：内田治「すぐわかる Excel による統計解析」東京図書