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