如何模拟类内相关数据？

Question

我正在为混合模型模拟随机数据集。为了使混合模型有意义，群集中的数据点需要关联（ICC）。如何生成此类数据？

这些是我的变量：

• ID
• 基线（基线，首次测量）
• 一半（第二次测量）
• 最终（最终评估）
• 对照组（对照组还是实验组？）
• 集群（显示人所属集群的因素变量）

我需要嵌套在每个群集中的三个度量值相互关联。

 N <- 200  
 ID <- seq(from = 1, to = N, by = 1)
 control <- rbinom(n=N, size = 1, prob = .5) 
 cluster <- as.factor(rbinom(n=N, size = 2, prob = .5) + 1)
 base <- rnorm(n=N,mean = 100,sd=15)
 half <- base + rnorm(n=N, mean = 0, sd = 3) + control*tau_1
 final <-  half + control*(tau_2-tau_1) + rnorm(n=N, mean = 0, sd = 3)
 df.wide <- data.frame(ID,control,age,cluster,base,half,final)

代码基于我在上面编写的内容+添加了一些噪音，因此

rnorm(n=N, mean = 0, sd = 3)

在第二和第三小节中。

如何使三个测量值在群集内相互关联？

谢谢。