密度图

作者

[编辑] 郑虎;

[审核] .

核密度图是一种用来观察连续型变量分布的图形。

环境配置

系统: Cross-platform (Linux/MacOS/Windows)
编程语言: R
依赖包: ggplot2; ggthemes

# 安装包
if (!requireNamespace("ggplot2", quietly = TRUE)) {
  install.packages("ggplot2")
}
if (!requireNamespace("ggthemes", quietly = TRUE)) {
  install.packages("ggthemes")
}

# 加载包
library(ggplot2)
library(ggthemes)

数据准备

# 加载数据
data <- read.delim("files/Hiplot/040-density-data.txt", header = T)

# 整理数据格式
data[,2] <- factor(data[,2], levels = unique(data[,2]))

# 查看数据
head(data)

  Value Group2
1   4.2 treat1
2  11.5 treat1
3   7.3 treat1
4   5.8 treat1
5   6.4 treat1
6  10.0 treat1

可视化

# 密度图
data["group_add_by_code"] <- "g1"

p <- ggplot(data, aes_(as.name(colnames(data[1])))) +
  geom_density(col = "white", alpha = 0.85,
               aes_(fill = as.name(colnames(data[2])))) +
  ggtitle("") +
  scale_fill_manual(values = c("#e04d39","#5bbad6")) +
  theme_stata() +
  theme(text = element_text(family = "Arial"),
        plot.title = element_text(size = 12,hjust = 0.5),
        axis.title = element_text(size = 12),
        axis.text = element_text(size = 10),
        axis.text.x = element_text(angle = 0, hjust = 0.5,vjust = 1),
        legend.position = "right",
        legend.direction = "vertical",
        legend.title = element_text(size = 10),
        legend.text = element_text(size = 10))

p