目的

Morisita の類似度指数 Cλ を計算する

使用法

Morisita(a, b)

引数

a     第一群の観察度数ベクトル
b     第二群の観察度数ベクトル

ソース

インストールは，以下の 1 行をコピーし，R コンソールにペーストする
source("http://aoki2.si.gunma-u.ac.jp/R/src/Morisita.R", encoding="euc-jp")

# Morisita の類似度指数 Cλ を計算する
Morisita <- function(        a,                      # 第一群の観察度数ベクトル
                        b)                      # 第二群の観察度数ベクトル
{
        stopifnot(length(a) == length(b))
        na <- sum(a)
        nb <- sum(b)
        return(2*sum(a*b)/((sum(a*(a-1)/na/(na-1))+sum(b*(b-1)/nb/(nb-1)))*na*nb))
}


使用例

> Morisita(c(3,2,3,4,3,4), c(2,6,4,3,2,1))
[1] 1.018433