Чи можливо, що 3 вектори мають усі негативні попарні кореляції?


16

Враховуючи три вектори , і , чи можливо, що кореляції між і , і , і є негативними? Тобто це можливо?abcabacbc

corr(a,b)<0corr(a,c)<0corr(b,c)<0

3
Негативні кореляції означають, що геометрично, що центрировані вектори взаємно складають тупі кути. У вас не повинно виникнути проблем з малюванням конфігурації трьох векторів у площині, які мають це властивість.
whuber

Вони не можуть бути повністю негативно корельованими ( ), але в цілому може бути деяка негативна кореляція, знову ж таки межі, встановлені іншими кореляціями. ρ=-1
карафа

2
@whuber Ваш коментар, здається, суперечить відповіді Хейкі Пулкінкіна, яка стверджує, що вектори в літаку неможливі. Якщо ви стоїте поруч, слід перетворити свій коментар у відповідь.
РМ

2
@RM Не існує суперечності між Уабером та Хейкі. Це питання задає матрицю даних розміром . Зазвичай ми б говорили про точок даних у 3 вимірах, але це Q говорить про три "вектори" в вимірах. Хейкі каже, що всі негативні кореляції не можуть відбутися, якщо (дійсно, дві точки після центрування завжди ідеально співвідносяться, тому кореляції повинні бути і не можуть бути всі ). Вюбер каже, що 3 вектори в розмірах можуть ефективно лежати у двовимірному підпросторі (тобто X - ранг 2) і пропонує уявити логотип Mercedes. Xn×3nnn=2±11nХ
Амеба каже, що повернеться до Моніки

Відповіді:


19

Це можливо, якщо розмір вектора становить 3 і більше. Наприклад

a=(1,1,1)b=(1,9,3)c=(2,3,1)

Кореляції -

cor(a,b)=0.80...cor(a,c)=0.27...cor(b,c)=0,34 ...

Ми можемо довести, що для векторів розміру 2 це неможливо:

cor(a,b)<02(iaibi)(iai)(ibi)<02(a1b1+a2b2)(a1+a2)(b1b2)<02(a1b1+a2b2)(a1+a2)(b1b2)<02(a1b1+a2b2)a1b1+a1b2+a2b1+a2b2<0a1b1+a2b2a1b2+a2b1<0a1(b1b2)+a2(b2b1)<0(a1a2)(b1b2)<0

Формула має сенс: якщо більший за , повинен бути більшим за щоб зробити кореляцію негативною.a 2 b 1 b 1a1a2b1b1

Аналогічно для кореляцій між (a, c) та (b, c) ми отримуємо

(a1a2)(c1c2)<0(b1b2)(c1c2)<0

Ясна річ, що всі ці три формули не можуть утримуватися в один і той же час.


3
Ще один приклад чогось несподіваного, що трапляється лише в розмірі три чи вище.
п -

1
З векторами розміру кореляції зазвичай (пряма через дві точки), і ви не можете мати три кореляції з трьома векторами будь-якого розміру± 1 - 12±11
Генрі

9

Так вони можуть.

Припустимо, у вас є багатовимірний нормальний розподіл . Єдине обмеження на - це те, що воно повинно бути позитивним напіввизначеним.ΣXR3,XN(0,Σ)Σ

Тож візьмемо такий прикладΣ=(10.20.20.210.20.20.21)

Усі його власні значення позитивні (1,2, 1,2, 0,6), і ви можете створювати вектори з негативною кореляцією.


7

почнемо з кореляційної матриці для 3 змінних

Σ=(1pqp1rqr1)

негативна визначеність створює обмеження для парних кореляцій які можна записати якp,q,r

pqrp2+q2+r212

Наприклад, якщо , значення обмежується , що примушує . З іншого боку, якщо , може бути в межах .p=q=1r2rr2+1r=1p=q=12r2±34

Відповідаючи на цікаве подальше запитання від @amoeba: "яка найменша кореляція, яку можуть мати одночасно всі три пари?"

Нехай , Знайдіть найменший корінь , який дасть вам . Можливо, для когось це не дивно.p=q=r=x<02x33x2+112

Більш сильний аргумент можна зробити, якщо одна з кореляцій, скажімо, . З того ж рівняння , ми можемо вивести, що . Тому, якщо два співвідношення дорівнюють , третім має бути .r=12pqp2+q2p=q11


2
Див. Stats.stackexchange.com/questions/72790/… , серед іншого.
whuber

2

Проста функція R для дослідження цього:

f <- function(n,trials = 10000){
  count <- 0
  for(i in 1:trials){
    a <- runif(n)
    b <- runif(n)
    c <- runif(n)
    if(cor(a,b) < 0 & cor(a,c) < 0 & cor(b,c) < 0){
      count <- count + 1
    }
  }
  count/trials
}

Як функція n, f(n)починається з 0, стає ненульовою при n = 3(з типовими значеннями близько 0,06), потім збільшується до приблизно 0,11 за n = 15, після чого, здається, стабілізується:

введіть тут опис зображення Таким чином, не тільки можливо, щоб усі три кореляції були негативними, це не здається надзвичайно рідкісним (принаймні, для рівномірного розподілу).

Використовуючи наш веб-сайт, ви визнаєте, що прочитали та зрозуміли наші Політику щодо файлів cookie та Політику конфіденційності.
Licensed under cc by-sa 3.0 with attribution required.