我正在使用plm软件包在R中运行固定效果的不平衡面板回归：
It works!我对未观察到的单个效果（即截距）的整体分布特别感兴趣。
我使用以下方法提取了各个效果：
model_fe = plm(yield_diff ~ Bid_ask, data = mydata.set, model = "within")
summary(model_fe)
这将返回每个人的截距（单独的固定效果）。但是，我对这些影响的总体分布感兴趣（四分位数，均值，中位数，最大值，最小值）
我似乎设法获得均值：
fixef(model_fe)
summary(fixef(model_fe))
是否可以获取其他摘要统计信息？