On Wed, 15 Sep 2010, Claudio wrote:
This is a stupid question: how to perform a Granger-Sims
causality in R?
I presume you mean, in gretl.
The GS causality test is like the Granger causality, but
includes not only lags, but also leads of the variable...
As Sims (1972) puts it, "If and only if causality runs one way
from current and past values of some list of exogenous variables
to a given endogenous variable, then in a regression of the
endogenous variable on past, current and future values of the
exogenous variables, the future values of the exogenous variables
should have zero coefficients."
You can't specify leads (negative lags) with gretl's VAR command,
but VARs are just estimated via OLS so you can easily do this in
gretl. Here, for example, is a Sims test for one-way causality of
US car sales (QNC) by US consumers' income (INCOME):
open data9-7
list Xlags = INCOME(-1 to -4)
list Xleads = INCOME(+1 to +4)
ols QNC 0 Xlags INCOME Xleads
omit Xleads --quiet
"omit" gives F(4,46) = 0.263 with p-value 0.90, so the hypothesis
of one-way causality (all the leads of INCOME have zero
coefficients) is not rejected.
Allin Cottrell