您如何使用交互式(非组合)预测器调用模型?

How do you call a model using interactive (not combined) predictors?

如何调用使用交互式(非组合)预测变量的模型?

我有一个非常有效的回归模型:

lm_comb1 <- lm(cbind(MD_EARN_WNE_P10, X40.year.NPV) ~ SAT_AVG*TUITIONFEE_OUT*AVGFACSAL* RET_FT4_POOLED*INEXPFTE*C100_4_POOLED,data = training_set_NORM)

当我调用summary(lm_comb1)时,我没有遇到任何问题。我想检查添加多项式会如何影响事物,所以我适合:

lm_comb1_POLY2 <- lm(formula = cbind(training_set_NORM$MD_EARN_WNE_P10, training_set_NORM$X40.year.NPV) ~ poly(training_set_NORM$SAT_AVG,2,raw=TRUE)*poly(training_set_NORM$TUITIONFEE_OUT,2,raw=TRUE)*poly(training_set_NORM$AVGFACSAL,2,raw=TRUE)*poly(training_set_NORM$RET_FT4_POOLED,2,raw=TRUE)*poly(training_set_NORM$INEXPFTE,2,raw=TRUE)*poly(training_set_NORM$C100_4_POOLED,2,raw=TRUE))

但是,当我添加 summary(lm_comb1_POLY2) 时,R-markdown 不会编织输出。它一直到 100%,但从未产生实际输出。我试过编织成 HTML 以及 PDF。我没有收到错误消息,但是当我删除 summary(lm_comb1_POLY2) 时,输出是编织的并且制作得很好。

有什么想法吗?我对 m_comb1_POLY2 的称呼有什么古怪之处吗? * 而不是 + 操作员是否以某种方式搞砸了?

注意:我没有提供任何数据,因为这个问题似乎不需要重现性。

干杯!

我突然想到要尝试

lm_comb1_POLY2 <- lm(cbind(MD_EARN_WNE_P10, X40.year.NPV) ~ poly(SAT_AVG*TUITIONFEE_OUT*AVGFACSAL* RET_FT4_POOLED*INEXPFTE*C100_4_POOLED,2,raw=TRUE), data = training_set_NORM)
summary(lm_comb1_POLY2.1)

这似乎有效。不过,如果有专家出现,请纠正我。希望这对其他人有帮助!