pairs用回帰直線付きpanel
pairsによる散布図にはpanelパラメータにより、panel.smoothなどを使って直線を書けますが、単純な一次回帰直線がなかったようなので作ってみました。 もし標準ライブラリにすでにある、とかpanelパラメータを使うべきでないといった情報がありましたらコメントいただければ幸いです。
関数 panel.lsfit() はpairs()のpanelパラメータに使うことを意図して作成しました。lsfit()を用いた一次回帰直線と個別データをあわせてプロットするものです。
panel.lsfit <- function(x,y,...) { f <- lsfit(x,y)$coef xx<- c ( min (x) , max(x) ) yy <- f["X"] * xx+ f["Intercept"] lines(xx,yy) points(x,y) }
関数 panel.lsfit() の引数は次の通りです。
> panel.lsfit(x, y)
x: 散布図の各点のx座標のベクタ
y: 散布図の各点のy座標のベクタ
data("swiss") pairs(swiss, panel=panel.lsfit)