# NOT RUN {
#### Example 6-1
## ------------------------------------------------------------------------
library("plm")
## ------------------------------------------------------------------------
y ~ x1 + x2 + x3 | x1 + x3 + z
y ~ x1 + x2 + x3 | . - x2 + z
## ------------------------------------------------------------------------
data("SeatBelt", package = "pder")
SeatBelt$occfat <- with(SeatBelt, log(farsocc / (vmtrural + vmturban)))
ols <- plm(occfat ~ log(usage) + log(percapin) + log(unemp) + log(meanage) +
log(precentb) + log(precenth)+ log(densrur) +
log(densurb) + log(viopcap) + log(proppcap) +
log(vmtrural) + log(vmturban) + log(fueltax) +
lim65 + lim70p + mlda21 + bac08, SeatBelt,
effect = "time")
fe <- update(ols, effect = "twoways")
ivfe <- update(fe, . ~ . | . - log(usage) + ds + dp +dsp)
rbind(ols = coef(summary(ols))[1,],
fe = coef(summary(fe))[1, ],
w2sls = coef(summary(ivfe))[1, ])
## ------------------------------------------------------------------------
SeatBelt$noccfat <- with(SeatBelt, log(farsnocc / (vmtrural + vmturban)))
nivfe <- update(ivfe, noccfat ~ . | .)
coef(summary(nivfe))[1, ]
# }
Run the code above in your browser using DataCamp Workspace