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

Question

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

假设我创建了三个正态分布变量 x、y 和 z（n=100，均值为 0，sd 为 1）。变量彼此独立（r~0）。然后我创建了两个新变量 xz 和 yz（即 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))
}