我想以对数比例绘制一些数据,但是R绘制的密度曲线并非从我拥有的第一个极限数据点开始。
以下是一些测试数据和我使用的命令:
d = c(1,1,1,2,4,7,10,20,40,100,200,600,1000,3000,15000,23000)
plot(density(d), log='x', xlim=c(1,23000))
不幸的是,我没有足够的声誉点来发布图像,但是如果您键入这两行,它应该会出现一个图,范围从12左右到1(第一个数据点)。有解决这个问题的主意吗?
谢谢大家!
如果你看
plot(density(d), xlim=c(-100,100), type="b")
您将看到它具有正负两个特定值。x的负值不会以对数刻度显示,最小的正值是
> min( density(d)$x[density(d)$x >0] )
[1] 13.72479
造成你所看到的。您可以使用以下方法防止密度估算值达到0或以下
plot(density(d, from=0.000001), log='x', xlim=c(1,23000))
看起来像这样
本文收集自互联网,转载请注明来源。
如有侵权,请联系[email protected] 删除。
我来说两句