# NOT RUN {
#### Example 10-14
## ------------------------------------------------------------------------
# }
# NOT RUN {
data("EvapoTransp", package = "pder")
library("splm")
data("etw", package = "pder")
evapo <- et ~ prec + meansmd + potet + infil + biomass + plantcover +
softforbs + tallgrass + diversity + matgram + dwarfshrubs + legumes
semsr.evapo <- spreml(evapo, data=EvapoTransp, w=etw,
lag=FALSE, errors="semsr")
summary(semsr.evapo)
## ------------------------------------------------------------------------
library("lmtest")
library("plm")
coeftest(plm(evapo, EvapoTransp, model="pooling"))
## ------------------------------------------------------------------------
coeftest(spreml(evapo, EvapoTransp, w=etw, errors="sem"))
#### Example 10-17
## ------------------------------------------------------------------------
saremsrre.evapo <- spreml(evapo, data = EvapoTransp,
w = etw, lag = TRUE, errors = "semsr")
summary(saremsrre.evapo)$ARCoefTable
round(summary(saremsrre.evapo)$ErrCompTable, 6)
# }
Run the code above in your browser using DataCamp Workspace