「統計学関連なんでもあり」の過去ログ--- 045

No.17058　式について　　【学生】　2012/06/12(Tue) 08:53

Rで以下の事象を解析しようと思っています。
式が正しいかどうか，確認させてください。

ある事象について，1個体につき，N回観察し，その中で反応した回数がα回，反応しなかった回数がN-α回でした。
ある処理をしたグループAと処理をしなかったグループBで各10個体ずつデータをとりました。
観察回数Nは個体によって異なります。
ある処理によって，反応確率が変わるのか（変わらないのか）について調べたいと思っております。以下の式を指定することで，正しく解析することができますでしょうか？

glm(formula=cbind(α, N-α）~group, family=binomial)

よろしくお願い致します。

No.17059　Re: 式について　　【青木繁伸】　2012/06/12(Tue) 11:42

まず，そのプログラムでは，エラーも結果も出ないようですよ。
> n <- sample(5:15, 20, replace=TRUE)
> m <- sapply(n, function(m) sample(1:(m-1), 1))
> group <- factor(rep(c("A", "B"), each=10))
> glm(formula=cbind(n, m）~group, family=binomial)