在模拟过程中如何保持一些变量固定而其他变量随机？

Question

我是R的新手，我无法找到解决问题的方法。

想象一下，我创建了三个正态分布的变量x，y和z（每个n = 100，平均值为0，sd为1）。变量彼此独立（r~0）。然后我创建了两个新变量x z和y z（即x，z和y，z之间的相互作用）。

现在假设我想在计算交互之前保持x和y固定并从z中随机采样（或创建新的z），并且在每个来自z的新随机样本之后，计算相互作用并计算相互之间的相关系数他们（为了争论，我想这样做100次）并将100个相关性打印到新的数据文件中以供进一步分析。我如何实现这个目标？

Answer 1

不确定这是你要求的。给它一张支票。

x<-rnorm(100)
y<-rnorm(100)
z<-rnorm(100)
xz<-x*z
yz<-y*z

co<-data.frame(nrow=100, ncol=2)
for (i in 1:100){
  z[i]<-rnorm(1)
  xz[i]<-z[i]*x[i]
  yz[i]<-z[i]*y[i]
 co[i,]<-c(cor(z,x), cor(z,y))
}