Создайте объект lm из данных/коэффициентов
Кто-либо знает о функции, которая может создать объект lm, учитывая набор данных и коэффициенты?
Я интересуюсь этим, потому что я начал играть с Усреднением байесовской модели (BMA), и я хотел бы смочь создать объект lm из результатов bicreg. Я хотел бы иметь доступ ко всем хорошим универсальным функциям lm как диагностическое графическое изображение, предсказать, cv.lm и т.д.
Если Вы вполне уверены, такая функция не существует, это также очень полезно для знания!
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)
2 ответа
Нет функции, что я знаю, что делает это. Можно было, конечно, быть сделано. Все, что вам нужно сделать вашу MagicFunction, это создать список с элементами:
> names(fakeModel)
[1] "coefficients" "residuals" "effects" "rank"
[5] "fitted.values" "assign" "qr" "df.residual"
[9] "xlevels" "call" "terms" "model"
Тогда сделайте его объектом LM
> class(fakeModel) <- c("lm")
, позвольте мне просто сказать, что я думаю, что это плохое представление. Чем сказать, что общая функция, которую вы применяете, будут применимы к объекту Bicreg. Например, как бы вы интерпретировали AIC (FAKEMODEL)?
Вы лучше создают свои собственные функции, чтобы сделать диагностику и прогноз.
кажется, что можно вычислить ваш объект lm, как обычно, и затем изменить коэффициенты впоследствии путем изменения $coefficients атрибут вашего лм () результат.
Дополнительную информацию см. в этом вопросе и результатах:
http://tolstoy.newcastle.edu.au/R/e2/help/07/08/24294.html
Не уверенный это соответствует тому, что вы хотите сделать, хотя...