R: Konstruowanie skorelowanych zmiennych
Mam zmienną o określonym rozkładzie (normale w moim przykładzie poniżej).
set.seed(32)
var1 = rnorm(100,mean=0,sd=1)
Chcę utworzyć zmienną (var2), która jest skorelowana z var1 ze współczynnikiem korelacji liniowej (w przybliżeniu lub dokładnie) równym „Corr”. Nachylenie regresji między var1 i var2 powinno być (szorstko lub dokładnie) równe 1.
Corr = 0.3
Jak mogę to osiągnąć?
Chciałem zrobić coś takiego:
decorelation = rnorm(100,mean=0,sd=1-Corr)
var2 = var1 + decorelation
Ale oczywiście podczas biegania:
cor(var1,var2)
, wynik nie jest zbliżony do Corr!