[英]How to add summary statistics in histogram plot using ggplot2?
我想在使用ggplot2
制作的直方图中添加汇总统计信息。 我正在使用以下代码
#Loading the required packages
library(dplyr)
library(ggplot2)
library(reshape2)
library(moments)
library(ggpmisc)
#Loading the data
df <- iris
df.m <- melt(df, id="Species")
#Calculating the summary statistics
summ <- df.m %>%
group_by(variable) %>%
summarize(min = min(value), max = max(value),
mean = mean(value), q1= quantile(value, probs = 0.25),
median = median(value), q3= quantile(value, probs = 0.75),
sd = sd(value), skewness=skewness(value), kurtosis=kurtosis(value))
#Histogram plotting
p1 <- ggplot(df.m) + geom_histogram(aes(x = value), fill = "grey", color = "black") +
facet_wrap(~variable, scales="free", ncol = 2)+ theme_bw()
p1+geom_table_npc(data = summ, label = list(summ),npcx = 0.00, npcy = 1, hjust = 0, vjust = 1)
每个方面都有所有变量的汇总统计。 我希望它应该只显示分面变量的汇总统计信息。 怎么做?
您需要拆分 data.frame:
p1+geom_table_npc(data=summ,label =split(summ,summ$variable),
npcx = 0.00, npcy = 1, hjust = 0, vjust = 1,size=2)
或嵌套您拥有的汇总表:
summ <- summ %>% nest(data=-c(variable))
# A tibble: 4 x 2
variable data
<fct> <list<df[,9]>>
1 Sepal.Length [1 × 9]
2 Sepal.Width [1 × 9]
3 Petal.Length [1 × 9]
4 Petal.Width [1 × 9]
p1+geom_table_npc(data = summ,label =summ$data,
,npcx = 0.00, npcy = 1, hjust = 0, vjust = 1)
声明:本站的技术帖子网页,遵循CC BY-SA 4.0协议,如果您需要转载,请注明本站网址或者原文地址。任何问题请咨询:yoyou2525@163.com.