¿Alguien sabe de una función que pueda crear un objeto lm dado un conjunto de datos y coeficientes?

Estoy interesado en esto porque empecé a jugar con el promedio de modelo bayesiano (BMA) y me gustaría poder crear un objeto lm a partir de los resultados de bicreg. Me gustaría tener acceso a todas las funciones genéricas agradables de lm, como el diagnóstico, la predicción, cv.lm, etc.

Si está bastante seguro de que tal función no existe, ¡también es muy útil saberlo!

library(BMA)
mtcars_y <- mtcars[, 1] #mpg
mtcars_x <- as.matrix(mtcars[,-1])
res <- bicreg(mtcars_x, mtcars_y)

summary(res)
res$postmean # bma coefficients

# The approximate form of the function
# I'm looking for
lmObject <- magicFunction(data=mtcars, coefficients=res$postmean)