1coef.ecoNP <- function(object, subset = NULL, obs = NULL, ...) {
2  mu <- object$mu
3  if (is.null(subset))
4    subset <- 1:nrow(mu)
5  else if (max(subset) > nrow(mu))
6    stop(paste("invalid input for `subset.' only", nrow(mu), "draws are stored."))
7
8  if (is.null(obs))
9    obs <- 1:dim(object$mu)[3]
10  else if (max(subset) > dim(object$mu)[3])
11    stop(paste("invalid input for `obs.' only", dim(object$mu)[3], "draws are stored."))
12
13  return(mu[subset,,obs])
14}
15