如何在每組中添加多個觀察值並在ggplot2 boxplot中使用組平均值？

Question

我做的哪裏y=age和x=Patient groups

age <- ggplot(data, aes(factor(group2), age)) + ylim(15, 80) 
age + geom_boxplot(fill = "grey80", colour = "#3366FF") 

我希望你能幫助我的幾件事情基本箱線圖：

1）是否有可能包括上述一些每組觀察每個組boxplot（但不是在我的組標籤的X軸上），而不必在油漆:)這樣做？ 我已經嘗試使用：

age + annotate("text", x = "CON", y = 60, label = "25") 

其中CON是第一組和y = 60是〜剛箱線圖這組以上。但是，該命令不起作用。我認爲它有一些事情要做，它把x看作一個連續的而不是一個分類變量。

2）雖然有很多關於使用均值而不是中值的問題，但我還沒有找到適合我的代碼？

3）同樣的問題，有沒有辦法可以在boxplot中包含平均羣體屬性？也許使用

age + stat_summary(fun.y=mean, colour="red", geom="point") 

然而，其中只包括一個點的意思所在。或者用

age + annotate("text", x = "CON", y = 30, label = "30") 

其中CON是第一組和y = 30是〜組齡平均。 知道如何靈活和豐富的ggplot2語法是我希望有一個更優雅的方式使用真正的統計輸出，而不是annotate。

任何建議/鏈接將不勝感激！

謝謝！

Answer 1

這是什麼你喜歡什麼？隨着stat_summary，如要求：

# function for number of observations 
give.n <- function(x){ 
    return(c(y = median(x)*1.05, label = length(x))) 
    # experiment with the multiplier to find the perfect position 
} 

# function for mean labels 
mean.n <- function(x){ 
    return(c(y = median(x)*0.97, label = round(mean(x),2))) 
    # experiment with the multiplier to find the perfect position 
} 

# plot 
ggplot(mtcars, aes(factor(cyl), mpg, label=rownames(mtcars))) + 
    geom_boxplot(fill = "grey80", colour = "#3366FF") + 
    stat_summary(fun.data = give.n, geom = "text", fun.y = median) + 
    stat_summary(fun.data = mean.n, geom = "text", fun.y = mean, colour = "red") 

黑數是觀測的數量，紅色數字是平均值。joran的回答顯示瞭如何把數字在框中 頂部

帽尖：https://stackoverflow.com/a/3483657/1036500

Answer 2

解答第一個問題。 要顯示框上方的值，您應該提供x值作爲數字而不是級別名稱。因此，繪製高於第一個值的值給x=1。

data(ToothGrowth) 
ggplot(ToothGrowth,aes(supp,len))+geom_boxplot()+ 
    annotate("text",x=1,y=32,label=30) 

Answer 3

我想這就是你正在尋找的？

myboxplot <- ddply(mtcars, 
        .(cyl), 
        summarise, 
        min = min(mpg), 
        q1 = quantile(mpg,0.25), 
        med = median(mpg), 
        q3 = quantile(mpg,0.75), 
        max= max(mpg), 
        lab = length(cyl)) 
ggplot(myboxplot, aes(x = factor(cyl))) + 
    geom_boxplot(aes(lower = q1, upper = q3, middle = med, ymin = min, ymax = max), stat = "identity") + 
    geom_text(aes(y = max,label = lab),vjust = 0) 

我才意識到我錯誤地使用了中位數當你問的意思，但你可以明顯地使用任何功能的middle審美請你。