1"summary.varest" <-
2function(object, equations = NULL, ...){
3  ynames <- colnames(object$y)
4  obs <- nrow(object$datamat)
5  if (is.null(equations)) {
6    names <- ynames
7  }
8  else {
9    names <- as.character(equations)
10    if (!(all(names %in% ynames))) {
11      warning("\nInvalid variable name(s) supplied, using first variable.\n")
12      names <- ynames[1]
13    }
14  }
15  eqest <- lapply(object$varresult[names], summary)
16  resids <- resid(object)
17  covres <- cov(resids) * (obs - 1) / (obs - (ncol(object$datamat) - object$K))
18  corres <- cor(resids)
19  logLik <- as.numeric(logLik(object))
20  roots <- roots(object)
21  result <- list(names = names, varresult = eqest, covres = covres, corres = corres, logLik = logLik, obs = obs, roots = roots, type = object$type, call = object$call)
22  class(result) <- "varsum"
23  return(result)
24}
25
26