r - 如何从分位数回归 rq() 中绘制交互效应

Question

我想在 R 中运行分位数回归，其中包括一个三向交互项，其中位数为三个分类预测变量（tau = 0.5）。为了解释效果，我想创建一个情节。对于线性模型，我使用“效果”包中的函数 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)

给