Я запускаю несбалансированную панель регрессии с фиксированными эффектами в R с использованием пакета plm:
model_fe = plm(yield_diff ~ Bid_ask, data = mydata.set, model = "within")
summary(model_fe)
Меня особенно интересует общее распределение ненаблюдаемых отдельных эффектов, то есть перехватов.
Я извлек отдельные эффекты, используя следующее:
fixef(model_fe)
summary(fixef(model_fe))
Возвращает перехват (отдельные фиксированные эффекты) для каждого из отдельных. Однако меня интересует общее распределение этих эффектов (квартили, среднее значение, медиана, максимум, мин)
Мне кажется, мне удается получить среднее значение как:
weighted.mean(fixef(model_fe), as.numeric(table(index(model_fe)[[1]])))
Есть ли способ получения другой сводной статистики?