比をプロットする

[2007-05-19-87]で見たように一方が他方の何倍かになっている、と結論づけたいときがある。しかし、経験的分位点分位点プロットでは正確に分からない可能性がある。そこで、y/xをプロットさせることを考える。仮説が正しければy/xは定数cになるはずである。

ykrs.qqozone<-qqplot(ykrs.ozone,stmf.ozone,plot=FALSE)$x
plot(ykrs.qqozone,(sort(stmf.ozone)/ykrs.qqozone))
lines(lowess(ykrs.qqozone,(sort(stmf.ozone)/ykrs.qqozone)f=1/2))

一行めのコマンドはqqplotを行ったときのx軸に並ぶよような数字を取り出したものと考えればよい。lowessなどについては[2007-05-19-97]を参照のこと。