# 第51回日本生態学会大会 (釧路) の自由集会での久保さんの話題提供を実行する手順
# 資料は下のページ
# http://hosho.ees.hokudai.ac.jp/~kubo/ce/2004/kubo/kuboJES2004.pdf

##	データの準備
#各グループの名称を設定
level <- c(rep("A",20), rep("B",10), rep("C",20))

#グループ分けを設定
level.mapped1 <- c(rep("A+B",30), rep("C",20))
level.mapped2 <- c(rep("A+C",20), rep("B",10), rep("A+C",20))
level.mapped3 <- c(rep("A",20), rep("B+C",30))
level.mapped4 <- c(rep("A",20), rep("B",10), rep("C",20))

#文字列を因子オブジェクトに変換
level <- factor(level)
level.mapped1 <- factor(level.mapped1)
level.mapped2 <- factor(level.mapped2)
level.mapped3 <- factor(level.mapped3)
level.mapped4 <- factor(level.mapped4)

#ポワソン分布の乱数発生
n.flower <- c(rpois(20,2.5), rpois(10,3.5), rpois(20,3.5))


##	データの解析
result0 <- glm(n.flower ~ 1, family = poisson (link = log))
result1 <- glm(n.flower ~ level.mapped1 - 1, family = poisson (link = log))
result2 <- glm(n.flower ~ level.mapped2 - 1, family = poisson (link = log))
result3 <- glm(n.flower ~ level.mapped3 - 1, family = poisson (link = log))
result4 <- glm(n.flower ~ level.mapped4 - 1, family = poisson (link = log))


##	結果の表示
result0$aic
result1$aic
result2$aic
result3$aic
result4$aic