vcov.WCE <- function(object, ...){ # obtain variance covariance matrix from WCE object
  object$vcovmat
}