1print.eco <- function(x, digits = max(3, getOption("digits") -3), 2 ...){ 3 cat("\nCall:\n", deparse(x$call), "\n\n", sep="") 4 5 if (is.null(x$N)) 6 N <- rep(1, nrow(x$X)) 7 else N <- x$N 8 9 W.mean <- cbind(mean(x$W[,1,] %*% (x$X*N/sum(x$X*N))), 10 mean(x$W[,2,] %*% ((1-x$X)*N/sum((1-x$X)*N)))) 11 colnames(W.mean) <- c("W1", "W2") 12 rownames(W.mean) <- "posterior mean" 13 14 cat("Aggregate In-sample Estimates:\n\n") 15 print.default(format(W.mean, digits = digits), print.gap = 2, quote = 16 FALSE) 17 cat("\n") 18 invisible(x) 19} 20