孟德尔随机化的R语言代码

以下是使用R语言实现孟德尔随机化的一个例子：

# 定义实验对象
subjects <- c("A", "B", "C", "D", "E", "F")

# 定义处理和对照组
treatment_group <- c()
control_group <- c()

# 随机分组
for (subject in subjects) {
  if (runif(1) < 0.5) {
    treatment_group <- c(treatment_group, subject)
  } else {
    control_group <- c(control_group, subject)
  }
}

# 打印结果
print(paste("Treatment Group:", treatment_group))
print(paste("Control Group:", control_group))

在这个例子中，我们使用`c()`函数定义实验对象、处理组和对照组。然后使用`runif()`函数生成一个随机数，在每个实验对象上进行随机分组，并将其添加到相应的组中。最后，打印出处理组和对照组的结果。这个例子仅供参考，实际应用中需要根据实际情况进行调整和优化。