一、前言

峰巒圖是核密度估計(jì)圖的變種，主要用于展示多數(shù)據(jù)系列的核密度估計(jì)圖。

1.1 示例文獻(xiàn)

文獻(xiàn)來源

原始圖片

二、R包

本期使用的R包：

• tidyverse：萬金油包；
• ggridges：繪制峰巒圖；
• ggsci：提供繪圖顏色;
• gWQS包：提供演示數(shù)據(jù)。

library(tidyverse)
library(ggridges)
library(ggsci)

三、演示數(shù)據(jù)

演示數(shù)據(jù)簡介：gWQS包中有一個(gè)內(nèi)置數(shù)據(jù)集，內(nèi)置數(shù)據(jù)集的名稱叫wqs_data，wqs_data數(shù)據(jù)集有34種多環(huán)芳烴暴露數(shù)據(jù)、25種鄰苯二甲酸酯暴露數(shù)據(jù)和其他類型數(shù)據(jù)。

本期僅使用wqs_data數(shù)據(jù)集的前5種多環(huán)芳烴暴露數(shù)據(jù)和性別。

# PCBs name
PCBs_name <- c("LBX074LA","LBX099LA","LBX105LA","LBX118LA","LBX138LA")
# get the first 5 PCBs exposure data and sex
PCBs <- wqs_data[c(PCBs_name,"sex")]
# get PCBs' absolute value
PCBs[PCBs_name] <- abs(PCBs[PCBs_name])
# view PCBs data
head(PCBs)

四、R語言實(shí)現(xiàn)

# draw plot
ggplot()+
  # geometric layer
  geom_density_ridges_gradient(data=PCBs,mapping=aes(x=LBX074LA,y=sex,fill=sex))+
  # visual mapping
  scale_fill_manual(name="Sex",labels=c("Male","Female"),values=pal_npg("nrc")(10))+
  # coordinate adjustment
  scale_y_discrete(name="Sex",labels=c("Male","Female"))+
  # theme adjustment
  theme_light()+
  theme(axis.text=element_text(color="black"),
        legend.position="top",
        legend.direction="horizontal")

五、結(jié)果解讀

NHANES數(shù)據(jù)庫中多環(huán)芳烴的編碼與對(duì)應(yīng)名稱。

PCB74在男性和女性尿液中的分布均屬于右偏型分布。

本文由mdnice多平臺(tái)發(fā)布