Существуют ли какие-либо функции, которые могут имитировать переменную отклика из надежной линейной модели или модели квантильной регрессии, как для линейных моделей (например, 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)