> fit <- glm(y ~ x, data = d, family = poisson) > fit # プリントでもいい Call: glm(formula = y ~ x, family = poisson, data = d) Coefficients: (Intercept) x 1.29172 0.07566 Degrees of Freedom: 99 Total (i.e. Null); 98 Residual Null Deviance…

> d <- read.csv("data3a.csv") > d y x f 1 6 8.31 C 2 6 9.44 C 3 6 9.50 C # 略 99 7 10.86 T 100 9 9.97 T > d$x # xだけ表示 [1] 8.31 9.44 9.50 9.07 10.16 8.32 10.61 10.06 9.93 10.43 10.36 10.15 [13] 10.92 8.85 9.42 11.11 8.02 11.93 8.55 7.19…

> load("data.RData") # データを読み込む > data # データの確認 [1] 2 2 4 6 4 5 2 3 1 2 0 4 3 3 3 3 4 2 7 2 4 3 3 3 4 3 7 5 3 1 7 6 4 6 5 2 [37] 4 7 2 2 6 2 4 5 4 5 1 3 2 3 > length(data) #データの個数確認 [1] 50 > summary(data) #データの基…

・線形モデルLMは等分散正規分布を仮定 そうでないものは一般線形モデルGLMデータ解析のための統計モデリング入門――一般化線形モデル・階層ベイズモデル・MCMC (確率と情報の科学)作者: 久保拓弥出版社/メーカー: 岩波書店発売日: 2012/05/19メディア: 単行…