r関数cor.test（）：ピアソン相関のp値はどのようにして計算されますか？

Question

これはおそらく非常に簡単な質問です。

私はここでR.

Answer 1

でcor.test（）関数内のp値の計算の背後にある方法論は、ピアソンの相関のp値を算出し、コードがある見つけることができません。

x <- c(44.4, 45.9, 41.9, 53.3, 44.7, 44.1, 50.7, 45.2, 60.1) 
y <- c(2.6, 3.1, 2.5, 5.0, 3.6, 4.0, 5.2, 2.8, 3.8) 
ct = cor.test(x, y, method = "pearson") 
ct$p.value ## this is what cor.test() gives 

n <- length(x) 
r <- cor(x, y) 
df <- n - 2 
t = sqrt(df) * r/sqrt(1 - r^2) 
pval = 2 * min(pt(t, df), pt(t, df, lower.tail = FALSE)) ## this is calculated manually 

ct$p.value == pval