R语言分析协变量之间的非线性关系
2020-04-07 本文已影响0人
拓端tecdat
最近我被问到我的R和Stata软件包是否能够适应协变量之间的非线性关系。答案是肯定的，在这篇文章中，我将说明如何做到这一点。
为了说明，我们将模拟具有两个协变量x1和x2以及连续结果y的非常大的数据集。
< - 10000 \nx1 < - rnorm（n）\nx2 < - x1 ^ 2 + rnorm（n）\ny < - x1 + x2 + rnorm（n）\n [（runif（n）< - NA \nmydata < - data.frame（x1， X2，Y）"}"> set.seed（123） n < - 10000 x1 < - rnorm（n） x2 < - x1 ^ 2 + rnorm（n） y < - x1 + x2 + rnorm（n） [（runif（n）< - NA mydata < - data.frame（x1， X2，Y）
need-to-insert-img
因此， 模型的真实系数是0（截距） 。注意，实体模型中没有非线性，但x2对x1的依赖性存在非线性。
< - （mydata，smtype =“lm” ，\n numit = 50，method = c（“”，“norm”，“”））"}"> imps1 < - （mydata，smtype =“lm” ， numit = 50，method = c（“”，“norm”，“”））
need-to-insert-img
< - imputationList（imps1 $ impDatasets）\n"}"> impobj < - imputationList（imps1 $ impDatasets）
need-to-insert-img
输出：
