理学療法士がまとめたノート

読者です 読者をやめる 読者になる 読者になる

統計学備忘録 since2016

理学療法、医療、福祉、小児

Rで簡単(二元配置分散分析:対応なし)

統計学

サンプルをコピペで2式を実行

方法は Rで簡単(一元配置分散分析:対応なし)を参照

 

検定は帰無仮説が大切です!

商品のの主効果:帰無仮説=カップ麺でもインスタント麺でも評価の母平均は等しい

スープの主効果:帰無仮説=スープが違っても評価の母平均は等しい

商品スープの主効果:帰無仮説=商品とスープの組合せと評価は関係ない

   f:id:yoshida931:20161118113459p:plain

 

商品 スープ 評価
カップ とんこつ 10
カップ とんこつ 11
カップ とんこつ 11
カップ とんこつ 9
カップ とんこつ 9
カップ 醤油 11
カップ 醤油 8
カップ 醤油 10
カップ 醤油 8
カップ 醤油 8
カップ みそ 7
カップ みそ 3
カップ みそ 5
カップ みそ 2
カップ みそ 3
インスタント とんこつ 11
インスタント とんこつ 11
インスタント とんこつ 10
インスタント とんこつ 9
インスタント とんこつ 10
インスタント 醤油 9
インスタント 醤油 7
インスタント 醤油 9
インスタント 醤油 8
インスタント 醤油 8
インスタント みそ 4
インスタント みそ 4
インスタント みそ 2
インスタント みそ 3
インスタント みそ 3

 

x<-read.table("clipboard",header=T)
summary(aov(x$評価~x$商品*x$スープ))

 

f:id:yoshida931:20161118121202p:plain

 

スープによる主効果あり!

 

次の式でも同じ結果になります

summary(aov(x$評価~x$商品+x$スープ+x$商品:x$スープ))

 

ベクトルを使用した方法

評価<-c(x$評価)

商品<-factor(c(rep("カップ",15),rep("インスタント",15)))

スープ<-factor(rep(c(rep("とんこつ",5),rep("醤油",5),rep("みそ",5)),2))

summary(aov(x$評価~x$商品*x$スープ))

 

検証するために図を描きましょう!

interraction.plot(商品,スープ,評価)

interraction.plot(スープ,商品,評価)

 

f:id:yoshida931:20161118121118p:plain

f:id:yoshida931:20161118121138p:plain