贝叶斯核机回归估计混合物健康效应 【BKMR包】——实操篇

 1. 下载及加载 BKMR 包

# 下载方式1
install.packages("bkmr")

# 下载方式2
install.packages("devtools")
devtools::install_github("jenfb/bkmr")

# 导入
library(bkmr)

2. 两种变量选择方法

2.1 分量变量选择（component-wise variable selection ）

2.1.1 生成一些数据，以备分析

set.seed(111)
dat <- SimData(n = 50, M = 4) # 生成一个包含50个样本(observations)和4个预测变量(predictor variables)的数据集。
y <- dat$y
X <- dat$X
Z <- dat$Z

# 查看生成的数据
head(y)
# [1]   6.325343   6.709022  11.462322 -12.538379   3.040817   8.985771
head(Z)
#    z1         z2         z3          z4
# [1,]  0.2352207  0.1916634  0.5996197 -0.09152026
# [2,] -0.330