一个人旅行*-*的博客当拟合二分类BKMR时,可能会得到如下报错信息: "Error in checkSymmetricPositiveDefinite(H, name = “H”) : H must be positive definite" 或 “Error in if (log(runif(1)) <= logalpha) { : missing ...
探路者Myra的博客基础BKMR模型首先,对于中介变量Mi的模型(12): [ Mi = β0 + ℎM Ai + ϵMi ]然后,对于健康结果Yi的模型(13): [ Yi = θ0 + ℎY Ai, Mi + ϵYi ]最后,对于总效应的混合物模型(14): [ Yi = γ0 + ℎTE Ai ...