是否有功能可以模拟来自稳健的线性回归或分位数回归模型的数据？

Question

是否有任何函数可以像线性模型（例如stats::simulate.lm）那样模拟健壮的线性模型或分位数回归模型的响应变量？

如果没有，是否有一种方法可以使代码适用于任何一个模型？

以下是我正在处理的数据和模型类型的示例：

#Data
df <- data.frame(Response = c(1:30 + rnorm(n = 30)), Covariate =  c(seq(from = 10, to = 1.3, by = -0.3)))

#Robust linear regression
fit.rlm <- MASS::rlm(Response ~ Covariate, data = df)

#Quantile regression
fit.qr <- quantreg::rq(Response ~ Covariate, data = df, tau = c(0.025,0.5,0.975)