Ramsey's RESET test for functional form.
resettest(formula, power = 2:3, type = c("fitted", "regressor",
"princomp"), data = list(), vcov = NULL, …)
a symbolic description for the model to be tested
(or a fitted "lm"
object).
integers. A vector of positive integers indicating the powers of the variables that should be included. By default, the test is for quadratic or cubic influence of the fitted response.
a string indicating whether powers of the fitted response, the regressor variables (factors are left out), or the first principal component of the regressor matrix should be included in the extended model.
an optional data frame containing the variables in the model.
By default the variables are taken from the environment which resettest
is called from.
optional arguments to be passed to waldtest
for carrying out the F test.
An object of class "htest"
containing:
the test statistic.
the corresponding p-value.
degrees of freedom.
a character string with the method used.
a character string with the data name.
The RESET test is a popular diagnostic for correctness of functional
form. The basic assumption is that under the alternative the model can be
written in the form
\( y = X\beta + Z\gamma + u\).
Z
is generated by taking powers either of the fitted response, the
regressor variables, or the first principal component of X
. A standard
F-Test is then applied to determine whether these additional variables have
significant influence. The test statistic under \(H_0\) follows an F
distribution with parameter
degrees of freedom.
This function was called reset
in previous versions of the package. Please
use resettest
instead.
Examples can not only be found on this page, but also on the help pages of the
data sets bondyield
, currencysubstitution
,
growthofmoney
, moneydemand
,
unemployment
, wages
.
J.B. Ramsey (1969), Tests for Specification Errors in Classical Linear Least-Squares Regression Analysis. Journal of the Royal Statistical Society, Series B 31, 350--371
W. Kr<e4>mer & H. Sonnberger (1986), The Linear Regression Model under Test. Heidelberg: Physica
# NOT RUN {
x <- c(1:30)
y1 <- 1 + x + x^2 + rnorm(30)
y2 <- 1 + x + rnorm(30)
resettest(y1 ~ x, power=2, type="regressor")
resettest(y2 ~ x, power=2, type="regressor")
# }
Run the code above in your browser using DataLab