4

Rでいくつかの密度プロットを実行しようとしています。元々は密度プロットを使用していましたが、視覚的にggplot2を好むため、ggplot2の密度プロットに変更しました。

そこで、密度プロット関数を使用して密度プロットを実行し、ggplot2で密度プロットを実行しました（以下を参照）が、プロットが同一ではないことがわかりました。ggplot2（右のプロット）でy値の一部が失われたかドロップされたようです。これには特別な理由がありますか？ggplotを運命プロット（左のプロット）と同一にするにはどうすればよいですか。

ここに画像の説明を入力してください

コード：

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))
4

2 に答える 2

6

試行錯誤の結果、あなたが望むように見えます

 + xlim(c(10300,11000))

それよりも

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

coord_cartesianプロットの制限を拡張しますが、プロット内に描画されるものはまったく変更されません...

于 2012-08-24T15:09:50.610 に答える
0

失われた値の問題ではありません。関数plot(density())は極値の平滑化に進みますが、小さなデータセットではあまり正確ではありません。より大きなデータセットの場合、2つのプロットは同じになります。

于 2012-08-24T15:13:14.320 に答える