R:近似曲線のひき方

撒布図なんかによくある近似曲線の書き方。

例えば貝殻の殻口と殻幅について

plot(ShellLength ~ ShellWidth)

とプロットしたとすると、

 result <- lm(ShellLength ~ ShellWidth)

 abline(result,col="red")

で赤い近似曲線が引ける。簡単。

ちなみに

 text(10,25,substitute(paste(R^2,"=",x),list(x=round(summary(result)$r.squared,digits=5))))

で赤い字でR^2値が出てくる。