在分析考试数据的时候，我们可能会需要分维度看听说读写的单项得分的概率分布。我们当然可以分四次画出4个图，不过ggplot的facet功能可以一次把四个图都画出来，而且是在同一个尺度上，更方便我们作对比。
啥也不说啦，直接上代码
library(readxl)
library(tidyr)
library(ggplot2)
es <- read_xlsx("~/School Level Data .xlsx")
es_long <- gather(es,nap_dim,nap_score,3:6) # 这句是重点，需要把宽数据变成长数据
ggplot(data=es_long,aes(x=nap_score,fill=nap_dim))+
geom_density(alpha=0.3)+
facet_wrap(~nap_dim) # 给听说读写不同的维度分别画概率分布曲线
## Warning: Removed 61 rows containing non-finite values (stat_density).
