所以我使用rpart创建了一个回归树,它输出到reg_tree
# show summary statistics of reg_tree
summary(reg_tree)
# store top variables as new values
topvars <- reg_tree$variable.importance
# output of topvars
topvars
q_21fb1900 q_2b3296a0 q_0 q_fde6a01e q_7fa850ed q_323d6cee q_c6ab3657 q_eb2ad90d q_5dcb2b57
5.303283e+15 5.196871e+15 4.002239e+15 4.412505e+14 2.616730e+14 2.162128e+14 2.035465e+14 1.354927e+14 5.095959e+13
q_af2830be q_caa61b2c q_a6828865 q_99f5a0bd q_be83fe28 q_efdc29dd q_9e86aa7f q_2ea0e2aa q_5049294d
2.176437e+13 1.210118e+13 1.126591e+13 8.387189e+12 4.951978e+12 4.115929e+12 3.864235e+12 1.449853e+12 5.436949e+11
q_5ae0f0cd q_518fba14
5.436949e+11 5.412242e+11
我想将这些名称中的每一个提取为 xvar1、xvar2 并自动将它们放入以下模型中,其中每个 xvar 都对应于列标题:lm(y_var ~ xvar1 + xvar2 + xvar3 + ... +,data)
.
IE
lm(y_var ~ q_21fb1900 + q_2b3296a0 + q_0 + ... +,data)
我该如何做到这一点,以便我可以放入新的数据集,而不必担心将来自定义更改每个 xvar?