R/boa.gewekePwr.R

In brian-j-smith/boa: Bayesian Output Analysis Program (BOA) for MCMC

"boa.gewekePwr" <- function(link)
#-------------------------------------------------------------------------------
#-------------------------------------------------------------------------------
{
   spec <- NULL
   pnames <- boa.pnames(link)
   if(nrow(link) > 10) {
      for(i in pnames) {
         spec <- c(spec, spectrum0(link[, i])$spec)
      }
   } else {
      spec <- rep(NA, length(pnames))
   }

   return(structure(spec, names = pnames))
}