R - plot(density)

画密度分布图的时候。

1）如果想只取某一范围(x, y)的数值，不要直接把范围外的数值删掉。可用 R 中参数: from=x, to=y.

2）画密度分布图时它会根据数值范围来划bin，这样当数值范围很大时，细小的波动便看不到了。于是我们可以调bin的宽度。width＝z （z可以为0.1或更小，根据自己需要的resolution）

example：

plot(density(EastAsia, from=0, to=2, width=0.1), col="red", xaxt="n", type="l", main="TMRCA of Four Ancestries", xlab='TMRCA(million)', ylab="density")

points(density(CentralAsiaSiberia, from=0, to=2, width=0.1), col="brown", type="l")

详见 R plot(density)