- （１つの飼育容器に10匹入れて選好性試験：AとBのどっちに行くか）×20反復
- 母比率=0.5を検定したい
- 切片だけのロジスティック回帰？（もちろん各個体は他の個体と独立に選ぶという前提で）
# ロジスティック関数・切片だけモデル・シミュレーション
# 関数定義：test(p=母比率) 省略時 test() は0.5
test <- function(p=0.5){
x <- rbinom(20, 10, p) # ２項分布に従うランダムを20反復取り出す
g <- glm ( cbind(x, 10 - x) ~ 1, binomial)
return (summary(g))
}test()
# パラメータ（切片）の分布
# （ひどいミス発見(汗）修正。2012/11/15）alpha <- NULL
for(i in 1:1000) alpha <- c(alpha, test()$coef[1])
hist(alpha)
後日、よく見たら係数をリストにして関数に吐き出させればいいのか。相変わらず。
