> d4 <- d7[d7$x == 4,] #葉数4のサブセット > table(d4$y) #生存数がyiだった個体をカウント 0 1 2 3 4 5 6 7 8 3 1 4 2 1 1 2 3 3 > c(mean(d4$y), var(d4$y)) #平均と分散 [1] 4.050000 8.365789 > library(glmmML) > glmmML(cbind(y, N - y) ~ x, data =…

いろいろ参照してます http://qh73xe.jimdo.com/r%E3%81%AE%E5%9F%BA%E6%9C%AC/r%E3%81%AE%E4%BE%BF%E5%88%A9%E3%81%AA%E9%96%A2%E9%80%A3%E3%82%BD%E3%83%95%E3%83%88/rtools/ http://tjo.hatenablog.com/entry/2013/10/04/190701 https://github.com/stan-…