В настоящее время я использую модель фиксированных эффектов с использованием функции plm
из одноименного пакета. У меня есть кластеризация данных, поэтому я использую функцию vcovCR
из пакета clubSandwich
, чтобы получить стандартные ошибки оцененных коэффициентов регрессии.
Я хочу сделать следующее: создать график, который показывает прогнозируемое значение Y при заданном X1 (когда другие переменные сохраняются постоянными при их средних значениях), а также добавить к графику диапазоны уверенности 95%.
Что я пробовал:
- Использование функции
ggpredict()
из пакетаggeffects
. Это не сработало, потому что пакетplm
, очевидно, был проблематичным дляggeffects
. - Вручную с использованием подхода оценки фиксированных эффектов «центрирования группы» с последующим запуском обычной
lm()
функции для данных, центрированных по группе. Запускggpredict()
на этих данных дает мне хорошие полосы уверенности, но значения X и Y затем центрируются по группе (сюрприз, сюрприз!), И я хочу, чтобы прогнозы и значения X были в исходной шкале, а не в центре среднего. шкала.
Любые предложения о том, как получить доверительные интервалы для оценки фиксированного эффекта, при использовании устойчивых стандартных ошибок кластера, были бы весьма признательны.