从r中随机正态变量生成相关的总和量表

Question

我想在r中生成5个随机正态变量，每个变量具有给定的均值和sd

x1 <- rnorm(100,20,5)
x2 <- rnorm(100,22,4)
y1 <- rnorm(100,21,3)
y2 <- rnorm(100,23,4)
y3 <- rnorm(100,22,5)

并且使用具有预定义相关性的总和标度x和y作为

x <- rowSums(cbind(x1,x2))
y <- rowSums(cbind(y1,y2,y3))
cor(x,y) # must be about 0.6

我认为使用循环应该是一个低效的解决方案。 y=r*x+sqrt(1-r*r)*rnorm(n)也不适用。我有点困惑如何开始。任何建议都将不胜感激。