R:ggplot2中的密度图与密度图

时间:2012-08-24 14:54:16

标签: r ggplot2

我正在尝试在R中做一些密度图。我最初使用密度图,但我改为ggplot2中的密度图,因为我在视觉上更喜欢ggplot2。

因此,我使用密度绘图函数绘制了密度图,并在ggplot2中进行了密度图(见下文),但我发现这些图并不相同。看起来ggplot2(右图)中的某些y值已丢失或丢失。这有什么特别的原因吗?如何使ggplot与命运图相同(左图)。

enter image description here

代码:

library(ggplot2)
library(grid)
par(mfrow=c(1,2))

# define function to create multi-plot setup (nrow, ncol)
vp.setup <- function(x,y){
 grid.newpage()
 pushViewport(viewport(layout = grid.layout(x,y)))
}

# define function to easily access layout (row, col)
vp.layout <- function(x,y){
  viewport(layout.pos.row=x, layout.pos.col=y)
}


vp.setup(1,2)

dat <- read.table(textConnection("
low high
10611.0 14195.0
10759.0 14437.0
10807.0 14574.0
10714.0 14380.0
10768.0 14448.0
10601.0 14239.0
10579.0 14218.0
10806.0 14510.0
"), header=TRUE, sep="\t")

plot(density(dat$low))

dat.low = data.frame(low2 = c(dat$low), lines = rep(c("low")))

low_plot_gg = (ggplot(dat.low, aes(x = low2, fill = lines)) +
 stat_density(aes(y = ..density..)) + 
 coord_cartesian(xlim = c(10300, 11000)) 
)
print(low_plot_gg, vp=vp.layout(1,2))

2 个答案:

答案 0 :(得分:6)

根据一些试验和错误,它看起来像你想要的

 + xlim(c(10300,11000))

而不是

 + coord_cartesian(xlim = c(10300, 11000)) 

coord_cartesian扩展了情节的限制,但并没有改变它们内部绘制的内容......

答案 1 :(得分:0)

这不是丢失价值的问题。函数plot(density())继续平滑极值,但对于您的小数据集来说并不是很准确。对于更大的数据集,两个图将是相同的。