アラカシ

アラカシ Quercus glauca 2005.4.17 京都府京都市

コバノミツバツツジ

コバノミツバツツジ Rhododendron reticulatum 2005.4.17 京都府京都市

Mac OS X 10.3.9

Mac OS X 10.3.9アップデートが配布されている。しかしあと13日でTiger発売。

クサイチゴ

クサイチゴ Rubus hirsutus 2005.4.14 京都府京都市

今日のjeconet

jeconetに会員名簿廃止と規程改定の結果を流す。

今日は、人工洞穴を生息地とするコウモリの保全が話題となっていた。

(セイヨウ)カラシナ

カラシナ(セイヨウカラシナ) Brassica juncea 2005.4.9 京都府京都市

Jeconet管理人の作業

Jeconetのページを更新した。会員名簿廃止というのは、管理人の日々の作業量を減らすという効果もある(はず)。もっとも、三中さんの指摘するような問題もあるわけだが。

アオサギと松本酒造

アオサギ(幼鳥か)。

松本酒造の酒蔵。昨日は、カメラを持った人が大勢おりました。そのうちの一人だったわけですが。

タグ：鳥

類似度指数・重複度指数を計算するR関数(改訂)

前回のは、forでループを回したりするのがRらしくなかったので、書き直してみる。

MorisitaのCλ

clambda <- function(c1, c2) {
	if (length(c1) != length(c2)) stop("different vector lengths")
	n1 <- sum(c1)
	n2 <- sum(c2)
	l1 <- sum(c1 * (c1 - 1)) / (n1 * (n1 - 1))
	l2 <- sum(c2 * (c2 - 1)) / (n2 * (n2 - 1))
	cc <- sum(c1 * c2)
	return (2 * cc / ((l1 + l2) * n1 * n2))
}

Piankaのα

pianka.alpha <- function(c1, c2) {
	if (length(c1) != length(c2)) stop("different vector lengths")
	p1 <- sum(c1 * c1)
	p2 <- sum(c2 * c2)
	pp <- sum(c1 * c2)
	return (pp / (sqrt(p1) * sqrt(p2)))
}

サクラ

2005-04-08 職場のソメイヨシノ