знаходження p-значення в грунтовій кореляції в R

26

Чи можна знайти значення р у грушовій кореляції в R?

Щоб знайти кореляцію груші, я зазвичай роблю це

col1 = c(1,2,3,4)
col2 = c(1,4,3,5)
cor(col1,col2)
# [1] 0.8315218

Але як я можу знайти p-значення цього?

r correlation p-value pearson-r

— тубі
джерело

1

cor?corcor.test

— Довідка

32

Ви можете використовувати cor.test:

col1 = c(1,2,3,4) 
col2 = c(1,4,3,5)
cor.test(col1,col2)

що дає:

# Pearson's product-moment correlation   
# data:  col1 and col2   
# t = 2.117, df = 2, p-value = 0.1685   
# alternative hypothesis: true correlation is not equal to 0   
# 95 percent confidence interval:   
#  -0.6451325  0.9963561   
# sample estimates:   
#       cor    
# 0.8315218

Більше інформації про статистику та додаткові параметри на офіційній сторінці:

https://stat.ethz.ch/R-manual/R-patched/library/stats/html/cor.test.html

— брумар
джерело

дякую, не вистачає балів, щоб використати його, хоча

— tubby

9

Якщо потрібно лише значення P:

> cor.test(col1,col2)$p.value
[1] 0.1684782

— rnso
джерело

7

Виконайте наступне:

 library(Hmisc) # You need to download it first.
 rcorr(x, type="pearson") # type can be pearson or spearman

Тут x - це кадр даних, і rcorr повертає кожну кореляцію, яку можна сформувати з кадру даних "x".

Або ви могли самі порахувати статистику:

т = \frac{\hat{ρ}}{\sqrt{\frac{1 - {\hat{ρ}}^{2}}{н - 2}}}

$t = \frac{\hat{\rho}}{\sqrt{\frac{1-\hat{\rho}^2}{n-2}}}$

$\hat{\rho}$

— Репмат
джерело

Дякую, але що таке х? Я думаю, що це деяка конкатенація col1 і col2, тому що нам потрібні два вектори для обчислення кореляції груші. Але ви можете мені сказати, що таке х?

— тубі

Це кадр даних, дивіться моє оновлення.

— Репмат