相関係数の差の検定（同一集団が1変数を共有している場合） Rの関数

参考：Weaver, B., & Wuensch, K. L. (2013). SPSS and SAS programs for comparing Pearson correlations and OLS regression coefficients. Behavior research methods, 45(3), 880-895. のp. 892

目的：同じ集団において，変数Aと変数C，変数Bと変数Cの相関係数が有意に異なるか知りたい

用意するもの：r12（変数Aと変数Cの相関），r13（変数Bと変数Cの相関），r23（変数Aと変数Bの相関），n（サンプルサイズ）

以下関数

cor_significance_test <- function(n, r12, r13, r23){

R = 1 - r12^2 - r13^2 - r23^2 + 2*r12*r13*r23

t=(r12-r13)*sqrt((n-1)*(1+r23)/(R*2*((n-1)/(n-3))+((r12+r13)^2/4)*(1-r23)^3))

p=(1-pt(abs(t), df=n-3))*2

print(t)

print(p)

}

t値とp値を返す