| fxlm {fxregime} | R Documentation |
Interface to lm for fitting exchange rate regression
models (Frankel-Wei models).
fxlm(formula, data, ...)
formula |
a "formula" describing the linear model to be fit.
For details see below. |
data |
a "zoo" time series. |
... |
arguments passed to lm. |
fxlm is a function for fitting exchange rate regression models also
known as Frankel-Wei models. It is a simple convenience interface to lm:
data is assumed to be a "zoo" series in which, by default, the
first column is the dependent variable. If formula is omitted, the first
column is regressed on the remaining columns in data. The main difference
compared to plain lm models is that the error variance is reported as
a full parameter (estimated by maximum likelihood) in the coef method
and the estfun method (but currently not in the vcov method).
Furthermore, the index (also known as the time stamps) of the underlying data set
can be extracted by the time/index method.
An object of class "fxlm" inheriting from "lm".
Shah A., Zeileis A., Patnaik I. (2005), What is the New Chinese Currency Regime?, Report 23, Department of Statistics and Mathematics, Wirtschaftsuniversitaet Wien, Research Report Series, November 2005.
Zeileis A., Shah A., Patnaik I. (2008), Testing, Monitoring, and Dating Structural Changes in Maximum Likelihood Models, Report 70, Department of Statistics and Mathematics, Wirtschaftsuniversitaet Wien, Research Report Series, August 2008.
## load package and data
library("fxregime")
data("FXRatesCHF", package = "fxregime")
## compute returns for CNY (and explanatory currencies)
## for one year after abolishing fixed USD regime
cny <- fxreturns("CNY", frequency = "daily",
start = as.Date("2005-07-25"), end = as.Date("2006-07-24"),
other = c("USD", "JPY", "EUR", "GBP"))
## estimate full-sample exchange rate regression
fm <- fxlm(CNY ~ USD + JPY + EUR + GBP, data = cny)
coef(fm)
summary(fm)
## test parameter stability (with double max test)
scus <- gefp(fm, fit = NULL)
plot(scus, aggregate = FALSE)
## which shows a clear increase in the variance in March 2006
## alternative tests: Andrews' supLM ...
plot(scus, functional = supLM(0.1))
## ... or Nyblom-Hansen test (Cramer-von Mises type test)
plot(scus, functional = meanL2BB)