Какова функция, которая предоставит вам нижнюю и верхнюю границы коэффициента корреляции отдельно?

Question

m=c(1,2,5,4,6,8) 
h=c(1,2,9,8,7,3) 
cor(m,h) 
#[1] 0.4093729 

Если оценить коэффициент корреляции (R), то вы можете также оценить 95% доверительный интервал для коэффициента корреляции (R), в результате чего, например, что-то вроде

R = 0.40 [0.33 0.56] 

где «лучше» оценка для R равна 0.40, и есть вероятность, что истина R находится между 0.3 и 0.56. (Обратите внимание, что эти числа полностью составлены.)

Я ищу функцию, которая будет обеспечивать нижнюю и верхнюю границы R отдельно. Чтобы иметь что-то вроде:

R = 0.40 
upper [0.33] 
lower [0.56] 

что-то simiilar этому в MATLAB:

  [R,P,RLO,RUP]=corrcoef(...) also returns matrices RLO and RUP, of the same size as R,    
     containing lower and upper bounds for a 95% confidence interval for each coefficient. 

Answer 1

В "Смотри также" раздел справочной странице cor «s, это говорит

кор .test для доверительных интервалов (и испытаний)

> cor.test(m, h) 

    Pearson's product-moment correlation 

data: m and h 
t = 0.8974, df = 4, p-value = 0.4202 
alternative hypothesis: true correlation is not equal to 0 
95 percent confidence interval: 
-0.6022868 0.9164582 
sample estimates: 
     cor 
0.4093729 

Или получить в интервале более непосредственно:

> x = cor.test(m, h) 
> x$conf.int 
[1] -0.6022868 0.9164582 
attr(,"conf.level") 
[1] 0.95