1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
|
\name{bgtest}
\alias{bgtest}
\alias{vcov.bgtest}
\alias{df.residual.bgtest}
\title{Breusch-Godfrey Test}
\description{
\code{bgtest} performs the Breusch-Godfrey test for higher-order
serial correlation.
}
\usage{
bgtest(formula, order = 1, order.by = NULL, type = c("Chisq", "F"),
data = list(), fill = 0)
}
\arguments{
\item{formula}{a symbolic description for the model to be tested
(or a fitted \code{"lm"} object).}
\item{order}{integer. maximal order of serial correlation to be tested.}
\item{order.by}{Either a vector \code{z} or a formula with a single explanatory
variable like \code{~ z}. The observations in the model
are ordered by the size of \code{z}. If set to \code{NULL} (the
default) the observations are assumed to be ordered (e.g., a
time series).}
\item{type}{the type of test statistic to be returned. Either
\code{"Chisq"} for the Chi-squared test statistic or
\code{"F"} for the F test statistic.}
\item{data}{an optional data frame containing the variables in the
model. By default the variables are taken from the environment
which \code{bgtest} is called from.}
\item{fill}{starting values for the lagged residuals in the auxiliary
regression. By default \code{0} but can also be set to \code{NA}.}
}
\details{
Under \eqn{H_0} the test statistic is asymptotically Chi-squared with
degrees of freedom as given in \code{parameter}.
If \code{type} is set to \code{"F"} the function returns
a finite sample version of the test statistic, employing an \eqn{F}
distribution with degrees of freedom as given in \code{parameter}.
By default, the starting values for the lagged residuals in the auxiliary
regression are chosen to be 0 (as in Godfrey 1978) but could also be
set to \code{NA} to omit them.
\code{bgtest} also returns the coefficients and estimated covariance
matrix from the auxiliary regression that includes the lagged residuals.
Hence, \code{\link{coeftest}} can be used to inspect the results. (Note,
however, that standard theory does not always apply to the standard errors
and t-statistics in this regression.)
}
\value{
A list with class \code{"bgtest"} inheriting from \code{"htest"} containing the
following components:
\item{statistic}{the value of the test statistic.}
\item{p.value}{the p-value of the test.}
\item{parameter}{degrees of freedom.}
\item{method}{a character string indicating what type of test was
performed.}
\item{data.name}{a character string giving the name(s) of the data.}
\item{coefficients}{coefficient estimates from the auxiliary regression.}
\item{vcov}{corresponding covariance matrix estimate.}
}
\references{
Breusch, T.S. (1978): Testing for Autocorrelation in Dynamic Linear
Models, \emph{Australian Economic Papers}, 17, 334-355.
Godfrey, L.G. (1978): Testing Against General Autoregressive and
Moving Average Error Models when the Regressors Include Lagged
Dependent Variables', \emph{Econometrica}, 46, 1293-1301.
Wooldridge, J.M. (2013): \emph{Introductory Econometrics: A Modern Approach},
5th edition, South-Western College.
}
\author{David Mitchell <david.mitchell@dotars.gov.au>, Achim Zeileis}
\seealso{\code{\link{dwtest}}}
\examples{
## Generate a stationary and an AR(1) series
x <- rep(c(1, -1), 50)
y1 <- 1 + x + rnorm(100)
## Perform Breusch-Godfrey test for first-order serial correlation:
bgtest(y1 ~ x)
## or for fourth-order serial correlation
bgtest(y1 ~ x, order = 4)
## Compare with Durbin-Watson test results:
dwtest(y1 ~ x)
y2 <- filter(y1, 0.5, method = "recursive")
bgtest(y2 ~ x)
bg4 <- bgtest(y2 ~ x, order = 4)
bg4
coeftest(bg4)
}
\keyword{htest}
|
