国固定効果

Question

国の固定効果をカントリーダミーで見積もりしようとしています。

fe1b <- plm( 
    bond_GDP_local ~ real_r + equity_volatility, model = 'within', data = panel_eme_filtered 
) 

これは私の1以下と同じ係数を与える：

fe1bc <- plm( 
    bond_GDP_local ~ real_r + equity_volatility +country_code, model = 'within', data = panel_eme_filtered 
) 

私は私の式に国ダミーを入力していますが、私は結果でそれを見ることはできません。 最初のモデルにはすでに組み込まれていますか？

Oneway (individual) effect Within Model 

Call: 
plm(formula = bond_GDP_local ~ real_r + equity_volatility, data = panel_eme_filtered, 
    model = "within") 

Balanced Panel: n=8, T=60, N=480 

Residuals : 
    Min. 1st Qu. Median 3rd Qu. Max. 
-2.7200 -0.3450 -0.0927 0.2200 5.6200 

Coefficients : 
         Estimate Std. Error t-value Pr(>|t|) 
real_r   -0.0331088985 0.0171886368 -1.926 0.0547 . 
equity_volatility -0.0000003838 0.0000006396 -0.600 0.5488 
--- 
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 

Total Sum of Squares: 345.7 
Residual Sum of Squares: 342.7 
R-Squared:  0.008731 
Adj. R-Squared: -0.01025 
F-statistic: 2.06979 on 2 and 470 DF, p-value: 0.1274 

別の質問：

あなたは

それらの両方が私にこれを与えている感謝どのようにこのモデルで堅牢なパネルの時系列データの標準誤差を推定することができますか？

Answer 1

おそらくpanel_eme_filteredは、country_codeで索引付けされたpdata.frameですか？その場合、回帰式にcountry_codeを含めることは重要ではありません。これを行う別の方法はlfeです。

library(lfe) 

fe2 <- felm(
    bond_GDP_local ~ real_r + equity_volatility | country_code, 
    data = panel_eme_filtered 
) 
summary(fe2, robust = T) # heteroskedastic robust SE's 

また

fe3 <- felm(
    bond_GDP_local ~ real_r + equity_volatility | country_code | 0 | country_code, 
    data = panel_eme_filtered 
) 
summary(fe3)