Ｒでロジスティック回帰は割と楽にできる。
x1 y1 y2
3 4 6
10 9 4
8 7 10
3 1 4
6 8 10
9 5 4
g<-glm(cbind(y1,y2)~x,binomial,data)
summary(g)
独立変数が離散でmulti-categoriesのときもほぼ同様(...~x-1で切片なしモデルにするのが普通なのかしら。そうしているのをよく見る）。なんだが、その後にカテゴリー間の比較をしようとするとAICを使ってモデル比較ということになるようだ。これがけっこうメンドい。
参考サイト：久保さんのところの「多重検定と多重比較をめぐって」あたり。
http://hosho.ees.hokudai.ac.jp/~kubo/ce/index.html
まさにそのもののデータなんだが・・・。とりあえず「比率」→「変数変換」→「ANOVA」というconventionalな方法でごまかそうという気になっている。みんなそうしてるしね。いや、それじゃいけないのかな。まあともかくもう眠い。
