新・医学と統計(23)

前回の Log-linear Regression を次の「R」プログラムで実行してみましょう。

 

前回用いた、図１の Excel ファイル「Frequency.csv」 の 「A1：C5」を選択しグリップボードにコピーしておき、以下を実行します。

 

***
# R console に記入し実行
dat<- read.delim("clipboard", header=T) 

 

# 以下、Rエディタに記入し実行
dat
library(MASS)

# （１）相互作用項を含むモデルの予測
TableLM1<-loglm(Count~ Gender + TC + Gender*TC, data=dat)
TableLM1
fitted(TableLM1)

 

# （２）相互作用項を含まないモデルの予測
TableLM2<-loglm(Count~ Gender + TC, data=dat)
TableLM2
fitted(TableLM2)

 

出力結果：
図３　予測の結果　

 

 

以上から、相互作用項を含むモデルの当てはめの良さがわかる。

なお、

「関数 loglin」を用いるときは度数による分割表を用います。

 

次回に続く！

 

情報統計研究所をご気軽にご利用くださいませ。

 

情報統計研究所刊：「すぐに役立つ統計のコツ」立読み

↓
https://e-hon.cloudpages.jp/viewer/asp/9784_274_218163

 

スタットの医学看護統計-blog