GLMの基礎

R による階層モデルと混合効果モデル

Richard Erickson

Data Scientist

正規性の仮定

  • 残差は「正規」である
  • 変換
  • 「アークサインは無意味:...」
  • 代替分布
R による階層モデルと混合効果モデル

RにおけるGLMの構文

glm(y ~ x, family = "gaussian")
?family
R による階層モデルと混合効果モデル

λ=3のポアソン分布の回帰プロット例

R による階層モデルと混合効果モデル

ポアソン回帰の例

glm(y ~ x, family = "poisson")
R による階層モデルと混合効果モデル

ロジスティック回帰のプロット例。

R による階層モデルと混合効果モデル

ロジスティック回帰の例

# binary: y = 0 or 1
glm(y ~ x, family = "binomial")

# Wilkinson-Rogers: cbind(success, failure) glm(cbind(success, failure) ~ x, family = "binomial")
# Weighted format: y = 0.3, weights = 10 glm(y ~ x, weights = weights, family = "binomial")
R による階層モデルと混合効果モデル

練習しましょう!

R による階層モデルと混合効果モデル

Preparing Video For Download...