portes (version 1.04)

Hosking: The Modified Multivariate Portmanteau Test, Hosking (1980)

Description

The modified multivariate portmanteau test suggested by Hosking (1980).

Usage

Hosking(obj,lags=seq(5,30,5),order=0,SquaredQ=FALSE)

Arguments

obj

a fitted time-series model with any of the classes "ar", "Arima", "arima0", "FitAR", or "FitFGN" or a residuals with class "ts", "numeric", "matrix"<

lags

vector of lag values.

order

$p+q$, where ($p$,$q$) is the ARMA or the VARMA order.

SquaredQ

TRUE, apply the test to the squared values. This is needed to check the Autoregressive Conditional Heteroscedastic, ARCH, effects. When SquaredQ = FALSE, apply for the usual residuals.

Value

• The multivariate test statistic suggested by Hosking (1980) and its associated p-values for different lags based on asymptotic chi-square distribution with k^2(lags-order) degrees of freedom.

References

Hosking, J. R. M. (1980). "The Multivariate Portmanteau Statistic". Journal of American Statistical Association, 75, 602-608.

See Also

acf, Box.test, BoxPierce, LjungBox, LiMcLeod, GVStat, portes, Get.Resid

Examples

##############################################################
## Quarterly, west German investment, income, and consumption 
## from first quarter of 1960 to fourth quarter of 1982: 
##############################################################
data(WestGerman)
DiffData <- matrix(numeric(3 * 91), ncol = 3)
  for (i in 1:3) 
    DiffData[, i] <- diff(log(WestGerman[, i]), lag = 1)
fit <- ar.ols(DiffData, intercept = TRUE, order.max = 2)
lags <- c(5,10)
## Apply the test statistic on the fitted model 
Hosking(fit,lags,order = 2)  
## Apply the test statistic on the residuals
res <- ts((fit$resid)[-(1:2), ])
Hosking(res,lags,order = 2)