2014-08-01から1ヶ月間の記事一覧
d$y.rnd <- rpois(100, lambda = mean(d$y)) # 新のモデルから生成 mean(d$y) # 平均 fit1 <- glm(y.rnd ~ 1, data = d, family = poisson) fit2 <- glm(y.rnd ~ x, data = d, family = poisson) >fit1$deviance - fit2$deviance #逸脱度の差 [1] 1.273795 g…
d$y.rnd <- rpois(100, lambda = mean(d$y)) # 新のモデルから生成 mean(d$y) # 平均 fit1 <- glm(y.rnd ~ 1, data = d, family = poisson) fit2 <- glm(y.rnd ~ x, data = d, family = poisson) >fit1$deviance - fit2$deviance #逸脱度の差 [1] 1.273795 g…