共分散分析を使えば、2つの回帰直線の傾きや切片が異なるかどうか検討できる。
直接的な傾きの差の検定や切片の差の検定ではないが、同様な結果が簡単に得られるので、共分散分析でよいと思っている。
データの準備
こちらを題材として、EZRで再解析してみる。
https://toukeier.hatenablog.com/entry/2019/01/06/120845toukeier.hatenablog.com
データはこのように準備する。
Rスクリプトはこちら。
NO2.c <- c(0.025,0.030,0.016,0.028,0.021,0.041,0.012,0.015,0.030,0.018,0.023) pre.c <- c(4.1,6.8,2.6,11.5,2.8,10.7,2.7,2.8,3.5,5.0,3.6) are.c <- rep("c",length(NO2.c)) NO2.o <- c(0.024,0.033,0.030,0.017,0.022,0.027,0.016,0.017,0.020,0.022) pre.o <- c(6.8,8.7,11.4,4.5,5.6,10.5,5.2,6.8,6.0,8.0) are.o <- rep("o",length(NO2.o)) cough.prev.c <- data.frame(NO2=NO2.c,pre=pre.c,are=are.c) cough.prev.o <- data.frame(NO2=NO2.o,pre=pre.o,are=are.o) cough.prev <- rbind(cough.prev.c,cough.prev.o)
散布図と回帰直線を描く
層別の散布図を描かせてみる。
「グラフと表」から「散布図」を選択。
XにNO2、Yにせき・たんの有症割合 pre (prevalence) を選択する。
設定のポイントは、層別因子の設定。
層別因子は、are (area)である。
これなら傾きは異ならなそうであることが想像される。
共分散分析を実施する
「統計解析」→「連続変数の解析」→「共分散分析ANCOVA」を選択して、共分散分析を実施する。
目的変数がpre、比較する群はare、補正に用いる変数がNO2とする。
交互作用はp=0.783で統計学的有意ではなく、are別には傾きが違わないことが分かった。
areの効果は統計学的有意なので、2群で切片は異なることがわかった。
グラフを見れば一目瞭然ではある。
NO2も有意なので、傾きも意味がある。
つまりは、NO2が高くなるにつれて、せき・たんの有症割合も高くなると言える。
それを表しているのが先ほどの2本の回帰直線がひかれた散布図である。
共分散分析を実施すると自動で描かれた。
まとめ
2つの回帰直線の比較は、共分散分析を使うと簡単に実施できる。
EZRで共分散分析を実施し、グラフを描かせるのはとても簡単だ。
参考になれば。
おすすめ書籍
EZR公式マニュアル
コメント