plm package
パネルモデルの計量経済学を行うために、Stata から R ( ) に移行しています。Stata では、変量効果などのパネル モデルは通常、内、間、および全体の R 二乗を報告します。
私はplm
、ランダム効果モデルで報告されたR 二乗が R 二乗内に対応することを発見しました。では、R を使用して R-squared 全体と R-squared を取得する方法はありますplm package
か?
R と Stata で同じ例を参照してください。
library(plm)
library(foreign) # read Stata files
download.file('http://fmwww.bc.edu/ec-p/data/wooldridge/wagepan.dta','wagepan.dta',mode="wb")
wagepan <- read.dta('wagepan.dta')
# Random effects
plm.re <- plm(lwage ~ educ + black + hisp + exper + expersq + married + union + d81 + d82 + d83 + d84 + d85 + d86 + d87,
data=wagepan,
model='random',
index=c('nr','year'))
summary(plm.re)
スタタで:
use http://fmwww.bc.edu/ec-p/data/wooldridge/wagepan.dta
xtset nr year
xtreg lwage educ black hisp exper expersq married union d81 d82 d83 d84 d85 d86 d87, re
R で報告された R-sq (0.18062) は、少なくともこの場合、Stata で報告された R-sq Within (0.1799) に似ています。Stata で報告された (0.1860) と全体 (0.1830) の間の R 平方を R で取得する方法はありますか?