Displaying 1 result from an estimated 1 matches for "gr_lme".
Did you mean:
gr_lm
2010 Apr 09
0
panel regression with twoways random effects, on unbalanced data?
...ple,
library("AER")
library("plm")
data("Grunfeld", package = "AER")
pgr <- plm.data(Grunfeld, index = c("firm", "year"))
gr_re <- plm(invest ~ value + capital, data = pgr,
model = "random", effect="individual")
gr_lme <- lme(invest ~ value + capital, Grunfeld, random=~1|firm)
summary(gr_re)
summary(gr_lme)
For this "individual" RE model, both plm() and lme() yield apparently
similar results. While for the "twoways" one not quite:
gr_re <- plm(invest ~ value + capital, data = pgr,
mod...