> 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…