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