R/m_effective.R

#'  meff.
#' @param R R
#'
m_effective <- function(R) {

  evs <- eigen(R)$values
    abs.evs <- abs(evs)
  k <- length(evs)
    m <- 1 + (k - 1) * (1 - var(evs) / k)

  m <- floor(m)

  return(m)

}
snembrini/poolVIM documentation built on May 4, 2019, 3:10 a.m.