我正在模拟混合模型的随机数据集。为了使混合模型有意义,需要关联集群内的数据点 (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)
在第 2 和第 3 措施中。
如何使三个测量值在集群内相互关联?
谢谢你。