У меня есть фрейм данных со 126 столбцами, и я хочу запустить линейные регрессии по всем столбцам.Я сделал это с помощью функции lapply ():
my_lms <- lapply(1:126, function(x) df[,x] ~ df$x1))
Снова с помощью функции lapply () я получаю сводную статистику:
lapply(my_lms, summary)
Мой вопрос заключается в том, как получить такое же резюместатистика с использованием ошибок Ньюи-Уэста?
Я пытался:
coeftest(my_lms, vcov. = NeweyWest)
, но это выдает мне сообщение об ошибке:
Ошибка в UseMethod ("estfun"): не применимометод для 'estfun', примененный к объекту класса "list"
Спасибо