Preciso gerar séries de dados que possuem correlações definidas usando o R. Utilizei um método que encontrei aqui no SO (Como gerar variáveis correlacionados no R?) e consegui criar as variáveis com a correlação desejada, no entanto, ao tentar automatizar este processo para a criação de 1000 estimações e para correlações diferentes, o resultado obtido é uma matriz 1000x5 com todos os valores idênticos. O código que estou utilizando é o seguinte:
set.seed(2423049)
corr = matrix(,1000,5)
for(k in 1:5){
for (i in 1:5){
for(j in 1:1000){
rho = c(-0.7,-0.3,0,0.1,0.5) # correlações que preciso utilizar
xstar=rnorm(1000,2,2) # x* com distribuicao normal N(2,2)
a2=rnorm(1000,2,2) # parametro criado para obter w a partir da correlacao rho com x*
w = rho[k]*xstar+sqrt(1-rho[k]^2)*a2 # w calculado a partir de uma correlacao definida com x
corr[j,i]=cor(xstar,w) # matriz de correlacoes entre x* e w
}
}
}
Através deste processo, o resultado obtido foi uma matriz 1000x5 onde todos os valores eram 0.5499732
O que estou fazendo de errado?