如何从分位数回归rq（）中绘制交互效果

Question

我想在R中运行分位数回归，其中包括一个三向交互项，其中三个分类预测因子位于中位数（tau = 0.5）。为了解释效果，我想创建一个情节。对于线性模型，我使用包“effect”中的函数effect（）。不幸的是，这不适用于“quantreg”包中的rq（）。您对如何显示交互效果有什么建议吗？

library(quantreg)
data(mtcars)
summary(rq(formula = mpg ~ vs * am * carb, tau = 0.5, data = mtcars))

Answer 1

您可以使用visreg套餐，但必须简化（由您）三方互动，一次绘制两个条款，例如：

mod <- rq(formula = mpg ~ vs * am * carb, tau = 0.5, data = mtcars)

library(visreg)
visreg(mod, "carb", by="am", overlay=TRUE)

给出