R/morisita.horn.R

Defines functions `morisita.horn`

`morisita.horn` <-
function(x,y) {
  aN <- sum(x)
  bN <- sum(y)
  da <- sum(x^2)/aN^2
  db <- sum(y^2)/bN^2
  top <- sum(x*y)
  mor <- 2*top/((da+db)*aN*bN)
  return(mor)
}

Try the fossil package in your browser

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

fossil documentation built on March 23, 2020, 5:06 p.m.