如何在ggplot2中獲得正確的「densityplot」？

Question

我有一個數據框，其中有兩個不同的數據組的Recall值。我使用下面的代碼製作了密度圖，但我對它的外觀感到困惑。對於一個組（在這種情況下爲GE），召回非常高（始終爲1），在圖中我預計在x軸= 1時尖峯和相對較高的峯值（wrt Epi組），而在輸出中不是這種情況圖表（見附圖）。

命令我使用：

ggplot(epiGE, aes(x=Rec)) + geom_density(aes(group=Class, colour=Class, fill=Class), alpha=0.3) + theme(axis.title.x=element_blank(), axis.title.y=element_blank()) 

我也曾嘗試GGPLOT2直方圖，但有作爲在x軸= 1我不明白GE組的右側分佈，作爲barplot顯示與GE組相比，Epi組中Recall=1的值更多。

ggplot(epiGE, aes(x=Rec)) + geom_histogram(aes(group=Class, colour=Class, fill=Class), alpha=0.3) + theme(axis.title.x=element_blank(), axis.title.y=element_blank()) 

有人可以請指導我什麼，我做錯了什麼？謝謝。

這裏是我的數據dput：

> dput(epiGE) 
structure(list(Class = structure(c(1L, 2L, 1L, 2L, 1L, 2L, 1L, 
2L, 1L, 2L, 1L, 2L, 1L, 2L, 1L, 2L, 1L, 2L, 1L, 2L, 1L, 2L, 1L, 
2L, 1L, 2L, 1L, 2L, 1L, 2L, 1L, 2L, 1L, 2L, 1L, 2L, 1L, 2L, 1L, 
2L, 1L, 2L, 1L, 2L, 1L, 2L, 1L, 2L), .Label = c("Epi", "GE"), class = "factor"), 
    Rec = c(1, 1, 0.928571428571429, 1, 1, 1, 0.625, 1, 0.935897435897436, 
    1, 0.911764705882353, 1, 0.790697674418605, 1, 0.891891891891892, 
    1, 1, 1, 0.98019801980198, 1, 0.949367088607595, 1, 0.877551020408163, 
    1, 0.944444444444444, 1, 1, 1, 0.955445544554455, 1, 0.5, 
    1, 0.170731707317073, 1, 0.513513513513513, 1, 0, 1, 0.9, 
    1, 0.875, 1, 0.884615384615385, 1, 0.5, 1, 0.9, 1)), .Names = c("Class", 
"Rec"), class = "data.frame", row.names = c(NA, -48L)) 

Answer 1

這與平滑的帶寬。您可以使用參數adjust來調整帶寬。這是一個例子。

ggplot(epiGE, aes(x=Rec)) + 
    geom_density(aes(group=Class, colour=Class, fill=Class), alpha=0.3, adjust = 1/10) + 
    theme(axis.title.x=element_blank(), axis.title.y=element_blank()) 

設置adjust到1/10意味着使用默認帶寬1/10。