以前にも同じような質問をしたことがあるので、申し訳ありませんが、回答は今までうまくいっていました。この方法を使用して見栄えの良いプロットを 6 つ作成しましたが、奇妙なプロットが 2 つ得られました。次の例を使用して、この「不適合」を確認できます。
x=c(9222,187720,42162,7005,3121,7534,21957,272901,109667,1394312,12230,69607471,79183,6389,64859,32479,3535,9414098,2464,67917,59178,2278,33064,357535,11876,21036,11018,12499632,5160,84574)
y=c(0,4,1,0,1,0,0,1,5,13,0,322,0,0,1,1,1,32,0,0,0,0,0,0,0,0,0,33,1,1)
lin=lm(y~x)
plot(x, y, log="xy")
abline(lin, col="blue", untf=TRUE)
これは、実際のデータを使用して作成したプロットです (左がログ-ログ、右が通常):
lin が 0 の値を考慮に入れていると想定していたので、欠落している 0 の値についてはあまり心配していませんでしたが、ログ プロットでわかるように、線は (1,1) の近くでも開始されていません。現在の様子から、ポイントが (1000,10) 付近にあると予想されます。
何が起こっているか知っている人はいますか？lin の係数を手動でプロットすることは役に立ちますか? もしそうなら、誰かが私にこれを行う方法を説明できますか?