我正在从Stata迁移到R(plm包
),以进行面板模型计量经济学。在Stata中,面板模型(如随机效应)通常报告内部、外部和总体的R平方。
我发现,plm
随机效应模型中报告的R平方对应于内部R平方。那么,在R中使用plm包
是否有办法获得总体和外部的R平方?
以下是使用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)
在Stata中:
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平方值(0.18062)与Stata报告的R平方值(0.1799)类似。是否有办法在R中获得Stata报告的R平方值之间(0.1860)和总体(0.1830)的R平方值?