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)

}

Try the poolVIM package in your browser

Any scripts or data that you put into this service are public.

poolVIM documentation built on May 1, 2019, 9:13 p.m.