是否有任何函数可以模拟来自稳健线性模型或分位数回归模型的响应变量,例如线性模型(即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)