| cointRegFM {cointReg} | R Documentation |
Fully Modified OLS
Description
Computes the Phillips and Hansen (1990) Fully Modified OLS estimator.
Usage
cointRegFM(x, y, deter, kernel = c("ba", "pa", "qs", "tr"),
bandwidth = c("and", "nw"), demeaning = FALSE, check = TRUE, ...)
Arguments
x |
[ |
y |
[ |
deter |
[ |
kernel |
[ |
bandwidth |
[ |
demeaning |
[ |
check |
[ |
... |
Arguments passed to |
Details
The equation for which the FM-OLS estimator is calculated:
y = \delta \cdot D + \beta \cdot x + u
with D as the deterministics matrix.
Then \theta = (\delta', \beta')' is the full parameter vector.
The calculation of t-values and the variance-covariance matrix is only
possible, if y is one-dimensional.
Value
[cointReg]. List with components:
delta[numeric|matrix]-
coefficients as vector / matrix
beta[numeric|matrix]-
coefficients as vector / matrix
theta[numeric|matrix]-
combined coefficients of
betaanddeltaas vector / matrix sd.theta[numeric]-
standard errors for
theta t.theta[numeric]-
t-values for
theta p.theta[numeric]-
p-values for
theta residuals[numeric]-
FM-OLS residuals (first value is always missing)
omega.u.v[numeric]-
conditional long-run variance based on OLS residuals.
varmat[matrix]-
variance-covariance matrix
Omega[list]-
the whole long-run variance matrix and parts of it
beta.OLS[numeric|matrix]-
OLS coefficients as vector / matrix
delta.OLS[numeric|matrix]-
OLS coefficients as vector / matrix
u.OLS[numeric]-
OLS residuals
bandwidth[list]-
numberandnameof bandwidth kernel[character]-
abbr. name of kernel type
References
Phillips, P.C.B. and B. Hansen (1990): "Statistical Inference in Instrumental Variables Regression with I(1) Processes," Review of Economic Studies, 57, 99–125, DOI:10.2307/2297545.
See Also
Other cointReg: cointRegD,
cointRegIM, cointReg,
plot.cointReg, print.cointReg
Examples
set.seed(1909)
x1 = cumsum(rnorm(100, mean = 0.05, sd = 0.1))
x2 = cumsum(rnorm(100, sd = 0.1)) + 1
x3 = cumsum(rnorm(100, sd = 0.2)) + 2
x = cbind(x1, x2, x3)
y = x1 + x2 + x3 + rnorm(100, sd = 0.2) + 1
deter = cbind(level = 1, trend = 1:100)
test = cointRegFM(x, y, deter, kernel = "ba", bandwidth = "and")
print(test)