R中的多重图形与正常

问题描述:

我试图用正常曲线绘制实验数据(用于比较)。我可以绘制其中一条曲线,但正常曲线无法正确绘制,我只能看到一条线在屏幕上向下。我得到相同的结果为dcauchy。唯一正确绘制的曲线是“pdf_cauchy”变量。R中的多重图形与正常

这就是我要做的:

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 

考虑你的x轴是在这些图什么。您的初始绘图没有给出x轴,因此仅对100柯西值使用1:100。现在,您添加一条从0开始的正态曲线,以100为均值为0.您只会得到曲线的右半部分,并且会得到极值。平均值为50,平均值为25,你会得到更合理的东西,你也可以用dcauchy做类似的事情。

您可能确实需要在初始绘图中放置一些x值,但我不知道它们会是什么。您是否希望展示位置的峰值为0?还有其他值?只有你知道x轴的位置和范围应该是什么。

+0

我想为均值(中心)为大约0.1 – Josh

+0

我有一种感觉,它可能是与X域的一个问题,我怎么可能改变所以我只观察-5到5的值? – Josh

+0

您制作x值,只需将其粘贴到第一个绘图命令中即可。我会把它改成'x John

这是画出一个好办法:

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)