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