0

実験データを通常の曲線でグラフ化しようとしていました(比較用)。曲線の 1 つを問題なくグラフ化できましたが、通常の曲線を適切にプロットできず、画面に線が走っているだけでした。dcachy でも同じ結果が得られます。正しくプロットされる唯一の曲線は、「pdf_cauchy」変数です。

これは私がやろうとしていることです:

plot(pdf_cauchy, type = "l", col = "red", lwd = 4)
curve(dnorm(x, 0, 1), add = T, col = "blue", lwd = 4)
x <- seq(-5, 5, by = .001)
lines(dcauchy(x, 0, 1, FALSE), type = "l", col = "green", lwd = 4)

これが「pdf_cauchy」のデータです

c(0.00127277, 0.00132412, 0.0013786, 0.00143646, 0.00149797, 
0.00156345, 0.00163324, 0.00170772, 0.00178731, 0.00187247, 0.00196372, 
0.00206165, 0.00216691, 0.00228022, 0.00240241, 0.00253439, 0.0026772, 
0.00283202, 0.00300018, 0.00318319, 0.00338278, 0.0036009, 0.0038398, 
0.00410205, 0.00439061, 0.00470887, 0.00506072, 0.00545067, 0.00588391, 
0.00636639, 0.00690497, 0.00750753, 0.00818301, 0.00894155, 0.00979444, 
0.010754, 0.0118334, 0.0130459, 0.0144036, 0.015916, 0.0175867, 
0.0194097, 0.0213637, 0.0234058, 0.0254655, 0.0274413, 0.0292036, 
0.0306076, 0.0315168, 0.0318319, 0.0315168, 0.0306076, 0.0292036, 
0.0274413, 0.0254655, 0.0234058, 0.0213637, 0.0194097, 0.0175867, 
0.015916, 0.0144036, 0.0130459, 0.0118334, 0.010754, 0.00979444, 
0.00894155, 0.00818301, 0.00750753, 0.00690497, 0.00636639, 0.00588391, 
0.00545067, 0.00506072, 0.00470887, 0.00439061, 0.00410205, 0.0038398, 
0.0036009, 0.00338278, 0.00318319, 0.00300018, 0.00283202, 0.0026772, 
0.00253439, 0.00240241, 0.00228022, 0.00216691, 0.00206165, 0.00196372, 
0.00187247, 0.00178731, 0.00170772, 0.00163324, 0.00156345, 0.00149797, 
0.00143646, 0.0013786, 0.00132412, 0.0012
4

2 に答える 2

3

これらのグラフの x 軸を考えてみましょう。最初のプロットには x 軸が指定されていないため、100 個のコーシー値に 1:100 を使用します。ここで、0 から始まり、平均 0 で 100 になる通常の曲線を追加します。曲線の右半分と極値のみを取得します。平均が 50 で sd が 25 の場合、より合理的な結果が得られ、dcauchy で同様のことを行うこともできます。

おそらく最初のプロットにいくつかの x 値を入れる必要があるかもしれませんが、それらがどうなるかわかりません。配置のピークを 0 にしたかったのですか? ... 他の値を指定しましたか? x 軸の位置と範囲を知っているのはあなただけです。

于 2012-11-14T04:23:58.210 に答える
0

これはプロットする良い方法です:

plot(x, dcauchy(x, 0, 1, FALSE), type = "l", col = "red", lwd = 4)
par(new = TRUE)
plot(x, pdf_cauchy, type = "l", col = "blue", lwd = 2)
par(new = TRUE)
plot(x, dnorm(x, 0, 1), type = "l", col = "green", lwd = 3)
于 2012-11-14T13:40:48.030 に答える