RでのGLM備忘録

エクセルファイルのデータをRで読む込み、GLMで分析する。Stepwise regression法でモデル選択する。
1. エクセルファイルをcsvファイルにする。
2. Rでcsvファイルを読み込む。
> z<-read.csv("〜.csv")
3. カテゴリ変数がある場合は、GLMで分析する前に、as.factor()でカテゴリ変数とする。
> z$xn<-as.factor(z$xn)
4. GLMで分析する。今回は0/1データを連続変数で説明する。二項分布でモデルがlogitの場合は、
> result1=glm(y~x1+x2+x3+...xn, family=binomial(link="logit"), data=z)
logitではなくprobitにする場合は、linkのところをprobitにするだけ。
5. GLMMの場合は、clusterでrandom effectを指定する。
> result2=glmmML(y~x1+x2+x3+...xn, family=binomial(link="logit"), data=z, cluster=xm)
6. Stepwise regression法でモデル選択をする場合は、
>library(MASS)
>stepAIC(result1, direction="both")
デフォルトは後ろ向き。directionで方向を決められる。