Dear R users,
I have a question here
library(AER)
library(plm)
library(sandwich)
## take the following data
data("Gasoline", package="plm")
Gasoline$f.year=as.factor(Gasoline$year)
Now I run the following regression
rhs <- "-1 + f.year + lincomep+lrpmg+lcarpcap"
m1<- lm(as.formula(paste("lgaspcar ~", rhs)), data=Gasoline)
###Now I want to find the autocorrelation,heteroskedasticity adjusted
standard errors as a part of coeftest
### Basically I would like to take care of the within country serial
correlaion
###that is I want to do
coeftest(m1, vcov=function(x) vcovHAC(x,order.by=...))
Please suggest what should be the argument of order.by and whether that will
give me the desired result
[[alternative HTML version deleted]]
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
